声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2869|回复: 8

[FFT] 时域信号需要满足什么条件,才能利用IFFT进行恢复?

[复制链接]
发表于 2013-8-16 13:52 | 显示全部楼层 |阅读模式

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

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

x
下面是两段程序,其中红色字体为唯一的不同之处,程序一能恢复而程序二不能恢复的原因是什么?
程序一
clear;
N=1000;
fs=10000000;
n=0:N-1;
t=n/fs;
fc=200000;
y=cos(2*pi*fc*t);
figure
plot(t,y);
grid;
df=fs/N;
n2=1:N/2+1;
ff=(n2-1)*df;
x=fft(y,N);
z=abs(x(n2));
pyy=angle(x(n2));
Pyy=180/pi*pyy;
z(2:end)=z(2:end).*exp(j*Pyy(2:end));
z=[z conj(z(end-1:-1:2))];
xr=real(ifft(z,'symmetric'));
figure;
plot(t,xr,'r');
grid;

程序二
clear;
N=1000;
fs=10000000;
n=0:N-1;
t=n/fs;
fc=200000;
t0=2/fc;
y=exp(-(t-0.0000250).^2/t0^2).*cos(2*pi*fc*t);
figure
plot(t,y);
grid;
df=fs/N;
n2=1:N/2+1;
ff=(n2-1)*df;
x=fft(y,N);
z=abs(x(n2));
pyy=angle(x(n2));
Pyy=180/pi*pyy;
z(2:end)=z(2:end).*exp(j*Pyy(2:end));
z=[z conj(z(end-1:-1:2))];
xr=real(ifft(z,'symmetric'));
figure
plot(t,xr,'r');
grid;

点评

能进行FFT变换的时域信号就能进行IFFT变换。 其实没有什么FFT变换,就是傅里叶变换!  发表于 2013-8-16 16:00
回复
分享到:

使用道具 举报

 楼主| 发表于 2013-8-16 19:56 | 显示全部楼层
程序二中经过IFFT恢复的时域信号与原始信号不一样,为什么呢?
发表于 2013-8-16 20:11 | 显示全部楼层

贴图上来看!
不过我也不太懂!
检查自己的程序是否正确。
不可能。对于随机信号,逆变换后时域上完全匹配貌似不可能(牛人给解释);对于一般的正弦信号几乎完全重合(理所当然认为的)
发表于 2013-8-16 22:39 | 显示全部楼层
路过,学习一下
 楼主| 发表于 2013-8-18 10:22 | 显示全部楼层
yghit08 发表于 2013-8-16 20:11
贴图上来看!
不过我也不太懂!
检查自己的程序是否正确。

程序一

程序一

程序二

程序二

这就是两个程序的图片,原因好像是程序二中的时域信号为非周期信号,不知道对不对,谢谢yghit08的支持!

点评

第一个正确逆变换。 第二个从图上看有泄露。非周期段出现在延时段和迅速衰减段,具体原因有待论坛内其他人给你解释,程序貌似没问题。不妨试试,变换后不做其他操作直接做逆变换看看(没做工,不知道逆变换的操作)  发表于 2013-8-18 10:52
发表于 2014-7-11 15:44 | 显示全部楼层
程序错了,这句exp(j*Pyy(2:end)),里面的角度因该是弧度不是。改成pyy就对了。

评分

1

查看全部评分

发表于 2014-7-11 18:07 | 显示全部楼层
傅里叶变换只能处理周期平稳信号,正逆变换,前后基本不会变化。

fft之后的频域信号,只记录原信号的幅值与相位。

当信号为非平稳时,ifft会把信号当做全时间信号恢复出来。

你可以试一下,前段是周期信号,后段是0的信号。ifft恢复出来的依然是全时间段的周期信号。

傅里叶变换是一种全局变换,处理非平稳信号很捉襟见肘。

你可以用短时傅里叶变换和逆短时傅里叶变换,对非平稳信号进行处理。恢复效果会很好。

评分

1

查看全部评分

发表于 2014-7-21 18:55 | 显示全部楼层
采用短时傅立叶变换,试试看效果会不会好点
发表于 2014-8-5 20:15 | 显示全部楼层
yghit08 发表于 2013-8-16 20:11
贴图上来看!
不过我也不太懂!
检查自己的程序是否正确。

能不能邀请我去清影PT 我的邮箱liubao12398@126.com
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-25 06:36 , Processed in 0.078317 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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