基于matlab的语音信号采集与处理课程设计报告

基于matlab的语音信号采集与处理课程设计报告内容摘要：

%绘制原始语音信号的频率响应图 title(39。 频率响应图 39。 ) figure(3) subplot(2,1,1)。 %创建 两行一列绘图区间的第 1个绘图区间 plot(abs(y1(1:512))) %做原始语音信号的 FFT频谱图 4 title(39。 原始语音信号 FFT频谱 39。 ) subplot(2,1,2)。 plot(f,abs(y1(1:512)))。 %abs是 绝对值 ， plot是 直角坐标下线性刻度曲线 title(39。 原始语音信号频谱 39。 ) xlabel(39。 Hz39。 )。 ylabel(39。 幅值 39。 )。 时域波形和频谱 ： 图 1 原始语音信号 5 图 2 语音信号频率响应图 图 3 原始语音信号 FFT 与信号频谱 2. 给原始的 语音 信号加上一个高频余弦 噪声 ，频率为 5500hz。 对 加噪后的 语音进行分析，并画出其 信号时域和频谱图。 6 程序： fs=22050。 x1=wavread(39。 C:\Documents and Settings\Administrator\桌面 \声音文件 .wav39。 )。 %读取语音信号的数据，赋给变量 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=。 %噪声幅值 d=[Au*sin(2*pi*5500*t)]39。 %所加的噪声是正弦信号 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(39。 加噪后的信号 39。 )。 xlabel(39。 time n39。 )。 %x轴的名字是“ time n” ylabel(39。 fuzhi n39。 )。 %y轴的名字是“ fuzhi n” figure(2) %创建图形窗 2 subplot(2,1,1)。 % 创建两行一列绘图区间的第 1个绘图区间 plot(f,abs(y1(1:512)))。 %做原始语音信号的频谱图 title(39。 原始语音信号频谱 39。 )。 xlabel(39。 Hz39。 )。 %x轴的名字是“ time n” ylabel(39。 fuzhi39。 )。 % y轴的名字是“ fuzhi” subplot(2,1,2)。 %创建两行一列绘图区间的第 2个绘图区间 plot(f,abs(y2(1:512)))。