论坛 › K01：MatLab › 精细积分MATLAB程序

zzwjd 发表于 2005-6-16 06:32:27

精细积分MATLAB程序

在钟万勰院士最新书上的原代码

% Precise Integration Method
clear;
A=zeros(2);
C=A;
D=;
B=[-6,2;2,-4];
f0=;
f1=zeros(size(f0));
H=;
I=eye(size(H));
iH=inv(H);

tf=20;
step=; % 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=;
    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

silt 发表于 2005-6-16 19:48:03

Re:精细积分MATLAB程序

钟先生的精细积分仅仅几行啊

:)

wwh8888 发表于 2005-6-20 08:22:44

Re:精细积分MATLAB程序

用matlab就是这么几行了

tto1986 发表于 2013-4-28 10:09:43

新手报道，小弟才开始看这个程序，请教楼主，jj为什么要从1:3循环呢

查看完整版本: 精细积分MATLAB程序