zfx620 发表于 2012-2-16 00:30

求助fouier逆变换

在进行数值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 发表于 2012-2-16 05:34

回复 1 # zfx620 的帖子

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

zfx620 发表于 2012-2-16 10:16

回复 2 # 321forever 的帖子

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

zfx620 发表于 2012-2-16 10:42


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












321forever 发表于 2012-2-17 00:23


fs=100;%设定采样频率
N=128;
n=0:N-1;
t=n/fs;
f0=10;%设定正弦信号频率
%生成正弦信号
x=sin(2*pi*f0*t);
figure(1);
subplot(311);
plot(t,x);%作正弦信号的时域波形
%进行FFT变换并做频谱图
y=fft(x,N)/N;%进行fft变换
mag=abs(y);%求幅值
f=(0:N-1)'*100/N;%进行对应的频率转换
subplot(312);
plot(f,mag);%做频谱图
%用IFFT恢复原始信号
xifft=ifft(y)*N;
magx=real(xifft);
subplot(313);
plot(t,magx);
以前写的从sin经过fft，再ifft出来的结果一样的。希望lz参考

zfx620 发表于 2012-2-17 10:53

回复 5 # 321forever 的帖子

多谢{:{39}:}
页: [1]
查看完整版本: 求助fouier逆变换