matlab实习报告.doc

实验报告 题目一 : A sinusoid of frequency ω0=0.1πand duration of 300 samples, that is ,0≤n300, is input to a (causal) filter with transfer function ,where a=0.97。Adjust the scale factor b such that the filter’s 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 filter’s magnitude response. 实验目的 给定一个系统的传输函数，输入一个正弦序列，经过系统的传输函数后求其差分方程，同时考虑零点，极点对幅频特性曲线的影响。 实验原理 一般时域离散系统或网络可以用差分方程，单位脉冲响应以及系统函数进行描述。 实验步骤简述 求系数b的值，利用其模值求出b的值。求出其差分方程y(n)，利用滤波函数将其传输函数进行逆变换。求出传输函数的幅频特性曲线，画出输入正弦函数的图形，画出差分方程的图形 实验程序框图 实验源程序 当a=0.97时 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); x1=[xn zeros(1,150)]; a=[1 -a1]; yn=filter(b,a,x1); figure(name,实验一,color,r); subplot(211); plot(xn); title( x(n)); xlabel(n); ylabel(幅度); t1=toc subplot(212); plot(yn); title( y(n)); xlabel(n); ylabel(幅度); figure(name,实验一,color,r); freqz(b,a); pause close all; 6.实验结果 当a=0.35时 close all; clear all; clc tic a1=0.35; w=0.1*pi; b=abs(1-a1*exp(-j*w)); N=300; n=1:N; xn=sin(w*n); x1=[xn zeros(1,150)]; a=[1 -a1]; yn=filter(b,a,x1); figure(name,实验一,color,r); subplot(211); plot(xn); title( x(n)); xlabel(n); ylabel(幅度); t1=toc subplot(212); plot(yn); title( y(n)); xlabel(n); ylabel(幅度); figure(name,实验一,color,r); freqz(b,a); pause close all; 实验结果 ] 7. 实验总结: 当a值变小后，它的相频响应的曲线会变的平缓 题目 （二） 一：实验题目： It is desired to extract signal x(n)=s(n)+v(n)=s+v(n),where v(n) is zero-mean white Caussian noise of variance б.(1)First-order IIR smoother; To this end,the following IIR lowpass filter is used:H(Z)=b/(1-az-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+h1z-1+h2z-2+h3z-3.(提示：当四个系数相等时，其噪声抑制效果最佳。) 二：实验目的： 让确定信号分别通过有限长脉冲响应（FIR）数字滤波器和无限长脉冲响应（IIR）数字滤波器，通过分析输出图像，掌握滤波器各参数对性能的影响。 三：实验原理: 确定信号为一个常数与高斯白噪声相加。通过滤波器，对于