%实验一第一题A=[1234];B=[3456];C=A+B;D=A-B;E=A.*B;F=A./B;G=A.^B;%a.^b,是a中每个元素b次方，形成的新矩阵subplot(2,4,1)%两行四列第一个stem(A);title('A');%标题subplot(2,4,2)stem(B);title('B');subplot(2,4,3)stem(C);title('C');subplot(2,4,4)stem(D);title('D');subplot(2,4,5)stem(E);title('E');subplot(2,4,6)stem(F);title('F');subplot(2,4,7)stem(G);title('G');%实验一第二题n=0:1:15;x1=0.8.^n;%求n次方x2=exp((0.2+3j)*n);%虚数jx3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);N=0:1:63;x4=3*cos(0.125*pi*rem(N,16)+0.2*pi)+2*sin(0.25*pi*rem(N,16)+0.1*pi);%以16为周期进行延拓N2=0:1:39;x5=3*cos(0.125*pi*rem(N2,10)+0.2*pi)+2*sin(0.25*pi*rem(N2,10)+0.1*pi);%以10为周期进行延拓figure(1);stem(n,x3);title('x(n)');%标题figure(2);subplot(2,1,1);stem(N,x4);title('x(n)以16为周期进行周期延拓');subplot(2,1,2);stem(N2,x5);title('x(n)以10为周期进行周期延拓');第三题n=0:1:3;x=[1-135];%circshift矩阵循环平移,[a,b],a是列向移动,b是行向移动x1=2*circshift(x,[0-2]);x2=x1-circshift(x,[01]);X1=x2+2*x;y1=circshift(x,[01]);y2=circshift(x,[02]);y3=circshift(x,[03]);y4=circshift(x,[04]);y5=circshift(x,[05]);X2=1*y1+2*y2+3*y3+4*y4+5*y5;subplot(2,1,1);stem(n,X1);title('x1(n)=2*x(n+2)-x(n-1)-2*x(n)');xlabel('n');ylabel('x1(n)');subplot(2,1,2);stem(n,X2);title('x2(n)=1*x(n-1)+2*x(n-2)+3*x(n-3)+4*x(n-4)+5*x(n-5)');xlabel('n');ylabel('x2(n)');第四题functionstepshift(n0,n1,n2)n1=input('请输入起点:');n2=input('请输入终点:');n=n1:n2;n0=input('请输入原点:');y=[(n-n0)>=0];%当逻辑值(n-n0)>=0为真的时候，y=1stem(n,y);xlabel('序列值n');ylabel('函数值');title('函数stepshift(n0,n1,n2)');第五题B=[1sqrt(2)1];A=[1-0.670.9];[H,w]=freqz(B,A);%离散系统频率响应的函数fregz()fp=20*log(abs(H));%幅频响应xp=angle(H);%相频响应subplot(2,1,1);plot(w,fp);title('幅频响应');xlabel('时间/t');ylabel('幅频响应/dB');subplot(2,1,2);plot(w,fp);title('相频响应');xlabel('时间/t');ylabel('相频响应');第六题h=[8-2-123];x=[23-1-3];n=0:7;y=conv(h,x);%h和x的卷积，结果放在y里stem(n,y);title('[8-2-123]与[23-1-3]的离散卷积');第七题n=0:1:30;u=[n>=0];A=[10.1-0.06];B=[1-2];h=impz(B,A,n);subplot(2,1,1);stem(n,h);title('单位脉冲响应');subplot(2,1,2);y=zplane(B,A);title('零极点图');实验二第一题f1=50;f2=500;f3=1000;f=8000;n=0:299;y=sin(f1/f*2*pi*n)+sin(f2/f*2*pi*n)+sin(f3/f*2*pi*n);stem(n,y);title('抽样信号');xlabel('n');ylabel('幅值');[h,w]=freqz(y);figure(2);plot(w(1:299),abs(h(1:299)));title('幅频响应');xlabel('角度');ylabel('幅度');第二题f1=50;f2=500;f3=1000;f=800;n=0:99;y=sin(f1/f*2*pi*n)+sin(f2/f*2*pi*n)+sin(f3/f*2*pi*n);stem(n,y);title('抽样信号');xlabel('n');ylabel('幅值');[h,w]=freqz(y);figure(2);plot(w(1:499),abs(h(1:499)));title('幅频响应');xlabel('角度');ylabel('幅度');第三题一f=1;fs=16;N=52;n=0:N-1;x=cos(2*pi*(f/fs)*n);subplot(1,2,1);stem(n,x(1:N));title('原信号');xlabel('n');ylabel('幅度');[h,w]=freqz(x);subplot(1,2,2);plot(w(1:499),abs(h(1:499)));title('原信号幅频响应');xlabel('角度');ylabel('幅度');M=4;y=decimate(x,M);%抽取m=0:N/M-1;figure(2);subplot(1,2,1);stem(m,y(1:N/M));title('抽取信号');xlabel('n');ylabel('幅度');[h1,w1]=freqz(y);subplot(1,2,2);plot(w1(1:499),abs(h1(1:499)));title('抽取信号幅频响应');xlabel('角度');ylabel('幅度');第三题二f=1;fs=16;N=52;n=0:N-1;x=cos(2*pi*(f/fs)*n);subplot(1,2,1);stem(n,x(1:N));title('原信号');xlabel('n');ylabel('幅度');[h,w]=freqz(x);subplot(1,2,2);plot(w(1:499),abs(h(1:499)));title('原信号幅频响应');xlabel('角度');yla...
