数字信号处理陈后金版作业.docVIP

数字信号处理作业： 基于MATLAB的有噪声的语音信号分析与处理设计 一、题目要求： 选择一个语音信号作为分析对象，或录制一段语音信号； 对语音信号进行采样，画出采样后语音信号的时域波形和频谱图； 利用MATLAB中的正弦函数产生噪声加入到语音信号中，使语音信号被污染，然后进行频谱分析； 设计FIR和IIR数字滤波器，并对被噪声污染的语音信号进行滤波，画出滤波前后信号的时域波形和频谱，并对滤波前后的信号进行比较，分析信号的变化； 进行信号频谱的搬移，进行变声实验 回放语音信号 二、题目的实现方法 语音信号的采集 利用PC 机上的声卡和WINDOWS 操作系统可以进行数字信号的采集。将话筒输入计算机的语音输入插口上,启动录音机。按下录音按钮,接着对话筒说话“语音信号处理”,说完后停止录音,屏幕左侧将显示所录声音的长度。点击放音按钮,可以实现所录音的重现。以文件名“sound”保存入g :\the sound文件夹中。可以看到,文件存储器的后缀默认为. wav ,这是WINDOWS 操作系统规定的声音文件存的标准。 语音信号的时频分析 Matlab软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数 Wavread 函数调用格式 y=wavreadfile)%读取file所规定的wav文件，返回采样值放在向量y中。 [y,Ft,nbits]=wavread(file) %采样值放在向量y中，fs表示采样频率（z），nbits表示采样位数。 y=wavreadfile，N%读取钱N点的采样值放在向量y中。 y=wavreadfile，[N1,N2]%读取从N1到N2点的采样值放在向量y中。 对语音信号shengyin.wav进行采样其程序如下： [y,Ft,nbits]=wavered (shengyin); %把语音信号进行加载入Matlab 仿真软件平台中 Ft = 22050 nbits = 16 首先画出语音信号的时域波形，然后对语音信号进行频谱分析。在matlab中利用fft对信号进行快速傅里叶变换，得到信号的频谱特性。 其程序如下： [x,Ft,bits]=wavread(shengyin);%读入语音信号,返回采样频率Fs x_length=length(x);%计算语音信号的长度 sound(x,Ft,bits);%以Fs的频率播放语音信号,Fs的改变可以改变声调，以16位的精度重新播放语音信号播放 subplot(311); plot(x);%x是采集得到的离散声音信号 xlabel(抽样时刻t); title(原始信号波形); grid; fz=fftshift(fft(x,x_length));%计算语音离散信号x的x_length点的DFT，并重新排列,把DFT移位到对称位置； W=linspace(-pi,pi,x_length); f=W*Fs/(2*pi); subplot(312); plot(f,abs(fz));%画出幅频图 xlabel(频率); title(原始信号的频谱); grid; subplot(313); plot(f,angle(fz)); xlabel(频率); title(原始信号的相位); grid; 程序运行结果如下图1： 图1 原始信号的波形 语音信号加噪与频谱分析 用正弦函数sin(wt)对原始语音信号进行加噪声处理。 程序如下： clf;%清除原先曲线波形 t=0:length(x)-1; Fn=1500;%干扰正弦信号的频率 omiga=2*pi*Fn/Ft; noise1=0.1*sin(omiga*t);%加入正弦干扰信号，频率1.5khz x1=x+noise1; subplot(311); plot(t,x1); grid; xlabel(时间（t）); title(加噪声的信号波形); fz1=fftshift(fft(x1,length(x))); subplot(312); plot(f,abs(fz1));%画出加入高频干扰之后的语音信号幅频图 xlabel(频率); title(加噪声的信号频谱); grid; subplot(313); plot(f,angle(fz)); xlabel(频率); title(加噪声的信号相位); grid on; sound(x1,Ft);%播放加入干扰信号后的声音信号 程序运行结果如下图2： 说明：由程序运行结果可以看到，加入噪声后信号的时域信号波形里面已经有噪声成分，可以看到频谱图中两个幅值很大的尖峰，其频率是1500Hz，说明程序加噪声已经成功，播放加噪声的语音信号时听到“滴滴”的噪声信号。 根据语音信号的特点给出有关滤波器的新能指标： 设计FIR数