实验2-MATLAB的数值计算一 实验记录.doc

西华大学实验报告第 组 西华大学实验报告（理工类）开课学院及实验室： 实验时间 ： 年 月 日学 生 姓 名学 号成 绩学生所在学院年级/专业/班课 程 名 称课 程 代 码实验项目名称实验二：MATLAB的数值计算（一）项 目 代 码指 导 教 师项 目 学 分一、实验目的1、学会如何创建、使用和保存MATLAB的矩阵和数组。2、学会矩阵和数组的基本操作。二、实验原理1、 MATLAB用于产生常用特殊矩阵的函数见教材。实验中用到的产生单位阵的函数为eye(m,n),ones(m,n),rand(m,n), m, n 为矩阵的行和列。用函数diag(v)可建立以向量v为主对角线元素的对角阵。2、 程序生成向量:生成有规律元素的向量，可用冒号运算符；线性等分向量，可用函数y=linspace(x1,x2,n)；对数等分向量，可用函数y=logspace(x1,x2,n)。3、 MATLAB通过确认矩阵下标，可以对矩阵进行插入子块、提取子块和重排子块的操作。A(m,n)：提取第m行，第n列元素；A(:,n)：提取第n列元素；A(m,:)：提取第m行元素；A(m1:m2,n1:n2)：提取第m1行到第m2行和第n1列到第n2列的所有元素（提取子块）。4、矩阵扩展：如果在原矩阵中一个不存在的地址位置上设定一个数（赋值），则该矩阵会自动扩展行列数，并在该位置上添加这个数，而且在其他没有指定的位置补零。三、实验设备、仪器及材料高性能计算机，MATLAB软件。四、实验步骤（按照实际操作过程）1、 输入元素为从0到20，步距为2的一个向量P。 p=0:2:20p = 0 2 4 6 8 10 12 14 16 18 202、 生成一个元素从1到100，线性等分的(1*10)维的行向量。 a=1:100a = Columns 1 through 16 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Columns 17 through 32 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 Columns 33 through 48 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 Columns 49 through 64 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 Columns 65 through 80 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 Columns 81 through 96 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 Columns 97 through 100 97 98 99 100 a=linspace(1,100,10)a = 1 12 23 34 45 56 67 78 89 1003、 输入矩阵，并将矩阵的2行3列元素改为10。 M=1 2 3 4;3 2 1 0M = 1 2 3 4 3 2 1 0 M6=10? M6=10 |Error: Unbalanced or misused parentheses or brackets. M(2,3)=10M = 1 2 3 4 3 2 10 0 M(6)=10M = 1 2 3 4 3 2 10 04、 产生一个对角阵，主对角线元素分别为1、2、3、4。 b=1 0 0 0;0 2 0 0;0 0 3 0;0 0 0 4b = 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0 4 diag(b)ans = 1 2 3 45、 产生一个4阶随机矩阵R，令矩阵R1等于矩阵R的第3行，令矩阵R2等于矩阵R的第1、3列。 R=rand(4,4)R = 0.8132 0.1987 0.0153 0.4660 0.0099 0.6038 0.7468 0.4186 0.1389 0.2722 0.4451 0.8462 0.2028 0.1988 0.9318 0.5252 R1=R(3,:)R1 =0.1389 0.2722 0.4451 0.8462 R2=R(1,:)R2 = 0.8132 0.0099 0.1389 0.2028 R2=R(3,:)R2 =0.0153 0.7468 0.4451 0.9318 6、 输入矩阵，使用全下标方式取出元素“3”；使用单下标方式取出元素“8”；提取A矩阵的第2行和第3行元素形成一个新的子矩阵B。A1 2 3;4 5 6 ;7 8 9A = 1 2 3 4 5 6 7 8 9 C=A(1,3)C = 3 D=A(6)D =8 q=A(:,2) q = 2 5 8 a=A(:,2)a = 2 5 8 B=(a,q)? B=(a,q) |Error: Incomplete or misformed expression or statement. B=a;qB = 2 5 8 2 5 87、 输入A为3*3的魔方阵，B为3*3的单位阵，由小矩阵组成3*6的大矩阵C和6*3的大矩阵D。 A=magic(3)A = 8 1 6 3 5 7 4 9 2 B=eye(3,3)B = 1 0 0 0 1 0 0 0 1 A(3,6)=0A = 8 1 6 0 0 0 3 5 7 0 0 0 4 9 2 0 0 0 B=(3,6)? B=(3,6) |Error: Incomplete or misformed expression or statement. B(3,6)=0B = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 C=A-BC = 7 1 6 0 0 0 3 4 7 0 0 0 4 9 1 0 0 0 D=transpose(C)D = 7 3 4 1 4 9 6 7 1 0 0 0 0 0 0 0 0 0法二： A=magic(3)A = 8 1 6 3 5 7 4 9 2 B=eye(3,3)B = 1 0 0 0 1 0 0 0 1 A=(:,4:6)=B? A=(:,4:6)=B |Error: Missing variable or function. A(:,4:6)=BA = 8 1 6 1 0 0 3 5 7 0 1 0 4 9 2 0 0 1 C=AC = 8 1 6 1 0 0 3 5 7 0 1 0 4 9 2 0 0 1 transpose(C)ans = 8 3 4 1 5 9 6 7 2 1 0 0 0 1 0 0 0 18、 计算函数的值，其中，t为，z为0.707，步长为。 t=0:0.1*pi:2*pi; z=0.707; fx=10*exp(-2*t)*sin(4*t)/sqrt(1-z2)? Error using = mtimesInner matrix dimensions must agree. fx=10*exp(-2*t).*sin(4*t)/sqrt(1-z2)fx = Columns 1 through 10 0 7.1743 2.3655 -1.2620 -1.0893 -0.0000 0.3100 0.1022 -0.0545 -0.0471 Columns 11 through 20 -0.0000 0.0134