freq_analyse滤波前后频谱分析图.m

来自「RLS算法」· M 代码 · 共 41 行

%频谱图绘制
%dxt
%EI，HUST
%2008.12
clear all
clc
[xn,Fsy,bitsy] = wavread('signalnoise.wav');%读入信号噪声文件
y_single=xn(:,1)';                          %读取一个声道
X=fft(y_single');                           %进行fft变换
Mag=abs(X);                                 %求幅值
f=(0:length(X)-1)*Fsy/length(X);            %进行对应的频率转换
subplot(3,1,1);
plot(f,Mag)                                 %画频谱图
xlabel('f(Hz)');
ylabel('幅值');
title('混有噪声的信号频谱');

[n,Fsy,bitsy] = wavread('noise');           %读入噪声文件
n_single=n(:,1)';                           %读取一个声道
X=fft(n_single');                           %进行fft变换

Mag=abs(X);                                 %求幅值
f=(0:length(X)-1)*Fsy/length(X);            %进行对应的频率转换
subplot(3,1,2);
plot(f,Mag,'r');                            %做频谱图
xlabel('f(Hz)');
ylabel('幅值');
title('噪声频谱');

[s,Fsy,bitsy] = wavread('filter');         %读入滤波后的语音文件
s_single=s(:,1)';                          %单声道
X=fft(5*s_single');                        %进行fft变换
Mag=abs(X);                                %求幅值
f=(0:length(X)-1)*Fsy/length(X);           %进行对应的频率转换
subplot(3,1,3);
plot(f,Mag)                                %做频谱图
xlabel('频率(Hz)');
ylabel('幅值');
title('滤波后的信号频谱');