MATLAB程序设计入门.doc

实验一、MATLAB程序设计入门实验目的熟悉matlab界面;掌握matlab的变量的命名规则;掌握常用数学运算符号和标点符号的使用方法;掌握常用数学函数的使用方法;掌握建立和使用M文件的方法;学会创建简单的数组并对数组元素进行访问,掌握数组的各种运算;学会创建数组并对数组元素进行访问和操做,掌握矩阵的各种运算规则;掌握关系与逻辑运算符的使用规则；实验原理与方法一、变量与函数1、变量MATLAB中变量的命名规则是：（1）变量名必须是不含空格的单个词；（2）变量名区分大小写；（3）变量名最多不超过19个字符；（4）变量名必须以字母打头，之后可以是任意字母、数字或下划线，变量名中不允许使用标点符号.（5）关键字(如if、while等)以及固有函数(如max,sum等)不能作为变量名。特殊变量表2、数学运算符号及标点符号（1）MATLAB的每条命令后，若为逗号或无标点符号，则显示命令的结果；若命令后为分号，则禁止显示结果. （2）“%” 后面所有文字为注释. （3） “.”表示续行.对于点乘与点乘幂运算,主要针对的是向量和数组的运算（数组的概念在下面有）.如:x=1 1;1 1;x.2 表示对每个元素取幂为2的运算 结果：x=1 1;1 1x2 表示x矩阵的2次幂 结果：x=2 2;2 2再有 y=2 2;2 2;x*y 表示普通的矩阵乘法 结果：4 4;4 4x.*y 表示两个矩阵对应元素相乘 结果：2 2;2 2 当然若y是常数,则*与.*没有区别/与./的区别 与* .* 类似3、数学函数二、数组与矩阵1、创建简单的数组x=a b c d e f 创建包含指定元素的行向量x=first：last 创建从first开始，加1计数，到last结束的行向量x=first：increment：last 创建从first开始，加increment计数，last结束的行向量x=linspace(first，last，n） 创建从first开始，到last结束，有n个元素的行向量x=logspace(first，last，n） 创建从first开始，到last结束，有n个元素的对数分隔行向量. 例：x=linspace(1,2,10)y=10.(x)z=logspace(1,2,10) 注:z的输出是: z =10.0000 12.9155 16.6810 21.5443 27.8256 35.9381 46.4159 59.9484 77.4264 100.0000即 z(k)=10(k*t) t=(last-first)/(n-1)2、 数组元素的访问（1）访问一个元素： x(i)表示访问数组x的第i个元素.x=1 2;3 4; x(1)=1 x(2)=3 x(3)=2 x(4)=4 即先访问列也可以用x(i,j), 即第i行第j列的元素（2）访问一块元素： x(a ：b ：c)表示访问数组x的从第a个元素开始，以步长为b到第c个元素（但不超过c），b可以为负数，b缺损时为1. 如：x(2:2:6)访问第2，4，6个元素x(2:2:5)访问第2，4个元素 没有第5个元素x(2:5)访问第2，3，4，5个元素（3）直接使用元素编址序号. x(a b c d) 表示提取数组x的第a、b、c、d个元素构成一个新的数组x(a) x(b) x(c) x(d). 3、数组的方向前面例子中的数组都是一行数列，是行方向分布的. 称之为行向量. 数组也可以是列向量，它的数组操作和运算与行向量是一样的，唯一的区别是结果以列形式显示. 产生列向量有两种方法： 直接产生 例 c=1；2；3；4 转置产生 例 b=1 2 3 4; c=b 单引号是转置运算 说明：以空格或逗号分隔的元素指定的是不同列的元素，而以分号分隔的元素指定了不同行的元素. 4、数组的运算（1）标量-数组运算数组对标量的加、减、乘、除、乘方是数组的每个元素对该标量施加相应的加、减、乘、除、乘方运算. 设：a=a1,a2,an, c=标量则：a+c=a1+c,a2+c,an+c a.*c=a1*c,a2*c,an*ca./c= a1/c,a2/c,an/c(右除） a.c= c/a1,c/a2,c/an (左除） a.c= a1c,a2c,anc c.a= ca1,ca2,can （2）数组-数组运算当两个数组有相同维数时，加、减、乘、除、幂运算可按元素对元素方式进行的，不同大小或维数的数组是不能进行运算的. 设：a=a1,a2,an, b=b1,b2,bn则：a+b= a1+b1,a2+b2,an+bn a.*b= a1*b1,a2*b2,an*bn a./b= a1/b1,a2/b2,an/bn a.b=b1/a1,b2/a2,bn/ana.b=a1b1,a2b2,anbn（3）数组维数：length(a) 即数组的列数5、矩阵的建立逗号或空格用于分隔某一行的元素，分号用于区分不同的行；除了分号，在输入矩阵时，按Enter键也表示开始一新行； 输入矩阵时，严格要求所有行有相同的列.例：m=1 2 3 4 ；5 6 7 8；9 10 11 12 p=1 1 1 1 2 2 2 2 3 3 3 3特殊矩阵的建立：a= 产生一个空矩阵，当对一项操作无结果时，返回空矩阵，空矩阵的大小为零. b=zeros(m，n) 产生一个m行、n列的零矩阵c=ones(m，n) 产生一个m行、n列的元素全为1的矩阵d=eye(m，n) 产生一个m行、n列的单位矩阵6、矩阵中元素的操作（1）矩阵A的第r行：A（r，：）（2）矩阵A的第r列：A（：，r）（3）依次提取矩阵A的每一列，将A拉伸为一个列向量：A（：）（4）取矩阵A的第i1i2行、第j1j2列构成新矩阵:A(i1:i2, j1:j2)（5）以逆序提取矩阵A的第i1i2行，构成新矩阵:A(i2:-1：i1，：) -1表示逆序（6）以逆序提取矩阵A的第j1j2列，构成新矩阵:A(:, j2:-1：j1 ） 以顺序提取矩阵A的第i1i2行，构成新矩阵:A(i1:i2,:)以顺序提取矩阵A的第j1j2列，构成新矩阵:A(:,j1:j2 )（7）删除A的第i1i2行，构成新矩阵:A(i1:i2，：)= （8）删除A的第j1j2列，构成新矩阵:A(：， j1:j2)= （9）将矩阵A和B拼接成新矩阵：A，B把B接在A右面 A；B把B接在A下面 A，B 要求行数相同 A；B 要求列数相同7、矩阵的运算（1）标量-矩阵运算：同标量-数组运算。（2）矩阵加法：A+B（3）矩阵乘法：A*B（4）方阵的行列式：det（A）（5）方阵的逆：inv（A）（6）方阵的特征值与特征向量：V，D=eig(A)（7）求矩阵维数：m,n=size(A)（8）求数组或矩阵元素之和：sum(A) 是每一列求和,得到一个行向量（9）求数组或矩阵元素的均值：mean(A) 是每一列求一个平均值（10）求矩阵或数组的最大最小值：max,min 同样是每一列求一个最大(小)值（11）求数组或矩阵的非零元素的下标： find（A） 默认A大于0find（AC）例：、X = 1 0 4 -3 0 0 0 8 6;indices = find(X)结果：indices = 1 3 4 8 9例：X = 1 0 4 -3 0 0 0 8 6;find(X 2)结果：ans = 3 8 9例：X = 3 2 0; -5 0 7; 0 0 1;r,c,v = find(X)则结果为：r为矩阵X的非零元素的行指标所构成的向量c为矩阵X的非零元素的列指标所构成的向量v为矩阵X的非零元素的所构成的向量。例：X = 3 2 0; -5 0 7; 0 0 1;r,c,v = find(X2)（11）求矩阵的范数:norm(A,1) 计算矩阵A的1范数norm(A,2) 计算矩阵A的2范数norm(A,inf) 计算矩阵A的无穷范数可参考: /doc-view-501.html 三、关系与逻辑运算1、关系操作符2、逻辑运算符四、程序设计Matlab有两种工作方式：1）人机交互的命令行指令操作方式，即在命令窗口每输入一条命令，则立即运行该命令得到结果。2）进行控制流的程序设计，即编制一种可存储的以M为扩展名的文件（简称M文件），M文件有两种形式，命令式（Script）和函数式(Function).两者相同之处在于它们都是以m作为扩展名的文本文件，不进入命令窗口，而是由文本编辑器来创建的外部文本文件。M文件分两种：(1)命令式M文件：也称脚本文件（Script file）， 就是将Matlab的一系列命令按顺序编制成一个文本文件，文件名后缀为M，然后在command window 下运行文件名，则按顺序执行文件中的命令。文件建立方法：1. 在Matlab中，点:File-New-M-file2. 在编辑窗口中输入程序内容3. 点：File-Save，输入文件名，后缀为M，存盘例：建立命令式脚本文件qwe.m在编辑窗口输入如下命令：a=1b=2; %行尾加分号，注意运行结果。c=a+b然后存盘。在command window 下输入qwe并回车。注： 在运行次文件之前，需要把它所在目录加到MATLAB的搜索路径上去，或将文件所在目录设为当前目录。 （2）函数式M文件：（function file） MATLAB的内部函数是有限的，有时为了研究某一个函数的各种性态，需要为MATLAB定义新函数，为此必须编写函数文件. 一个函数M文件与脚本文件类似之处在于它们都是一个有.m 扩展名的文本文件。如同脚本M文件一样，函数M文件不进入命令窗口，而是由文本编辑器所创建的外部文本文件。一个函数的M文件与脚本文件在通信方面是不同的。函数与MATLAB工作空间之间的通信，只通过传递给它的变量和通过它所创建的输出变量。在函数内中间变量不出现在MATLAB工作空间，或与MATLAB工作空间不交互。一个函数的M文件的第一行把M文件定义为一个函数，并指定它的名字。它与文件名相同，但没有.m 扩展名。它也定义了它的输入和输出变量。这类文件的第一行必须是一特殊字符function开始，格式为： function Y1,Y2Ym=函数名（X1,X2,Xn）其中X1,X2,Xn为输入变量，Y1,Y2Ym为输出变量，均可，X为数、数组或者矩阵。M文件建立方法：1. 在Matlab中，点:File-New-M-file2. 在编辑窗口中输入程序内容3. 点：File-Save，存盘，M文件名必须与函数名一致。注: 函数式文件执行之后，只保留最后结果(输出变量)，不保留中间过程，所定义的变量也仅在函数内部起作用，并随调用的结束而被清除。例：定义函数 f(x1,x2)=100(x2-x12)2+(1-x1)21.建立M文件：fun.mfunction f=fun(x)f=100*(x(2)-x(1)2)2+(1-x(1)22. 可以直接使用函数fun.m例如：计算 f(1,2), 只需在Matlab命令窗口键入命令：x=1 2fun(x)注：上述所定义的函数fun的输入变量只有一个X,X是一个二维数组。也可使输入变量为两个，但是调用该函数时，输入参数必须也是两个数。1.建立M文件：fun.mfunction f=fun(x1,x2)f=100*(x2-x12)2+(1-x1)22. 可以直接使用函数fun.m例如：计算 f(1,2), 只需在Matlab命令窗口键入命令：fun(1,2)10)控制语句MATLAB提供三种决策或控制流结构：for循环、while循环、if-else-end结构. 这些结构经常包含大量的MATLAB命令，故经常出现在MATLAB程序中，而不是直接加在MATLAB提示符下1、for循环：允许一组命令以固定的和预定的次数重复for x=arraycommandsend在for和end语句之间的命令串commands按数组（array）中的每一列执行一次. 在每一次迭代中，x被指定为数组的下一列，即在第n次循环中，x=array(：，n)例：求1+2+100先建立脚本文件：myfun.mmysum=0;for i=1:1:100 mysum=mysum+i;end存盘，然后在命令窗口输入myfun，并回车。2、While循环与for循环以固定次数求一组命令相反，while循环以不定的次数求一组语句的值. while expression命令语句1命令语句2命令语句nend只要在表达式(expression)里的所有元素为真，就执行while和end语句之间的命令串例 求1+2+100建立脚本文件：mysum.mmysum=0;i=1;while (i=100) %此处的括号可加可不加（但要与前面的关键字隔开） %下面的if-else同样 mysum=mysum+i; i=i+1;end 例：Fibonacci数组的元素满足Fibonacci 规则： ，；且。现要求该数组中第一个大于10000的元素。a(1)=1;a(2)=1;i=2;while a(i)=10000 a(i), %输出数组中第一个大于10000的元素 break; %跳出所在的一级循环。 end;end,i %输出元编号 ans = 10946i =21 5return：终止后面语句的执行，并返回结果。6、If-Else-End结构（1）有一个选择的一般形式是：if expression命令语句1命令语句2命令语句nend如果在表达式(expression)里的所有元素为真，就执行if和end语句之间的命令串 （2）有两个选择的一般形式是：if expression1elseend(3) 有三个或更多的选择的一般形式是：if expression1elseif expression2 elseif expression3elseend例：建立符号函数（fhfun.m）function f=fhfun(x)if x0 f=1;elseif x=0 f=0;else f=-1;end可以从命令窗口输入fhfun(3), fhfun(-2)分别计算函数