作业帮matlab求解微分方程组问题.采用m文件形式编写一个function函数,
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/07/03 15:53:05
matlab求解微分方程组问题.采用m文件形式编写一个function函数,
采用m文件形式编写一个function函数,求解该方程组
方程组为：
三个角度的初值为[1,-1,2],w0=0.001
Wx,Wy和Wz作为function函数的变量输入,输出为三个角度的数值解和三个角度对时间求导的瞬时值
采用m文件形式编写一个function函数,求解该方程组
方程组为：
三个角度的初值为[1,-1,2],w0=0.001
Wx,Wy和Wz作为function函数的变量输入,输出为三个角度的数值解和三个角度对时间求导的瞬时值
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])
再问： 主函数使用simulink编写的，这只是其中一部分，用MATLAB function模块来完成，所以要采用m文件编写function函数的，不是在command windows里直接输入的
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里直接输入的