MATLAB上机习题一
MATLAB上机模拟考试题
M A T L A B上机模拟考试题(一)(总8页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--(1)在MATLAB的命令窗口中执行_____命令，将命令窗口的显示内容清空。
（）off(2)在MATLAB的命令窗口中执行_____命令，使数据输出显示为十六进制表示。
（） long rat hex short e(3)下列变量名中_____是合法的。
（） *y,a,1 \y,a1234 ,1 bcx ,i,j(4)已知x=0:5,则x有_____个元素。
（）(5)一下运算符中哪个的优先级最高_____。
（） A./ B.^ C.~=D.&(6)使用检测函数isnumeric(10)的结果是_____。
（）(7)三维图形中默认视角是_____度。
（） A.方位角=0 俯仰角=90B.方位角=90 俯仰角=0C.方位角= 仰俯角=30D.方位角=0 仰俯角=180(8)将符号表达式化简为因式分解因式分解因式分解因式分解形式，使用_____函数。
（）(9)运行以下命令，则_____描述是正确的。
（）>>syms a b cd >>A=[a b;c d] 占用的内存小于100B B.创建了5个符号变量占用的内存是a b c d的总和 D.不存在（10）已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。
（11）运行命令bitor（8,7）的结果是_____。
（12）运行以下命令：>>x=0:10; >>y1=sin(x); >>y2=5*sin(x);>>y3=[10*sin(x);20*sin(x)]; >>plot(x,y1,x,y2,x,y3) 则在一个图形窗口中，可以看到_____条曲线。
（13）符号表达式“g=sym(sin(a*z)+cos(w*v))”中的自由符号变量是_____。
MATLAB 上机 习题及答案
15、今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5，要求先求得P(x)=P1(x)+P2(x)，然后计算xi=0.2*i各点上的P(xi)（i=0,1,2,…,5）值。
p1=[1.0 0.0 0.0 -2.0 1.0];>> p2=[0.0 0.0 1.0 4.0 -0.5];>> p1x=poly2sym(p1);p2x=poly2sym(p2);>> p=p1x+p2xp =x^4+2*x+1/2+x^2>> x=0:5;>> x.^4+2*x+1/2+x.^2ans =0.5000 4.5000 24.5000 96.5000 280.5000 660.50001、试个MATLAB的工作空间中建立以下2个矩阵：A＝[1 2]1234B⎡⎤=⎢⎥⎣⎦，求出矩阵A和B的乘积，并将结果赋给变量C。
>> A=[1 2]A =1 2>> B=[1 23 4]B =1 23 4>> C=A*BC =7 102、利用MATLAB提供的帮助信息，了解inv命令的调用格式，并作简要说明。
help invINV Matrix inverse.INV(X) is the inverse of the square matrix X.A warning message is printed if X is badly scaled ornearly singular.See also SLASH, PINV, COND, CONDEST, LSQNONNEG, LSCOV. Overloaded methodshelp gf/inv.mhelp zpk/inv.mhelp tf/inv.mhelp ss/inv.mhelp lti/inv.mhelp frd/inv.mhelp sym/inv.mhelp idmodel/inv.m3、使用help命令查询函数plot的功能以及调用方法，然后利用plot命令绘制函数y=sin(x)的图形，其中0xπ≤≤。
Matlab上机练习题及答案
Matlab 上机练习题及答案---------------------------------------------------------------------1、矩阵Y=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3472123100451150425，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。
取出子矩阵⎥⎦⎤⎢⎣⎡21301，并求该矩阵的维数。
解：命令为：Y=[5,2,4;0,15,1;45,100,23;21,47,3]Y(2,3)Y(10)sub2ind([43],2,3)[i,j]=ind2sub([43],10)find(Y==100)sub2ind([43],3,2)B=Y(2:2:4,3:-2:1)或B=Y([24],[31])[m n]=size(Y)---------------------------------------------------------------------2、已知矩阵A=[10-1；241；-205]，B=[0-10；213；112]求2A+B 、A 2-3B 、A*B 、B*A 、A.*B ，A/B 、A\B 解：命令为：A=[10-1；241；-205]B=[0-10；213；112]E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B---------------------------------------------------------------------3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵，并计算两者的乘积。
解：命令为：A=eye(3,4)B=8*ones(4)C=A*B---------------------------------------------------------------------4、创建矩阵a=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------7023021.5003.120498601,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b构成矩阵d ，计算a*b 、c<d ,c&d,c|d ,~c|~d 解：命令为：a=[-1,0,-6,8;-9,4,0,12.3;0,0,5.1,-2;0,-23,0,-7]b=a(:,[12])c=a([12],:)d=b’e=a*b f=c<d g=c&d h=c|d i=~c|~d---------------------------------------------------------------------5、求!201∑=n n 解：命令文件为sum=0;s=1;for n=1:20s=n*s;sum=sum+s;end sum---------------------------------------------------------------------6、求a aa aaa aa a S n ++++=得值，其中a 是一个数字，由键盘输入，表达式中位数最多项a 的个数，也由键盘输入。
MATLAB上机练习题(一)
MATLAB 上机练习题（一）
1. 计算下式
7
2009/42008/41949222?+-⨯-⨯=
2. 试计算d
3461.1,30,5b d a b c
a b c +=-+===其中
要求只输出d 的值，其他变量的值不输出。
**************************************************************
3. 计算下式：
181(2)sin()cos()6,3
y
x x y π-+-==（，其中 0.1104
2log 805(3)2009
6(4)1,5;
2e
e x x x -⨯+-=其中
4. 完成以上题目后清除工作空间中的所有变量。
**************************************************************
5. 求下列方程的根：
（1）325x x -=
（2）32340x x x -+=
6. 根据课本任意找出4个矩阵在MA TLAB 中输入和输出。
附1：增加的一些常用函数pi：常数
round(x)：四舍五入至最近整数
fix(x) ：无论正负，舍去小数至最近整数min(x): 向量x的元素的最小值
max(x): 向量x的元素的最大值
diff(x): 向量x的相邻元素的差
length(x): 向量x的元素个数
sin(x)：正弦函数
cos(x)：余弦函数
tan(x)：正切函数
asin(x)：反正弦函数
acos(x)：反余弦函数
atan(x)：反正切函数。
MATLAB上机题目
1、熟悉MATLAB 运行环境，找出命令历史窗口、工作空间窗口，记下MATLAB 当前浏览路径。
2、求下列联立方程的解
8102569
583
24754
12743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x
3、MATLAB 环境下，给出下列值的大小
eps=？ Realmax=? Realmin=?
4、生成矩阵
a=102030405060710⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦
，取出第三行第二列元素；一、二行，二至五列元素（用一个命令）；
a(17)=?
若a(1,2)=[],则a=？
若a(20)=4.5,则a=？
5、i=inf,j=NaN,则i/j=? i*j=?
6、生成矩阵
b=123435672i i i i +-⎡⎤⎢⎥-⎢⎥⎢⎥+⎣⎦
,求其共轭、转置、共轭转置矩阵 7、检查工作空间中的变量，以及这些变量的详细信息
8、列出2×2阶的单位矩阵I ，4×4阶的魔方矩阵M 和4×2阶的全1矩阵A 、全零矩阵B 。
(2)将这些矩阵拼接为6×6阶的矩阵C ：
⎥⎦
⎤⎢⎣⎡=M B A I C ' (3)取出C 的第2，4，6行，组成3×6阶的矩阵1C ，取出第2，4，6列，组成6×3阶的矩阵2C
9、产生8×6阶的正态分布随机数矩阵1R ,求其各列的平均值和标准差，并求其全体的平均值和标准差。
10、产生4×6阶的均匀分布随机数矩阵R ，要求其元素在1到16之间，并求此矩阵前四列组成的方正的逆阵。
matlab习题(1)
1.输入20个数，求其中最大数和最小数。
要求分别用循环结构和调用MA TLAB的max函数、min函数来实现。
a=rand(1,20);max1=max(a) %用max函数求最大值min1=min(a) %用max函数求最大值%以下用循环方法求解：max2=a(1);min2=a(1);for i=2:20if a(i)>max2max2=a(i);endif a(i)<min2min2=a(i);endendmax2min2调用MATLAB的max函数、min函数来实现:>> a=rand(1,20)zuidazhi=max(a)zuixiaozhi=min(a)a =Columns 1 through 80.7513 0.2551 0.5060 0.6991 0.8909 0.9593 0.5472 0.1386 Columns 9 through 160.1493 0.2575 0.8407 0.2543 0.8143 0.2435 0.9293 0.3500 Columns 17 through 200.1966 0.2511 0.6160 0.4733zuidazhi =0.9593zuixiaozhi =0.13862. 求Fibonacci数列(1)大于4000的最小项。
(2)5000之内的项数。
(1)function t=fibonacc1f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendt=f(n);return(2)function n=fibonacc2f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendn=n-1;return3. 写出下列程序的输出结果:s=0;a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=afor j=1:4if rem(k(j),2)~=0s=s+k(j);endendendSAns=1081、解方程组Ax ＝b ，分别用求逆解法与直接解法求其解。
matlab上机练习（附答案）
1.以下两种说法对吗？（1）MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。
（2）MATLAB指令窗中显示的数据有效位数不超过七位。
2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同？DIARY filename causes a copy of all subsequent command window inputand most of the resulting command window output to be appended to thenamed file. If no file is specified, the file 'diary' is used.DIARY OFF suspends it.DIARY ON turns it back on.DIARY, by itself, toggles the diary state.Use the functional form of DIARY, such as DIARY('file'),when the file name is stored in a string.3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上？“位于搜索路径上的目录”与“当前目录”在MATLAB中的功用相同吗？4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组？用直接键入法？用数组编辑器？用M文件编辑器？5.运用数组算术运算符去掉下面程序里的for/end循环：x=11:15for k=1:length(x)z(k)=x(k)^2+2.3*x(k)^0.5;endx=11:15 z1=x.^2+2.3*x.^0.56.不使用数组算术运算符，重写下面的程序代码：x=[2 1 4]z=1./(1+x.^2)x=2;k=1;while i<=4,z2(k)=1/(1+i^2);i=i+1;x=x+1;end7．某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
matlab上机练习答案
实验一MA TLAB的基本命令与基本函数１已知矩阵a =11 12 13 1421 22 23 2431 32 33 3441 42 43 44求(1) A(:,1) (2) A(2,:)(3) A(:,2:3) (4) A(2:3,2:3)(5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:)(9) ones(2,2) (10) eye(2)(11) [A,[ones(2,2);eye(2)]](12) diag(A) (13) diag(A,1)(14) diag(A,-1) (15) diag(A,2)2(1)输入如下矩阵Ａ０π/3Ａ＝π/6 π/2(2) 求矩阵B1,B1中每一元素为对应矩阵Ａ中每一元素的正弦函数(3) 求矩阵B2, B2中每一元素为对应矩阵Ａ中每一元素的余弦函数(4) 求B12+B22(5) 求矩阵Ａ的特征值与特征矢量：称特征矢量为Ｍ，而特征值矩阵为Ｌ(6) 求Msin(L)M-13已知水的黏度随温度的变化公式为μ=μ0/(1+at+bt2)其中μ0=1.785×10-3，a=0.03368,b=0.000221,求水在0，20，40，80℃时的黏度。
程序如下：miu0=1.785e-3;a=0.03368;b=0.000221;t=0:20:80miu=miu0./ (1+a*t+b*t.^2)（2）一个长管，其内表面半径为a，温度为Ta ；外表面半径为Tb；则其径向和切向应力可分别表示为：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛+----=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛----=r b a b r b a b a a b v T T E r b a b r b a b a a b v T T E b a t b a r ln ln 11)/ln()1(2)(ln ln 1)/ln()1(2)(2222222222ασασ式中r 为管子的径向坐标，E 为管子材料的弹性模量，ɑ为热膨胀系数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB上机习题一
请按以下步骤完成上机实验：
1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内；
2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面；
3）如果习题是问答题，请将答案写在题目后；
4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后；
5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”；
6）上传该文件到FTP的相关目录。
1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。
2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。
3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。
4. 通过使用帮助确定内置变量ispc 的含义。
6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。
5. 创建double类型的变量，并进行计算
1）a=87，b=190，计算a+b、a-b、a*b；
2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致。
6. 计算如下表达式：
1）()sin 60
2）3e
3）3cos 4π⎛⎫ ⎪⎝⎭
4）2
7562323336
+⨯-⨯⨯
7. 设u=2，v=3，计算：
1）4log uv v
2）()
2
2u e v v u +-
3
8. 计算如下表达式：
1）()()3542i i -+
2）()sin 28i -
9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较：
1）420<
2）420≤
3）420==
4）4~20=
5）''''b B <
10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b >
2）a b <
3）&&a b b c >>
4）a d ==
5）|a b c >
6）~~d
11. 判断下列变量名哪些是合法的MATLAB 变量：
1）fred
7）fred! 2）book_1
8）book-1 3）2ndplace
9）Second_Place 4）#1
10）No_1 5）vel_5
11）vel.5 6）tan
12）while 提示，可以使用isvarname 对上述变量名进行检验。
12. 编写脚本程序，命名为barbell.m ，完成如下计算：
如图1所示的杠铃，
每个球的半径是10cm ，两个求直接链接杆的长度是15cm ，杆的直径为1cm ，计算杠铃的体积和表面积；使用disp 函数显示计算结果。
（提示：disp 函数的用法参考教材的例2-31，面积计算时相交部分球冠可近似为圆，体积计算时可忽略重叠的球冠体积）
解答：
barbell.m （拷贝文件内容到如下方框中）：
运行屏幕截图：
13. 在高中课本中大家学过描述气体压强（p ）、温度（T ）、体积（V ）和气体摩尔数（n ）之间的理想气体状态方程：
pV nRT =
其中，R 为理想气体常数。
上述状态方程描述的是低压强和高温度时气体的特性。
在1873年，范德瓦尔斯对这个方程进行了修正，使其更好地描述气体在不同压强和温度条件下的状态，该方程被称为范德瓦尔斯方程：
()22n a p V nb nRT V ⎛⎫+-= ⎪⎝
⎭ 该方程新增了a 和b 两个变量，用来表示纯净气体的属性。
1）编写脚本程序，命名为state_equation.m ，利用如下数据：
分别使用理想气体状态方程和范德瓦尔斯方程计算水蒸气的温度，并使用disp 函数显示计算结果。
2）根据范德瓦尔斯方程及第1问的数据，编写getTemperature(p,V,n)函数实现温度计算，其中a、b、R作为常量。
调用该函数测试计算结果。
（提示，注意单位制，1 bar=105 Pa）
解答：
state_equation.m（拷贝到如下方框中）：
运行屏幕截图：
getTemperature.m（拷贝到如下方框中）：
运行屏幕截图：
14. 种群繁殖符合指数规律：
0rt
P P e
其中，P为目前的种群数，P0为原始种群数，r为持续增长率，t为时间。
现有100只兔子，持续增长率是每年90%（r=0.9），编写rabbit(n)函数用来计算n年后兔子的数量，并测试10年结果。
（提示，注意利用2.2.1节中的数学函数）
解答：
rabbit.m（拷贝到如下方框中）：
运行屏幕截图：
15. 一般来说，站得越高会看的越远，那么到底能看多远呢？主要取决于山的高度和地球的半径，如下图所示：
由于月球和地球的半径相差很大，所以能看到的地平线距离相差也很远，现已知地球半径是6378km，月球半径是1737km，编写脚本文件horizon.m计算0到8000m高的山上（珠穆朗玛峰海拔8850m）在地球和月球上看到地平线的距离分别是多少。
解答：
horizon.m（拷贝文件内容到如下方框中）：
运行屏幕截图：。