MATLAB实验离散信号的基本运算.doc

MATLAB实验：离散信号的基本运算 一、实验目的 利用MATLAB进行离散时间序列的基本运算，掌握基本的MATLAB函数的编写和调试方法，预习离散序列的傅立叶变换。 二、实验内容 1、学习和掌握离散信号的基本运算。 2、综合运用上述命令，完成下列任务。 （1） 已知：，调用定义的函数，实现： a. b. 要求：分别画出a和b的图像。 函数定义： 移位： function [y,n]=sigshift(x,m,n0) n=m+n0;y=x; 反转： function[y,n]=sigfold(x,n) y=fliplr(x);n=-fliplr(n); 相加： function[y,n]=sigadd(x1,n1,x2,n2) n=min(min(n1),min(n2)):max(max(n1),max(n2)); y1=zeros(1,length(n));y2=y1; y1(find((n=min(n1))(n=max(n1))==1))=x1; y2(find((n=min(n2))(n=max(n2))==1))=x2; y=y1+y2; 相乘： function [y,n]=sigmult(x1,n1,x2,n2) n=min(min(n1),min(n2)):max(max(n1),max(n2)); y1=zeros(1,length(n));y2=y1; y1(find((n=min(n1))(n=max(n1))==1))=x1; y2(find((n=min(n2))(n=max(n2))==1))=x2; y=y1.*y2; Matlab程序代码： n=-4:3; x=[1,2,5,6,8,7,4,5]; [x11,n11]=sigshift(x,n,3); [x12,n12]=sigshift(x,n,-5); [x13,n13]=sigfold(x,n); [x13,n13]=sigshift(x,n,2); [x01,n01]=sigadd(2*x11,n11,-4*x12,n12); [x1,n1]=sigadd(x01,n01,x13,n13); subplot(2,1,1),stem(n1,x1);title(x1(n)=2x(n-3)-4x(n+5)+x(2-n)); [x21,n21]=sigfold(x,n); [x21,n21]=sigshift(x,n,3); [x22,n22]=sigshift(x,n,4); [x22,n22]=sigmult(x,n,x22,n22); [x2,n2]=sigadd(x21,n21,x22,n22); subplot(2,1,2),stem(n2,x2);title(x2(n)=x(3-n)-x(n)x(n-4)); (2)已知：分别在4张图画出它的幅度、相位、实部和虚部。 n=-10:10; x=exp((-0.2+j*0.5)*n); magx=abs(x); angx=angle(x); realx=real(x); imagx=imag(x); subplot(2,2,1);plot(magx);title(x(n)的幅度); subplot(2,2,2);plot(angx);title(x(n)的相位); subplot(2,2,3);plot(realx);title(x(n)的实部); subplot(2,2,4);plot(imagx);title(x(n)的虚部); （3）已经：，；，，求卷积，并求出的能量信号 要求：分别用stem作出，要求在图中有起始位置和结束位置信息。 x=[-2,4,5,3,0,6,8,9];nx=-3:4; h=[4,5,9,5,-3,-1,5,7];nh=-4:3; [y,ny]=conv_m(x,nx,h,nh) subplot(3,1,1);stem(nx,x);xlabel(n);xlabel(x(n));title(x(n)的图像); subplot(3,1,2);stem(nh,h);xlabel(n);xlabel(h(n));title(h(n)的图像); subplot(3,1,3);stem(ny,y);xlabel(n);xlabel(y(n));title(y(n)的图像); E=sum(abs(y).^2) y = -8 6 22 63 86 66 48 122 197 153 36 -5 73 101 63 ny = -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6