[hrt]

首先建立函数M文件如下：
function dx=fmat(t,x) % x=[x1,x2,x3,x4],其中x3=dx1、x4=dx2；
dx=[x(3);x(4);x(2)-x(1);30+x(1)-2*x(2)];
end

再在命令行中输入：
>> t=0:100;
>> x0=[12,23,0,0];
>> [t,x]=ode45(@fmat,t,x0);
>> [t,x]

可得到结果。
~ ~