求助fouier逆变换

zfx620

在进行数值ifft时出现了困难，希望各位帮忙！
例如：
时间域函数：a=exp(-x^2);
起fourier变换的频率域表达式为：A=pi/2*exp(-1/4*w^2);其中w是对应频域中的变量。
现在用ifft进行验算：
x=1:10;
w=(1:10)/10;
a=exp(-x.^2);
A=pi/2*exp(-1/4*w.^2);
B=fft(a);
结果a不等于B??
以及
w=1:10;
x=(1:10)/10;
a=exp(-x.^2);
A=pi/2*exp(-1/4*w.^2);
B1=ifft(A);
结果A不等于B1??

如果能简单代替的话，则应该相同，问题出在哪里呢，希望达人告知

321forever

回复 1 # zfx620 的帖子

1. lz好像大小写的a没有写清楚。
2.Y = fft(X) returns the discrete Fourier transform (DFT) of vector X，lz的A是不是连续傅立叶变换的结果呢

zfx620

回复 2 # 321forever 的帖子

A是a的连续傅立叶变换的结果
那么离散的FFT HE IFFT结果应该是一致的呀？

zfx620

换个问题，对于一个微分方程，进行fourier变换后，最终表达式中含有变换参数w，在进行ifft时，w如何选取，有的时候w为分母是w不能从零开始，这种情况下该怎么办？？？？

321forever

1. 
   
2. fs=100;%设定采样频率
   
3. N=128;
   
4. n=0:N-1;
   
5. t=n/fs;
   
6. f0=10;%设定正弦信号频率
   
7. %生成正弦信号
   
8. x=sin(2*pi*f0*t);
   
9. figure(1);
   
10. subplot(311);
    
11. plot(t,x);%作正弦信号的时域波形
    
12. %进行FFT变换并做频谱图
    
13. y=fft(x,N)/N;%进行fft变换
    
14. mag=abs(y);%求幅值
    
15. f=(0:N-1)'*100/N;%进行对应的频率转换
    
16. subplot(312);
    
17. plot(f,mag);%做频谱图
    
18. %用IFFT恢复原始信号
    
19. xifft=ifft(y)*N;
    
20. magx=real(xifft);
    
21. subplot(313);
    
22. plot(t,magx);

复制代码

以前写的从sin经过fft，再ifft出来的结果一样的。希望lz参考

zfx620

回复 5 # 321forever 的帖子

多谢