常微分方程数值求解实验（实验五）.docVIP

常微分方程数值求解实验（实验五） 求解如下洛仑兹方程（初值问题），并将解在相空间中动态显示方程的解。 程序：function dx=lorenz(t,x) dx=[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+35*x(2)-x(3)]; t_final=100; x0=[31;-1;-1]; [t x]=ode45(lorenz,[0,t_final],x0); plot(t,x) figure; comet3(x(:,1),x(:,2),x(:,3)) 求微分方程边值问题的数值解，并在相平面中绘制解的图形。 程序：function dy=bq(t,y) dy=[y(1)*(1-sin(y(2)));y(1)+y(2)]; [t y]=ode45(bq,[0,100],[1;1]); plot(t,y) 蜘蛛与螺线： 四只蜘蛛从一个边长为6米的正方形四个角开始爬行．每只蜘蛛都以每秒1 厘米的速度朝它右边的一个蜘蛛爬行（如图所示），结果它们都以一定的速率朝中心移动，问多少分钟后它们会在中心相遇？并绘出轨迹线。不妨改变速度，再计算相遇时间，从计算出的结果，可以得到怎样的理论结果？另外，对其它形状的正多边形（如正六边形）考虑同样的问题． 注：蜘蛛所走的路线