【求助】钟院士新书上的精细积分程序来源！

zhousidun

% 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;

  foriter=1:N

   Ta=2*Ta+Ta*Ta;

  end

T=I+Ta;

vk=[0;0;0;0];

  foriter=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

很多人都说这是钟万勰院士新书上的程序，谁能告知一下到底是哪本书？

牛小贱

《应用力学的辛数学方法》——钟万勰院士

zhousidun

牛小贱 发表于 2014-3-16 23:42
《应用力学的辛数学方法》——钟万勰院士

多谢！！！！！！

zhousidun

牛小贱 发表于 2014-3-16 23:42
《应用力学的辛数学方法》——钟万勰院士

哥们，你对精细积分有了解么？

牛小贱

zhousidun 发表于 2014-3-18 16:08
哥们，你对精细积分有了解么？

不是很了解，只是知道原理，懂点皮毛，但是这个是比较简单的，你可以参考一下这个帖子：http://forum.vibunion.com/forum.php?mod=viewthread&tid=105856

zhousidun

牛小贱 发表于 2014-3-18 16:14
不是很了解，只是知道原理，懂点皮毛，但是这个是比较简单的，你可以参考一下这个帖子：http://forum.chi ...

《应用力学的辛数学方法》0.2.8公式里面的r0和r1是什么啊？

沧海一声笑

zhousidun 发表于 2014-3-18 17:23
《应用力学的辛数学方法》0.2.8公式里面的r0和r1是什么啊？

楼主，这种问题，你可以自己琢磨一下，自己看一下嘛……虽然没读过这本书，但是感觉应该不难的样子啊

zhousidun

找到源程序的出处了，现在上传。

xinglong-liu

zhousidun 发表于 2014-3-20 23:47
找到源程序的出处了，现在上传。

thank yo very much!