matlab求解微分方程组问题.采用m文件形式编写一个function函数,采用m文件形式编写一个function函数,求解该方程组方程组为:三个角度的初值为[1,-1,2],w0=0.001Wx,Wy和Wz作为function函数的变量输入,输

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/18 09:55:21
matlab求解微分方程组问题.采用m文件形式编写一个function函数,采用m文件形式编写一个function函数,求解该方程组方程组为:三个角度的初值为[1,-1,2],w0=0.001Wx,W

matlab求解微分方程组问题.采用m文件形式编写一个function函数,采用m文件形式编写一个function函数,求解该方程组方程组为:三个角度的初值为[1,-1,2],w0=0.001Wx,Wy和Wz作为function函数的变量输入,输
matlab求解微分方程组问题.采用m文件形式编写一个function函数,
采用m文件形式编写一个function函数,求解该方程组
方程组为:

三个角度的初值为[1,-1,2],w0=0.001
Wx,Wy和Wz作为function函数的变量输入,输出为三个角度的数值解和三个角度对时间求导的瞬时值

matlab求解微分方程组问题.采用m文件形式编写一个function函数,采用m文件形式编写一个function函数,求解该方程组方程组为:三个角度的初值为[1,-1,2],w0=0.001Wx,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])