[求助]向高手求救：Gaver-stehfest 算法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问：如何利用matlab实现Gaver-stehfest 算法的拉普拉斯变换？

我写的程序，可是效果很差，请高手指正

format long
n=9;
t=1;
sum1=0;
p=0;
for k=1:n
    d1=(-1)^(n-k)*k^n/(factorial(k)*factorial(n-k))*log(2)/t...
        *factorial(2*k+4)/(factorial(k+2)*factorial(k+1));
    for m=0:k+2
        f=@(x)x*0.02+0.5*x^2*0.05^2+2*(0.5*30/(30-x)+0.5*30/(30+x)-1)-(m+k+2)*log(2)/t;
        x1=fzero(f,-60);
        x2=fzero(f,-20);
        a1=(30+x2)/30*(-x1)/(x2-x1);
        a2=(-x1-30)/30*(-x2)/(x2-x1);
        d3=t/((m+k+2)*log(2))*(a1*exp(0.9*(-x2))+a2*exp(0.9*(-x1)));
        d2=(-1)^m*factorial(k+2)/(factorial(m)*factorial(k+2-m));
        sum1=sum1+d2*d3;
    end
    p=p+d1*sum1;
end
p

[此贴子已经被作者于2008-7-5 22:34:54编辑过]

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：gave Teh Est AVE Ver 求助 高手 算法 求救

自己先顶一个

虽然没人来顶贴，但是还是很厉害。顶你！！应该有很多人；来看过