matlab求解微分方程组问题。采用m文件形式编写一个function函数，求解该方程组!!!
采用m文件形式编写一个function函数，求解该方程组
方程组为：
三个角度的初值为[1,-1,2],w0=0.001
Wx，Wy和Wz作为function函数的变量输入，输出为三个角度的数值解和三个角度对时间求导的瞬时值
第1个回答 2014-06-15
clc
wx=;
wy=;
wz=;
w0=;
f=@(t,y)([wx*cos(y(2))+wz*sin(y(2))+w0*sin(y(3));
wy+w0*cos(y(1))*cos(y(3))+tan(y(1))*(wx*sin(y(2))-wz*cos(y(2))+w0*sin(y(1))*cos(y(3)));
-(wx*sin(y(2))-wz*cos(y(2))+w0*sin(y(1))*cos(y(3)))/cos(y(1))]);
[t,Y]=ode45(f,[0 1],[1 -1 2])追问
wx=;
wy=;
wz=;
w0=;
f=@(t,y)([wx*cos(y(2))+wz*sin(y(2))+w0*sin(y(3));
wy+w0*cos(y(1))*cos(y(3))+tan(y(1))*(wx*sin(y(2))-wz*cos(y(2))+w0*sin(y(1))*cos(y(3)));
-(wx*sin(y(2))-wz*cos(y(2))+w0*sin(y(1))*cos(y(3)))/cos(y(1))]);
[t,Y]=ode45(f,[0 1],[1 -1 2])追问
主函数使用simulink编写的，这只是其中一部分，用MATLAB function模块来完成，所以要采用m文件编写function函数的，不是在command windows里直接输入的
本回答被提问者采纳
相似回答