fs=22050; %语音信号的采样频率为 22050 x1=wavread(f:轻音乐-古筝.wav); %读取语音信号的数据，赋给变量 x1 sound(x1,22050); %播放语音信号 y1=fft(x1,1024); %对信号做 1024 点 FFT 变换 figure(1) plot(x1) %做原始语音信号的时域图形 title(原始语音信号); xlabel(time n); ylabel(fuzhi n); figure(2) freqz(x1) %绘制原始语音信号的频率响应图 title(频率响应图) figure(3) subplot(2,1,1); plot(abs(y1(1:512); %做原始语音信号的 FFT 频谱图 title(原始语音信号 FFT 频谱) subplot(2,1,2); plot(f,abs(y1(1:512); title(原始语音信号频谱) xlabel(Hz); ylabel(fuzhi);clear; fs=22050; x1=wavread(f:轻音乐-古筝.wav); f=fs*(0:511)/1024; t=0:1/22050:(length(x1)-1)/22050; %将所加噪声信号的点数调整到与原始信 号相同 Au=1 d=0.5*(2*pi*1000*t); %噪声为 1KHz 的余弦信号cos x2=x1+d; sound(x1,8000); pause(40); sound(x2,22050); %播放加噪声后得语音信号 y2=fft(x2,1024); figure(1) plot(t,x2) title(加噪后的语音信号); xlabel(time n);ylabel(fuzhi n); figure(2) subplot(2,1,1); plot(f,abs(x1(1:512); title(原始语音信号频谱); xlabel(Hz); ylabel(fuzhi);