离散系统的差分方程,冲击响应和卷积分析.docVIP

实 验 报 告 专业： 班级： 指导老师： 姓名： 学号： 实验室： 实验名称：离散系统的差分方程，冲击响应和卷积分析 时间： 一，实验目的及要求 （一）实验目的 加深对离散系统的差分方程，冲击响应和卷积分析的理解 二，实验要求 给出理论计算结果和程序计算结果并讨论 （二）实验原理 在matlab中，可以用函数y=Filter(p,d,x)求解差分方程，也可以用函数y=Conv(x,h)计算卷积。 三，实验内容 （1）卷积的计算 n=1:50; hb=zeros(1,50); hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1; close all subplot(3,1,1);stem(hb);title(25FF.系统hb[n]); m=1:50; A=444.128; a=50*sqrt(2.0)*pi; T=0.001; w0=50*sqrt(2.0)*pi; x=A*exp(-a*m*T).*sin(w0*m*T); subplot(3,1,2);stem(x);title(25FF.输入信号x[n]); y=conv(x,hb); subplot(3,1,3);stem(y);title(25FF.输出信号y[n]); （2）卷积定律验证 k=-25:25; X=x*(exp(-j*pi/12.5)).^(n*k); magX=abs(X); subplot(3,2,1);stem(magX);title(25FF.输入信号的幅度谱); angX=angle(X); subplot(3,2,2);stem(angX);title(25FF.输入信号的相位谱); Hb=hb*(exp(-j*pi/12.5)).^(n*k); magHb=abs(Hb); subplot(3,2,3);stem(magHb);title(25FF.系统响应的幅度谱); angHb=angle(Hb); subplot(3,2,4);stem(angHb);title(25FF.系统响应的相位谱); n=1:99; k=1:99; Y=y*(exp(-j*pi/12.5)).^(n*k); magY=abs(Y); subplot(3,2,5);stem(magY);title(25FF.输出信号的幅度谱); angY=angle(Y); subplot(3,2,6);stem(angY);title(25FF.输出信号的相位谱); XHb=X.*Hb; Subplot(2,1,1);stem(abs(XHb));title(25FF.x(n)的幅度谱与hb(n)幅度谱相乘); Subplot(2,1,2);stem(abs(Y));title(25FF.y(n)的幅度谱);axis([0,60,0,8000]) 编制程序求解下列两个系统的单位冲击响应和阶跃响应，并绘出其图形。 y[n]+0.75y[n-1]+0.125y[n-2]=x[n]-x[n-1] y[n]=0.25{x[n-1]+x[n-2]+x[n-3]+x[n-4]} （1）程序如下： N=21; a=[1 -1]; b=[1 0.75 0.125]; x1=[1 zeros(1,N-1)]; x2=ones(1,N); k=0:1:N-1; y1=filter(a,b,x1); y2=filter(a,b,x2); subplot(2，1，1);stem(k,y1); xlabel(n);ylabel(幅度);title(25FF.单位冲击响应); subplot(2，1，2);stem(k,y2); xlabel(n);ylabel(幅度);title(25FF.阶跃响应); （2）程序如下： N=21; a=0.25*[1 1 1 1]; b=[1]; x1=[1 zeros(1,N-1)]; x2=ones(1,N); k=0:1:N-1; y1=filter(a,b,x1); y2=filter(a,b,x2); subplot(2，1，1);stem(k,y1);title(25FF.单位冲击响应); xlabel(n);ylabel(幅度); subplot(2，1，2);stem(k,y2);title(25FF.阶跃响应); xlabel(n);ylabel(幅度); （1） （2） 结 论 教 师 批 阅 年 月 日