数字信号处理原理及应用的matlab程序

复指数的matlab程序：

n=[0:40];

x=exp((0.02+0.3j)*n);

subplot(2,2,1);

stem(n,real(x));

xlabel('n');ylabel('Re(x)');

subplot(2,2,2);

stem(n,imag(x));

xlabel('n');ylabel('Im(x)');

subplot(2,2,3);

stem(n,abs(x));

xlabel('n');ylabel('Magnitude');

subplot(2,2,4);

stem(n,(180/pi)*angle(x));

xlabel('n');ylabel('Phase');

移位的matlab程序：

m=[-1:5];

x=(1/3)*(1/3).^m;

n=m-1;

x1=x;

n1=m+1;

x2=x;

subplot(3,1,1);stem(m,x);ylabel('x(n)');xlabel('n');

subplot(3,1,2);stem(n,x1);ylabel('x(n+1)');xlabel('n');

subplot(3,1,3);stem(n1,x2);ylabel('x(n-1)');xlabel('n');

反褶的matlab程序：

n=[-1:5];

x=(1/3)*(1/3).^n;

y=fliplr(x);

m=-fliplr(n);

subplot(2,1,1);stem(n,x);ylabel('x(n)');xlabel('n');

subplot(2,1,2);stem(m,y);ylabel('x(-n)');xlabel('n');

加减乘积的matlab程序：

m1=[-2:-1];y1=2..^m1;

m2=[0:5];y2=m2+1;

m=min(min(m1),min(m2)):max(max(m1),max(m2));

z1=zeros(1,length(m));z2=z1;

z1(find((m>=min(m1))&(m<=max(m1))==1))=y1;

z2(find((m>=min(m2))&(m<=max(m2))==1))=y2;

y=z1+z2;

n=[-1:5];

x=(1/2)*(1/2).^n;

nn=min(min(m),min(m)):max(max(n),max(n));

z11=zeros(1,length(nn));z22=z11;

z11(find((nn>=min(m))&(nn<=max(m))==1))=y;

z22(find((nn>=min(n))&(nn<=max(n))==1))=x;

z=z11+z22;

zz=z11.*z22;

subplot(4,1,1);stem(n,x);title('x(n)');

subplot(4,1,2);stem(m,y);title('y(n)');

subplot(4,1,3);stem(nn,z);title('z1(n)=x(n)+y(n)'); subplot(4,1,4);stem(nn,zz);title('z2(n)=x(n)*y(n)'); xlabel('n');