【求助】钟院士新书上的精细积分程序来源!
本帖最后由 牛小贱 于 2014-3-18 16:08 编辑% 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; foriter=1:N Ta=2*Ta+Ta*Ta; end T=I+Ta; vk=; 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很多人都说这是钟万勰院士新书上的程序,谁能告知一下到底是哪本书?
《应用力学的辛数学方法》——钟万勰院士 牛小贱 发表于 2014-3-16 23:42
《应用力学的辛数学方法》——钟万勰院士
多谢!!!!!!{:3_49:} 牛小贱 发表于 2014-3-16 23:42
《应用力学的辛数学方法》——钟万勰院士
哥们,你对精细积分有了解么? zhousidun 发表于 2014-3-18 16:08
哥们,你对精细积分有了解么?
不是很了解,只是知道原理,懂点皮毛,但是这个是比较简单的,你可以参考一下这个帖子:http://forum.vibunion.com/forum.php?mod=viewthread&tid=105856 牛小贱 发表于 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 发表于 2014-3-20 23:47
找到源程序的出处了,现在上传。
thank yo very much!
页:
[1]