MATLAB 求卷积

4567

收藏 2009-10-05

请高手指点，我想用matlab求两个函数的卷积，X~N(0.3, 0.3^2), T~0.7t(4)，这里t(4)是自由度为4的t分布。
X~f(x)=1/(sqrt(2pi)*0.3)*exp(-(x-0.3)^2/(2*0.3^2) ); 0.7T~h(t)=Gamma(5/2)/Gamma(2)*(4pi)^(-1/2) * [1+t^2/(4*0.7^2) ]^(-5/2) /0.7. 求和Y=X+T的分布，X 和 T 独立。
我用matlab算他们的卷积，但是matlab只能求离散的，这样会不太精确。还有，这样的傅里叶变换怎么变的。

扫码加我拉你入群

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

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

全部回复

liujun197987

2009-10-6 10:11:37

上MATLAB论坛去找

扫码加我拉你入群

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

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

ltx1215

2009-10-9 00:16:24

这是我写的代码
syms x t;%输入两个符号数
f=1/(sqrt(2pi)*0.3)*exp(-(x-0.3)^2/(2*0.3^2) ); %输入f（x）
h=gamma(5/2)/gamma(2)*(4*pi)^(-1/2)*[1+t^2/(4*0.7^2)]^(-5/2)/0.7;%输入 h（t）
Y=fourier(f)*fourier(h);%进行傅里叶变换后相乘，因为时域卷积等于频域相乘
ifourier(Y)%对结果进行反傅里叶变换就得到了你想要的东东

扫码加我拉你入群

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

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

xunmeng_7

2009-10-10 23:44:11

3# ltx1215 但是还是有积分号，这得到的是个密度函数，我要求它在(y1,y2,.....,yn)的值，求不出来呀。下面是算完逆fourier变换后的结果：
ifourier(Y)=129414594274953184449/144115188075855872/pi^(3/2)*2^(1/2)*int(20/3*pi^(3/2)/(49+25*_U1^2)^(5/2)*2^(1/2)*exp(-1/18*(10*x+10*_U1-3)^2),_U1 = -Inf .. Inf)

扫码加我拉你入群

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

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

hucs

2010-5-15 15:07:06

4# xunmeng_7
后来问题解决没?

扫码加我拉你入群

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

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

hucs

2010-5-15 15:16:09

我也碰到类似的问题，是关于laplace逆变换的......

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群