能不能由振幅和相角重建时域序列？
s=1:8;x=fft(s);
amp=abs(x);%此处取模丧失了正负性
phs=atan2(imag(x),real(x));
a=tan(phs);%由振幅和相角还原复数的实部和虚部
rl=(amp.*amp/(1+a.*a))^0.5;
im=rl.*a;
y=abs(ifft(rl+im*i));%ifft
plot(1:length(s),s,y)
ifft结果y不等于原始信号s，因为amp=abs(x)丢失了正负，有没有什么办法可以还原？ 本帖最后由 wdhd 于 2016-9-18 14:08 编辑
原帖由 classcc 于 2009-5-31 23:49 发表
ifft结果y不等于原始信号s，因为amp=abs(x)丢失了正负，有没有什么办法可以还原？
整个程序计算中不是由“amp=abs(x)丢失了正负”，而是另有原因。在下程序中还是用了amp=abs(x)语句，能恢复原数据。
s=1:8;
x=fft(s);
amp=abs(x);%此处取模丧失了正负性
phs=atan2(imag(x),real(x));
%a=tan(phs);%由振幅和相角还原复数的实部和虚部
%rl=(amp.*amp/(1+a.*a))^0.5;
%im=rl.*a;
rl=amp.*cos(phs);
im=amp.*sin(phs);
y=real(ifft(rl+im*i));%ifft
plot(1:length(s),s,'r','linewidth',2); hold on;
plot(1:length(s),y);
原程序的错误是：
rl=(amp.*amp/(1+a.*a))^0.5;
rl是由平方根求得，永远为正值，丢失了负数。同时
y=abs(ifft(rl+im*i));
也会丢失负数（本数据都是正数，如果有负数就被丢失）。 原来是这样，非常感谢。我的基本功不好，简单推导都出错了...
页:
[1]