由幅度谱反求信号

zhuguanwen · 发表于 2012-6-8 20:26

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

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

x

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 · 发表于 2012-6-8 20:36

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

zhuguanwen · 发表于 2012-6-30 11:03

哎，没人指点啊

翔雨 · 发表于 2012-7-13 14:11

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

SandVNo2 · 发表于 2012-7-13 14:39

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

登录/注册后可看大图

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

也这么认为

翔雨 · 发表于 2012-7-13 14:57

回复 5 # SandVNo2 的帖子

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

SandVNo2 · 发表于 2012-7-14 01:46

回复 6 # 翔雨的帖子

哈哈哈，:lol

账号		自动登录	找回密码
密码			我要加入

[其他] 由幅度谱反求信号

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