MATLAB的基本操作与编程基础.doc
实验一MATLAB的基本操作、编程基础姓名： 学号： 网选班级： 1.简述实验目的及实验原理答：本次实验的目的是为了了解matlab软件和matlab语言的用法。简要介绍了一些常用的运算符和函数，还有一些信号系统在matlab软件里怎样实现的。实验原理是通过matlab软件实现信号系统的一些计算。2.记录实验内容3、4的命令或程序及图形内容3.画出以下各序列在给定区间的图形。1) x(n)=2(n+2)-(n-4)+(n) ,-5n52) x(n)=nu(n)-u(n-10)+10e-0.3(n-10)u(n-10)-u(n-20),0n203) x(n)=2-cos(0.4n)- cos(0.8n） ,0n50答：（1）程序为：n=-5:5;x1=(n+2)=0;x2=(n-4)=0;x3=(n+0)=0;x=2.*x1-x2+x3;stem(n,x);xlabel('n');ylabel('x');(2)程序为：%x(n)=nu(n)-u(n-10)+10e-0.3(n-10)u(n-10)-u(n-20),n=0:20;x1=n>=0;x2=n>=10;x3=n>=20;x4=exp(-0.3).(n-10);x=x1-x2+10.*x4.*(x2-x3);stem(n,x);xlabel('n 作者：陈哲');ylabel('x');（3）程序为：%x(n)=2-cos(0.4¦Ðn)- cos(0.8¦Ðn£© ,0¡Ün¡Ü50%x(n)=2-cos(0.4¦Ðn)- cos(0.8¦Ðn£© ,0¡Ün¡Ü50n=0:50;x1=cos(0.4*n.*pi);x2=cos(0.8*n.*pi);x=2-x1-x2;subplot(211);stem(n,x);xlabel('n 作者：陈哲');ylabel('x');title('采样点很少的情况');n1=0:0.01:50;x11=cos(0.4*n1.*pi);x21=cos(0.8*n1.*pi);x0=2-x11-x21;subplot(212);stem(n1,x0);xlabel('n 作者：陈哲');ylabel('x');title('采样点很多的情况');4.令x(n)=1,2,3,4,5,6,7,6,5,4,3,2,1,画出下列序列的图形1) x1(n)=2x(n-5)-3x(n-4)2) x2(n)=x(3-n)-x(n)x(n-2)3) xe(n)和xo(n) -通过编写求xe(n)、xo(n)的M-file来实现(需要给出测试函数) xe(n) 和xo(n)分别为x(n)的偶部和奇部解：(1)程序为：%x1(n)=2x(n-5)-3x(n-4),x(n)=1,2,3,4,5,6,7,6,5,4,3,2,1%x1=x(n-5) x2=x(n-4)x=1,2,3,4,5,6,7,6,5,4,3,2,1;n=0:12;x1,n1=sigshift(x,n,5);x2,n2=sigshift(x,n,4);y,n3 = sigadd(2.*x1),n1,(-3.*x2),n2)subplot(211);stem(n,x);xlabel('n');ylabel('x(n)');subplot(212);stem(n3,y);xlabel('n');ylabel('x1(n)=2x(n-5)-3x(n-4)');自定义函数有：%实现信号加的M-file :y(n)=x1(n)+x2(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; %移位的M-file :y(n)=x(n-n0)function y,n = sigshift(x,m,n0)n = m+n0; y = x;（2）程序为：%x2(n)=x(3-n)-x(n)x(n-2)%x0=x(3-n);x1=x(n)x(n-2);x=1,2,3,4,5,6,7,6,5,4,3,2,1;n0=0:12;x00,n10 = sigshift(x,3,n0);x0,n1 = sigfold(x00,n10);x11,n11 = sigshift(x,2,n0);x1,n2 = sigmult(x,n0,x11,n11);x2,n = sigadd(x0,n1,-x1,n2);stem(n,x2);xlabel('n');ylabel('x2(n)');title('x2(n)=x(3-n)-x(n)x(n-2)');自定义函数有：%实现信号乘的M-file : y(n)=x1(n) x2(n)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; %折叠的M-file :y(n)=x(-n)function y,n = sigfold(x,n)y = fliplr(x); n = -fliplr(n);（3）程序为：奇部函数：function y1,m1 =jishu( x,n ) for i=1:2:n if i<=n m1(i+1)/2)=(i+1)/2; y1(i+1)/2)=x(i); end end偶部函数：function y1,m1 =oushu( x,n ) for i=2:2:n if i<=n m1(i)/2)=(i)/2; y1(i)/2)=x(i); end end测试程序：x=1,2,3,4,5,6,7,6,5,4,3,2,1;n0=1:13;n=13;y1,m1=jishu(x,n);y2,m2=oushu(x,n);subplot(311);stem(n0,x);xlabel('n');ylabel('x(n)');title('x(n)');subplot(312);stem(m1,y1);xlabel('n');ylabel('x(n)Æ沿');title('x(n)Æ沿');subplot(313);stem(m2,y2);xlabel('n');ylabel('x(n)ż²¿');title('x(n)ż²¿');
