由幅度谱反求信号

zhuguanwen

F=[0:1/(Fs/2):1-1/(Fs/2)];%设置抽样点的频率，抽样频率必须包含0和1，单位为pi
F1=(1:500);
F2=(501:length(F));
A1=120*ones(1,length(F1));
A2=120-(120-93.4248)*log(F2/500)/log(Fs/2/500);
A=[A1 A2];
Rx=ifft(A);%由功率谱A求自相关
y3=levinson(Rx);%Levinson-Durbin算法求信号
Len=max(size(y3));
Y3=fft(y3)/(Len/2);
Y3=abs(Y3);
f3=(0:Len-1)*Fs/Len;
figure
plot(f3(1:Len/2),A(1:Len/2));
figure
plot(f3(1:Len/2),Y3(1:Len/2));%y频谱图

这是幅值为A的图，


这是由A先求自相关再求信号，然后求信号的傅里叶变换的频域图，为啥本来应该衰减求出来的信号是增加的呢？

zhuguanwen

哦，图片的横坐标单位都是Hz

zhuguanwen

哎，没人指点啊

翔雨

ifft不是反傅里叶变换嘛？怎么是求自相关呢？

SandVNo2

翔雨 发表于 2012-7-13 00:11

                               
登录/注册后可看大图

ifft不是反傅里叶变换嘛？怎么是求自相关呢？

也这么认为

翔雨

回复 5 # SandVNo2 的帖子

呵呵。。你和我一样是菜鸟！你再看看其他贴，就知道FFT和IFFT滤波怎么一回事了

SandVNo2

回复 6 # 翔雨 的帖子

哈哈哈，:lol