| 查看: 1220 | 回复: 1 | |||
[求助]
matlab傅里叶变换和积分问题求教
|
| 我有一个函数B(rho)，他是以rho为自变量的一个已知函数，现在需要对其进行傅里叶变换，使得S(eta)=fft(B(rho))，变换后变成了以eta为自变量的函数，现在需要这个以eta为自变量的函数与一个已知的函数MTF(eta)相乘后，在eta为（-20：20）的范围内积分，如何用matlab实现。 |
回复此楼» 猜你喜欢
大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。
已经有11人回复
-
两块石头
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
|
具体问题是这样：B(rho)是个时域的以rho为变量的函数，是是需要对其傅里叶变换后，变换到频域espolin的函数与另一个频域的函数进行相乘后积分，也就是说程序的问题在这里fft(B(eposilon)),应该怎么改？ function y=B(rho) lamda=532e-9; C_2n=1e-14; R=3000; k=2*pi/lamda; aaa=k*rho.^2/(2*R); Omega_21=1.23*C_2n*k^(7/6)*(2*R)^(11/6); eta_y=3*(1+0.69*Omega_21^(6/5)); eta_x=2.61/(1+1.11*Omega_21^(6/5)); y=exp(0.15*Omega_21*eta_x^(7/6)*hypergeom(7/6,1,-k*rho.^2/(8*R)*eta_x)... +1.265*Omega_21*(k*rho.^2/(2*R)/eta_y).^(5/12).*besselk(5/6,(k*rho.^2/(2*R)*eta_y).^0.5))-1; end v=0:50; dim=size(v,2); y=zeros(size(v)); for i=1:dim fun=@(eposilon) abs(fft(B(eposilon)).*(abs(exp(-2.2*(log(eposilon./v(i))).^2))).^2; y(i)=integral(fun,0,250); end plot(v,y); |
2楼2017-03-07 10:48:13
100