急 用MATLAB编写质点的运动轨迹图像
来源:学生作业帮 编辑:百度作业网作业帮 分类:综合作业 时间:2024/07/23 04:15:23
急 用MATLAB编写质点的运动轨迹图像
函数文件编写时,要提供输入输出变量检测功能,包括变量的个数与类型的检测,若不符合相关要求给出相应错误信息;函数中应提供适当的帮助信息,内容包括函数的功能,用法等信息(中英文均可);程序中要有必要的注释信息.
![](http://img.wesiedu.com/upload/0/ad/0ad1b937fc9d96a6181d531c43a6c317.jpg)
函数文件编写时,要提供输入输出变量检测功能,包括变量的个数与类型的检测,若不符合相关要求给出相应错误信息;函数中应提供适当的帮助信息,内容包括函数的功能,用法等信息(中英文均可);程序中要有必要的注释信息.
![](http://img.wesiedu.com/upload/0/ad/0ad1b937fc9d96a6181d531c43a6c317.jpg)
![急 用MATLAB编写质点的运动轨迹图像](/uploads/image/z/8017024-40-4.jpg?t=%E6%80%A5+%E7%94%A8MATLAB%E7%BC%96%E5%86%99%E8%B4%A8%E7%82%B9%E7%9A%84%E8%BF%90%E5%8A%A8%E8%BD%A8%E8%BF%B9%E5%9B%BE%E5%83%8F)
function dx=appl(t,x,flag,v)dx=[-t/(x(1)^2+x(2)^2)*x(1)/sqrt(x(1)^2+x(2)^2);v-t/(x(1)^2+x(2)^2)*x(2)/sqrt(x(1)^2+x(2)^2)]; clear;clc;x0=[1,0]; %初值 [t,y1]=ode45('appl',[0,40],x0,[],0.8); %利用ode45求解微分方程[t,y2]=ode45('appl',[0,40],x0,[],1);[t,y3]=ode45('appl',[0,40],x0,[],1.3);[t,y4]=ode45('appl',[0,40],x0,[],sqrt(2));[t,y5]=ode45('appl',[0,40],x0,[],1.5);figure(1)plot(y1(:,1),y1(:,2),'b',y2(:,1),y2(:,2),'k',y3(:,1),y3(:,2),'r',y4(:,1),y4(:,2),'g',y5(:,1),y5(:,2),'m')%画图axis([-1 2 0 1])title('质点运动轨迹');xlabel('X坐标'); ylabel('Y坐标');legend('v=0.8','v=1','v=1.3','v=sqrt(2)','v=1.5')![](http://img.wesiedu.com/upload/1/e6/1e636ca100a94374ade9f5d358cc65b2.jpg)
![](http://img.wesiedu.com/upload/1/e6/1e636ca100a94374ade9f5d358cc65b2.jpg)