数字信号实验报告1.docVIP

数字信号处理导论实验报告 姓名：金涛 学号：201005090209 实验一 信号、系统及系统响应 当系统的输入输出差分方程为： Y(n)-0.8y(n-1)-0.5y(n-2)=0.7x(n)+0.3x（n-1） 通过 MATLAB 编程实现并考虑如下问题： （1）当系统的输入为单位冲激函数时，分别利用filter 函数和impz 函数得到的系统单位冲激响应曲线。 （2）当系统的输入为单位阶跃函数时，分别利用filter 函数和impz 函数得到的系统单位阶跃响应曲线。 （3）对于不同的输入，系统的输出有什么变化，试讨论之。 +1)] a=[a(1),a(2),...,a()] 并要求a（1）=1，如果a(1)≠1,则程序自动将其规划为1 实验内容 源程序 x=ones(100);t=1:100; %产生单位阶跃序列 b=[.7,.3]; %b向量 a=[1,.8,.5]; %a向量 y=filter(b,a,x); %实现实验1 ，图（1） plot(t,x,g.,t,y,k-); [h,t]=impz(b,a,40); %求出单位抽样响应。图（2） stem(t,h,.);grid on; t=0:20; x=[1,zeros(1,20)]; %产生单位采样序列 b=[.7,.3]; %形成b向量 a=[1,.8,.5]; %形成a向量 y=filter(b,a,x); %filter函数 图（3） plot(t,x,g.,t,y,k-); t=0:20; x=[1,zeros(1,20)]; b=[.7,.3]; a=[1,.8,.5]; [h,t]=impz(b,a,40); %impz函数 图（4） stem(t,h,.);grid on; 实验结果 （1） （2） （3） （4） 分析结果 输入为离散是，输出为连续。相反，输入为连续，输出为离散。 对于单位阶跃和单位抽样输入来说，输出没有变化 实验二 使用FFT作频谱分析 （1） 使用FFT 对MATLAB 中的数据noissin.dat 进行谱分析。 （2） 使用FFT 对语音数据noisyspeech.wav 进行谱分析。 实验原理 （1）离散傅里叶变换(DFT)公式为： X(k)=∑x(n)W^nk; x(n)=∑X(k)W^-nk; 其中w=e^(2∏nk/N),N为离散序列的长度。 （2）快速傅里叶变换(FFT)是利用w因子的取值特点，减少DFT的复数乘法的次数。其中一种是时间抽取基2算法，它将时间按奇偶逐级分开，直到两点的DFT。MATLAB提供了fft函数可用于计算FFT，器调用形式为;X=fft（x）或X=fft（x，N），N为2的整数次幂，若x的长度小于N，则补零，若超过则舍去N后的数据。 （3）函数形式[y,fs,bits]=wavread(Blip,[N1 N2]);用于读取语音，采样值放在向量y中，fs表示采样频率(Hz)，bits表示采样位数。[N1 N2]表示读取从N1点到N2点的值（若只有一个N的点则表示读取前N点的采样值）。 实验内容 xx=wavread(noisyspeech.wav); fs=100;N=128; x=xx(1:N); n=1:N; X=abs(fft(x,N)); subplot(221);plot(n,x); xlabel(n);ylabel(x(n)); grid on; subplot(222);plot(n,X); xlabel(k);ylabel(|X(k)|); grid on; load noissin; s=noissin(1:20); S=fft(s); subplot(223);plot(abs(s)); xlabel(n);ylabel(|s(n)|); grid on; subplot(224);plot(abs(S)); xlabel(k);ylabel(|S(k)|); grid on; 实验结果及分析 实验三 使用双线性Z变换设计IIR滤波器 使用双线性Z 变换法设计一个低通数字IIR 滤波器，给定的数字滤波器的技术指标为fp=100Hz,fs=300Hz,ap=3dB,as=20dB,抽样频率Fs=1000Hz 实验原理 1）设计滤波器就是要设计一个系统是其能让一定频率