C3_class.ppt

原则： ? 用[ ]将矩阵内容括起来 ? 每行的元素用空格或逗号分开 ? 行间用分号或回车键表结束 用matlab函数创建特殊矩阵: zeros m，n ：全部元素都为0的矩阵 ones m，n ：全部元素都为1的矩阵 eye m，n ：单位矩阵 rand m，n ：随机矩阵 magic n ：n维魔方阵 pascal n ： n维pascal矩阵 ?启动M文件编辑器，输入待创建矩阵； ?将该文件保存为*.m文件； ?在MATLAB命令窗口运行该文件，即可生成文件中包含的矩阵。 a.\b b./a 第3章练习题 （3）分别建立3×2、与矩阵A同样大小的零矩阵。 （4）利用M文件建立矩阵B，其变量名为B，M文件名为mydata.m。 （2）创建如下矩阵 A 和 B，对其作关系运算，标识出两矩阵中元素相等的位置，元素值不等的位置。 （1）求当 x 1, y 2 时如下表达式的值。 （5）求解如下两个线性方程组： （7）已知多项式 求多项式根。 （6）已知多项式 ，求二者乘积。 （8）已知多项式 求其导数。 （9）已知数据： （0, -0.447）， （0.1, 1.978）， （0.2, 3.11）， （0.3, 5.25）， （0.4, 5.02），（0.5, 4.66）， （0.6, 4.01）， （0.7, 4.58）， （0.8, 3.45）， 0.9, 5.35 ， 1, 9.22 ，拟合其三次多项式。 The End 数组运算符 使用原则： 数组必须同维，或其中的某个数组为标量； 数组运算的结果仍是数组，与原数组同维。 a [1,2;3,4;5,6]; b [2,3;4,5;6,7]; a*b Error using mtimes Inner matrix dimensions must agree. a.*b ans 2 6 12 20 30 42 a./b b.\a A.^p A.^B p.^A A、B为数组，p为标量 数组的乘方 3.2.3 矩阵与数组函数 矩阵函数 矩阵的秩 rank 矩阵的逆 矩阵的特征值及特征向量 inv eig 矩阵伪逆 矩阵正交化运算 pinv orth 方阵的行列式值 det 函数功能 函数名 矩阵函数 [R,p] chol x [Q,R] qr A 矩阵分解函数 [L,U] lu A 数组函数 三角函数 指数函数 对数函数 复数运算函数 多项式的表达 3.3.1 多项式的创建和表达 把多项式表达成一个行向量 : ①向量中的元素按多项式降幂排列 ②若多项式中缺某幂次项，用0填补该系数 多项式的行向量表达原则: 直接输入向量构造多项式 P [2,5,0,4,1,4]; P_1 poly2sym P P_1 2*x^5+5*x^4+4*x^2+x+4 P_2 poly2str P,x P_2 2 x^5 + 5 x^4 + 4 x^2 + x + 4 whos P_1 P_2 Name Size Bytes Class P_1 1x1 166 sym object P_2 1x33 66 char array 用函数poly创建多项式 poly A 创建的是多项式 的系数向量。 a [2,-4,7]; P poly a P 1 -5 -22 56 Ps poly2str P,t Ps t^3 - 5 t^2 - 22 t + 56 3.3.2 多项式运算 多项式四则运算 加、减运算符： 加减运算的多项式向量大小要相等；多项式阶次不同时，低阶多项式首零 填补到高阶多项式的阶次。 多项式四则运算 乘运算： C conv A , B 执行两个向量A、B的卷积； 若向量A、B的长度分别为m、n，则C的长度为m+n-1 多项式四则运算 除运算： [ Q , R ] deconv B , A 执行两个向量B、A的解卷积 B conv A , Q + R 多项式导数运算 polyder P P为多项式的行向量表示 多项式求根 roots P P为多项式的行向量表示 根以列向量形式保存 多项式求值 polyval P, x P为多项式的行向量表示 X为矩阵或向量 3.3.3 多项式拟合和插值 多项式拟合 P polyfit X, Y, N X,Y为已知数据 N为要拟合的多项式阶次 练习： 已知5个点： （1, 3.2）， （2, 46.7）， （3, 138.3），（4, 278.6）， （5, 467.