《MATLAB变量和数据操作》由会员分享，可在线阅读，更多相关《MATLAB变量和数据操作（38页珍藏版）》请在金锄头文库上搜索。
1、第2讲 变量和数据操作 2.1 变量与赋值 2.2 内存变量的管理 2.3 MATLAB常用数学函数 2.4 MATLAB矩阵 2.5 M文件,2.1 变量与赋值,2.1.1 变量命名,2.1.2 赋值语句,(1) 变量=表达式 (2) 表达式其中表达式是用运算符将有关运算量连接起来的式子，其结果是一个矩阵。,例2-1 计算表达式的值，并显示计算结果。在MATLAB命令窗口输入命令：x=1+2i;y=3-sqrt(17);z=(cos(abs(x+y)-sin(78*pi/180)/(x+abs(y)其中pi和i都是MATLAB预先定义的变量，分别代表代表圆周率和虚数单位。输出结果是： z =
2、 -0.3488 + 0.3286i,2.1.3 预定义变量,在MATLAB工作空间中，还驻留几个由系统本身定义的变量。例如，用 pi 表示圆周率的近似值，用 i，j 表示虚数单位。预定义变量有特定的含义，在使用时，应尽量避免对这些变量重新赋值。,2.2 内存变量的管理,2.2.1 内存变量的删除与修改MATLAB 工作空间窗口专门用于内存变量的管理。在工作空间窗口中可以显示所有内存变量的属性。当选中某些变量后，再单击 Delete按钮，就能删除这些变量。当选中某些变量后，再单击 Open 按钮，将进入变量编辑器。通过变量编辑器可以直接观察变量中的具体元素，也可修改变量中的具体元素。,clea
3、r命令用于删除MATLAB工作空间中的变量。who和whos这两个命令用于显示在MATLAB工作空间中已经驻留的变量名清单。who命令只显示出驻留变量的名称。whos在给出变量名的同时，还给出它们的大小、所占字节数及数据类型等信息。,2.2.2 内存变量文件利用 MAT文件可以把当前MATLAB工作空间中的一些有用变量长久地保留下来，扩展名是 .mat。MAT文件的生成和装入由save和load命令来完成。常用格式为：save 文件名 变量名表 -append-asciiload 文件名 变量名表 -ascii,其中，文件名可以带路径，但不需带扩展名.mat，命令隐含一定对.mat文件进行操作
4、。变量名表中的变量个数不限，只要内存或文件中存在即可，变量名之间以空格分隔。当变量名表省略时，保存或装入全部变量。-ascii选项使文件以ASCII格式处理，省略该选项时文件将以二进制格式处理。save命令中的 append 选项控制将变量追加到MAT文件中。,2.2.3 数据的输出格式,MATLAB 用十进制数表示一个常数，具体可采用日常记数法和科学记数法两种表示方法。在一般情况下，MATLAB 内部每一个数据元素都是用双精度数来表示和存储的。数据输出时用户可以用format命令设置或改变数据输出格式。format命令的格式为：format 格式符 其中格式符决定数据的输出格式,各种 for
5、mat 格式,2. 3 MATLAB常用数学函数,MATLAB 提供了许多数学函数，函数的自变量规定为矩阵变量，运算法则是将函数逐项作用于矩阵的元素上，因而运算的结果是一个与自变量同维数的矩阵。,Matlab中常见数学函数,函数使用说明：(1) 三角函数以弧度为单位计算。(2) abs 函数可以求实数的绝对值、复数的模、字符串的ASCII码值。(3) 用于取整的函数有 fix、floor、ceil、round，要注意它们的区别。(4) rem与 mod函数的区别。rem(x,y)和mod(x,y) 要求 x, y 必须为相同大小的实矩阵或为标量。,系统预定义变量,pi 圆周率 ，其值为 ima
6、g(log(-1),inf/Inf 无穷大,nan/NaN Not-a-Number，一个不定值，如 0/0,eps 浮点运算相对精度,特殊变量 ans,i/j 虚部单位，即,应尽量避免给系统预定义变量重新赋值！,2.3.2 MATLAB特殊变量表,2.4.1 向量与矩阵的生成,从矩阵中抽取行或列,2.4 MATLAB矩阵,向量与矩阵的生成（续）,向量与矩阵运算,矩阵的生成,直接输入: A=1, 2, 3; 4, 5, 6; 7, 8, 9,由向量生成,由函数生成,通过编写m文件生成,例： C=magic(3),常见矩阵生成函数,2.4.2 矩阵的拆分,1. 矩阵元素通过下标引用矩阵的元素，例
7、A(3,2)=200采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在 MATLAB 中，矩阵元素按列存储，先第一列，再第二列，依次类推。例如 A=1,2,3;4,5,6;A(3) ans = 2,显然，序号(Index)与下标(Subscript )是一一对应的，以mn矩阵A为例，矩阵元素A(i,j) 的序号为 (j-1)*m+i。其相互转换关系也可利用sub2ind和ind2sub函数求得。,2. 矩阵拆分(1) 利用冒号表达式获得子矩阵 A(:,j)表示取A矩阵的第j列全部元素；A(i,:)表示A矩阵第i行的全部元素；A(i,j)表示取A矩阵第i行、第j列
8、的元素。 A(i:i+m,:)表示取A矩阵第ii+m行的全部元素；A(:,k:k+m)表示取A矩阵第kk+m列的全部元素，A(i:i+m,k:k+m)表示取A矩阵第 ii+m行内，并在第kk+m列中的所有元素。此外，还可利用一般向量和end运算符来表示矩阵下标，从而获得子矩阵。end表示某一维的末尾元素下标。,(2) 利用空矩阵删除矩阵的元素在MATLAB中，定义 为空矩阵。给变量X赋空矩阵的语句为X= 。注意，X= 与clear X不同，clear是将X从工作空间中删除，而空矩阵则存在于工作空间中，只是维数为0。,2.4.3 特殊矩阵,1. 通用的特殊矩阵常用的产生通用特殊矩阵的函数有：ze
9、ros：产生全0矩阵(零矩阵)。ones：产生全1矩阵(幺矩阵)。eye：产生单位矩阵。rand：产生01间均匀分布的随机矩阵。randn：产生均值为 0，方差为 1的标准正态分布随机矩阵。,例2-3 分别建立33、32和与矩阵A同样大小的零矩阵。(1) 建立一个33零矩阵：zeros(3)(2) 建立一个32零矩阵：zeros(3,2)(3) 设A为23矩阵，则可以用zeros(size(A)建立一个与矩阵A同样大小零矩阵。A=1 2 3;4 5 6; %产生一个23阶矩阵Azeros(size(A) %产生一个与矩阵A同样大小的零矩阵,例2-4 建立随机矩阵：(1) 在区间20,50内均匀
10、分布的5阶随机矩阵。(2) 均值为 0.6、方差为 0.1的5阶正态分布随机矩阵。命令如下：x=20+(50-20)*rand(5)y=0.6+sqrt(0.1)*randn(5)此外，常用的函数还有reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成mn的二维矩阵。,2. 用于专门学科的特殊矩阵(1) 魔方矩阵 魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于 n 阶魔方阵，其元素由 1,2,3,n2 共 n2 个整数组成。MATLAB提供了求魔方矩阵的函数magic(n)，其功能是生成一个n阶魔方阵。,例2-5 将101125等25个数填入
11、一个5行5列的表格中，使其每行每列及对角线的和均为565。 M=100+magic(5),(2) 范得蒙矩阵范得蒙(Vandermonde)矩阵最后一列全为1，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在MATLAB中，函数vander(V)生成以向量V为基础向量的范得蒙矩阵。例如，A=vander(1;2;3;5)即可得到上述范得蒙矩阵。,(3) 希尔伯特矩阵在MATLAB中，生成希尔伯特矩阵的函数是hilb(n)。使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。MATLAB中，有一个专门求希尔伯特矩阵的逆的函数i
12、nvhilb(n)，其功能是求 n阶的希尔伯特矩阵的逆矩阵。,例2-6 求4阶希尔伯特矩阵及其逆矩阵。命令如下：format rat %以有理形式输出H=hilb(4)H=invhilb(4),(4) 托普利兹矩阵托普利兹(Toeplitz)矩阵除第一行第一列外，其他每个元素都与左上角的元素相同。生成托普利兹矩阵的函数是toeplitz(x,y)，它生成一个以x为第一列，y为第一行的托普利兹矩阵。这里x, y均为向量，两者不必等长。toeplitz(x)用向量x生成一个对称的托普利兹矩阵。例如 T=toeplitz(1:6),(5) 伴随矩阵MATLAB生成伴随矩阵的函数是compan(p)，
13、其中p是一个多项式的系数向量，高次幂系数排在前，低次幂排在后。例如，为了求多项式的x3-7x+6的伴随矩阵，可使用命令：p=1,0,-7,6;compan(p),(6) 帕斯卡矩阵我们知道，二次项(x+y)n展开后的系数随n的增大组成一个三角形表，称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡(Pascal)矩阵。函数pascal(n)生成一个n阶帕斯卡矩阵。,例2-7 求(x+y)5的展开式。在MATLAB命令窗口，输入命令： pascal(6)矩阵次对角线上的元素1,5,10,10,5,1即为展开式的系数。,用 Matlab 语言编写的程序称为 M 文件。,M 文件是由若干 Matla
14、b 命令组合在一起构成的，它可以完成某些操作，也可以实现某种算法。事实上，Matlab 提供的内部函数以及各种工具箱，都是利用 Matlab 语言开发的 M 文件。用户也可以结合自己的工作需要，开发自己的程序或工具箱。,2.5 M 文件,M 文件以 .m 为扩展名。,M 文件是一个文本文件，可以用任何文本编辑器来建立和编辑，通常使用 Matlab 自带的 M 文件编辑器。,M 文件的建立与打开,在命令窗口中输入 f2cs，即可执行该 M 文件。,例：编写一个脚本文件将华氏温度转化为摄氏温度,脚本文件举例,脚本文件运行后，文件中所产生的所有变量都驻留在 Matlab 当前工作空间，可以用 whos 查看。,程序2： 首先建立函数文件f2c.m。 function c=f2c(f) c=5*(f-32)/9 然后在MATLAB的命令窗口调用该函数文件。 clear; y=input(Input Fahrenheit temperature：); x=f2c(y) 输出情况为： Input Fahrenheit temperature：70 c =21.1111 x =21.1111,
