Question

matlab中怎样写微分方程组

Answer 1

以下列例子，来说明怎样写微分方程组。

解  1、建立m-文件rigid.m如下：

       function dy=rigid(t,y)

       dy=zeros(3,1);

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

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

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

2、取t0=0，tf=12，输入命令：

    [T,Y]=ode45('rigid',[0 12],[0 1 1]);

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

3、结果如图

图中，y1的图形为实线，y2的图形为“*”线，y3的图形为“+”线。

Answer 2

需要将微分方程离散化