这里以三元二次常微分方程组做一个例子,更多元更高次的都类似。
比如下列方程组:
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来解方程组,系统图如下:
结果如下:
结果是完全一致的。
请发表评论