oulaqiancha.m

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

function [yt,th]=oulaqiancha(f,y0,t0,td,h)
y=y0;
t=t0;
yt(:,1)=double(y0);
th(1)=double(t0);
for i=2:fix((td-t0)/h)+1
fn=double(eval(f));
yt(:,i)=yt(:,i-1)+fn'.*h;
th(i)=th(i-1)+h;
y=yt(:,i);t=th(i);
end
%n=length(fn);
%C=zeros(1,n);C(1)=1;
%out=C*yt;