实验二Matlab矩阵的初等运算与其答案.docx
- 文档编号:11317351
- 上传时间:2023-02-26
- 格式:DOCX
- 页数:21
- 大小:43.82KB
实验二Matlab矩阵的初等运算与其答案.docx
《实验二Matlab矩阵的初等运算与其答案.docx》由会员分享,可在线阅读,更多相关《实验二Matlab矩阵的初等运算与其答案.docx(21页珍藏版)》请在冰豆网上搜索。
实验二Matlab矩阵的初等运算与其答案
实验二Matlab矩阵的初等运算
实验目的:
掌握Matlab的运算方法
实验内容:
在Matlab命令窗口输入:
H1=ones(3,2)H2=zeros(2,3)H3=eye(4)
观察以上各输入结果,并在每式的后面标注其含义。
>>formatcompact
>>H1=ones(3,2),disp('3行2列的全1矩阵')
H1=
11
11
11
3行2列的全1矩阵
>>H2=zeros(2,3),disp('2行3列的全零矩阵')
H2=
000
000
2行3列的全零矩阵
>>H3=eye(4),disp('4阶的单位矩阵')
H3=
1000
0100
0010
0001
4阶的单位矩阵
已知
,
,
,
试把这四个矩阵组合为一个大矩阵,看看有几种组合方式?
8
>>formatcompact
>>Q=[123;456];P=[789];R=[1;0];S=3;
>>[Q,R;P,S]
ans=
1231
4560
7893
>>[R,Q;P,S]
ans=
1123
0456
7893
>>[Q,R;S,P]
ans=
1231
4560
3789
>>[R,Q;S,P]
ans=
1123
0456
3789
>>[S,P;R,Q]
ans=
3789
1123
0456
>>[S,P;Q,R]
ans=
3789
1231
4560
>>[P,S;R,Q]
ans=
7893
1123
0456
>>[P,S;Q,R]
ans=
7893
1231
4560
建立一个字符串向量,删除其中的大写字母。
提示:
利用find函数和空矩阵。
>>a='ABCDefgijKLMN123'
a=
ABCDefgijKLMN123
>>k=find(a>='A'&a<='Z')
k=
123410111213
>>a(k)=[]
a=
efgij123
在命令窗中分别输入who和whos,观察检查结果是否与所得结果相符。
>>who
Yourvariablesare:
H1H2H3PQRSaansk
>>whos
NameSizeBytesClassAttributes
H13x248double
H22x348double
H34x4128double
P1x324double
Q2x348double
R2x116double
S1x18double
a1x816char
ans3x496double
k1x864double
已知矩阵
,计算以下表达式的结果,体会*,^,sqrtm,expm与.*,.^,sqrt,exp的区别。
(1)r1=S^2,r2=2.^S,r3=S.^2,
(2)u1=sqrtm(S),v1=u1*u1
(3)u2=sqrt(S),v2=u2.*u2
(4)u3=expm(S),v3=logm(u3)
(5)u4=exp(S),v4=log(u4)
>>S=[13125;4709;7162;82113]
S=
13125
4709
7162
82113
>>r1=S^2
r1=
1374613971
10479147110
693814262
1175519589
>>r2=2.^S
r2=
28409632
161281512
1282644
256420488
>>r3=S.^2
r3=
1914425
1649081
491364
6441219
>>u1=sqrtm(S)
u1=
+---
+++-
-+++
-+-+
>>v1=u1*u1
v1=
-+-+
+--+
+-+-
+--
>>u2=sqrt(S)
u2=
0
>>u3=expm(S)
u3=
+008*
>>v3=logm(u3)
v3=
>>u4=exp(S)
u4=
+005*
>>v4=log(u4)
v4=
13125
4709
7162
82113
对中定义的矩阵S完成以下变换,输出变换后的矩阵:
(1)将矩阵S上下翻转
(2)将矩阵S左右翻转
(3)将矩阵S重组为一个2行8列的矩阵
(4)将矩阵S整体逆时针旋转90°
(5)提取矩阵S对角线上的元素
(6)建立一个对角阵T,对角线上的元素为S对角线上的元素,其余元素为0
(7)取出矩阵S的左下三角部分
(8)取出矩阵S的右上三角部分
(9)把矩阵S的元素按列取出排成一行
>>S=[13125;4709;7162;82113]
S=
13125
4709
7162
82113
>>flipud(S)
ans=
82113
7162
4709
13125
>>fliplr(S)
ans=
51231
9074
2617
31128
>>reshape(S,2,8)
ans=
173112652
487201193
>>rot90(S)
ans=
5923
120611
3712
1478
>>diag(S)
ans=
1
7
6
3
>>T=diag(S)
T=
1
7
6
3
>>T=diag(T)
T=
1000
0700
0060
0003
>>tril(S)
ans=
1000
4700
7160
82113
>>triu(S)
ans=
13125
0709
0062
0003
>>S(:
)'
ans=
147837121206115923
已知矩阵A=[135],B=[246]
(1)求C=A+B,D=A-2,E=B-A
>>formatcompact
>>A=[135],B=[246]
A=
135
B=
246
>>C=A+B,D=A-2,E=B-A
C=
3711
D=
-113
E=
111
>>%都是按元素群运算
(2)求F1=A*3,F2=A.*B,F3=A./B,F4=A.\B,F5=B.\A,F6=B.^A,F7=2./B,F8=B.\2.
>>F1=A*3
F1=
3915
>>%按矩阵运算
>>F2=A.*B
F2=
21230
>>%按元素群运算A点乘B
>>F3=A./B
F3=
>>%按元素群运算A各个元素右除B
>>F4=A.\B
F4=
>>%按元素群运算A各个元素左除B
>>F5=B.\A
F5=
>>%按元素群运算B各个元素左除A
>>F6=B.^A
F6=
2647776
>>%按元素群运算B各个元素的A次幂
>>F7=2./B
F7=
>>%按元素群运算2右除B
>>F8=B.\2
F8=
>>%按元素群运算B各个元素左除2
(3)求Z1=A*B’,Z2=B’*A观察以上各输出结果,比较各种运算的区别,并在每式的后面标注其含义。
>>Z1=A*B'
Z1=
44
>>%按矩阵运算A乘以B的转置
>>Z2=B'*A
Z2=
2610
41220
61830
>>%按元素群运算B的转置乘以A
已知矩阵
,
,求H1=I*J’,H2=I’*J,H3=I.*J并求它们的逆阵。
(1)H1=I*J’
>>formatcompact
>>I=[14813;-36-59;2-7-12-8]
I=
14813
-36-59
2-7-12-8
>>J=[543-2;6-23-8;-13-97]
J=
543-2
6-23-8
-13-97
>>H1=I*J'
H1=
19-8230
-24-117129
-385429
>>det(H1)
ans=
+004
>>inv(H1)
ans=
(2)H2=I’*J
>>H2=I'*J
H2=
-1516-2436
63-1793-105
226117-60
12710138-154
>>det(H2)
ans=
因为H2的行列式接近于零,H2没有逆矩阵
(3)H3=I.*J
>>H3=I.*J
H3=
51624-26
-18-12-15-72
-2-21108-56
>>pinv(H3)
ans=
因为H3为一个非满秩的矩阵,所以H3没有逆矩阵,而有一个伪逆矩阵.
已知矩阵
(1)用指令检查C的阶数;
>>formatcompact
>>C=[;;]
C=
0
>>[m,n]=size(C)
m=
3
n=
4
矩阵C是一个3行4列的矩阵.
(2)找出C(2,3)的值;
>>C(2,3)
ans=
(3)找出值为的元素的下标。
>>[a,b]=find(C==
a=
2
3
1
b=
1
2
4
a为行下标,b为列下标。
求下列表达式的值。
(1)
,注意:
Matlab当中三角函数的运算按弧度进行。
>>Z1=2*sin(85*pi/180)/(1+exp
(2))
Z1=
(2)
,其中
>>formatcompact
>>x=[21+2i;5]
x=
+
>>Z2=1/2*log(x+sqrt(1+x))
Z2=
+
(3)
,
>>formatcompact
>>a=:
:
a=
Columns1through9
Columns10through18
Columns19through27
Columns28through36
0
Columns37through45
Columns46through54
Columns55through61
>>Z3=(exp.*a)-exp.*a))/2.*sin(a+
Z3=
Columns1through9
Columns10through18
Columns19through27
Columns28through36
0
Columns37through45
Columns46through54
Columns55through61
求下列联立方程的解
3x+4y-7z-12w=4
5x-7y+4z+12w=-3
x+8z-5w=9
-6x+5y-2z+10w=-8
解:
令A=[34-7-12;5-7412;108-5;-65-210]
X=[x;y;z;w]B=[4;-3;9;-8]
X=A\B
>>formatcompact
>>A=[34-7-12;5-7412;108-5;-65-210];B=[4;-3;9;-8];X=A\B
X=
(1)列写2×2阶的单位矩阵I,4×4阶的魔方矩阵M和4×2阶的全幺矩阵A,全零矩阵B。
>>formatcompact
>>I=eye
(2),A=ones(4,2),B=zeros(4,2),M=magic(4)
I=
10
01
A=
11
11
11
11
B=
00
00
00
00
M=
162313
511108
97612
414151
(2)将这些矩阵拼接为6×6阶的矩阵C:
>>formatcompact
>>I=eye
(2);A=ones(4,2);B=zeros(4,2);M=magic(4);C=[I,A';B,M]
C=
101111
011111
00162313
00511108
0097612
00414151
(3)取出C的第2,4,6行,组成3×6阶的矩阵C1,取出第2,,4,6列,组成6×3阶的矩阵C2。
>>formatcompact
>>C1=C([2,4,6,],:
),C2=C(:
[2,4,6,])
C1=
011111
00511108
00414151
C2=
011
111
0213
0118
0712
0141
(4)求D=C1C2及D1=C2C1。
>>formatcompact
>>D=C1*C2,D1=C2*C1
D=
13535
0313281
0281345
D1=
00925259
0110262610
006220421529
008723323096
008324525068
0074168155113
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 Matlab 矩阵 初等 运算 与其 答案