急,matlab 这是一个球体在空中飞行受空气阻力和气压升力的模型,输出的图形应该是抛物线,但是我拿MATLAB怎么也做不出来这个图,求大神帮我看下我的程序,或者帮我写写程序(最好有注释)Y1
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/04 23:08:55
![急,matlab 这是一个球体在空中飞行受空气阻力和气压升力的模型,输出的图形应该是抛物线,但是我拿MATLAB怎么也做不出来这个图,求大神帮我看下我的程序,或者帮我写写程序(最好有注释)Y1](/uploads/image/z/399351-39-1.jpg?t=%E6%80%A5%2Cmatlab+%E8%BF%99%E6%98%AF%E4%B8%80%E4%B8%AA%E7%90%83%E4%BD%93%E5%9C%A8%E7%A9%BA%E4%B8%AD%E9%A3%9E%E8%A1%8C%E5%8F%97%E7%A9%BA%E6%B0%94%E9%98%BB%E5%8A%9B%E5%92%8C%E6%B0%94%E5%8E%8B%E5%8D%87%E5%8A%9B%E7%9A%84%E6%A8%A1%E5%9E%8B%2C%E8%BE%93%E5%87%BA%E7%9A%84%E5%9B%BE%E5%BD%A2%E5%BA%94%E8%AF%A5%E6%98%AF%E6%8A%9B%E7%89%A9%E7%BA%BF%2C%E4%BD%86%E6%98%AF%E6%88%91%E6%8B%BFMATLAB%E6%80%8E%E4%B9%88%E4%B9%9F%E5%81%9A%E4%B8%8D%E5%87%BA%E6%9D%A5%E8%BF%99%E4%B8%AA%E5%9B%BE%2C%E6%B1%82%E5%A4%A7%E7%A5%9E%E5%B8%AE%E6%88%91%E7%9C%8B%E4%B8%8B%E6%88%91%E7%9A%84%E7%A8%8B%E5%BA%8F%2C%E6%88%96%E8%80%85%E5%B8%AE%E6%88%91%E5%86%99%E5%86%99%E7%A8%8B%E5%BA%8F%EF%BC%88%E6%9C%80%E5%A5%BD%E6%9C%89%E6%B3%A8%E9%87%8A%EF%BC%89Y1)
急,matlab 这是一个球体在空中飞行受空气阻力和气压升力的模型,输出的图形应该是抛物线,但是我拿MATLAB怎么也做不出来这个图,求大神帮我看下我的程序,或者帮我写写程序(最好有注释)Y1
急,matlab
这是一个球体在空中飞行受空气阻力和气压升力的模型,输出的图形应该是抛物线,但是我拿MATLAB怎么也做不出来这个图,求大神帮我看下我的程序,或者帮我写写程序(最好有注释)
Y1是Y轴,Y2是Y轴的速度,Y3是Z轴,Y4是Z轴的速度
m,b,c是参变量,y2和y4是初始速度
我的程序如下:
定义函数
function factor = route (t,y,m,b,c)
factor = [y(2);(-b)/m*sqrt(y(2)^2+y(4)^2)*y(2)+c/m*y(4);y(4);(-9.8)-b/m*sqrt(y(2)^2+y(4)^2)*y(4)-c/m*y(2)];
控制台:
m=0.43;
>> b=0.0114;
>> c=0.714;
[t,y] = ode45(@(t,y) route(t,y,m,b,c),[0 2.05],[0 10 0 10]);
这是微分方程组
急,matlab 这是一个球体在空中飞行受空气阻力和气压升力的模型,输出的图形应该是抛物线,但是我拿MATLAB怎么也做不出来这个图,求大神帮我看下我的程序,或者帮我写写程序(最好有注释)Y1
将下面所有的内容(当然除了这一行之外)复制到“新建的m文件”,F5运行即可
function matlab100000198914
% 这是一个球体在空中飞行受空气阻力和气压升力的模型,输出的图形应该是抛物线,但是我拿MATLAB怎么也做不出来这个图,求大神帮我看下我的程序,或者帮我写写程序(最好有注释)
% Y1是Y轴,Y2是Y轴的速度,Y3是Z轴,Y4是Z轴的速度
% m,b,c是参变量,y2和y4是初始速度
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5 1e-5]);
[T,Y] = ode45(@factor,[0 2],[0 0 0 0],options);
plot(T,Y(:,1),'-r',T,Y(:,2),'-.g',T,Y(:,3),'.b',T,Y(:,4),'-om');
h = legend('y1','y2','y3','y4',3);
function dy = factor(t,y)
dy = zeros(4,1); % a column vector
m=0.43; b=0.0114; c=0.714;
dy = [y(2);
-b/m*sqrt(y(2)^2+y(4)^2)*y(2)+c/m*y(4);
y(4);
(-9.8)-b/m*sqrt(y(2)^2+y(4)^2)*y(4)-c/m*y(2)];