信号系统及系统响应Matlab1解读

实验一信号系统及系统响应一、实验目的1、了解常用的时域离散信号及其特点；2、掌握Matlab产生常用时域离散信号的方法；3、掌握时域离散信号简单的基本运算方法；4、掌握求解离散时间系统冲激响应和阶跃响应的方法；5、进一步理解卷积定理，掌握应用线性卷积求解离散时间系统响应的基本方法；6、掌握离散系统的响应特点。二、实验内容1、自己设定参数，分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、实指数序列、随机序列；2、自己设定参数，分别表示并绘制信号移位、信号相加、信号相乘、信号翻转、信号和、信号积、信号能量。3、已知信号其他040614452)(nnnnx（1）描绘)(nx序列的波形；（2）用延迟的单位脉冲序列及其加权和表示)(nx序列；（3）描绘一下序列的波形)2()(),2(2)(),2(2)(321nxnxnxnxnxnx4、请分别用impz和dstep函数求解下面离散时间系统的冲激响应和阶跃响应。（1）系统的差分方程为：)(866.0)2(64.0)1(8.0)(nxnynyny（2）系统的系统函数为：21115.01)(zzzzH5、运行例题2.3，理解卷积过程和程序中每一句的意义。6、利用第1题求得的系统冲激响应求解系统在激励)3()(nunx下的响应。三、实现步骤1、自己设定参数，分别表示并绘制单位抽样序列、单位阶跃序列、正弦序列、实指数序列、随机序列。输出图形如图1所示。x=zeros(1,10);x(2)=1;subplot(3,2,1);stem(x,'filled')axis([0,10,-0.2,1]);title('单位抽样序列');N=10;u=ones(1,N);subplot(3,2,2);stem(u,'filled')axis([-10,10,0,1]);title('单位阶跃序列');x=-20:1:20;y=sin(0.2*pi.*x+0.5*pi);subplot(3,2,3);stem(x,y,'filled');axis([-20,20,-2,2]);title('正弦序列');n=0:10;a1=1/2;y1=a1.^n;subplot(3,2,4);stem(n,y1,'filled');axis([0,10,0,1]);title('实指数序列，a=1/2');n=0:10;a2=2;y2=a2.^n;subplot(3,2,5);stem(n,y2,'filled');title('实指数序列,a=2');y=rand(1,20);subplot(3,2,6);stem(y,'filled');title('随机序列');图12、自己设定参数，分别表示并绘制信号移位、信号相加、信号相乘、信号翻转、信号和、信号积、信号能量。信号的移位：n=-3:10;k0=3;k1=-3;%实现信号的移位x=cos(2*pi*n/10);x1=cos(2*pi*(n-k0)/10);x2=cos(2*pi*(n-k1)/10);subplot(3,1,1),stem(n,x,'filled');ylabel('x(n)');subplot(3,1,2),stem(n,x1,'filled');ylabel('x(n-2)');subplot(3,1,3),stem(n,x2,'filled');ylabel('x(n+2)');信号的相加相乘：n=-3:20;x1=cos(2*pi*n/10);subplot(2,2,1);stem(n,x1,'filled');title('x(1)');axis([-4,20,-2,2]);x2=cos(2*pi*n/10);subplot(2,2,2);stem(n,x2,'filled');title('x(2)');axis([-4,20,-2,2]);y=x1+x2;subplot(2,2,3);stem(n,y,'filled');title('信号相加');axis([-4,20,-2,2]);y=x1.*x2;subplot(2,2,4);图2stem(n,y,'filled');title('信号相乘');axis([-4,20,-2,2]);图3信号的翻转：n=-5:5;x=exp(-0.4*n);x1=fliplr(x);n1=-fliplr(n);subplot(2,1,1),stem(n,x,'filled');title('x(n)');subplot(2,1,2),stem(n1,x1,'filled');title('x(-n)');%见下图4信号和、信号积、信号能量：x=[1,2,3,4,5,6,7,8,9];%和y1=sum(x)%积y2=prod(x)%能量E1=sum(x.*conj(x))得到：y1=45y2=362880E1=285图43、已知信号其他040614452)(nnnnx（1）描绘)(nx序列的波形；（2）用延迟的单位脉冲序列及其加权和表示)(nx序列；（3）描绘一下序列的波形)2()(),2(2)(),2(2)(321nxnxnxnxnxnxfunctionf=u(t)f=(t=0);subplot(2,1,1)n=-10:10;y1=(2*n+5).*(u(n+4)-u(n))+6.*(u(n)-u(n-5));stem(n,y1,'filled')axis([-10,10,-3,6]);title('序列波形');t=-10:10;subplot(2,1,2)y=(-3)*(u(t+4)-u(t+3))+(-1)*(u(t+3)-u(t+2))+(u(t+2)-u(t+1))+3*(u(t+1)-u(t))+6*(u(t)-u(t-1))+6*(u(t-1)-u(t-2))+6*(u(t-2)-u(t-3))+6*(u(t-3)-u(t-4))+6*(u(t-4)-u(t-5));stem(t,y,'filled')axis([-10,10,-3,6]);title('用单位脉冲序列及其加权和表示序列波形');图5subplot(2,2,1)t=-10:10;y=(-3)*(u(t+4)-u(t+3))+(-1)*(u(t+3)-u(t+2))+(u(t+2)-u(t+1))+3*(u(t+1)-u(t))+6*(u(t)-u(t-1))+6*(u(t-1)-u(t-2))+6*(u(t-2)-u(t-3))+6*(u(t-3)-u(t-4))+6*(u(t-4)-u(t-5));stem(t,y,'filled')axis([-10,10,-6,12]);title('x(n)');subplot(2,2,2)t=-10:10;y=(-3)*(u(t+4)-u(t+3))+(-1)*(u(t+3)-u(t+2))+(u(t+2)-u(t+1))+3*(u(t+1)-u(t))+6*(u(t)-u(t-1))+6*(u(t-1)-u(t-2))+6*(u(t-2)-u(t-3))+6*(u(t-3)-u(t-4))+6*(u(t-4)-u(t-5));stem(t+2,2*y,'filled')axis([-10,10,-6,12]);title('2x(n-2)');subplot(2,2,3)t=-10:10;y=(-3)*(u(t+4)-u(t+3))+(-1)*(u(t+3)-u(t+2))+(u(t+2)-u(t+1))+3*(u(t+1)-u(t))+6*(u(t)-u(t-1))+6*(u(t-1)-u(t-2))+6*(u(t-2)-u(t-3))+6*(u(t-3)-u(t-4))+6*(u(t-4)-u(t-5));stem(t-2,2*y,'filled')axis([-10,10,-6,12]);title('2x(n+2)');subplot(2,2,4)t=-10:10;y=(-3)*(u(t+4)-u(t+3))+(-1)*(u(t+3)-u(t+2))+(u(t+2)-u(t+1))+3*(u(t+1)-u(t))+6*(u(t)-u(t-1))+6*(u(t-1)-u(t-2))+6*(u(t-2)-u(t-3))+6*(u(t-3)-u(t-4))+6*(u(t-4)-u(t-5));stem(2-t,y,'filled')axis([-10,10,-6,12]);title('x(2-n)');图64、请分别用impz和dstep函数求解下面离散时间系统的冲激响应和阶跃响应。（1）系统的差分方程为：)(866.0)2(64.0)1(8.0)(nxnynyny%y(n)=0.8y(n-1)-0.64y(n-2)+0.866x(n)%y(n)-0.8y(n-1)+0.64y(n-2)=0.866x(n)%a0=1,a1=-0.8,a2=0.64%b0=0.866,b1=0,b2=0a=[1,-0.8,0.64];b=[0.866,0,0];n=20;hn=impz(b,a,n);%冲激响应gn=dstep(b,a,n);%阶跃响应subplot(2,1,1),stem(hn,'filled');%显示冲激响应曲线title('系统的单位冲激响应');ylabel('h(n)');xlabel('n');axis([0,n,1.1*min(hn),1.1*max(hn)]);subplot(2,1,2),stem(gn,'filled');%显示阶跃响应曲线title('系统的单位阶跃响应');ylabel('g(n)');xlabel('n');axis([0,n,1.1*min(gn),1.1*max(gn)]);（2）系统的系统函数为：21115.01)(zzzzH%a0=1,a1=-1,a2=1%b0=1,b1=-0.5,b2=0a=[1,-1,1];b=[1,-0.5,0];n=20;hn=impz(b,a,n);%冲激响应gn=dstep(b,a,n);%阶跃响应subplot(2,1,1),stem(hn,'filled');%显示冲激响应曲线title('系统的单位冲激响应');ylabel('h(n)');xlabel('n');axis([0,n,1.1*min(hn),1.1*max(hn)]);subplot(2,1,2),stem(gn,'filled');%显示阶跃响应曲线title('系统的单位阶跃响应');ylabel('g(n)');xlabel('n');axis([0,n,1.1*min(gn),1.1*max(gn)]);图7图85、运行例题2.3，理解卷积过程和程序中每一句的意义。clf;nf1=0:20;%f1的时间向量f1=0.8.^nf1;lf1=length(f1);nf2=0:10;%f2的时间向量lf2=length(nf2);%取f2的时间向量的长度f2=ones(1,lf2);lmax=max(lf2,lf1);iflf2lf1nf2=0;nf1=lf2-lf1;elseiflf2lf1nf1=0;nf2=lf1-lf2;elsenf2=0;lf1=0;endlt=lmax;u=[zeros(1,lt),f2,zeros(1,nf2),zeros(1,lt)];图9t1=(-lt+1:2*lt);f1=[zeros(1,2*lt),f1,zeros(1,nf1)];hf1=fliplr(f1);N=length(hf1);y=zeros(1,3*lt);fork=0:2*ltp=[zeros(1,k),hf1(1:N-k)];y1=u.*pyk=sum(y1);y(k+lt+1)=yk;subplot(4,1,1);stem(t1,u);subplot(4,1,2);stem(t1,p);subplot(4,1,3);stem(t1,y1);subplot(4