matlab经典习题及解答

第1章MATLAB概论1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。1.2MATLAB系统由那些部分组成？MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。1.4MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock……菜单项就可以将独立的窗口重新防止的桌面上。1.5如何启动M文件编辑/调试器？在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。1.6存储在工作空间中的数组能编辑吗？如何操作？存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。1.7命令历史窗口除了可以观察前面键入的命令外，还有什么用途？命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。1.8如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的SetPath菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。1.9在MATLAB中有几种获得帮助的途径？在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLABHelp菜单项可以打开帮助浏览器；（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help函数名”可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗口键入“lookfor关键词”可以搜索出一系列与给定关键词相关的命令和函数（4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数。注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。第2章MATLAB矩阵运算基础2.1在MATLAB中如何建立矩阵194375，并将其赋予变量a？a=[573;491]2.2有几种建立矩阵的方法？各有什么优点？可以用四种方法建立矩阵：①直接输入法，如a=[2573]，优点是输入方法方便简捷；②通过M文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改；③由函数建立，如y=sin(x)，可以由MATLAB的内部函数建立一些特殊矩阵；④通过数据文件建立，该方法可以调用由其他软件产生数据。2.3在进行算术运算时，数组运算和矩阵运算各有什么要求？进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a与b相乘（a*b）时必须满足a的列数等于b的行数。2.4数组运算和矩阵运算的运算符有什么区别？在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b为矩阵乘，a.*b为数组乘。2.5计算矩阵897473535与638976242之和。a=[535;374;798];b=[242;679;836];a+bans=777914131512142.6求i44i93i49i67i23i57i41i72i53i84x的共轭转置。x=[4+8i3+5i2-7i1+4i7-5i;3+2i7-6i9+4i3-9i4+4i];x’ans=4.0000-8.0000i3.0000-2.0000i3.0000-5.0000i7.0000+6.0000i2.0000+7.0000i9.0000-4.0000i1.0000-4.0000i3.0000+9.0000i7.0000+5.0000i4.0000-4.0000i2.7计算572396a与864142b的数组乘积。a=[693;275];b=[241;468];a.*bans=12363842402.8“左除”与“右除”有什么区别？在通常情况下，左除x=a\b是a*x=b的解，右除x=b/a是x*a=b的解，一般情况下，a\bb/a。2.9对于BAX，如果753467294A，282637B，求解X。A=[492;764;357];B=[372628]’;X=A\BX=-0.51184.04271.33182.10已知：987654321a，分别计算a的数组平方和矩阵平方，并观察其结果。a=[123;456;789];a.^2ans=149162536496481a^2ans=3036426681961021261502.11463521a，263478b，观察a与b之间的六种关系运算的结果。a=[123;456];b=[8–74;362];abans=010101a=bans=010101abans=101010a=bans=101010a==bans=000000a~=bans=1111112.127.0802.05a，在进行逻辑运算时，a相当于什么样的逻辑量。相当于a=[11011]。2.13在sin(x)运算中，x是角度还是弧度？在sin(x)运算中，x是弧度，MATLAB规定所有的三角函数运算都是按弧度进行运算。2.14角度604530x，求x的正弦、余弦、正切和余切。x=[304560];x1=x/180*pi;sin(x1)ans=0.50000.70710.8660cos(x1)ans=0.86600.70710.5000tan(x1)ans=0.57741.00001.7321cot(x1)ans=1.73211.00000.57742.15用四舍五入的方法将数组[2.45686.39823.93758.5042]取整。b=[2.45686.39823.93758.5042];round(b)ans=26492.16矩阵728365219a，分别对a进行特征值分解、奇异值分解、LU分解、QR分解及Chollesky分解。[v,d]=eig(a,b)v=-0.4330-0.2543-0.1744-0.56570.9660-0.6091-0.70180.04720.7736d=13.54820004.83030003.6216a=[912;563;827];[u,s,v]=svd(a)u=-0.56010.5320-0.6350-0.4762-0.8340-0.2788-0.67790.14620.7204s=15.52340004.56480003.3446v=-0.82750.3917-0.4023-0.3075-0.9156-0.2592-0.4699-0.09070.8781[l,u]=lu(a)l=1.0000000.55561.000000.88890.20411.0000u=9.00001.00002.000005.44441.8889004.8367[q,r]=qr(a)q=-0.69030.3969-0.6050-0.3835-0.9097-0.1592-0.61360.12210.7801r=-13.0384-4.2183-6.82600-4.8172-1.0807003.7733c=chol(a)c=3.00000.33330.666702.42671.1447002.29032.17将矩阵4257a、3817b和2695c组合成两个新矩阵：（1）组合成一个43的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按列顺序排列的b矩阵元素，第三列为按列顺序排列的c矩阵元素，即237912685574（2）按照a、b、c的列顺序组合成一个行矢量，即296531877254a=[42;57];b=[71;83];c=[59;62];%(1)d=[a(:)b(:)c(:)]d=475586219732%(2)e=[a(:);b(:);c(:)]'e=452778135692或利用（1）中产生的de=reshape(d,1,12)ans=452778135692第3章数值计算基础3.2求解多项式x3-7x2+2x+40=0的根。r=[1-7240];p=roots(r);-0.21510.44590.79490.27073.4计算多项式乘法(x2+2x+2)(x2+5x+4)。c=conv([122],[154])c=17161883.5计算多项式除法(3x3+13x2+6x+8)/(x+4)。d=deconv([31368],[14])d=3123.7计算多项式9514124234xxxx的微分和积分。p=[4–12–145];pder=polyder(p);pders=poly2sym(pder)pint=polyint(p);pints=poly2sym(pint)pders=12*x^2-24*x-14pints=x^4-4*x^3-7*x^2+5*x3.8解方程组66136221143092x。a=[290;3411;226];b=[1366]';x=a\bx=7.4000-0.2000-1.40003.11矩阵943457624a，计算a的行列式和逆矩阵。a=[42-6;754;349];ad=det(a)ai=inv(a)ad=-64ai=-0.45310.6562-0.59370.7969-0.84370.9062-0.20310.1562-0.09373.12y=sin(x)，x从0到2，x=0.02，求y的最大值、最小值、均值和标准差。x=0:0.02*pi:2*pi;y=sin(x);ymax=max(y)ymin=min(y)ymean=mean(y)ystd=std(y)ymax=1ymin=-1ymean=2.2995e-017ystd=0.70713.1354321x，108642y，计算x的协方差、y的协方差、x与y的互协方差。x=[12345];y=[246810];cx=cov(x)cy=cov(y)cxy=cov(x,y)cx=2.5000cy=10cxy=2.50005.00005.000010.00003.15有一正弦衰减数据y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi，用三次样条法进行插值。x0=0:pi/5:4*pi;y0=sin(x0).*exp(-x0/10);x=0:pi/20:4*pi;y=spline(x0,y0,x);plot