求解微分方程组matlab

 

function dy=rigid(t, y)  %其中y包含为三个未知函数对应的三个未知量，t为自变量

dy=zeros(3,1);

dy(1)=y(2)*y(3);

dy(2)=-y(1)*y(3);

dy(3)=-0.51*y(1)*y(2);



运行命名：

[T,Y]=ode45('rigid',[0 12],[0 1 1]);   %其中[0,12]为自变量取值范围，[0 1 1]为初始值



plot(T,Y(:,1),'k-',T,Y(:,2),'r*',T,Y(:,3),'b+')

也可以利用dsolve函数将他求解

下面的代码运行有错误等后面我会慢慢改正

syms y t
[x,y,z]=dsolve('Dy(1)=y(2)*y(3)','Dy(2)=-y(1)*y(3)','Dy(3)=-0.51*y(1)*y(2)','y(1)(0)=0','y(2)(0)=0','y(3)(0)=1','t')

求解解析解的时候利用dsolve函数求解，再定义自变量的取值范围，代入到刚求解的函数使用plot画图即可

具体操作如下：

syms y t
dsolve('Dy=-y+t+1','y(0)=1','t')
t=1:0.1:2
z =t + exp(-t)
plot(t,z)

 

__EOF__

本文作者：userName
本文链接：https://www.cnblogs.com/20kkk/p/16537563.html
关于博主：评论和私信会在第一时间回复。或者直接私信我。
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！