MATLAB课程论文推荐.docVIP

基于MATLAB的语音信号采集与处理 1．用MATLAB对原始语音信号进行分析，画出它的时域波形和频谱 fs=22050; %语音信号采样频率为22050 x1=wavread(C:\Documents aSettings\Administrator\桌面\声音文件.wav); %读取语音信号的数据，赋给变量x1， sound(x1,22050); %播放语音信号 y1=fft(x1,1024); %对信号做1024点FFT变换 f=fs*(0:511)/1024; %将0到511，步长为1的序列的值与fs相乘并除以1024的值，赋值给f figure(1) %创建图形窗xlabel(time n); %x轴的名字是“time n” ylabel(fuzhi n); %y轴的名字是“fuzhi n” figure(2) freqz(x1) %绘制原始语音信号的频率响应图 title(频率响应图) figure(3) subplot(2,1,1); %创建subplot(2,1,2); plot(f,abs(y1(1:512))); %abs是绝对值直角坐标下线性刻度曲线xlabel(Hz); ylabel(幅值); 时域波形和频谱 给原始的语音信号加上一个高频余弦噪声，频率为加噪后的语音信号时域和频谱图fs=22050; x1=wavread(C:\Documents and Settings\Administrator\桌面\声音文件.wav); %读取语音信号的数据，赋给变量x1 f=fs*(0:511)/1024; %将0到511，步长为1的序列的值与fs相乘并除以1024的值，赋值给f t=0:1/fs:(length(x1)-1)/fs; %将0到x1的长度减1后的值除以fs的值，且步长为1/fs的值，的序列的值，赋予t Au=0.03; %噪声幅值 d=[Au*sin(2*pi*5500*t)]; %所加的噪声是正弦信号 x2=x1+d; %将正弦信号噪声加在语音信号上 sound(x2,22050); %播放语音信号 y1=fft(x1,1024); %对信号y1做1024点FFT变换 y2=fft(x2,1024); %对信号y2做1024点FFT变换 figure(1); %创建图形窗1 plot(t,x2); %做加噪后的信号时域图形（绘出t对x2的线性图） title(加噪后的信号); xlabel(time n); %x轴的名字是“time n” ylabel(fuzhi n); %y轴的名字是“fuzhi n” figure(2) %创建图形窗2 subplot(2,1,1); % 创建两行一列绘图区间的第1个绘图区间 plot(f,abs(y1(1:512))); %做原始语音信号的频谱图 title(原始语音信号频谱); xlabel(Hz); %x轴的名字是“time n” ylabel(fuzhi); % y轴的名字是“fuzhi” subplot(2,1,2); %创建两行一列绘图区间的第2个绘图区间 plot(f,abs(y2(1:512))); %做加噪后的语音信号的频谱图（abs是绝对值，plot是直角坐标下线性刻度曲线） title(加噪后的信号频谱); xlabel(Hz); %x轴的名字是“time n” ylabel(fuzhi); % y轴的名字是“fuzhi” 时域波形和频谱与原始信号对比，区别wp=0.25*pi; %通带截止频率 ws=0.3*pi; %阻带截止频率 Rp=1;