matlab15练习.docx
《matlab15练习.docx》由会员分享，可在线阅读，更多相关《matlab15练习.docx（13页珍藏版）》请在冰豆网上搜索。
matlab15练习
1.请生成从0到10间隔为0.1的列向量。
A=0:
0.1:
10
2.矩阵A为5×5的魔方矩阵，请从A中提取8和16产生新的列向量B。
a=magic（5）
b=（a[12],[45]）
已知A=
请通过矩阵变换，由A生成B=
A=[1234;5678;9123;4567];
B=A（:
[4231]）
3.写出产生序列[
…
]的程序。
A=pi/100:
pi/100:
pi
4.下列变量名中,属于合法变量名的是____A_________。
A.flower2B.2flowerC._whatD.who's_it
5.矩阵M=[112733;295712;734537]
则M（2,3）=____A_________。
A.12B.29C.37D.45
6.表达式'
'的执行结果是_____B________。
A.a=38.18B.a=293.18C.a=292.18D.a=37.56
7.某同学设计了一个程序文件myprogram.m，并将其保存到了c:
\mydocument中，但在命令窗口中输入文件名>>myprogram后，MATLAB系统提示：
?
?
?
Undefinedfunctionorvariable‘myprogram’
试分析产生错误的原因并给出解决办法。
没有修改路径；file-setpath-addFolder
8.MATLAB语句后，加上___A______，则运行时不显示中间结果。
A.分号B.逗号C.句号D.空格
9.使用语句t=0:
110生成_____C________个元素的向量。
A.1100B.910C.111D.110
10.已知两个矩阵：
，
将矩阵A中所有等于-1的元素改为-3，将矩阵B中等于0的元素值改为新的矩
阵A中相应位置元素的值，最后将新的矩阵A的第三、四行的元素全改为1，
编写相应的程序。
A=[-3-76-5;-213-1;0013;26-10];
A（[24],[34]）=[-3]
B=[1468;2370;-2-157;11-10];
A（[24],4）=[0]
11.struct函数创建空结构体person，结构体的域包括name、score、rank和class，给name域赋值字符串变量MATLAB，然后删除（参考rmfield命令）class域，编写程序实现。
person=struct（'name','','score','','rank','','class',''）；
person.name='MATLAB'
rmfield（person,'class'）
12.已知多项式
和
，编写程序计算：
1）多项式a（x）+b（x）
p1=[1201];
>>p2=[0101];
>>p1+p2
2）多项式a（x）的根。
roots（p1）
13.已知多项式
和
，编写程序计算：
1）多项式a（x）*b（x）
p1=[02035];
>>p2=[10301];
>>p=conv（p1,p2）
2）多项式a（x）/b（x）的部分分式展开。
[r,p,k]=residue（p1,p2）
14.求解下列方程：
a=[34-7-12;5-742;108-5;-65-210]
b=[4;-3;9;-8]
a\b
15.一天中各个时段温度观测结果如下：
表1室内外温度观测结果（℃）
时间/h
5
7
9
11
13
15
17
室内t1
15.0
18.0
20.0
23.0
27.0
28.0
26.0
室外t2
12.0
15.0
22.0
28.0
32.0
33.0
30.0
试用三次样条（spline）插值法分别求出室内外5：
30-16：
30之间每隔2h，各点的近似温度（℃）。
t=5:
2:
17;
>>y1=[15182023272826];
>>y2=[12152228323330];
>>tt=5.5:
2:
16.5;
>>yy1=spline（t,y1,tt）
yy2=spline（t,y2,tt）
16.求解在x=8时多项式（x-1）（x-2）（x-3）（x-4）的值。
p1=[1-1];
p2=[1,-2];
p3=[1-3];
p4=[1-4];
a=conv（p1,p2）
b=conv（a,p3）
c=conv（b,p4）
polyval（c,8）
17.关于稀疏矩阵，下列哪句话是正确的______C____________。
A.用sparse函数可将一个稀疏矩阵转化为一个满矩阵
B.若A为满矩阵，B为稀疏矩阵，那么A+B仍然为一个稀疏矩阵
C.满矩阵能够参与的运算稀疏矩阵同样能够参与
D.稀疏矩阵是指元素值为0的元素个数占矩阵元素总个数的90%以上
18.当x=－3.2，使用取整函数得出－4，则该取整函数名为______D_________。
A.fixB.roundC.ceilD.floor
19.使用命令LOGSPACE（X1,X2）创建X1=0，X2=4
的行向量，有20个元素，查看其元素分布情况。
（1）写出程序代码；logspace（0,4*pi,20）
logspace（0,4*pi,20）
（2）分析结果中“Columns1through12”为什么都是0.0000？
如何显示正常的结果？
没有选择合适的数值显示格式
File-preference-numberformat-long
20.已知函数
，求该函数在区间
的最小值。
x=-1:
0.01:
3;
p=[10-3];
y=polyval（p,x）
min（y）
21.符号表达式“f=sym（'sin（a*b）+cos（c*d）'）”中的自由符号变量是d___。
22.创建符号表达式
的命令是_______。
symsabtx
a=[2+3,a*b;sin（t）,sqrt（x+1）]
23.用两种方法，求解下列方程：
（1）a=[34-7-12;5-742;108-5;-65-210]
b=[4;-3;9;-8]
a\b
（2）eq1=sym（'3*x+4*y-7*z-12*w=4'）;
eq2=sym（'5*x-7*y+4*z+2*w=-3'）;
eq3=sym（'x+8*z-5*w=9'）;
eq4=sym（'-6*x+5*y-2*z+10*w=-8'）;
s=solve（eq1,eq2,eq3,eq4）
24.请用62位数值解的精度表示‘pi’。
vpa（pi,62）
25.求方程组
关于
的解。
eq1=sym（'u*y^2+vz+w=0'）;
eq2=sym（'y+z+w=0'）;
[y,z]=solve（eq1,eq2）
26.利用符号积分法求
，并得到数值解。
27.symsx
f=1/sqrt（2*pi）*exp（（-x^2）/2）
I=int（f,0,1）
28.计算定积分
的命令，给出计算结果。
f=（abs（x）+sin（x））*x^2
I=int（f,-1,1）
29.求
的3阶导数的命令。
y=atan（（x+3）/（x-2））-log（1+exp（-2*x））
diff（y,3）
30.求符号极限
的命令。
symsxnm
f=tan（m*x）/n*x
limit（f）
31.用符号计算验证三角等式：
sin
（1）cos
（2）-cos
（1）sin
（2）=sin（1-2）
symsphi1phi2
y=sin（phi1）*cos（phi2）-cos（phi1）*sin（phi2）
y1=simple（y）
true=y1==sin（phi1-phi2）
32.在[0,2pi]范围内，用红线画sin（x）,用绿圈画cos（x）。
并加标注x轴：
“自变量X”、y轴：
“函数Y”、图名：
“示意图”，并加格栅。
Grid
x=0:
0.1:
2*pi;
plot（x,sin（x）,'r',x,cos（x）,'g'）
xlabel（'自变量x'）
ylabel（'函数Y'）
title（'示意图'）
grid
33.请用程序画出下图
x=0:
0.1:
2*pi;
plot（x,sin（x）,'b',x,cos（x）,'g'）
text（sin（pi）,'t=pi'）
text（1.5,cos（1.5）,'cos（x）'）
legend（'sin（x）','cos（x）',4）
34.将屏幕分割成四块，并分别画出y=sin（x）,z=cos（x）,a=sin（x）cos（x）,b=sin（x）/cos（x）,并加标注，如下图所示。
x=0:
0.1:
2*pi;
y1=sin（x）;
subplot（2,2,1）
plot（y1）;
title（'sin（x）'）;
y2=cos（x）;
subplot（2,2,2）
plot（y2）;
title（'cos（x）'）;
y3=sin（x）.*cos（x）
subplot（2,2,3）
plot（y3）;
title（'sin（x）*cos（x）'）;
y4=sin（x）./cos（x）
subplot（2,2,4）
plot（y4）;
title（'sin（x）/cos（x）'）;
35.利用plotyy指令在同一张图上绘制
和
在
上的曲线。
x=0:
0.1:
4;
plotyy（x,sin（x）,x,10.^x）
36.用曲面图表示函数
（如下图）。
x=-2:
0.1:
2;
y=x;
[xx,yy]=meshgrid（x,y）;
z=xx.^2+yy.^2;
surf（xx,yy,z）
37.请画出如下3D网线图和3D曲面图，其中
，
x=-4:
0.1:
4;
y=x;
[xx,yy]=meshgrid（x,y）;
z=9*（1-xx）.^2*exp（-xx.^2/2-（yy+1）.^2）;
surf（xx,yy,z）
38.请用符号函数绘图法绘制函数x=sin（3t）cos（t），y=sin（3t）sin（t）的图形，t的变化范围为[0,2]。
ezplot（'sin（3*t）*cos（t）','sin（3*t）*sin（t）',[0,2*pi]）
39.在同一图形窗口绘制
其中，y1以黑色实线标识，y2以蓝色虚线标识，
，给横坐标加上标注“x轴”
t=0:
0.1:
3*pi;
y1=2*exp（-t）.*sin（2*pi*t）;
y2=exp（-2*t）.*cos（2*pi*t+pi/3）;
plot（t,y1,'k-'）
hold
plot（t,y2,'b--'）
40.请找出下列程序中出现错误的行____A_______。
A.x=0:
0.2:
10;y=sqrt（x^2+sqrt（x））;
B.h=plot（x,y,'o',x,y）;
C.set（h,'color','r','linestyle','--'）;
D.xlabel（'时间'）;ylabel（'温度'）;
41.命令plot（x,y,'m:
*'）,其中字符串“m:
*”中,第一个字符“m”____B________。
A.表示曲线为m符号构成的线
B.表示曲线颜色为紫色
C.表示曲线x坐标的最大值为m变量的值
D.表示曲线y坐标的最大值为m变量的值
42.编写实现比较a，b两个值大小功能的函数：
result=ex01（a,b），a,b中的大值赋值给result；若a=b，则输出result“两值相等”并同时屏幕提示“请输入两个大小不同的数值”，要求写出完整的函数体。
functionresult=ex01（a,b）
if（a>b）
reult=a
elseifaresult=b
else
reult='两值相等'
end
43.Fibonacci数列各元素为：
1、1、2、3、5、8、…，满足以下关系：
F1=1，F2=1，Fn=Fn-1+Fn-2，用M函数文件实现：
找到Fibonacci数组中第一个大于10000的元素，即n=？
时，Fn>10000.
%ex01
F1=1;F2=1;
whileFn<10000
Fn=Fn-1+Fn-2
n=n+1
end
44.
45.A=randn（4，3）>0.5,编程实现查找矩阵A的每行中第一个0元素所在的列。
将结果存放着一个列向量中。
升级会员 