关于音频信号的处理，谢谢

小月点刀

小弟最近对DSp技术比较着魔，接了一份作业回来，但是自己做很艰难，我已经翻遍了程老师的《数字信号处理》这本书，无奈还是有很多看不懂，请大家帮帮忙，帮我写一下代码，最好注明一下，让我能看懂哈，谢谢了

音频信号的谱分析及去噪

要求：

l
录制一段音频（如歌曲，说话声等），采用Matlab工具对此音频信号用FFT作谱分析。
l
录制一段加入噪声的音频（如在歌声中加入尖锐的口哨声或者其他噪声），采用Matlab工具对此音频信号用FFT作谱分析。
l
选择合适的指标，设计FIR数字滤波器，将音频中加入的噪声信号减弱或滤除。

我目前只是做到能进行简单的谱分析，但是FFT加入还是会错
大家帮我看看

%进行FFT频谱分析
N=2048;
[x1,fs,NBITS]=wavread('2.wav');
[m,n]=size(x1);
y1=fft(x1,N);           
%y1 = fftshift(y1);
f=fs*(1:1024)/2048;
figure(1)
plot(linspace(0,m/fs,m),x1);  
title('原始语音信号');
xlabel('时间');
ylabel('幅值');
grid on
figure(2)
subplot(2,1,1);
plot(abs(y1(1:1024)));     
subplot(2,1,2);
plot(f,abs(y1(1:1024)));
title('原始语音信号频谱')
xlabel('Hz');
ylabel('幅值');

拜托大家了～！

[ 本帖最后由 eight 于 2007-7-12 10:52 编辑 ]

zhlong

请到信号处理版看看频谱分析、滤波相关的帖子。