summaryrefslogtreecommitdiffstats
path: root/lectures/lec07_intro2_eulers
diff options
context:
space:
mode:
authorbnewbold <bnewbold@robocracy.org>2014-04-22 22:01:56 -0400
committerbnewbold <bnewbold@robocracy.org>2014-04-22 22:01:56 -0400
commit822bd0b04d1cce54c2fff57c4b206c51d1dcb940 (patch)
treee3cd895f9b54b7d44ff249dea3c50dda927304cc /lectures/lec07_intro2_eulers
parent222292b3aa9d967557c3165cce5cc7f5773c20ab (diff)
downloaddmmsb2014-822bd0b04d1cce54c2fff57c4b206c51d1dcb940.tar.gz
dmmsb2014-822bd0b04d1cce54c2fff57c4b206c51d1dcb940.zip
rename notes -> lectures
Diffstat (limited to 'lectures/lec07_intro2_eulers')
-rw-r--r--lectures/lec07_intro2_eulers46
1 files changed, 46 insertions, 0 deletions
diff --git a/lectures/lec07_intro2_eulers b/lectures/lec07_intro2_eulers
new file mode 100644
index 0000000..17630d2
--- /dev/null
+++ b/lectures/lec07_intro2_eulers
@@ -0,0 +1,46 @@
+
+euler's method:
+
+dx/dt =~ ( x(t + Dt) - x(t) ) / Dt, for very small Dt
+
+so, x(t + Dt) = x(t) + f(x) * Dt, which is how to integrate the system
+
+"if Dt is too large, becomes highly unstable" (duh)
+
+scary MATLAB advice:
+- watch out of totally crazy values (very high)
+- non-negative values go negative
+
+MATLAB built-in ODE solvers: ode23, ode15s
+
+runge-kutta (use dxdt from between t_n and t_(n+1) )
+variable time-step methods
+
+summary: Euler's method sucks, news at 11.
+
+### misc julia notes
+
+the only place to find history (?!?!!) is in ~/.julia_history
+
+### codes
+
+using Winston
+
+a=20
+b=2
+c=5
+dt = 0.05
+tlast = 2
+
+iterations = int( round(tlast/dt) )
+xall = zeros(iterations, 1)
+x = c
+
+for i = 1:iterations
+ xall[i] = x
+ dxdt = a - b*x
+ x = x + dxdt*dt
+end
+
+time = dt * [0:iterations-1]'
+plot(time,xall)