急急急!毕业设计用!MATLAB中ode45有循环怎么办?

发布网友发布时间：2022-10-05 22:50

我来回答

共2个回答

热心网友时间：2023-10-21 01:56

% ode45没有像你那么用的！具体doc ode45查看帮助！

% 你就是想在一幅图上画出k分别为0和0.1的

% plot3(y(:,1),y(:,3),y(:,5))图像是不？可以这样如下，

% 而且你的g未定义！

function dy=rigid(t,y)

global m g k w r ;

dy=zeros(6,1);

dy(1)=y(2);

dy(2)=2*w*y(4)*sin(r)-k*y(2);

dy(3)=y(4);

dy(4)=-2*w*(y(6)*cos(r)+y(2)*sin(r))-k*y(4);

dy(5)=y(6);

dy(6)=2*w*y(4)*cos(r)-k*y(6)-g;

%以上的微分方程没有变

global w r k g;

w=0;

r=pi/4;

g=10;%我自己随意定义的

t0=linspace(0,21);

a=200*cos(pi/4);

b=200*cos(pi/4);

c=200*cos(pi/4);

kk=[0,0.1];

color=['r','b'];%我定义了颜色，要不全是蓝色，k=0对应红色k=0.1对应蓝色

for i=1:2

k=kk(i);

[t,y]=ode45('rigid',t0,[0 a 0 b 0 c]);

plot3(y(:,1),y(:,3),y(:,5),color(i))

hold on

end

title('考虑和不考虑空气阻力的运动轨迹');

xlabel('x')

ylabel('y')

zlabel('z')

grid on

legend('考虑空气','不考虑空气')

热心网友时间：2023-10-21 01:56

错误
1.k=[0,0.1]在子函数里不能用，k必须是一个数，或者跟t同样长的数组；
2.子函数最后一行用到了g，而g是多少没有定义；
3.画图的hold on用在了循环外面，感觉是不是你搞错了
5..主函数里面多了个end
针对此几个错误我给改了下
global k w r g i;%循环子函数用到g，i,加了两个全局变量
k=[0,0.1];
w=0;%w=0，所以子函数里面乘以w的项我给删掉了，你可以加上
r=pi/4;
g=5;%给g赋了个值
t=linspace(0,21)';
a=200*cos(pi/4);
b=200*cos(pi/4);
c=200*cos(pi/4);
% figure
for i=1:2
[t,y]=ode45('rigid',t,[0 a 0 b 0 c],[],k(i),w,r);
plot3(y(:,1),y(:,3),y(:,5));%这一句和下一句放在循环里面
hold on
end

title('考虑和不考虑空气阻力的运动轨迹');
xlabel('x')
ylabel('y')
zlabel('z')
hold off
grid on
%%%%%%%%%%%%%%%%%%%%%%%%%
function dy=rigid(t,y)
global m g k w r i ;%响应的子函数全局变量也要改变
dy=zeros(6,1);
dy(1)=y(2);
dy(2)=-k(i)*y(2);%这个用k(i)
dy(3)=y(4);
dy(4)=-k(i)*y(4);%这个也要用k(i)
dy(5)=y(6);
dy(6)=-g;