数字信号处理实验-程序全集

程序代码、绝对正确无误

实验二：10小题图形一方法一 x1=[0,1,2,3,4,3,2,1,0];n1=-2:6; [y1,n2]=sigfold(x1,n1); [y2,n3]=sigshift(x1,n2,-2); y3=2*y2; subplot(2,2,1);stem(n1,x1);title('序列x1') xlabel('n');ylabel('x1(n)'); subplot(2,2,2);stem(n2,y1);title('序列折叠') xlabel('n');ylabel('y1(n)'); subplot(2,2,3);stem(n3,y2);title('序列移位') xlabel('n');ylabel('y2(n)'); subplot(2,2,4);stem(n3,y3);title('序列y3') xlabel('n');ylabel('y3(n)'); 方法二 x1=[0,1,2,3,4,3,2,1,0];n1=-2:6; [y1,n2]=sigshift(x1,n1,2); [y2,n3]=sigfold(x1,n2); y3=2*y2; subplot(2,2,1);stem(n1,x1);title('序列x1') xlabel('n');ylabel('x1(n)'); subplot(2,2,2);stem(n2,y1);title('序列移位') xlabel('n');ylabel('y1(n)'); subplot(2,2,3);stem(n3,y2);title('序列折叠') xlabel('n');ylabel('y2(n)'); subplot(2,2,4);stem(n3,y3);title('序列y3') xlabel('n');ylabel('y3(n)'); 图形二方法 x1=[0,1,2,3,4,3,2,1,0]; n1=-2:6; n2=n1/3; n3=fix(-2/3):fix(6/3); n3 m=find(n1/3==fix(n1/3)); m y4=x1(m) subplot(3,1,1);stem(n1,x1);title('序列x1') xlabel('n');ylabel('x1(n)'); subplot(3,1,2);stem(n2,x1);title('变化后') xlabel('n');ylabel('x1(n)'); subplot(3,1,3);stem(n3,y4);title('序列y4') xlabel('n');ylabel('y4(n)'); 共轭反对称，x2为复共轭x*(-n)，x3是公式 n=[0:1:20]; alpha=-0.1+0.5j; x=exp(alpha*n); [x1,n1]=sigfold(x,n); x2=conj(x1); x3=0.5*[x+x2]; subplot(2,2,1);stem(n,real(x3)); title('实部');xlabel('n') subplot(2,2,3);stem(n,imag(x3)); title('虚部');xlabel('n') subplot(2,2,2);stem(n,abs(x3)); title('振幅');xlabel('n') subplot(2,2,4);stem(n,(180/pi)*angle(x3)); title('相位');xlabel('n') 共轭对称，x2为复共轭x*(-n)，x4是公式 n=[0:1:20]; alpha=-0.1+0.5j; x=exp(alpha*n); [x1,n1]=sigfold(x,n); x2=conj(x1); x4=0.5*[x-x2]; subplot(2,2,1);stem(n,real(x4)); title('实部');xlabel('n') subplot(2,2,3);stem(n,imag(x4)); title('虚部');xlabel('n') subplot(2,2,2);stem(n,abs(x4)); title('振幅');xlabel('n') subplot(2,2,4);stem(n,(180/pi)*angle(x4)); title('相位');xlabel('n') 验证共轭对称与共轭反对称 x5=x3+x4; subplot(2,2,1);stem(n,real(x5)); title('实部');xlabel('n') subplot(2,2,3);stem(n,imag(x5)); title('虚部');xlabel('n') subplot(2,2,2);stem(n,abs(x5)); title('振幅');xlabel('n') subplot(2,2,4);stem(n,(180/pi)*angle(x5)); title('相位');xlabel('n') 将x5与x

程序代码、绝对正确无误

比较，除横坐标相反，其余均相同 n=[0:1:20]; alpha=-0.1+0.5j; x=exp(alpha*n); subplot(2,2,1);stem(n,real(x)); title('实部');xlabel('n') subplot(2,2,3);stem(n,imag(x)); title('虚部');xlabel('n') subplot(2,2,2);stem(n,abs(x)); title('振幅');xlabel('n') subplot(2,2,4);stem(n,(180/pi)*angle(x)); title('相位');xlabel('n')