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

由幅度谱反求信号

F=;%设置抽样点的频率,抽样频率必须包含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=;
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 static/image/common/back.gif
ifft不是反傅里叶变换嘛?怎么是求自相关呢?

也这么认为

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

回复 5 # SandVNo2 的帖子

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

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

回复 6 # 翔雨 的帖子

哈哈哈,:lol
页: [1]
查看完整版本: 由幅度谱反求信号