ex2_1_2.m

1．编写欧拉前差、后差、梯形公式。 2．编写二阶、三阶龙格库塔法通用程序。 3．编写汉明积分法通用程序. 4．编写用状态转移法对连续系统状态方程进行离散化的通用程序。

syms f
f(1)='y(2)';
f(2)='y(3)+24';
f(3)='-24*y(1)-26*y(2)-9*y(3)-192';

[y,t]=oulahoucha(f,[1 0.5 -24]',0,1,0.02)

%stairs(t,out)

subplot(2,2,1),stairs(t,y(1,:)),axis([0 1 1 1.17]),title('y')
subplot(2,2,2),stairs(t,y(1,:)),axis([0 1 1 1.17]),title('x1')
subplot(2,2,3),stairs(t,y(2,:)),axis([0 1 -0.3 0.5]),title('x2')
subplot(2,2,4),stairs(t,y(3,:)),axis([0 1 -25.5 -24]),title('x3')