数字信号处理仿真指导书.docVIP

题目一 1．实验题目 A sinusoid of frequenc y ω0=0.1π and duration of 300 samples,that is, x(n)=sin(ω0 n),0≤n300,is input to (causal) filter with transfer function H(z)=b/(1-a/z),where a=0.97.Adjust the scale factor b such that the filters gain at ω0 is unity. Determine and plot the input x(n) and the output y(n) of the filter over the interval 0≤n450,by iterating the difference equation of the filter. At the same time, plot the filters magnitude response. 2．实验目的 （1）熟悉连续信号理想采样方法，加深对时域采样定理的理解。 （2）掌握数字滤波器的计算机仿真方法。 （3）通过观察对输入信号的滤波作用，获得数字滤波器的感性知识。 3．实验内容 打印输入输出序列。 画出滤波器幅频响应。 4．实验步骤 已知a=0.97，ω0=0.1π，H(z)=b/(1-a/z),且|H(Z)|=1.从而得出b=0.3096。 写出差分方程为y(n)=0.97y(n-1)+bx(n)。 画出输入序列x(n)=sin(ω0 n)及输出序列y(n)= x(n)*h(n)的时域图形。 由差分方程画出系统函数的频率响应特性曲线。 5．实验程序 close all; clear all; clc tic a1=0.97; w=0.1*pi; b=abs(1-a1*exp(-j*w)); N=300; n=1:N; xn=sin(w*n); a=[1 -a1]; yn=filter(b,a,xn); subplot(211); plot(n,xn); title( x(n)); xlabel(n); ylabel(幅度); t1=toc subplot(212); plot(n,yn); title( y(n)); xlabel(n); ylabel(幅度); figure(2); freqz(b,a); pause close all; 6．实验结果 1．输入序列x(n)=sin(ω0 n)及输出序列y(n)= x(n)*h(n)的时域图形 2．由差分方程画出系统函数的频率响应特性曲线 题目二：滤波器的设计 1．实验题目 It is desired to extract a constant signal x(n)=s(n)+v(n),where v(n) is zero-mean while is Caussian noise of varianceδ^2.fist-order IIR smoother; To this end ,the following IIR lowpass filter is used: H(z)=b/(1-a*z^(-1)),where the parameter a is restrict to the range 0a1.This can be achieve in two ways.(2) FIR averaging filters :Consider ,for example,the third-order filter H(z)=h0+h1*z^-1+h2*z^-2+h3*z^-3.(提示：当h0,h1,h2,h3四个系数相等时，其噪声抑制效果最佳) 2．实验目的 通过实验以加强对MATLAB的了解和应用。 通过调整题目中的参数，来观察滤波器的滤波效果。并与理论结果相比较，看是否与理论结果相吻合。 3．实验分析 由于待提取的分量为直流分量，而干扰信号的频谱占据整个频带，本题将从零点、极点分布特性出发，利用所谓的零、极点累试方法进行滤波器的设计。滤波器的特性取决于零、极点的情况。极点位置主要影响系统的幅度特性峰值位置及尖锐程度，零点位置主要影响系统的幅度特性谷值位置及其凹下去的程度；且通过零、极点的反复调整，并将得到的幅度特性与期望的幅度特性惊醒比较，最终得到所期望系统的幅度特性。这种设计方法称为零、极点累试方法。需要注意的是：①极点必须位于z平面单位圆内，保证数字滤波器因果稳定；②复 数零极点必须共轭成对出现，从而保证系统函数有理式的系数是实的。 系统函数为：H(z)=b/(1-az^(-1)),其极点为z=a,其对应的频谱响应为H(e^(jw)