matlab练习程序（高阶常微分方程组数值解）

这里以三元二次常微分方程组做一个例子，更多元更高次的都类似。

比如下列方程组：

x\'\' = x\' - x + y\' -z\'

y\'\' = y\' - y - x\'

z\'\' = z\' - z + x\'

matlab代码如下：

main.m：

clear all;
close all;
clc;

[t,x]=ode45(\'testfun\',[0:0.1:5],[1 1 -1 1 0.5 0]);　　%[0:0.1:5]是t的范围,[1 1 -1 1 0.5 0]是[x(0) y(0) z(0) x\'(0) y\'(0) z\'(0)]

plot(t,x(:,1))
grid on;

testfun.m：

function dx=testfun(t,x) 　  
dx=[x(4);                   %x\'
    x(5);                   %y\'
    x(6);                   %z\'
    x(4)-x(1)+x(5)-x(6);    %x\'\'
    x(5)-x(2)-x(4);         %y\'\'
    x(6)-x(3)+x(4)];        %z\'\'

结果如下：

同样的，可以用simulink来解方程组，系统图如下：

结果如下：

结果是完全一致的。