信号与系统实验报告 7.docx

信号与系统实验实验七报告一、实验内容1-1k=-20:20;f=cos(k*pi/8);stem(k,f,filled);1-2k=-20:20;f=sin(k*pi/8);stem(k,f,filled);5subplot(2,2,1);k1=1:4;f1=ones(1,length(k1);k2=5:10;f2=zeros(1,length(k2);k3=-1:0;f3=zeros(1,length(k3);stem(k1,f1,filled);hold onstem(k2,f2,filled);hold onstem(k3,f3,filled);axis(-1,10,-0.5,1.5);subplot(2,2,2);k4=1:2;f4=ones(1,length(k4);k5=3:7;f5=zeros(1,length(k5);stem(k4,f4,filled);hold onstem(k5,f5,filled);hold onstem(k3,f3,filled);axis(-1,7,-0.5,1.5);subplot(2,2,3);k6=2 4f6=1 1k7=3f7=0k8=-1:1;f8=zeros(1,length(k8);k9=5:9;f9=zeros(1,length(k9);stem(k6,f6,filled);hold onstem(k7,f7,filled);hold onstem(k8,f8,filled);hold onstem(k9,f9,filled);axis(-1,9,-0.5,1.5);subplot(2,2,4);k10=2:3;f10=ones(1,length(k10);k11=-1:1;f11=zeros(1,length(k11);k12=4:10;f12=zeros(1,length(k12);stem(k10,f10,filled);hold onstem(k11,f11,filled);hold onstem(k12,f12,filled);axis(-1,10,-0.5,1.5)/*自动生成*/k6 = 2 4f6 = 1 1k7 = 3f7 = 07function y=recur(a,b,n,f,f0,y0);N=length(a);y=y0 zeros(1,length(n);M=length(b)-1;f=f0 f;a1=a(N:-1:1);b1=b(M+1:-1:1);for i=N+1:N+length(n),y(i)=-a1*y(i-N:i-1)+b1*f(i-N:i-N+M);endy=y(N+1:N+length(n);调用的函数（m文件）function y=recur(a,b,n,f,f0,y0)a=-1.5 1;b=0 0 2;y0=2 1;f0=0 0;k=0:30;subplot(2,1,1)f=ones(1,length(k);y=recur(a,b,k,f,f0,y0);stem(k,y,o),xlabel(k),ylabel(y(k)subplot(2,1,2)n=length(k);f=zeros(1,n);f(1,1)=1;y=recur(a,b,k,f,f0,y0);stem(k,y,o),xlabel(k),ylabel(y(k)a=1 -1.5 1b=0 0 2;t=0:30;x=ones(1,length(t);y=filter(b,a,x);stem(t,y);title(离散系统单位阶跃响应)xlabel(k);ylabel(g(k)/*自动生成*/a = 1.0000 -1.5000 1.0000a=1 -1.5 1b=0 0 2;impz(b,a);title(离散系统单位样值响应)/*自动生成*/a = 1.0000 -1.5000 1.00008-1a=1 -1 -2b=1;subplot(1,2,1);impz(b,a);title(离散系统单位样值响应)subplot(1,2,2);t=0:15;x=ones(1,length(t);y=filter(b,a,x);stem(t,y);title(离散系统单位阶跃响应)xlabel(k);ylabel(g(k)/*自动生成*/a = 1 -1 -29-1a=1 -0.25 0.5;b=1 1;t=0:20;x=(1/2).t;y=filter(b,a,x)subplot(2,1,1)stem(t,x)title(输入序列)subplot(2,1,2)stem(t,y)title(响应序列)/*自动生成*/y = Columns 1 through 5 1.0000 1.7500 0.6875 -0.3281 -0.2383 Columns 6 through 10 0.1982 0.2156 -0.0218 -0.1015 -0.0086 Columns 11 through 15 0.0515 0.0187 -0.0204 -0.0141 0.0069 Columns 16 through 20 0.0088 -0.0012 -0.0047 -0.0006 0.0022 Column 21 0.000811-1function f,k=dconv(f1,f2,k1,k2) %The function of compute f=f1*f2 % f: 卷积和序列f(k)对应的非零样值向量 % k: 序列f(k)的对应序号向量 % f1: 序列f1(k)非零样值向量 % f2: 序列f2(k)的非零样值向量 % k1: 序列f1(k)的对应序号向量 % k2: 序列f2(k)的对应序号向量 f=conv(f1,f2)%计算序列f1 与f2 的卷积和f k0=k1(1)+k2(1);%计算序列f 非零样值的起点位置 k3=length(f1)+length(f2)-2;%计算卷积和f 的非零样值的宽度k=k0:k0+k3 %确定卷积和f非零样值的序号向量subplot(2,2,1)stem(k1,f1)%在子图1 绘序列f1(k)时域波形图title(f1(k)xlabel(k)ylabel(f1(k)subplot(2,2,2)stem(k2,f2)%在图2绘序列f2(k)时波形图title(f1(k) xlabel(k)ylabel(f2(k)subplot(2,2,3)stem(k,f);%在子图3绘序列f(k)的波形图title(f(k)f1(k)与f2(k)的卷积和f(k) xlabel(k) ylabel(f(k)h=get(gca,position);h(3)=2.5*h(3);set(gca,position,h)%将第三个子图的横坐标范围扩为原来的2.5 倍调用函数k1=0:20;f1=sin(0.2*k1);k2=0:10;f2=sin(0.5*k2);f,k=dconv(f1,f2,k1,k2);三、MATLAB 计算离散卷积的方法，指出计算有限长序列卷积的dsconv()函数有什么特点。Matlab信号处理工具箱提供了一个计算两个离散序列卷积和的函数conv()，其调用格式为 y=conv(f,h) 式中，f、h分别为待卷积的两序列的向量表示，y是卷积的结果。 如 f1=2 2 2; f2=1 4 9; y=conv(f1,f2) y = 2 10 28 26 18 对于有限长序列，我们建立一个通用函数，它可以计算并画出两个有限长序列卷积的结果和波形。能使三个波形的横坐标统一，间隔相同。卷积结果显示在横坐标的中间位置。这