资源描述:
《西工大信号与系统_实验2.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、......西北工业大学《信号与系统》实验报告西北工业大学2016年9月一、实验目的专业word可编辑.......一、实验要求二、实验设备(环境)三、实验内容与步骤四、实验结果2.1MATLAB函数conva:0,其他;1,n=0,10;2,n=1,9y[n]=3,n=2,84,n=3,75,n=4,66,n=5b:代码如下:n=0:10;xn=[111111]y=conv(xn,xn)stem(n,y);运行结果如图,与图2.1一致专业word可编辑.......c:代码如下n=0:5;xn=[111111]ny=0:10hn=[012345]y=conv(x
2、n,hn)stem(ny,y);运行结果如下,与图2.2一致专业word可编辑.......d:因为h不同,经过了时移且序列长度增加了,因此卷积后的结果也不一样,由于卷积后序列长度等于被卷积的两序列长度之和减去1,比在3中导出的信号要长,且每个元素值不一样e:h=[0000012345];x=[111111];y=conv(x,h)len=length(y);ny=[0:10];%计算向量y的序号stem(ny,y);gridon;运行结果:y=[0000013610151514295]专业word可编辑.......2.4离散时间LTI系统的性质a:代码如下x1
3、=[1111100000]h1=[1-1310]h2=[254-10]fori=1:length(x1),nx1(i)=i-1;endfori=1:length(h1),nx2(i)=i-1;endsubplot(311)stem(nx1,x1);title('x1')subplot(312)stem(nx2,h1);title('h1')专业word可编辑.......subplot(313)stem(nx2,h2);title('h2')运行截图如下:b:b:专业word可编辑.......由上图结果可得conv的输出与卷积的顺序无关C:x1=[1111100
4、000];h1=[1-1310];h2=[254-10];y1=conv(x1,h1);y2=conv(x1,h2);y=y1+y2%先分别求卷积,然后求和yy=conv(x1,h1+h2)%求冲激响应求和,再卷积运行截图:可见专业word可编辑.......先分别求卷积,然后求和得出的结果,跟先求冲激响应求和在卷积得出的结果相同,即满足分配律D:x1=[1111100000];h1=[1-1310];h2=[254-10];y1=conv(x1,h1);y2=conv(h1,h2);y=conv(y1,h2)%先x1与h1卷积,所得结果再与h2卷积yy=conv
5、(x1,y2)%先h1与h2卷积,再x1与所得结果卷积运行结果:2.5线性和时不变性A:系统一的结果图专业word可编辑.......系统二的结果图系统三的结果图专业word可编辑.......代码如下:x1=[100000];x2=[010000];x3=[120000];w1=w(x1)w2=w(x2)w3=w(x3)forn=1:length(x1),ny(n)=n-1;endsubplot(221);stem(ny,w1);gridon;legend('w1');subplot(222);stem(ny,w2);gridon;legend('w2');su
6、bplot(223);stem(ny,w3);gridon;legend('w3');专业word可编辑.......subplot(224);stem(ny,w1+2*w2);gridon;legend('w1+2*w2');函数定义如下:function[y]=w(x)len=length(x);fori=1:lenifi==1,y(i)=x(i);elseifi==2,y(i)=x(i)+x(i-1);elsey(i)=x(i)+x(i-1)+x(i-2);endendendfunction[y1]=y(x)len=length(x);fori=1:leny
7、1(i)=cos(x(i));endfunction[y1]=z(x)len=length(x);fori=1:leny1(i)=i*x(i);end专业word可编辑.......2.6:非因果有限冲激响应滤波器A:满足2.16式的LTI系统的单位冲激响应为b[n];若系统非因果,则N1必须小于0B:N6=N2+N4,N5=N1+N3C:x=[1524-22];fori=-3:3h(i+4)=1-abs(i)/3end;nx=[0:5];nh=[-3:3];subplot(211);stem(nx,x);gridon;legend('x');subplot(21
8、2);st