MATLAB实验报告模板xin.docx
《MATLAB实验报告模板xin.docx》由会员分享,可在线阅读,更多相关《MATLAB实验报告模板xin.docx(17页珍藏版)》请在冰豆网上搜索。
MATLAB实验报告模板xin
《数学软件选讲》
——MATLAB
实验报告
班级:
信计2班
学号:
080350215
姓名:
陈龙
河北工程大学理学院
2011年11月28日
实验一MATLAB编程环境
一、实验内容
二、实验结果及说明
三、实验代码
实验二
一、实验内容
(1)设A和B是两个同维同大小的矩阵,问:
1)A*B和A.*B的值是否相等?
2)A./B和B.\A的值是否相等?
3)A/B和B\A的值是否相等?
4)A/B和B\A所代表的数学含义是什么?
(2)写出完成下列操作的命令。
1)将矩阵A第2—5行中第1,3,5列元素赋给矩阵B。
2)删除矩阵A的第7号元素。
3)将矩阵A的每个元素值加30。
4)求矩阵A的大小和维数。
5)将向量t的0元素用机器零来代替。
6)将含有12个元素的向量x转换成3*4矩阵。
7)求一个字符串的ASCII码。
8)求一个ASCII码所对应的字符。
(3)下列命令执行后,L1、L2、L3、L4的值分别是多少/
A=1:
9;B=10-A;
L1=A==B;
L2=A<=5;
L3=A>3&A<7;
L4=find(A>3&A<7);
二、实验结果及说明
1、1)不相等
2)不相等
3)不相等
4)矩阵左除A\B表示由Gauss消去法获得线性方程AX=B的解X。
矩阵右除B/A表示线性方程XA=B的解X。
对于矩阵来说,左除和右除表示两种不同的除数矩阵和被除数矩阵的关系。
2、A=
123101
456112
71009125
11116
45688
1)B=
462
795
116
468
2)A=
Columns1through16
14714210015369161011
Columns17through24
121812568
3)A=
3132334031
3435364132
37130394235
3131313136
3435363838
4)ans=5
ans=55
5)t=1203
r=3
t=1.00002.00000.00003.0000
6)x=24681012141618202224
X=
281420
4101622
6121824
7)d=happy
ans=10497112112121
8)h=9
ans=a
3、A=123456789
B=987654321
L1=000010000
L2=111110000
L3=000111000
L4=456
三、实验代码
a=[123;456;789]
a(6)=100
b=[101112]'
a=[a,b]
c=[1111]
a=[a;c]
A=a
B=magic(5)
e=[4568]
A=[A;e]
d=[12568]'
A=[A,d]
A*B
A.*B
A./B
B./A
A/B
B/A
B=A(2:
1:
5,1:
2:
5)
A(7)=[]
a=[123;456;789]
a(6)=100
b=[101112]'
a=[a,b]
c=[1111]
a=[a;c]
A=a
B=magic(5)
e=[4568]
A=[A;e]
d=[12568]'
A=[A,d]
A=A+30
a=[123;456;789]
a(6)=100
b=[101112]'
a=[a,b]
c=[1111]
a=[a;c]
A=a
B=magic(5)
e=[4568]
A=[A;e]
d=[12568]'
A=[A,d]
length(A)
size(A)
t=[1203]
r=find(t==0)
t(r)=eps
x=[2:
2:
24]
X=reshape(x,3,4)
d='happy'
abs(d)
h=97
char(h)
A=1:
9
B=10-A
L1=A==B
L2=A<=5
L3=A>3&A<7
L4=find(A>3&A<7)
、
实验三矩阵基本运算
(二)
一、实验内容
(1)已知
,
求下列表达式的值:
1)A+6B和A2-B+I(I为单位矩阵)
2)A*B,A.*B和B*A
3)A/B和B\A
4)[A,B]和[A([1,3],:
);B^2]
(2)已知
,取出其前三行构成矩阵B,其前两列构成矩阵C,其右下角3×2子矩阵构成矩阵D,B与C的乘积构成矩阵E,分别求E=10&A<25)
二、实验结果及说明
A=
-15-4
078
3617
B=
83-1
253
-320
ans=
4723-10
123726
-15737
ans=
-18-21717
22533109
21867526
ans=
141416
-105121
125328180
ans=
-8154
03524
-91220
ans=
-110-15
722853
3-128
ans=
1.2234-0.92552.9787
-0.94682.3511-0.9574
4.61703.872313.8936
ans=
-5.3173-3.56610.8942
-1.2768-0.47430.2411
2.12291.3173-0.2924
ans=-1012345678
ans=
-15-4
078
3617
73371
173713
-2019
A=
23.000010.0000-0.77800
41.0000-45.000065.00005.0000
32.00005.0000032.0000
6.0000-9.540054.00003.1400
B=
23.000010.0000-0.77800
41.0000-45.000065.00005.0000
32.00005.0000032.0000
C=
23.000010.0000
41.0000-45.0000
32.00005.0000
6.0000-9.5400
D=
65.00005.0000
032.0000
54.00003.1400
E=
1.0e+003*
0.9141-0.2239
1.20802.7123
1.1330-0.2103
ans=
01
00
01
ans=
11
01
11
ans=
11
11
11
ans=
00
10
00
ans=1
5
三、实验代码
A=[-15-4;078;3617]
B=[83-1;253;-320]
A+6*B
A*A-B+eye(3)
A*B
A.*B
B*A
A/B
B/A
[A:
B]
[A([1:
3],:
);B^2]
A=[2310-0.7780;41-45655;325032;6-9.54543.14]
B=A(1:
1:
3,1:
1:
4)
C=A(1:
1:
4,1:
1:
2)
D=A(2:
1:
4,3:
1:
4)
E=B*C
EE&D
E|D
~E|~D
find(A>=10&A<25)
实验四矩阵分析
一、实验内容
(1)使用函数,实现方阵左旋90°或右旋90°的功能。
例如,原矩阵为A,A左旋后得到B,右旋后得到C。
,
,
(2)建立一个方阵A,求A的逆矩阵和A的行列式的值,并验证A与A-1是互逆的。
(3)求下列矩阵的主对角线元素、上三角阵、下三角阵、秩和迹。
1)
2)
(4)求矩阵A的特征值和相应的特征向量。
、二、实验结果及说明
A=
14710
25811
36912
B=
101112
789
456
123
C=
321
654
987
121110
A=
1-123
51-42
3052
111509
B=
-0.17580.16410.2016-0.0227
-0.1055-0.1016-0.03910.0664
-0.0508-0.08590.15160.0023
0.3906-0.0313-0.18130.0281
ans=
1.0000-0.0000-0.00000.0000
0.00001.00000.0000-0.0000
-0.0000-0.00001.0000-0.0000
0.0000-0.00000.00001.0000
ans=
1.0000-0.00000.0000-0.0000
0.00001.00000.0000-0.0000
0.00000.00001.0000-0.0000
-0.0000-0.00000.00001.0000
C=
1280
A=
1-123
51-42
3052
111509
D=
1
1
5
9
ans=
1-123
01-42
0052
0009
ans=
1000
5100
3050
111509
ans=
4
ans=
16
B=
0.430043.00002.0000
-8.90004.000021.0000
ans=
0.4300
4.0000
ans=
0.430043.00002.0000
04.000021.0000
ans=
0.430000
-8.90004.00000
ans=
2
ans=
4.4300
A=
1.00001.00000.5000
1.00001.00000.2500
0.50000.25002.0000
V=
0.72120.44430.5315
-0.68630.56210.4615
-0.0937-0.69760.7103
D=
-0.016600
01.48010
002.5365
三、实验代码
A=[14710;25811;36912]
B=rot90(A)
C=rot90(A,3)
A=[1-123;51-42;3052;111509]
B=inv(A)
A*B
B*A
C=det(A)
A=[1-123;51-42;3052;111509]
D=diag(A)
triu(A)
tril(A)
rank(A)
trace(A)
B=[0.43432;-8.9421]
diag(B)
triu(B)
tril(B)
rank(B)
trace(B)
A=[110.5;110.25;0.50.252]
[V,D]=eig(A)
实验五程序控制结构设计