sszhuanyi.m

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

function [ft,fmt]=sszhuanyi(A,B,T)
syms s
ft=expm(A*T);                           %solve eAT

fmt=int(expm(A*(T-s))*B,s,0,T);         %solve int(0~T)eA(T-to)Bdto