计算物理课件1-3篇.pptVIP

⑵、主程序： m=1;b=0;p=0; [t,y]=ode45(＠kqzlptfun,[0:0.001:2],[0,5,0,8],[],b,p,m); subplot(3,1,1) plot(y(:,1),y(:,3)); title(无空气阻力抛体运动,color,b) %%加标题 hold on %%保持图形窗口继续画图 subplot(3,1,2) m=1;b=0.1;p=0; [t,y]=ode45(＠kqzlptfun,[0:0.001:2],[0,5,0,8],[],b,p,m); plot(y(:,1),y(:,3)); title(空气阻力与速度一次方成正比,color,b) subplot(3,1,3) m=1;b=0.1;p=1; [t,y]=ode45(＠kqzlptfun,[0:0.001:2],[0,5,0,8],[],b,p,m); plot(y(:,1),y(:,3)); title(空气阻力与速度二次方成正比,color,b) 例2：弹簧小球的非线性运动：质量为m=0.1kg的小球，挂在原长为l=1.0m，劲度系数为k=4.8N/m的轻弹簧的一端，弹簧的另一端固定。 1、建立如图坐标。设小球某时刻t位于p点，写出方程: x y o mg v Ox: Oy: 2、令y(1)