实验 1 矩阵的建立和基本运算一、实验目的熟悉 矩阵 (matrix laboratory) 初等变换 的方法以及 矩阵运算 的各种命令 . 二、实验内容与要求1. 启动与退出2. 数、数组、矩阵的输入(1)数的输入>>a=5 >>b=2-5i注意:在行尾加“; ”,该行结果不显示; 标点符号一定要在英文状态下输入!(2)数组的输入>>b=[1,3,5,7,9,11] >>c=1:2:11 >>d=linspace(1,11 ,6) 问题:若 b 为在 0~2π (π 用 pi 表示)之间均匀分布的22 个数据,c=(1.3,2.5,7.6,2,- 3),d=(23,20,17,14,11,8,5,2) ,各用何种方法输入较简单?(3)矩阵的输入>>A=[2,3,5;1,3,5;6,9,4]%行之间要用分号隔开A= 2 3 5 1 3 5 6 9 4 >>m=input('请输入初始量, m=') ;请输入初始量, m= 问题:输入A(2,3) ,结果如何?输入A(7)又如何?注意: 变量名开头必须是英文字母,变量名对字母大小写是区分的. 3. 矩阵大小的测试和定位numel( Number of elements in an array or subscripted array expression.)>>A=[3,5,6;2,5,8;3,5,9;3,7,9] >>d=numel(A)%测试定矩阵A 的元素数, 5.x版本没有此命令>>[n,m] = size(A) %测试 A 的行 (n) 、列 (m) 数>>[i,j] = find(A>3)%找出 A 中大于 3 的元素的行列数注意: 对一个数组可用n = length(A) ,A 若是矩阵, n 给出 A 的行、列数的最大值 . 4. 矩阵的块操作>>A(2,:)%取出 A 的第 2 行的所有元素>>A([1,3],:)%取出 A 的第 1,3 行的所有元素>>A(2:3,1:2)%取出 A 的 2,3 行与 1 ,2 列交叉的元素>>A([1,3],:) = A([3,1],:) %将 A 的 1 行和 3 行互换问题:如何将 A 的 2,3 列互换?>>A(2,:) = 4 %将 A 的第 2 行的所有元素用4 取代>>A(find(A==3))=-3 %将 A 中等于 3 所有元素换为 -3 >>A(2,:) = [] %删除 A 的第 2 行>>reshape (A,2,3 ) %返回以 A 的元素重新构造的2×3 维矩阵>>[A(1:3,2:3),A(2:4,1:2);A,A(:,2)]%由小矩阵构造大矩阵。5. 矩阵的翻转操作flip( 抛、弹、翻转 ) ,rot ation (旋转)>>flipud(A) %A进行上下翻转>>fliplr(A) %A进行左右翻转>>rot90(A)%A逆时针旋转 90 °6. 特殊矩阵的产生random(随意 , 任意),randn (Normally distributed random numbers.)randperm (Random permutation.(排列))>>A = eye(n) %产生 n 维单位矩阵>>A = ones...