由幅度谱反求信号
F=;%设置抽样点的频率,抽样频率必须包含0和1,单位为piF1=(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先求自相关再求信号,然后求信号的傅里叶变换的频域图,为啥本来应该衰减求出来的信号是增加的呢?
哦,图片的横坐标单位都是Hz 哎,没人指点啊 ifft不是反傅里叶变换嘛?怎么是求自相关呢?
翔雨 发表于 2012-7-13 00:11 static/image/common/back.gif
ifft不是反傅里叶变换嘛?怎么是求自相关呢?
也这么认为 回复 5 # SandVNo2 的帖子
呵呵。。你和我一样是菜鸟!你再看看其他贴,就知道FFT和IFFT滤波怎么一回事了
回复 6 # 翔雨 的帖子
哈哈哈,:lol
页:
[1]