Matlab数值向量和数组幻灯片.ppt

多项式的求值 * 第4章 数值向量和数组 * MATLAB 7提供了两个函数来对多项式进行求值，即polyval和polyvalm。前者以数组为计算单位，后者以矩阵为计算单位。 p=[1.0000 -20.0000 -16.0000 480.0000 98.0000] p = 1 -20 -16 480 98 x=4 x = 4 polyval(p,x) ans = 738 求多项式的根 * 第4章 数值向量和数组 * 在MATLAB 7语言里，多项式由一个行向量表示，设为p，它的系数按降序排列，使用roots函数可以求出该多项式的根。其使用格式为roots(p)。 p=[1 0 3 12 -7] p = 1 0 3 12 -7 roots(p) ans = 0.7876 + 2.4351i 0.7876 - 2.4351i -2.0872 0.5121 多项式的四则运算 * 第4章 数值向量和数组 * 加法和减法 如果两个多项式的向量阶数相同，标准的数组加法有效。当两个多项式的向量阶数不同时，需要在低阶多项式的前边补0，使得它与相加的高阶多项式有相同的阶数。 a=[8 2 2 8],b=[6 1 6 1] a = 8 2 2 8 b = 6 1 6 1 c=a+b c = 14 3 8 9 Y3=poly2sym(c) Y3 = 14*x^3+3*x^2+8*x+9 多项式的四则运算 * 第4章 数值向量和数组 * 乘法 使用conv函数对多项式进行乘法运算。 格式为c=conv(a,b)，其中a和b为两个多项式的系数向量，c为相乘所生成的多项式的系数向量。 a=[1 2 3 4],b=[5 6 7 8]; Y1=poly2sym(a) Y1 = x^3+2*x^2+3*x+4 Y2=poly2sym(b) Y2 = 5*x^3+6*x^2+7*x+8 c=conv(a,b); Y=poly2sym(c) Y = 5*x^6+16*x^5+34*x^4+60*x^3+61*x^2+52*x+32 多项式的四则运算 * 第4章 数值向量和数组 * 除法 在数值计算中，经常需要用一个多项式去除另一个多项式。在MATLAB 7语言中，使用decon函数来完成该项功能。 a=[1 2 3 4],b=[5 6 7 8]; c=conv(a,b); d=deconv(c,a); e=deconv(c,b) x=poly2sym(d),y=poly2sym(e) x = 5*x^3+6*x^2+7*x+8 y = x^3+2*x^2+3*x+4 多项式的四则运算 * 第4章 数值向量和数组 * 求导和积分 在MATLAB 7语言中，分别使用polyder函数和polyint函数来求多项式的导数与积分。 p=[3 1 8 8] p = 3 1 8 8 q=polyder(p) q = 9 2 8 p1=polyint(q) p1 = 3 1 8 0 关系操作符 逻辑操作符 关系与逻辑函数 NaNs和空矩阵 各种运算符的优先级 * 第4章 数值向量和数组 * 4.关系和逻辑运算 (1) 关系操作符 关系运算符 该运算符的功能 关系运算符 该运算符的功能 ? 小于 = 大于等于 ? 大于 == 等于 = 小于等于 ~= 约等于 * 第4章 数值向量和数组 * 关系操作符及其功能 (2) 逻辑操作符 逻辑操作符 功 能 与 | 或 ~ 非 * 第4章 数值向量和数组 * 逻辑操作符及其功能 (3) 关系与逻辑函数 关系与逻辑函数及其功能 关系和逻辑函数 使 用 功 能 xor(s,t) 异或运算，s或t非零(真)返回1，s和t都是零(假)或都是非零(真)返回0 any(x) 如果在一个向量x中，任何元素是非零，返回1；矩阵x中的每一列有非零元素，返回1 all(x) 如果在一个向量x中，所有元素非零，返回1；矩阵x中的每一列所有元素非零，返回1 * 第4章 数值向量和数组 * (4) NaNs和空矩阵 * 第4章 数值向量和数组 * NaNs和空矩阵([ ]) 在MATLAB 7中作特殊处理: 根据IEEE数学标准。对NaNs的几乎所有运算结果都得出NaNs。 空矩阵由MATLAB 7的