信号与系统实验报告3.docxVIP

信号与系统实验报告实验三信号的卷积小组成员：钟鲁超迈克、实验目的1、理解卷积的物理意义。2、掌握运用计算机进行卷积运算的原理和方法。3、熟悉卷积运算函数conv的应用。二、预习内容1、卷积的定义及物理意义；2、卷积计算的图解法。3、卷积的应用。三、实验原理1. 卷积的定义2. 卷积的计算3. 卷积的应用 3.1 求解系统响应 3.2 相关性分析4. 涉及的Matlab函数 4.1 conv函数四、实验内容给定如下因果线性时不变系统：y[n]+0.71y[n-1]-0.46y[n-2]-0.62y[n-3]=0.9x[n]-0.45x[n]+0.35[n-2]+0.002x[n-3](1)不用impz函数，使用filter命令，求出以上系统的单位冲激响应h[n]的前20个样本；num=[0.9 -0.45 0.35 0.002];den=[1 0.71 -0.46 -0.62];n=0:19;x=(n==0);h=filter(num,den,x);stem(n,h)xlabel(n);ylabel(h[n]);title(Filter函数求得的单位冲激响应);（2）得到h[n]后，给定x[n]，计算卷积输出y[n]；并用滤波器h[n]对输入x[n]滤波，求得y1[n]；x=[1 -2 3 -4 3 -2 1];y=conv(h,x);n=0:25;subplot(2,1,1);stem(n,y);xlabel(时间序号);ylabel(振幅);title(用卷积得到的输出);grid;x1=[x zeros(1,19)];y1=filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel(时间序号);ylabel(振幅);title(用滤波得到的输出);grid;（3）y[n]和y1[n]有差别吗？为什么要对x[n]进行补零得到的x[n]来作为输入来产生y1[n]?（4）思考：设计实验，证明下列结论1单位冲激信号卷积：δ(t)*f(t)=f(t), δ(t-t0)*f(t)=f(t-t0)n=0:40;m=(n==0);f=sin(n);f1=conv(m,f);subplot(3,1,1);f1=f1(1:41);stem(n,f1);title();grid;subplot(3,1,2);stem(n,f);title();grid;subplot(3,1,3);stem(n,f-f1);title();grid;n=0:50;t=10;m=(n-t==0);f=sin(n);f1=conv(m,f);subplot(3,1,1);f1=f1(1:51);stem(n,f1);title();grid;subplot(3,1,2);f=[zeros(1,t) f]f=f(1:51);stem(n,f);title();grid;subplot(3,1,3);stem(n,f-f1);title();grid;2卷积交换律：f(t)=f1(t)*f2(t)=f2(t)*f1(t)n=0:40;f1=sin(n);f2=cos(n);y1=conv(f1,f2);y1=y1(1:41);y2=conv(f2,f1);y2=y2(1:41);subplot(3,1,1);stem(n,y1);title(f1*f2);grid;subplot(3,1,2);stem(n,y2);grid;subplot(3,1,3);stem(n,y1-y2);grid;3卷积分配律：f1(t)*[f2(t)+f3(t)]=f1(t)*f2(t)+f1(t)*f3(t)n=0:40;f1=sin(n);f2=cos(n);y1=conv(f1,f2);y1=y1(1:41);y2=conv(f2,f1);y2=y2(1:41);subplot(3,1,1);stem(n,y1);title(f1*f2);grid;subplot(3,1,2);stem(n,y2);grid;subplot(3,1,3);stem(n,y1-y2);grid;