太原理工大学现代科技学院MATLAB 矩阵和数组的操作实验报告文档格式.docx

1、4,5,6;7,8,9 A = 3 2 1 4 5 6 7 8 9 （2）直接利用MATLAB提供的函数创建一个3*3矩阵。 如在命令区输入rand（3,3）即得到一个3*3的随机矩阵 B=rand（3,3） B = 0.9501 0.4860 0.4565 0.2311 0.8913 0.0185 0.6068 0.7621 0.8214 （3）利用MATLAB提供给的“Matrix Editor”完成输入。 步骤1 在命令区输入A=1. 步骤2 用鼠标单击工具栏的工作区浏览器，MATLAB弹出变量浏览器，选中变量A，鼠标左键双击A,打开矩阵浏览器。 步骤3 在左下脚的两个文本框中分别输入希

2、望得到的矩阵的行数和列数：3行3列，即得到一个3*3矩阵。 步骤4 要将上面矩阵改为一个45矩阵，只需改变矩阵的行数和列数即可。如将3行3列改为4行5列，即可得到一个4*5矩阵。若想修改其中的元素，只需用鼠标选中表格中我们想要修改的元素，将原来的元素修改为我需要的值。 步骤5 要命令区输入save data A（data 为我们变量文件起的名称，系统会自动沿设定好的路径以“.mat”格式存储文件），即可保存在上面例子中创建的矩阵A. 步骤6 在命令区 输入load data 即可保存在文件中的矩阵读到MATLAB的命令工作区的内存中来。 C=1 C = 1 C 1 0 0 0 0 0 1 0

3、2 0 3 0 4 0 5 0 6 0 7 0 8 0 0 0 0 0 save data C load data 2. 建立一个等差数列，然后由它产生一个对角阵。 步骤1 在命令区输入a=linspace（0,1.5,5）产生一个等差数列。 步骤2 在命令区输入B=diag（a）产生一个对角阵。 a=linspace（0,1.5,5） a = 0 0.3750 0.7500 1.1250 1.5000 B=diag（a） 0 0.3750 0 0 0 0 0 0.7500 0 0 0 0 0 1.1250 0 0 0 0 0 1.5000 3. 利用MATLAB的函数inv（A）求方阵A的逆

4、矩阵。 步骤1 在命令区输入A=1，2；5，6得到一个22 的方阵。 步骤2 再输入B=inv（A）求出A的逆矩阵。 A=1,2;5,6 1 2 5 6 B=inv（A） -1.5000 0.5000 1.2500 -0.2500 四、练习 1. 创建一个5*5矩阵，提取住对角线以上的部分。 A=rand（5,5） 0.4447 0.1763 0.8936 0.1389 0.1988 0.6154 0.4057 0.0579 0.2028 0.0153 0.7919 0.9355 0.3529 0.1987 0.7468 0.9218 0.9169 0.8132 0.6038 0.4451 0

5、.7382 0.4103 0.0099 0.2722 0.9318 B=triu（A） 0 0.4057 0.0579 0.2028 0.0153 0 0 0.3529 0.1987 0.7468 0 0 0 0.6038 0.4451 0 0 0 0 0.9318 2. A=rand（3），B=magic（3），C=rand（3,4），计算A*B*C A=rand（3）,B=magic（3）,C=rand（3,4）,D=A*B*C 0.6979 0.8537 0.8998 0.3784 0.5936 0.8216 0.8600 0.4966 0.6449 8 1 6 3 5 7 4 9 2

6、0.8180 0.2897 0.7271 0.5681 0.6602 0.3412 0.3093 0.3704 0.3420 0.5341 0.8385 0.7027 D = 22.3218 14.2486 22.6098 19.9167 16.4714 10.3190 15.9728 14.2467 18.3897 11.5943 19.1132 16.5834 3. 创建一个3*3矩阵，并求其转置，逆矩阵。 A=rand（3） 0.5466 0.6213 0.5226 0.4449 0.7948 0.8801 0.6946 0.9568 0.1730 B=A,C=inv（A） 0.5466

7、 0.4449 0.6946 0.6213 0.7948 0.9568 0.5226 0.8801 0.1730 5.9121 -3.2936 -1.1030 -4.4832 2.2521 2.0854 1.0603 0.7672 -1.3257 4. 用两种方法求Ax=b的解（A为4阶随机矩阵，b为4阶列向量）。 A=rand（4）,b=（1:4）,x=bA 0.9601 0.2679 0.2126 0.2071 0.7266 0.4399 0.8392 0.6072 0.4120 0.9334 0.6288 0.6299 0.7446 0.6833 0.1338 0.3705 b = 2

8、3 4 x = 0.2209 0.2227 0.1437 0.1598 ,linsolve（A,b） 0.8729 0.6649 0.8188 0.6873 0.2379 0.8704 0.4302 0.3461 0.6458 0.0099 0.8903 0.1660 0.9669 0.1370 0.7349 0.1556 ans = 141104154301010653655869722617759472290645667717336001371848048640/5039983259460903602964471145710749501860647252854259052213333883

9、9 188848975313652446796351161733576614585665265865775181025314340864/50399832594609036029644711457107495018606472528542590522133338839 153214171780514282408806451432253555713142265086291465303018700800/50399832594609036029644711457107495018606472528542590522133338839 -4710739571658066640612625989302

10、21063180088644804283425613542326272/50399832594609036029644711457107495018606472528542590522133338839 5. 创建一个4阶随机矩阵A，计算A3。 A=rand（4）,B=A3 0.4508 0.2548 0.9084 0.0784 0.7159 0.8656 0.2319 0.6408 0.8928 0.2324 0.2393 0.1909 0.2731 0.8049 0.0498 0.8439 1.7209 1.3508 1.4347 0.9875 2.8362 2.8175 1.8353 2

11、.2541 1.7809 1.3238 1.1502 1.0112 2.2622 2.6361 1.4024 2.1656 6. 求 100999 之间被21整除的数的个数。 （提示：先利用冒号表达式，再利用find和length函数。） A=100:999; B=find（mod（A,21）=0）; C=length（B） 43 7. 设有矩阵A和B （1）求它们的乘积C=AB.。 （2）将矩阵C的右下角32子矩阵赋给D。 A=1:5;6:10;11:15;16:20;21:25,B=3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 3 0 16 17 -6 9 0 23 -4 9 7 0 4 13 11 C=A*B 93 150 77 258 335 237 423 520 397 588 705 557 753 890 717 D=C（3:5,2:3） 520 397 705 557 890 717