信号与系统实验报告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 4]f6=[1 1]k7=[3]f7=[0]k8=-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 1]b=[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 1]b=[0 0 2];impz(b,a);title(离散系统单位样值响应)/*******自动生成*******/a = 1.0000 -1.5000 1.00008-1a=[1 -1 -2]b=[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=fil