matlab实验报告总结Word文档11. 求一份matlab的试验报告计算方法试验报告3【实验目的】检查各种数值计算方法的长期行为【内容】给定方程组x(t)二ay(t),y(t)二bx(t),x(0)=0,y(0)=b的解是x-y平面上的一个椭圆,利用你已经知道的算法,取足够小的步长,计算上述方程的轨道,看看那种算法能够保持椭圆轨道不变。(计算的时间步长要足够多)【实验设计】用一下四种方法来计算:1.Euler法2.梯形法3.4阶RK法4.多步法Adams公式【实验过程】1.Euler法具体的代码如下:clear;a=2;b=1;A=0a;-b0;U=;u(:,1)=0;b;n=1000000;h=6*pi/n;fori=l:ndelta(i)=(u(l,i)/a厂2+(u(2,i)/b厂2厂0.5;u(:,i+1)=u(:,i)+h*A*u(:,i);endt=1:n+1;subplot(1,2,1);plot(1:n,delta);gridon;subplot(1,2,2);plot(u(1,:),u(2,:);grid