爲什麼從'trapz'函數出來的矢量是NAN？

我想要計算矢量XRECW的逆傅里葉變換。出於某種原因，我得到了NAN的矢量。請幫忙！爲什麼從'trapz'函數出來的矢量是NAN？

t = -2:1/100:2;  
x = ((2/5)*sin(5*pi*t))./((1/25)-t.^2);  
w = -20*pi:0.01*pi:20*pi;  
Hw = (exp(j*pi.*(w./(10*pi)))./(sinc(w./(10*pi)))).*(heaviside(w+5*pi)-heaviside(w-5*pi));%low pass filter  
xzohw = 0;  
for q=1:20:400 
    xzohw = xzohw + x(q).*(2./w).*sin(0.1.*w).*exp(-j.*w*0.2*((q-1)/20)+0.5);%calculating fourier transform of xzoh  
end  
xzohw = abs(xzohw);  
xrecw = abs(xzohw.*Hw);%filtering the fourier transform high frequencies  

xrect=0; 
for q=1:401 

    xrect(q) = (1/(2*pi)).*trapz(xrecw.*exp(j*w*t(q))); %inverse fourier transform 
end 
xrect = abs(xrect); 
plot(t,xrect)

來源

2012-12-08 user1887990

爲什麼不使用fft（）和ifft（）函數來計算傅立葉變換？ http://www.mathworks.de/help/matlab/ref/fft.html;jsessionid=b512fda0853222a659094479e27f –

您是否試圖找出NaN來自哪裏（通過簡化代碼）？ –

這是直接回答你的問題，爲什麼有一個南。如果你運行你的代碼，那麼Nan在第7行除以零來計算xzohw。請注意，w包含零：

>> find(w==0) 
ans = 
    2001

，您可以在您通過的w與(2./w)因素的要素分割線7看看。

一個快速修復（儘管它不能保證你的代碼能做你想做的事）是避免在w中使用避免零的步驟來包含0。由於PI肯定是不整除100，你可以嘗試服用按照.01增量步驟：

w = -20*pi:0.01:20*pi;

利用這一點，你的代碼產生可能類似於你在找什麼陰謀。爲了做得更好，我們可能需要更多關於你想要做什麼或者這些變量代表什麼的細節。

希望這會有所幫助！

來源

2012-12-08 16:35:08 Sam

爲什麼從'trapz'函數出來的矢量是NAN？

回答

相關問題