精细积分MATLAB程序

2015-11-18 04:13| 发布者: aspen| 查看: 1232| 评论: 4|原作者: AaronSpark|来自: 声振论坛

摘要: 在钟万勰院士最新书上的原代码 % Precise Integration Method clear; A=zeros(2); C=A; D=; B=; f0=; f1=zeros(size(f0)); H=; I=eye(size(H)); iH=inv(H); tf=20; step=; % different step size N=20; figure; ...

在钟万勰院士最新书上的原代码
% Precise Integration Method
clear; A=zeros(2); C=A; D=[0.5,0;0,1]; B=[-6,2;2,-4]; f0=[0;0;0;10]; f1=zeros(size(f0));
H=[A,D;B,C]; I=eye(size(H)); iH=inv(H); tf=20; step=[2,0.5,0.1];
% different step size
N=20; figure; hold; str=['o','x','b-'];
for jj=1:3 %PIM begin
dt=step(jj)/2^N; Ta=H*dt+(H*dt)^2*(I+(H*dt)/3+(H*dt)^2/12)/2;
for iter=1:N, Ta=2*Ta+Ta*Ta; end
T=I+Ta; vk=[0;0;0;0];
for iter=1:tf/step(jj)
iter, t(:,iter)=step(jj)*(iter-1); v(:,iter)=vk(1);
vk=T*(vk+iH*(f0+iH*f1))-iH*(f0+iH*f1+f1*step(jj));
end
% PIM end
% figure(jj); plot(t(1:tf/step(jj)),v,str(jj));
end

[ 本帖最后由 ChaChing 于 2010-7-11 15:09 编辑 ]

收藏分享邀请

上一篇：maker系列--最速下降法求线性方程下一篇：[转帖]bp神经网络源程序

发表评论

账号		自动登录	找回密码
密码			我要加入

精细积分MATLAB程序

最新评论

相关分类