实验二数值数组.docVIP

温州大学瓯江学院 WENZHOU UNIVERSITY OUJIANG COLLEGE 《Matlab仿真及其应用》 实 验 报 告 专　 业 班　级 学生姓名 学　号 指导教师 温州大学瓯江学院教务部制 实验二 数值数组（矩阵）及其运算 一、实验目的 1、掌握MATLAB软件环境下进行数值数组（矩阵）的创建、访问的基本方法 2、掌握数值数组（矩阵）的算术运算、逻辑运算和关系运算的方法。 二、实验内容 1、关于“变量” 不需要事先定义，在遇到新的变量名时，自动创建并分配存储空间；当遇到已存在的变量时，更新其内容。 变量名由字母、数字和下划线构成，并且必须以字母开头。 区分大小写字母。 对变量的赋值可采用赋值语句：变量＝表达式[；] 如果行末加上分号，则在屏幕上不显示结果，这在设计M文件时非常有用。 MATLAB中提供了一些用户不能清除的固定变量： ans，eps，pi，i，j，NaN，Inf，realmax，realmin，nargin，nargout 2、 矩阵的创建和操作 （1）矩阵的创建 1）建立程序文件timu1_1.m， 输入元素列表 输入矩阵，记录结果。 2）利用MATLAB内部函数产生矩阵，建立程序文件timu1_2.m，完成下面功能，记录结果。 函数 功 能 函数 功 能 compan 伴随阵 zeros 元素全为0的矩阵 diag 对角阵 ones 元素全为1的矩阵 magic 魔方矩阵 rand 元素服从均匀分布的随机矩阵 linspace 线性空间向量 randn 元素服从正态分布的随机矩阵 logspace 对数空间向量 eye 对角线上元素为1的矩阵（单位矩阵） 产生一个3阶单位矩阵。 产生一个2×5的全1矩阵。 产生一个2×2的全0方阵。 产生一个[0，1]之间均匀分布的随机向量R(1×6)。 提示 R=rand(1,6) 产生一个[-5,5]之间均匀分布的随机阵列R1（2，6） 提示输入： R1=5-2*5*rand(2,6) 产生均值为0、方差为1的标准正态分布随机数4×4方阵S，输出结果。 提示：S=randn(4) （2）矩阵的下标——访问（存取）矩阵元素的依据 全下标：行下标+列下标 ，如a（i,j），几何概念清楚，引述简单 单下标（linear index）：矩阵元素在内存中是按“列”存储的 “一维编号”：把矩阵的所有列按先左后右的次序首尾相接排成“一维长列”后依次编号 两种下标的转换关系：A为（m×n）的二维数组（矩阵），其中某元素的“全下标”为：（r，c），则相应的“单下标”为：k=（c-1）×m+r MATLAB有两个指令可以实现以上标志方法的转换： sub2ind 据全下标换算出单下标，如IND = sub2ind(SIZ,I,J) ind2sub 据单下标换算出全下标，如[I,J] = ind2sub(SIZ,IND) 建立程序文件timu2_1.m 建立矩阵，然后： 访问元素。计算矩阵a中第1行、第2列元素与第2行、第3列元素之和，并将结果保存在变量aa中，输出结果。 2）修改元素。将a的第3行、第3列元素值改为19；第2行元素改为5 10 15，输出结果。 提示：可用a(2,1:3)=[5 10 15]批量修改 ☆ 扩大矩阵（增加行或列）。将a扩大为4×3的矩阵，并使第4行第2列元素为20 ，输出结果。 （3）矩阵行列删除 在前题的基础上，建立程序文件timu3_1.m， （4）矩阵连接 建立程序文件timu3_2.m进行如下操作，记录结果并说明操作的含义； （5）矩阵操作 diag 对角矩阵和矩阵的对角化 rot90 矩阵旋转90度 reshape 矩阵元素重新排列 tril 矩阵的下三角阵 triu 矩阵的上三角阵 fliplr 矩阵左右翻转 flipud 矩阵上下翻转 cat 矩阵连接 repmat 复制并平铺矩阵 建立程序文件timu3_3.m，输入矩阵，然后分别进行如下计算并观察结果。 3、矩阵运算 （1）算术运算 算术运算符：＋ － * / \ ^ 算术运算和对应的M函数 功 能 算术表达式 M函数 功 能 算术表达式 M函数 加 法 A+B plus（A，B） 矩阵左除 A \ B mldivide（A，B） 减 法 A-B minus（A，B） 数组左除 A． \ B ldivide（A，B） 矩阵乘法 A*B mtimes（A，B） 矩阵求幂 A^B mpower（A，B） 数组乘法 A．*B times（A，B） 数组求幂 A．^B power（A，B） 矩阵右除 A/B mrdivide（A，B） 矩阵共轭转置 A ct