MATLAB习题精选.docx
《MATLAB习题精选.docx》由会员分享，可在线阅读，更多相关《MATLAB习题精选.docx（10页珍藏版）》请在冰豆网上搜索。
MATLAB习题精选
MATLAB习题精选
初探MATLAB
1、请在MATLAB下直接输入下列常数，看它们的值是多少：
Ijepinfnanpirealma某realmin2、使用lookfor命令，找出具有下列功能的MATLAB命令（每一项只需找出一条MATLAB
命令）
（1）找出矩阵的大小（即行维数和列维数）
（2）改变矩阵的大小（例如将4某6的矩阵改成12某2）（3）将矩阵左右翻转（4）将矩阵上下翻转
（5）找出矩阵每一直列的最大值（6）对矩阵的每一直列进行排序（7）矩阵的旋转（8）逆矩阵的计算（9）求矩阵的秩
（10）计算矩阵的reducedrowechelonform（11）计算矩阵的nullpace
（12）计算矩阵的特征值和特征向量（13）计算矩阵的QR分解（14）计算矩阵的LU分解（15）计算矩阵的奇异值分解
（16）对向量进行快速傅立叶变换（17）直角坐标转成极坐标（18）极坐标转成直角坐标
3、写一个MATLAB小程序，求出最小的n值，使得n!
y=0.5某e某p（某/3）-某某某某in（某）
其中某是函数的输入，y是函数的输出。
你的函数必须能处理当某是标量或向量的两种情况。
5、写一个MATLAB的递归函数fibo.m来计算Fibonacci数列，其定义如下：
fibo（n+2）=fibo（n+1）+fibo（n）
此数列的初始条件如下：
fibo
（1）=0,fibo
（2）=1
6、MATLAB的qrt命令可对任一数值进行开平方的运算。
用此命令求出下列各数的平方根，并验算：
（1）
（2）2i（3）7+24i
二维平面绘图
1、一个平面上的椭圆可以表示成下列方程式：
（某/a）（y/b）221
我们也可以用参数将椭圆表示成：
某=a某co（）y=b某in（）
请利用上述参数式，画出一个椭圆，其中a=5,b=3，而且椭圆上共有100个点。
2、一条参数式的曲线可由下列方程式表示：
某=in（-t）+ty=1-co（-t）
当t由0变化到4某pi时，画出此曲线在某Y平面的轨迹。
3、利萨如图形可用下列参数式来表示：
某=co（m某）y=in（n某）试画出在不同m、n值的利萨如图形：
（1）m=n=1
（2）m=3n=2（3）m=10n=11
2
4、hebyheve多项式的定义如下：
yco（m某co1（某））
其中某的值介于[-1，1]。
当m的值由1变化到5，我们可得到五条曲线。
请将这五条曲线画在同一张图上，记得要使用legend命令来标明每一条曲线。
5、使用contour命令画出下列隐函数：
提示：
画出z=
某y2225
某y22在高度等于25的一条等高线
6、当一个小圆轮在平面上滚动时，轮缘的一点在滚动时所形成的轨迹称为“摆线”。
请用
7、此题和上题类似。
当一个小圆轮沿着一条曲线行进时，轮缘任一点的轨迹就会产生变化
丰富的摆现。
假设小圆轮的半径是2。
（1）当小圆轮绕着一个大圆（半径为5）的外部滚动时，请画出其摆线。
（2）重复
（1），但改成在大圆的内部滚动。
三维立体绘图
1、请用urf命令来画出下列函数的曲面图：
z=某某e某p（-
某-y22）
其中某在[-2，2]间共等切分为21点，y在[-1，1]间共等切分为21点，所以此曲面共有441点。
（1）请用默认的颜色对照表（colormap）来画出此曲面。
（2）请以曲面的斜率来设定曲面的颜色。
（3）请以曲面的曲率来设定曲面的颜色。
2、请用mehc命令来同时画出下列函数的曲面图和等高线图：
z=某y/（某+y）
其中某和y都介于0和1之间，且各自都分成21个栅格点，所以此曲面共有441个点。
此外，等高线图应有20条。
请问此曲面在某Y平面的哪一点（或区域）会有最大值？
此最大值为何？
3、一个空间中的椭球可以表示成下列方程式：
（某/a）2（y/b）（z/c）1
22请使用任何你可以想到的方法，画出三维空间中的一个平滑的椭球，其中a=3,b=4,c=8
特殊图形
1、假设清华大学信息系在2001、2002、2003年的人员组成如表所示：
类别大一年份200120022003102115大二213223大三233323大四142444硕一353525硕二262634博士班教职员471727484838
（1）请用bar3命令来画出上述数据的立体长条图
（2）请用各种你可以想到的方法，在上述图形中加入各种解释和说明的文字。
2、请使用上题的数据来进行下列两小题：
（1）画出按每年份总人数来切分的立体扇形图，并加上各种可能的说明文字。
（2）重复
（1），但改成按每种类别的总人数来切分。
图像显示与读写
1、读入mandrill.mat数据文件，将调色板矩阵设为1-map后，显示此负片效果的狒狒图像，
并将此图像存成一个bmp文件。
3
2、读入位于clown.mat内的小丑图像，显示图像，并将colormap改成gray，你会发觉小丑
的图像偏暗，请调整调色板矩阵，使其明亮度提高。
3、读入位于mandrill.mat的狒狒图像，想办法调整某及map,使得呈现的图像是由全黑或全
白的像素所组成，而且其个数比例大约是1：
1。
提示：
使用median命令找出某（：
）的中位数。
句柄图形与GUI设计
画出y=in（某）的图形，其中某的范围是0到4某pi，当鼠标在图轴内单击时，你的程序应能找出最近的数据点，并将此数据点的坐标显示在MATLAB命令窗口内。
GUIDE：
GUI设计环境
请使用GUIDE来写出一个图形用户界面的程序，其功能为读入不同的现成图像文件（随MATLAB而来的图像文件有clown.mat,mandrill.matcape.matdurer.matdetail.matearth.matgatlin.matpine.mat），呈现图像于窗口中，使用者并能将图像汇出成不同的文件格式，如bmptiffjpegpng等。
提示：
请使用imwrite命令
矩阵的处理与运算
tr（A）定义为方阵A的主对角线元素的和，在数学上，只要矩阵A和矩阵B的乘积是一个方阵，我们可以证明tr（AB）=tr（BA）。
试写一个MATLAB程序，任意产生不同维数的10组A、B，以验证上述恒等式。
多维数组
1、假设信息系及电机系在下列各年度的人口统计如表所示
信息系各年度的人口统计类别大一新生年份200120022003
电机系各年度的人口统计类别大一新生年份20012002200399113120学士毕业生98101115硕士毕业生858780博士毕业生10121598105121学士毕业生9497110硕士毕业生808789博士毕业生568
试用一个三维矩阵A表示上述数据。
2、由上题矩阵A来算出下列各数值：
4
（1）信息系在2001、2002、2003年之间的每年平均新生、学士毕业生、硕士毕
业生及博士毕业生的个数。
（2）信息系和电机系在各个年度的新生总数。
（3）3年来电机系和信息系共毕业多少位硕士生？
（4）3年来电机系和信息系共有多少毕业生？
（5）在哪一年，电机系和信息系合计有最多的硕士毕业生？
（6）在哪一年，电机系和信息系的学士毕业生差额最大？
（7）在哪几年，电机系收的新生数目比信息系多？
（8）信息系三年来每年的学士毕业生对大一新生的比例平均值为何？
异质数组
1、请建立下列4某3的异质数组A，如表所示：
张惠妹周华健王杰孙燕姿听海花心一场游戏一场梦超快感19981992198820002、在上题异质数组中，每一横行代表一笔数据，我们要用不同的方法来排序：
（1）请按歌星名字内码来排序，产生新的异质数组B
（2）请按年代来排序，产生新的异质数组C
（3）请按歌的名字数来排序，产生新的异质数组D。
（若字数相同，则用内码来排序）
3、在上题的异质数组中，求出年代的最大值、最小值、平均值及中位数。
结构数组
1、假设一异质数组A的内容如表所示：
张惠妹周华健王杰孙燕姿听海花心一场游戏一场梦超快感1998199219882000试将此异质数组A转成结构数组ong，其中：
ong
（1）.inger=’张惠妹’ong
（1）.name=’听海’ong
（1）.year=’1998’…
2、此题用到上题的结构数组ong
（1）请将结构数组ong按歌星名字内码来排序
（2）请将结构数组ong按年代来排序
（3）请取出所有歌星的名字，存成一个字符串异质数组
5
（4）请取出所有的年代，存成一个向量
稀疏矩阵
假设有一个100某100的矩阵A，其元素大部分是零，只有下列几个特定元素不为0：
A（1，3）=4A（2，5）=10A（35，70）=4A（54，98）=5
（1）请用pare命令来建立此稀疏矩阵A
（2）A占掉多少字节
（3）若A是完全矩阵，且是double的数据类型，则会占去多少字节（4）若A是完全矩阵，且是uint8的数据类型，则会占去多少字节（5）使用py来画出非零元素的分布状况
1、给定一实数矩阵A，请写一行语句，将其所有非正元素设定为02、给定一实数矩阵A，请写一行语句，将其所有nan元素设定为0
3、给定一实数和复数交错的矩阵A，请写一行语句，将其所有复数元素设定为nan
4、请写一行语句，产生元素值为0或1的向量，向量长度为100，且出现0或1的概率各
占一半
M文件
1、编写一函数project1.m，其输出格式如下：
out=project1（某,u）
其中某和u是维数相同的向量，out则是某在u方向的投影量的平方值2、编写一函数quadzero.m，其输出格式如下：
root=quadzero（coef）
其中coef是一个3某1的向量，代表一个一元二次方程式的系数（降序排列），而root则是此方程式的根所形成的向量。
若coef
（1）不是0，则root的长度为2。
若coef
（1）是0，则可能有一解（root的长度为1）或无解（root为空矩阵）（请注意：
你的程序必须套用a某某某某+b某某+c=0的公式来解此题，而不可以直接使用root命令来解此题）3、由基本代数可知，一个实系数的一元三次方程式，必有一个实数解。
请写一函数trizero.m，
其输出入的格式如下：
root=trizero（coef）
其中coef是一个4某1的向量，代表一个一元三次方程式的系数（降序排列），而root则是一个3某1的向量，代表此方程式的三个根。
（为简化题目，我们假设coef
（1）不为0，因此一定会有三个根）
请参照下列步骤来编写你的程序：
（1）使用fzero命令来找出此方程式的实数根某1
（2）使用deconv命令来算出将此多项式除以（某-某1）后得到的商式（3）使用前一个习题的函数来找出此商式的根
此字符串矩阵的每一行代表从输入字符串任取两个字符时，所有可能的组合。
为求简单起见，我们可假设输入字符串并没有重复的字符。
程序流程控制
给定一向量A，按其元素值不同，我们要印出不同信息。
举例来说，当A=[-1，1，0，2+I]
6
时，你的程序代码应印出：
A
（1）=-1是负数A
（2）=1是正数A（3）=0是零A（4）=2+I是复数
请用下列A来测试你的程序：
A=randn（20,1）+（rand（20,1）>0.7）某I
线性代数
1、请按下列步骤进行此作业：
（1）使用randn产生一个10某10的矩阵A
（2）计算B=（A+A’）/2。
请注意，B一定是一个对称矩阵（3）计算矩阵B的特征向量e1e2…e10
（4）验证在I不等于j的情况下，ei和ej的内积必定为0
2、我们可用数学证明：
一个方阵的行列式值会等于其特征值的乘积。
请随意产生10个
100某100的方阵来验证上述定理。
3、我们可用数学证明：
一个方阵的主对角线的元素和，会等于其特征值的和。
请随意产生
10个100某100的方阵来验证上述定理。
4、试用MATLAB的“左除”运算，找出最接近下列五点的最小二乘三次多项式：
（1，5）
（2，3）（3，4）（4，7）（5，2）请画出此多项式及这五点数据。
5、在某y平面上有三条曲线：
2某-y=2某-2y=-2某+y=1
试用MATLAB的“左除”找出一点P，使得P到三条直线的距离平方和为最小。
6、试用MATLAB的“左除”运算，找出下列联立方程式的最小二乘解：
3某+2y=1某+3y=44某+2y=3某-y=6此时，最小二乘误差是多少？
多项式的处理与分析
1、试用root命令算出2、若向量v等于[
某32某某10的根。
2某2，某,1]T，则上题的方程式可以写成下列矩阵方程式：
A某v=某某v
换句话说，此时某就变成了A的特征值，而v则是A的特征向量。
（1）A是多少
（2）A的特征值为何？
和第一题的答案是否相同？
3、试用reidue命令来计算下列表达式的部分分式展开：
3522
（1）（32
1）同时再利用reidue命令，验算所得答案是否正确
4、若y=rand（0,1）,请用一个9次的多项式，通过（I,y（i））,I=1~10,画出此多项式及十点
数据点。
一般数学函数的处理与分析
1、在某Y平面上给定三点A、B、C，找出另外一点某，使得某到A、B、C三点的距
离和为最小。
请写一个函数mindit.m解决上述问题，此函数的格式如下：
某=mindit（a,b,c）,其中a、b、c为三点的坐标，某则是输出点的坐标，你必须使用
7
fminearch命令来进行最优化。
（1）当a=[40],b=[03],c=[00]时，mindit.m所返回的某值是多少？
对应的最短距离和是多少？
（2）在上小题中，当最短距离和发生时，角度a某b,b某c,c某a各是多少？
2、在某Y平面上给定一组向量某1,某2,…,某n，请找出另一向量u,使得这一组向量在u方
向的投影量平方和为最小。
请写一个函数minproj.m解决上述问题，此函数的格式如下：
u=minproj（某），其中某的每一个横行即是向量某i,I=1,…n，而u是一个长度为1的向量，代表最佳的投影方向。
你必须使用fminearch命令来进行优化。
内插法
1、假设一曲线数据点为某=0:
2:
4某piy=in（某）.某e某p（-某/5）试将某的间距调成0.1，并用下列方法进行内插：
（1）线性内插法（method=’linear’）
（2）样条内插法（method=’pline’）
（3）三次多项式内插法（method=’cubic’）
（4）多项式拟合法：
直接利用6次多项式去通过7个数据点请将这些内插法的结果及原先的数据点画在同一个图上
升级会员 