matlab实验

实验一基础练习练习要求：1、练习数据和符号的输入方式，自定义数据和变量将前面的命令在命令窗口中执行通过sym_matrix=sym('[abc;JackHelp_MeNO_WAY]')sym_matrix=[a,b,c][Jack,Help_Me,NO_WAY]symsabc;M1=sym('Classical');M2=sym('Jazz');M3=sym('Blues');A=[abc;M1M2M3;sym([235])]A=[a,b,c][Classical,Jazz,Blues][2,3,5]2、输入A=[715;256;315],B=[111;222;333],在命令窗口中执行下列表达式，掌握其含义：A=[715;256;315]A=715256315B=[111;222;333]B=111222333A(2,3);%显示矩阵A第2行第3列的元素ansans=6A(:,2);%显示矩阵A第2列的元素ansans=151A(3,:);%显示矩阵A第3行的元素ansans=315A(:,1:2:3);%初值为1，增量为2，终值为3，显示矩阵A第1、3列的元素ansans=752635A(:,3).*B(:,2);%矩阵A的第3列元素与矩阵B第3列元素相乘ansans=51215A(:,3)*B(2,:);%矩阵A的第3列元素与矩阵B的第2行元素的矢量积ansans=101010121212101010A*B;%矩阵A与矩阵B的矢量积ansans=242424303030202020A.*B;%矩阵A、B的内积ansans=715410129315A^2;%等于A与A的矢量积ansans=661766423370381346A.^2;%等于A.*Aansans=49125425369125B/A;%矩阵右除，即1*ABansans=0.18420.2105-0.23680.36840.4211-0.47370.55260.6316-0.7105B./A;%即1*.ABansans=0.14291.00000.20001.00000.40000.33331.00003.00000.60003、输入C=1:2:20,则C（i）表示什么？其中i=1，2，3，……，10；C=1:2:20;%初值为1，增量为2，终值为20，C(i)表示C中第i个元素CC=135791113151719C(2)ans=3C(i)表示C中第i个元素4、查找已创建变量的信息，删除无用的变量symsabc;whosa%查看变量a的信息NameSizeBytesClassAttributesa1x1112symcleara%删除变量a的信息练习2编程练习要求：1、编写程序，计算1+3+5+7+……+（2n+1）的值（用input语句输入n值）在文本编辑器中输入：n=input('输入n的值：');s=0;fori=1:2:2*n+1s=s+i;enddisp('1+3+5+7+……+(2n+1)的值为');S输入n的值：41+3+5+7+……+(2n+1)的值为s=252、编写分段函数f(x)的函数文件，存放于文件ff.m中，计算出f(-3)f(2)f(∞)的值在文本编辑器中输入：functions=ff(x)if(x=0)&&(x1)s=x;elseif(x=1)&&(x=2)s=2-x;elses=0;end在MATLAB命令窗口中输入：s=ff(-3)s=0s=ff(sqrt(2))s=0.5858s=ff(inf)s=0练习3矩形计算练习要求：将前面的命令在命令窗口中执行通过，掌握矩阵运算的基本方法和操作。例3-1A=[1,1,1;1,2,3;1,3,6];B=[8,1,6;3,5,7;4,9,2];A+Bans=92747105128A-Bans=-70-5-2-3-4-3-64例3-2A=[1,1,1;1,2,3;1,3,6];B=[8,1,6;3,5,7;4,9,2];A+Bans=92747105128A-Bans=-70-5-2-3-4-3-64X=[23451221];Y=[011110001100];X*Yans=8563332*Xans=468102442X=[-102];Y=[-2-11];dot(X,Y)ans=4sum(X.*Y)ans=4例3-3a=[123];b=[456];c=cross(a,b)c=-36-3例3-4a=[123];b=[456];c=[-36-3];x=dot(a,cross(b,c))x=54a=[123;426;749];b=[4;1;2];x=a\bx=-1.50002.00000.5000例3-5A=[123;221;343];inv(A)ans=1.00003.0000-2.0000-1.5000-3.00002.50001.00001.0000-1.0000B=[1,2,3,1,0,0;2,2,1,0,1,0;3,4,3,0,0,1];C=rref(B)%化行最简形C=1.0000001.00003.0000-2.000001.00000-1.5000-3.00002.5000001.00001.00001.0000-1.0000X=C(:,4:6)X=1.00003.0000-2.0000-1.5000-3.00002.50001.00001.0000-1.0000例3-6A=[123;221;343];D=det(A)D=2.0000例3-7A=sym('[1/x,1/(x+1);1/(x+2),1/(x+3)]');B=sym('[x,1;x+2,0]');C=B-AC=[x-1/x,1-1/(x+1)][x-1/(x+2)+2,-1/(x+3)]D=A\BD=[-x*(2*x^2+7*x+6),(x*(x^2+3*x+2))/2][2*(x+1)^2*(x+3),-(x*(x+1)*(x+3))/2]例3-8symsxfactor(x^9-1)ans=(x-1)*(x^2+x+1)*(x^6+x^3+1)例3-10symsxyp=expand((x+1)^3)p=x^3+3*x^2+3*x+1q=expand(sin(x+y))q=cos(x)*sin(y)+cos(y)*sin(x)实验2离散时间信号分析1、利用matlab语言编程产生和汇制下列有限长序列：（1）、单位脉冲序列)(nn=0:15;x=[1zeros(1,15)];stem(n,x);grid;xlabel('n');ylabel('x(n)');title('\delta(n)的线图');axis([-116-0.11.1]);line([-1,16],[0,0]);（2）、单位阶跃序列u（n）n=-5:15;x=[zeros(1,5)ones(1,16)];stem(n,x);grid;xlabel('n');ylabel('x(n)');title('u(n)的线图');axis([-616-0.11.1]);（3）、矩形序列)(8nRn=-5:10;R8=[(n=0)&(n=7)];stem(n,R8);grid;xlabel('n');ylabel('x(n)');title('R8(n)的线图');axis([-611-0.11.1]);（4）、正弦型序列)35sin()(nAnxn=-21:21;x=2*sin(0.2*pi*n+1*pi/3);n1=-24:0.1:24;x1=2*sin(0.2*pi*n1+1*pi/3);stem(n,x,'.');holdon;plot(n1,x1,'--');xlabel('n');ylabel('x(n)');title('2sin(0.2\pin+1\pi/3)的线图');axis([-23.523.5-2.12.1]);（5）任意序列)4(5)3(4)2(3)1(2)()(nnnnnnx)3(2)2()1(2)()(nnnnnhn1=0:4;x1=[12345];n2=0:3;x2=[1212];subplot(211),stem(n1,x1,'.');xlabel('n');ylabel('x(n)');title('x(n)=\delta(n)+2\delta(n-1)+3\delta(n-2)+4\delta(n-3)+5\delta(n-4)的线图');axis([-1506]);grid;subplot(212),stem(n2,x2,'.');xlabel('n');ylabel('h(n)');title('h(n)=\delta(n)+2\delta(n-1)+\delta(n-2)+2\delta(n-3)的线图');axis([-1403]);grid;2、序列的运算利用matlab语言编程完成上诉两序列的移位、反褶、和、积、累加等运算。并绘制运算后序列的波形。（1）、移位clearall;closeall;clc;nx=0:4;x=[12345];y=x;ny=nx+3;%右移3位subplot(121),stem(nx,x,'.');axis([-19-16]);grid;xlabel('n');ylabel('x(n)');title('原序列');subplot(122),stem(ny,y,'.');axis([-19-16]);grid;xlabel('n');ylabel('y(n)');title('右移3位后的序列');（2）、反褶clearall;closeall;clc;nx=0:4;x=[12345];ny=-fliplr(nx);y=fliplr(x);subplot(211),stem(nx,x,'.');axis([-55-16]);grid;xlabel('n');ylabel('x(n)');title('原序列');subplot(212),stem(ny,y,'.');axis([-55-16]);grid;xlabel('n');ylabel('y(n)');title('翻转后的序列');（3）、和与积clearall;closeall;clc;n1=0:4;x1=[12345];n2=0:3;x2=[1212];n=min(min(n1),min(n2)):max(max(n1),max(max(n2)));x11=zeros(1,length(n));x22=zeros(1,length(n));x11(find(n=min(n1))&(n=max(n1)))=x1;x22(find(n=min(n2))&(n=max(n2)))=x2;xa=x11+x22;xm=x11.*x22;subplot(211),stem(n,xa);axis([-1507]);grid;xlabel('n');ylabel('x(n)+h(n)');title('x(n)和h(n)的和');subplot(212),stem(n,xm);axis([-1509]);grid;xlabel('n');ylabel('x(n)*h(n)');title('x(n)和h(n)的积');(5)、累加clearall;closeall;clc;n=0:4;x1=[12345];x2=cumsum(x1);subplot(211),stem(n,x1);axis([-15016]);grid;xlabel('n');ylabel('x(n))');title('原序列');subplot(212),stem(n,x2