MATLAB习题及答案.docx
《MATLAB习题及答案.docx》由会员分享,可在线阅读,更多相关《MATLAB习题及答案.docx(14页珍藏版)》请在冰豆网上搜索。
MATLAB习题及答案
习题:
一、填空
1.清除命令窗口内容的命令是(clc)。
2.删除工作空间中保存的变量x的命令是(clearx)。
3.将双精度实数的显示格式设置成15位定点小数格式的命令是(formatlong)。
4.要在命令窗口中查看对函数log的帮助信息,所用命令是(helplog)。
5.要查询在工作空间中定义的变量x的相关信息,可用命令(whosx)。
6.设x是复数,取x的实部存入变量y的语句是(y=real(x))。
7.设x是复数,取x的虚部存入变量y的语句是(y=imag(x))。
8.设x是复数,取x的模存入变量y的语句是(y=abs(x))。
9.设x是复数,取x的共轭存入变量y的语句是(y=conj(x))。
10.表达式(5<2)*120的值是(0)。
11.表达式(5>2)*(6~=5)的值是
(1)。
12.绘制极坐标图形的Matlab库函数名是(polar)。
13.将图题设置成“Fig.1示意图”的语句是(title('Fig.1示意图'))。
14.将横坐标轴标签设置成“时间(秒)”的语句是(xlabel('时间(秒)'))。
15.设置图例的Matlab库函数名是(legend)。
16.绘制三维线图的Matlab库函数名是(plot3)。
17.绘制三维网格图的Matlab库函数名是(mesh)。
18.绘制三维表面图的Matlab库函数名是(surf)。
19.关闭网格显示的Matlab语句是(gridoff)。
20.关闭坐标架显示的Matlab语句是(axisoff)。
21.使得各坐标具有相同刻度间隔的Matlab语(axis('equal'))。
22.将三维图形视角设置成方位角60度、仰角45度的语句是(view(60,45))。
23.绘制二维等高线图的Matlab库函数名是(surfc)。
24.在同一图形窗口中开多个子窗口的Matlab库函数名是(subplot)。
二、基本操作
1,计算
与
的数组乘积。
2,对于
,如果
,
,求解X。
3,已知:
,分别计算a的数组平方和矩阵平方,并观察其结果。
4,角度
,求x的正弦、余弦、正切和余切。
(应用sin,cos,tan.cot)
5,将矩阵
、
和
组合成两个新矩阵:
(1)组合成一个43的矩阵,第一列为按列顺序排列的a矩阵元素,第二列为按列顺序排列的b矩阵元素,第三列为按列顺序排列的c矩阵元素,即
(2)按照a、b、c的列顺序组合成一个行矢量,即
6,将(x-6)(x-3)(x-8)展开为系数多项式的形式。
(应用poly,polyvalm)
7,求解多项式x3-7x2+2x+40的根。
(应用roots)
8,求解在x=8时多项式(x-1)(x-2)(x-3)(x-4)的值。
(应用poly,polyvalm)
9,计算多项式
的微分和积分。
(应用polyder,polyint,poly2sym)
10,解方程组
。
(应用x=a\b)
11,矩阵
,计算a的行列式和逆矩阵。
(应用det,inv)
12,用符号计算验证三角等式:
(应用syms,simple)
sin
(1)cos
(2)-cos
(1)sin
(2)=sin(1-2)
13,因式分解:
(应用syms,factor)
14,
,用符号微分求df/dx。
(应用syms,diff)
15,符号函数绘图法绘制函数x=sin(3t)cos(t),y=sin(3t)sin(t)的图形,t的变化范围为[0,2]。
(应用syms,ezplot)
16,绘制曲线
,x的取值范围为[-5,5]。
(应用plot)
17,有一组测量数据满足
,t的变化范围为0~10,用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线,在图中添加标题
,并用箭头线标识出各曲线a的取值,并添加标题
和图例框。
(应用plot,title,text,legend)
18,表中列出了4个观测点的6次测量数据,将数据绘制成为分组形式和堆叠形式的条形图。
第1次
第2次
第3次
第4次
第5次
第6次
观测点1
3
6
7
4
2
8
观测点2
6
7
3
2
4
7
观测点3
9
7
2
5
8
4
观测点4
6
4
3
2
7
4
19,x=[6649715638],绘制饼图,并将第五个切块分离出来。
20,用sphere函数产生球表面坐标,绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。
(应用sphere,mesh,hiddenoff,surf,NaN)
答案:
1,计算
与
的数组乘积。
>>a=[693;275];
>>b=[241;468];
>>a.*b
ans=
12363
84240
2,对于
,如果
,
,求解X。
>>A=[492;764;357];
>>B=[372628]’;
>>X=A\B
X=
-0.5118
4.0427
1.3318
3,已知:
,分别计算a的数组平方和矩阵平方,并观察其结果。
>>a=[123;456;789];
>>a.^2
ans=
149
162536
496481
>>a^2
ans=
303642
668196
102126150
4,角度
,求x的正弦、余弦、正切和余切。
>>x=[304560];
>>x1=x/180*pi;
>>sin(x1)
ans=
0.50000.70710.8660
>>cos(x1)
ans=
0.86600.70710.5000
>>tan(x1)
ans=
0.57741.00001.7321
>>cot(x1)
ans=
1.73211.00000.5774
5,将矩阵
、
和
组合成两个新矩阵:
(1)组合成一个43的矩阵,第一列为按列顺序排列的a矩阵元素,第二列为按列顺序排列的b矩阵元素,第三列为按列顺序排列的c矩阵元素,即
(2)按照a、b、c的列顺序组合成一个行矢量,即
>>a=[42;57];
>>b=[71;83];
>>c=[59;62];
%
(1)
>>d=[a(:
)b(:
)c(:
)]
d=
475
586
219
732
%
(2)
>>e=[a(:
);b(:
);c(:
)]'
e=
452778135692
或利用
(1)中产生的d
>>e=reshape(d,1,12)
ans=
452778135692
6,将(x-6)(x-3)(x-8)展开为系数多项式的形式。
>>a=[638];
>>pa=poly(a);
>>ppa=poly2sym(pa)
ppa=
x^3-17*x^2+90*x-144
7,求解多项式x3-7x2+2x+40的根。
>>r=[1-7240];
>>p=roots(r);
-0.2151
0.4459
0.7949
0.2707
8,求解在x=8时多项式(x-1)(x-2)(x-3)(x-4)的值。
>>p=poly([1234]);
>>polyvalm(p,8)
ans=
840
9,计算多项式
的微分和积分。
>>p=[4–12–145];
>>pder=polyder(p);
>>pders=poly2sym(pder)
>>pint=polyint(p);
>>pints=poly2sym(pint)
pders=
12*x^2-24*x-14
pints=
x^4-4*x^3-7*x^2+5*x
10,解方程组
。
>>a=[290;3411;226];
>>b=[1366]';
>>x=a\b
x=
7.4000
-0.2000
-1.4000
11,矩阵
,计算a的行列式和逆矩阵。
>>a=[42-6;754;349];
>>ad=det(a)
>>ai=inv(a)
ad=
-64
ai=
-0.45310.6562-0.5937
0.7969-0.84370.9062
-0.20310.1562-0.0937
12,用符号计算验证三角等式:
(应用syms,simple)
sin
(1)cos
(2)-cos
(1)sin
(2)=sin(1-2)
>>symsphi1phi2;
>>y=simple(sin(phi1)*cos(phi2)-cos(phi1)*sin(phi2))
y=
sin(phi1-phi2)
13,因式分解:
>>symsx;
>>f=x^4-5*x^3+5*x^2+5*x-6;
>>factor(f)
ans=
(x-1)*(x-2)*(x-3)*(x+1)
14,
,用符号微分求df/dx。
(应用syms,diff)
>>symsax;
>>f=[a,x^2,1/x;exp(a*x),log(x),sin(x)];
>>df=diff(f)
df=
[0,2*x,-1/x^2]
[a*exp(a*x),1/x,cos(x)]
15,符号函数绘图法绘制函数x=sin(3t)cos(t),y=sin(3t)sin(t)的图形,t的变化范围为[0,2]。
>>symst
>>ezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),[0,pi])
16,绘制曲线
,x的取值范围为[-5,5]。
>>x=-5:
0.2:
5;
>>y=x.^3+x+1;
>>plot(x,y)
17,有一组测量数据满足
,t的变化范围为0~10,用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线,在图中添加标题
,并用箭头线标识出各曲线a的取值,并添加标题
和图例框。
>>t=0:
0.5:
10;
>>y1=exp(-0.1*t);
>>y2=exp(-0.2*t);
>>y3=exp(-0.5*t);
>>plot(t,y1,'-ob',t,y2,':
*r',t,y3,'-.^g')
>>title('\ity\rm=e^{-\itat}')
>>title('\ity\rm=e^{-\itat}','FontSize',12)
>>text(t(6),y1(6),'\leftarrow\ita\rm=0.1','FontSize',11)
>>text(t(6),y2(6),'\leftarrow\ita\rm=0.2','FontSize',11)
>>text(t(6),y3(6),'\leftarrow\ita\rm=0.5','FontSize',11)
>>title('\ity\rm=e^{-\itat}','FontSize',12)
>>legend('a=0.1','a=0.2','a=0.5')
18,表中列出了4个观测点的6次测量数据,将数据绘制成为分组形式和堆叠形式的条形图。
第1次
第2次
第3次
第4次
第5次
第6次
观测点1
3
6
7
4
2
8
观测点2
6
7
3
2
4
7
观测点3
9
7
2
5
8
4
观测点4
6
4
3
2
7
4
>>y=[3696;6774;7323;4252;2487;8744];
>>bar(y)
19,x=[6649715638],绘制饼图,并将第五个切块分离出来。
>>x=[6649715638];
>>L=[00001];
>>pie(x,L)
20用sphere函数产生球表面坐标,绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。
>>[x,y,z]=sphere(30);
>>mesh(x,y,z)
>>mesh(x,y,z),hiddenoff
>>surf(x,y,z)
>>z(18:
30,1:
5)=NaN*ones(13,5);
>>surf(x,y,z)