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

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

噪声 ，频率为 5500hz。 对 加噪后的 语音进行分析，并画出其 信号 时域和频谱图。 程序： 淮南师范学院电气信息工程学院 2020 届电子信息工程专业课程设计报告 fs=22050。 x1=wavread(39。 39。 )。 %读取语音信号的数据，赋给变量 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)))。 %做加噪后的语音信号的频谱图（ abs是绝对值，plot 是直角坐标下线性刻度曲线） title(39。 加噪后的信号频谱 39。 )。 xlabel(39。 Hz39。 )。 %x 轴的名字是“ time n” ylabel(39。 fuzhi39。 )。 % y 轴的名字是“ fuzhi” 淮南师范学院电气信息工程学院 2020 届电子信息工程专业课程设计报告 时域波形和频谱 ： 图 4 加噪后的信号时域图 图 5 原始信号与加噪后信号频谱对比 淮南师范学院电气信息工程学院 2020 届电子信息工程专业课程设计报告 与原始信号对比，区别 ： 先原始信号没 加噪音之前 0到 2020 有幅值，在 4000 到 6000 之间没有幅值，但是在加了噪音之后 4000 到 6000 之间出现最大幅值 12，超出正常值。 如图 5 原始信号与加噪后信号频谱对比图所示。 3. 设计合适的滤波器，滤除高频噪声，绘出滤波后的信号频域和时域波形： 程序： fs=22050。 [x,fs,Nbit]=readwav(39。 39。 )。 sound(x,fs)。 t=0:1/22050:(length(x)1)/22050。 figure(1) plot(x) title(39。 原始语音信号 39。 )。 xlabel(39。 time n39。 )。 ylabel(39。 fuzhi39。 )。 figure(2) freqz(x)。 title(39。 频率响应图 39。 )。 n=length(x)。 y1=fft(x,n)。 y2=fftshift(x)。 f=0:fs/n:fs*(n1)/n。 figure(3) subplot(2,1,1)。 plot(abs(y2))。 title(39。 原始语音信号 FFT 频谱 39。 )。