大家看看对吗？如有错误给改改

chy19840208

本设计步骤如下：
1.录制一段自己的语音信号，并对录制的信号进行采样；画出采样后的语音信号的频谱图；
   2. 插入噪音；画出含有噪音的频谱图；
   3.给定滤波器的性能指标，采用窗函数法，划出滤波器的频域响应
   4.用该滤波器对含噪音的信号进行滤波，画出滤波后信号的频谱；
   5.回放语音信号。
程序如下：
fs=8000;
x1=wavread('E:\我的论文\论文\luyin.wav');
y1=fft(x1);
figure(1);
plot(abs(y1));
title('原始信号频谱');
xlabel('频率');
ylabel('频谱幅度');
t=0:1/8000:(length(x1)-1)/8000;
d=sin(2*pi*2000*t);
x2=x1+d;
sound(x2,fs);
y2=fft(x2);
figure(2);
plot(abs(y2));
title(‘含噪音的信号频谱’);
xlabel(‘频率’);
ylabel(‘频谱幅度’);
wp=0.25*pi;
ws=0.3*pi;
wdelta=ws-wp;
N=ceil(6.6*pi/wdelta);              %取整
wn=(0.2+0.3)*pi/2;
b=fir1(N,wn/pi,hamming(N+1));       %选择窗函数，并归一化截止频率
f2=filter(bz,az,x2)
sound(f2,8000);                    %播放滤波后的语音信号
y3=fft(f2,);
figure(3)
y2=fft(x2);
subplot(3,1,1);
plot(abs(y2));
title('滤波前的频谱')
xlabel('频率');
ylabel('频谱幅度');
subplot(3,1,2)
f2=plot(abs(y3));
title('滤波后的频谱')
xlabel('频率');
ylabel('频谱幅度');