声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1310|回复: 5

[综合讨论] 求助fouier逆变换

[复制链接]
发表于 2012-2-16 00:30 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
在进行数值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??

如果能简单代替的话,则应该相同,问题出在哪里呢,希望达人告知
回复
分享到:

使用道具 举报

发表于 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是不是连续傅立叶变换的结果呢
 楼主| 发表于 2012-2-16 10:16 | 显示全部楼层
回复 2 # 321forever 的帖子

A是a的连续傅立叶变换的结果
那么离散的FFT HE IFFT结果应该是一致的呀?
 楼主| 发表于 2012-2-16 10:42 | 显示全部楼层

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












发表于 2012-2-17 00:23 | 显示全部楼层

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

以前写的从sin经过fft,再ifft出来的结果一样的。希望lz参考
 楼主| 发表于 2012-2-17 10:53 | 显示全部楼层
回复 5 # 321forever 的帖子

多谢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 08:40 , Processed in 0.098171 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表