MATLAB基本运算.ppt

《MATLAB基本运算.ppt》由会员分享，可在线阅读，更多相关《MATLAB基本运算.ppt（75页珍藏版）》请在三一办公上搜索。

1、1,第2章 MATLAB基本运算,2.1 数据类型 2.2 数组2.3矩阵2.4 字符串2.5 日期和时间2.6 结构体和元胞数组2.7 多维数组2.8 关系运算和逻辑运算2.9 数组的信息获取2.10 多项式,2,2.1数据类型,MATLAB 定义了15种基本的数据类型，包括整型、浮点型、字符型和逻辑型等，用户可以定义自己的数据类型。任何类型的数据均按数组的形式进行存储和运算。,3,常数和变量,1.常数：采用十进制表示，可用带小数点的形式直接表示或用科学计数法表示。2.变量MATLAB的变量命名规则如下：变量名区分字母的大小写。例如，“a”和“A”是不同的变量。变量名不能超过63个字符，第6

2、3个字符后的字符被忽略。变量名必须以字母开头，变量名的组成可以是任意字母、数字或者下划线，但不能含有空格和标点符号（如，。%等）。关键字（如if、while等）不能作为变量名。,练习,4,赋值语句(1)变量=表达式(2)表达式其中表达式是用运算符将有关运算量连接起来的式子，其结果是一个矩阵。预定义的特殊变量 在MATLAB工作空间中，还驻留几个由系统本身定义的变量。例如，用pi表示圆周率的近似值，用i，j表示虚数单位。预定义变量有特定的含义，在使用时，应尽量避免对这些变量重新赋值。,5,特殊变量表,6,例2-1 计算表达式的值，并显示计算结果。在MATLAB命令窗口输入命令：x=1+2i;y=

3、3-sqrt(17);z=(cos(abs(x+y)-sin(78*pi/180)/(x+abs(y)其中pi和i都是MATLAB预先定义的变量，分别代表代表圆周率和虚数单位。输出结果是：z=-0.3488+0.3286i,7,练习a=5+3ib=2+4jc=2*pid=2/02*d1/d,8,整数和浮点数,1.整数,练习,9,2.浮点数,10,2.1.3 复数,Matlab 把复数作为一个整体处理（而不必像其他程序语言，将实部和虚部分开处理）。虚数单位用预定义变量i或j表示。复数z=a+bi 直角坐标表示和极坐标表示之间转换的matlab指令如下：Real(z)给出复数z的实部 Imag(z

4、)给出复数z的虚部Abs(z)给出复数z的模Angle(z)给出复数z的相角,11,2.2 数组,空数组（empty array）：没有元素的数组；标量（scalar）：是指11的矩阵，即为只含一个数的矩阵；向量（vector）：是指1n或n1的矩阵，即只有一行或者一列的矩阵；矩阵（matrix）：是一个矩形的mn数组，即二维数组；数组（array）：是指多维数组mnk，其中矩阵和向量都是数组的特例。,12,2.2.1 创建简单的数组,x=a b c d e f 创建包含指定元素的行向量,x=first：last 创建从first开始，加1计数，到last结束的行向量,x=first：incr

5、ement：last 创建从first开始，加increment计数，last结束的行向量,x=linspace(first，last，n）创建从first开始，到last结束，有n个元素的行向量 y=1:7y=1 2 3 4 5 6 7,13,例如：,z=3:2:9z=3 5 7 9 v=y,zv=1 2 3 4 5 6 7 3 5 7 9 u=linspace(2,9,11)u=Columns 1 through 6 2.0000 2.7000 3.4000 4.1000 4.8000 5.5000 Columns 7 through 11 6.2000 6.9000 7.6000 8.3

6、000 9.0000,14,2.2.2 对数组元素的访问,（3）直接使用元素编址序号.x(a b c d)表示提取数组x的第a、b、c、d个元素构成一个新的数组 x(a)x(b)x(c)x(d).,（2）访问一块元素：x(a：b：c)表示访问数组x的从第a个元素开始，以步长为b到第c个元素（但不超过c），b可以为负数，b缺损时为1.,（1）访问一个元素：x(i)表示访问数组x的第i个元素.,例如：x=1:2:9x=1 3 5 7 9 y=x(2:1:5)y=3 5 7 9 z=x(1 3 5)z=1 5 9,15,前面例子中的数组都是一行数列，是行方向分布的.称之为行向量.数组也可以是列向量，

7、它的数组操作和运算与行向量是一样的，唯一的区别是结果以列形式显示.,产生列向量有两种方法：直接产生 例 c=1;2;3;4 转置产生 例 b=1 2 3 4;c=b 说明：以空格或逗号分隔的元素指定的是不同列的元素，而以分号分隔的元素指定了不同行的元素.,2.2.3 数组的方向,16,2.2.4 数组的运算,（1）标量-数组运算 数组对标量的加、减、乘、除、乘方是数组的每个元素对该标量施加相应的加、减、乘、除、乘方运算.,设：a=a1,a2,an,c=标量,则：a+c=a1+c,a2+c,an+c a.*c=a1*c,a2*c,an*c a./c=a1/c,a2/c,an/c(右除）a.c=c

8、/a1,c/a2,c/an(左除）a.c=a1c,a2c,anc c.a=ca1,ca2,can,17,（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/an a.b=a1b1,a2b2,anbn,18,例2-14 使用数组算术运算法则进行向量的运算。t=0:pi/3:2*pi;%t为行向

9、量 x=sin(t)*cos(t)?Error using=mtimesInner matrix dimensions must agree.x=sin(t).*cos(t)x=0 0.4330-0.4330-0.0000 0.4330-0.4330-0.0000 y=sin(t)./cos(t)y=0 1.7321-1.7321-0.0000 1.7321-1.7321-0.0000,19,2.2.5 数组元素的删除 在MATLAB中可以对数组中的单个元素、子矩阵和所有元素进行删除操作，删除就是将其赋值为空矩阵（用表示）。,20,2.3 矩 阵,逗号或空格用于分隔某一行的元素，分号用于区分不

10、同的行.除了分号，在输入矩阵时，按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,2.3.1 矩阵的建立,21,特殊矩阵的建立:,d=eye(m，n)产生一个m行、n列的单位矩阵,c=ones(m，n)产生一个m行、n列的元素全为1的矩阵,b=zeros(m，n)产生一个m行、n列的零矩阵,a=产生一个空矩阵，当对一项操作无结果时，返回空矩阵，空矩阵的大小为零.,22,2.3.2 矩阵中元素的操作,（1）矩阵A的第r行：A（r，：）,（2）矩阵A的第r列：A（：

11、，r）,（4）取矩阵A的第i1i2行、第j1j2列构成新矩阵:A(i1:i2,j1:j2),（5）以逆序提取矩阵A的第i1i2行，构成新矩阵:A(i2:-1：i1，：）,（6）以逆序提取矩阵A的第j1j2列，构成新矩阵:A(:,j2:-1：j1）,（7）删除A的第i1i2行，构成新矩阵:A(i1:i2，：)=,（8）删除A的第j1j2列，构成新矩阵:A(：，j1:j2)=,（9）将矩阵A和B拼接成新矩阵：A B；A；B,（3）依次提取矩阵A的每一列，将A拉伸为一个列向量：A（：）,23,（2）矩阵-矩阵运算 a）元素对元素的运算:与数组-数组运算相同.*./.,2.3.3 矩阵的运算,（1）标

12、量-矩阵运算 与标量-数组运算规则相同。,b）矩阵对矩阵运算：矩阵加减法：A+B A-B 矩阵乘法：A*B,矩阵除法：矩阵左除：AX=B，求 X。MATLAB 求解：X=AB若A为非奇异方阵，则 X=A-1B最小二乘解（若A不是方阵）矩阵右除：XA=B，求 X。MATLAB求解：X=B/A若A为非奇异方阵，则 X=BA-1最小二乘解（若A不是方阵）,24,例2-12 用矩阵除法求方程组的解，已知方程组：解：X=AB是方程A*X=B的解，将该方程变换成A*X=B的形式。其中：A=2-3 1;1-1 1;1 3 1;B=8;7;6;X=AB,练习,25,2.3.4 矩阵分析 2.3.4.1 对角阵

13、与三角阵1对角阵只有对角线上有非0元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为1的对角矩阵称为单位矩阵。,26,提取矩阵的对角线元素设A为mn矩阵，diag(A)函数用于提取矩阵A主对角线元素，产生一个具有min(m,n)个元素的列向量。diag(A)函数还有一种形式diag(A,k)，其功能是提取第k条对角线的元素。构造对角矩阵设V为具有m个元素的向量，diag(V)将产生一个mm对角矩阵，其主对角线元素即为向量V的元素。diag(V)函数也有另一种形式diag(V,k)，其功能是产生一个nn(n=m+k)对角阵，其第k条对角线的元素即为向量V的元素。

14、,27,例2-10 先建立55矩阵A，然后将A的第一行元素乘以1，第二行乘以2，第五行乘以5。,A=17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;.11,18,25,2,19;D=diag(1:5);D*A%用D左乘A，对A的每行 乘以一个指定常数,28,2三角阵三角阵又进一步分为上三角阵和下三角阵，所谓上三角阵，即矩阵的对角线以下的元素全为0的一种矩阵，而下三角阵则是对角线以上的元素全为0的一种矩阵。,29,上三角矩阵求矩阵A的上三角阵的MATLAB函数是triu(A)。triu(A)函数也有另一种形式triu(A,k)，其功能是求矩阵

15、A的第k条对角线以上的元素。例如，提取矩阵A的第2条对角线以上的元素，形成新的矩阵B。下三角矩阵在MATLAB中，提取矩阵A的下三角矩阵的函数是tril(A)和tril(A,k)，其用法与提取上三角矩阵的函数triu(A)和triu(A,k)完全相同。,30,2.3.4.2 矩阵的转置与旋转1矩阵的转置转置运算符是单撇号()。2矩阵的旋转利用函数rot90(A,k)将矩阵A旋转90的k倍，当k为1时可省略。,31,3矩阵的左右翻转对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第二列调换，依次类推。MATLAB对矩阵A实施左右翻转的函数是fliplr(A)4矩阵的上下翻转MAT

16、LAB对矩阵A实施上下翻转的函数是flipud(A)。,32,2.3.4.3 矩阵的逆对于一个方阵A，如果存在一个与其同阶的方阵B，使得：AB=BA=I(I为单位矩阵)则称B为A的逆矩阵，当然，A也是B的逆矩阵。求一个矩阵的逆是一件非常烦琐的工作，容易出错，但在MATLAB中，求一个矩阵的逆非常容易。求方阵A的逆矩阵可调用函数inv(A)。例2-11 用求逆矩阵的方法解线性方程组。Ax=b其解为：x=A-1b,33,2.3.4.4 方阵的行列式把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为矩阵所对应的行列式的值。在MATLAB中，求方阵A所对应的行列式的值的函数是det(A)

17、。,34,2.3.4.5 矩阵的秩与迹1矩阵的秩矩阵线性无关的行数与列数称为矩阵的秩。在MATLAB中，求矩阵秩的函数是rank(A)。2矩阵的迹矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。在MATLAB中，求矩阵的迹的函数是trace(A)。,35,2.3.4.6 矩阵的合并矩阵的合并就是把两个以上的矩阵连接起来得到一个新矩阵，“”符号可以作为矩阵合并操作符，命令格式如下：c=a b%将矩阵a和b水平方向合并为cc=a;b%将矩阵a和b垂直方向合并为c,36,2.4 字符串2.4.1 创建字符串,字符串由多个字符组成，是1n的字符数组；每一个字符都是字符数组的一个元素，以ASC

18、II码的形式存放并区分大小，而显示的形式则是可读的字符。创建字符串 s1=matlab 7.3s1=matlab 7.3,s2=字符串s2=字符串,37,（1）直接赋值用单引号()括起字符来直接赋值创建字符串。使用两个单引号()输入字符串中的单引号 s3=显示matlabs3=显示matlab（2）多个字符串组合 str2=s1;s2,;s3str2=matlab 7.3字符串 显示matlab 2.字符数组的存储空间MATLAB在存储字符串时，每一个字符以ASCII码的形式存放，占用两个字节。,练习,38,2.4.2 字符串函数,1.字符串合并strcat函数用于将字符串水平连接合并成一个新

19、字符串，合并的同时会将字符串尾的空格删除。语法格式如下：strcat(s1,s2,)%将s1，s2合并成一个长字符串char(s1,s2,)%将s1，s2合并成一个字符矩阵strvcat(s1,s2,)%将s1，s2合并成一个字符矩阵,39,2.字符串与数值的转换abs将字符串转换为ASCII码数值str2num将字符串转换为数值3.字符串的其他操作MATLAB 7.3还可以对字符串进行比较、查找、运行等操作,其函数见书57页表2-13。,40,例2-19 使用字符串函数进行运算。str=a+b,c+d,str=a+b,c+d,n=findstr(str,)%查找字符串中,的位置n=4 8 s

20、tr1=str(1:n(1)%取第一个,前的字符串str1=a+b,str1=str(n(1)+1:n(2)%取第二个,前的字符串str1=c+d,练习,41,2.5 日期和时间 日期时间的表示格式,MATLAB 没有专门的日期时间类型，以三种格式表示：日期字符串、连续的日期数值和日期向量，不同的日期格式可以相互转换。1.日期格式（1）日期字符串日期字符串是最常用的，有多种输出格式。例，“2007年1月1日”可以表示为：01-Jan-2007 08:50:10、01-Jan-2007、01/01/2007、等。（2）连续的日期数值 表示当前时间到起始时间的距离。（起始时间为公元元年1月1日）（

21、3）日期向量 year month day hour minute second,42,2.4.2 日期时间函数,1.获取系统时间date：按照日期字符串格式获取当前系统时间；now：按照连续的日期数值格式获取当前系统时间；clock：按照日期向量格式获取当前系统时间。2.提取日期时间信息分别使用year、month、day、hour、minute、second函数。3.日期时间的显示格式日期时间的显示可以使用datestr函数显示为字符串的样式。datestr函数的格式如下：datestr(d,f)%将日期按指定格式显示,43,例2-21 实用日期函数按指定格式显示日期时间。y=num2st

22、r(year(now)y=2009 m=num2str(month(now)m=2 d=num2str(day(now)d=27 s=今天是,y,年,m,月,d,日,datestr(now,HH:MM:SS PM)s=今天是2009年2月27日 3:48:49 PM,44,4.计时函数（1）cputime方法cputime是返回MATLAB启动以来的CPU时间：程序执行的时间程序代码执行结束后的cputime在程序代码执行前的cputime（2）tic/toc方法tic在程序代码开始用于启动的一个计时器；toc放在程序代码的最后，用于终止计时器的运行，并返回计时时间就是程序运行时间。（3）et

23、ime方法etime方法使用etime函数来获得程序运行时间，etime函数的命令格式如下：etime(t1,t0)%返回t1-t0的值,45,2.7关系运算和逻辑运算2.7.1 逻辑运算,MATLAB 中逻辑型（logical）数据只有“1”和“0”，分别表示true和false两种状态，逻辑型变量只占1个字节。函数logical可以用来将数值型转换为逻辑型，任何非零的数值都转换为逻辑1，数值0转换为逻辑0。,46,2.7.2 关系运算,MATLAB 常用的关系操作符有、=、=（等于）、=（不等于）。关系运算规则：如果比较的两个变量都是标量，则结果为1（true）或0（false）；如果比较

24、的两个变量都是数组，则必须尺寸大小相同，结果也是同样大小的数组；如果比较的是一个数组和一个标量，则把数组的每个元素分别与标量比较，结果为与数组大小相同的数组。,47,2.7.3 逻辑运算,1.元素的逻辑运算元素的逻辑运算是将数组中的元素一一进行逻辑运算，常用的逻辑运算符：&（与）、|（或）、（非）和xor（异或）。在逻辑运算中，非0元素表示true，0元素表示false。2.先决逻辑运算先决逻辑运算符有：&（先决与）:A&B|（先决或）:A|B,48,比较两个语句的不同：x=0 y=(x=0)&(100/x10)x=0 y=(x=0)&(100/x10)3.位逻辑运算位逻辑运算函数有：bita

25、nd（位与）、bitor（位或）、bitcmp（位非）和bitxor（位异或）。,练习,49,2.7.4 运算符优先级,各类运算符的优先级为：括号算术运算符关系运算符逻辑运算符各符号优先顺序为：括号()转置.幂.一元加减+-逻辑非 乘*.*除/./.加减+-冒号:关系运算=pi)&(xpi&x4*pi,50,2.8数组的信息获取,1.数组的尺寸,51,例：使用获取数组尺寸的函数获得数组信息 a=rand(2,3)*10a=9.5013 6.0684 8.9130 2.3114 4.8598 7.6210 m,n=size(a)m=2n=3 l=numel(a)l=6 mean=sum(a(:)

26、/lmean=6.5458,52,2.数组的检测函数以“is”开头，函数返回的结果为逻辑型，如果检测符合条件则返回1，不符合条件就返回0。例 利用数组检测函数获取数组信息。a=0 1 2 3 4 5 isempty(a)ans=0,53,2.9 多项式,一个多项式按降幂排列为：p(x)=anxn+an-1xn-1+a1x+a0在MATLAB中用行向量来表示多项式的各项系数，使用长度为n+1的行向量按降幂排列，用0表示多项式中某次幂的缺项，则表示为：p=an an-1 a1 a0例如，p(x)=x3-4x2+3x+1可表示为p=1-4 3 1；p(x)=x2+5x+2x 可表示为p=1 5 2

27、0。,54,2.9.1 多项式求根和求值,1.多项式求根使用roots函数来计算多项式的根，多项式的根以列向量的形式表示；反过来，也可以根据多项式的根使用poly函数获得多项式。2.多项式求值函数polyval和polyvalm可以用来计算多项式在给定变量时的值。,55,例 计算多项式的根并由根得出多项式。p1=1-6 11-6 0 p1=1-6 11-6 0 r1=roots(p1)r1=0 3.0000 2.0000 1.0000 p2=poly(r1(2)r1(3)p2=1.0000-5.0000 6.0000 多项式:p1=x4-6x3+11x2-6x=x(x-3)(x-2)(x-1)

28、p2=x2-5x+6,56,例 计算多项式的值。polyvalm(p1,5)ans=120 x=1 2;3 4x=1 2 3 4 polyvalm(p1,x)ans=48 64 96 144多项式:p1=x4-6x3+11x2-6x=x(x-3)(x-2)(x-1),57,2.9.2 多项式的算术运算,1.多项式的乘法和除法多项式的乘法和除法运算分别使用函数conv和deconv来实现。p=conv(pl,p2)%计算多项式p1和p2的乘积q,r=deconv(pl,p2)%计算多项式p1与p2的商,练习,58,2.部分分式展开B,A=residue(r,p,k)%将部分分式和形式转化为两个多

29、项式除法r,p,k=residue(B,A)%将分母多项式A和分子多项式B进行部分分式展开,59,例2-33 将两个表达式G1进行部分分式展开，a1=1-6 11-6 0;b1=10;r1,p1,k1=residue(b1,a1)%将G1部分分式展开r1=1.6667-5.0000 5.0000-1.6667p1=3.0000 2.0000 1.0000 0k1=,60,3.多项式的微积分使用polyder函数来计算多项式的微分：polyder(p)%计算p的导数polyder(a,b)%计算a和b乘积的导数q,d=polyder(b,a)%计算a和b商的导数MATLAB没有专门的多项式积分函

30、数，但可以通过以下的公式计算完成积分：p./(length(p):-1:1),k%计算多项式p的积分,61,2.9.3 多项式的拟合与插值,1.多项式的拟合多项式拟合是用一个多项式来逼近一组给定的数据，是数据分析上的常用方法。（1）拟合函数p=polyfit(x,y,n)%由x和y得出多项式p说明：x、y向量分别为数据点的横、纵坐标；n是拟合的多项式阶次；p为拟合的多项式，p是n+1个系数构成的行向量。,62,例2-35 使用多项式拟合的方法对曲线 的数据进行拟合，拟合后根据多项式绘制的曲线如图所示。x=0:0.5:10;y=2*sin(x)+sqrt(x);p1=polyfit(x,y,5)

31、%5阶拟合 p2=polyfit(x,y,8)%8阶拟合,63,MATLAB 7.3提供了方便的曲线拟合图形用户接口，可以直接在图形窗口中使用菜单进行多项式拟合，在图形窗口中选择菜单“Tools”“Basic Fitting”，则出现“Basic Fitting”窗口。,64,2.插值运算插值（interpolation）是在两个原始数据点之间根据一定的运算关系插入新的数据点，以便更准确地得出数据的变化规律。一维插值是指对一个自变量的插值，interp1函数用来进行一维插值：yi=interp1(x,y,xi,method),65,例2-36 使用插值运算计算曲线中横坐标为9.5的对应纵坐标值

32、，曲线表达式为。x=1:10;y=2*sin(x)+sqrt(x);y01=interp1(x,y,9.5)%采用线性插值方法计算y01=2.9492 y02=interp1(x,y,9.5,spline)%采用三次样条插值方法计算y02=2.9558 y03=interp1(x,y,9.5,cubic)%采用三次多项式插值方法计算y03=3.0586 y0=2*sin(9.5)+sqrt(9.5)%实际值y0=2.9319,66,练习：,x=uint8(2.3e10)，则x所占的字节是_个。A.1 B.2 C.4 D.8,答案：A,67,练习：,下列_是合法常量。A.3*e10 B.1e50

33、0 C.-1.85e-56 D.10-2,答案：C,68,练习：,已知x=0:10，则x有_个元素。A.10 B.11 C.9 D.12,答案：B,69,练习：,已知数组a=，则a(:,3)是指_ A.所有元素 B.第一行元素C.第三列元素 D.第三行元素,答案：C,70,练习,求解方程组,71,练习：,已知s=显示hello，则s的元素个数是_ A.12 B.9 C.7 D.18,答案：B,72,练习：,运行字符串函数strncmp(s1,s2,2)，则结果是_ A.1 B.0 C.false D.true,答案：B,73,练习：,计算x从0到20间隔为1的行向量，y=sin(x)中，x4范围中y的所有值。,x=0:20;x1=(xpi)&(x4*pi)x=x1.*xy=sin(x),74,练习：,计算三个多项式s1、s2和s3的乘积，则算式为_ A.conv(s1,s2,s3)B.s1*s2*s3C.conv(conv(s1,s2),s3)D.conv(s1*s2*s3),答案：C,75,将传递函数进行部分分式展开：,b=10a=conv(1 3,conv(1 1,1 2)r,p,k=residue(b,a),