tlab：一维数组及其应用.pptVIP

数学物理方法（陆全康 赵惠芬，第二版，高等教育出版社，2003）P87 * 8、多项式的除法 [q,r] = deconv(b,a ) q：用来存储b/a的商多项式， r：用来存储b/a的余数多项式 余数通常都与除数和被除数中较长的那个向量等长，因此r的前面一般有若干个0 此函数表示多项式b除以多项式a得到商多项式q和余数多项式r，如果r的元素全部为0，则表示多项式b可以整除多项式a。 多项式的除法就是多项式系数向量之间进行解卷积运算 * 练习 * 9、多项式的微分（求导） polyder( ) 多项式微分 k = polyder(p) 求多项式p的微分多项式 k = polyder(a,b) 求a x b的微分 [q,d] = poyder(b,a)：求有理分式b/a的微分 q：分子多项式 d：分母多项式 * 10、多项式的积分 p = [1,6,20,48,69,72,44] h = polyder(g) q = polyint(h,44) q = polyint(p,k) 对多项式p积分，k为积分常数 q = polyint(p) 对多项式p积分，积分常数为0 * 11、部分分式展开（留数计算） 找出两个多项式b(x)和a(x)之比b(x)/a(x)的留数、极点及直项向量。 将部分分式展开的形式还原为两个多项式b(x)和a(x)相除的形式 [r,p,k] = residue(b,a) [b,a] = residue(r,p,k) * * 求函数在指定点处的留数 * 数学物理方法（陆全康 赵惠芬，第二版，高等教育出版社，2003）P130 * 数学物理方法（陆全康 赵惠芬，第二版，高等教育出版社，2003）P157 * 数学物理方法（陆全康 赵惠芬，第二版，高等教育出版社，2003）P160 * 数学物理方法（陆全康 赵惠芬，第二版，高等教育出版社，2003）P161 * 12、多项式的矩阵运算 p是多项式的系数向量，相当于用矩阵X代替多项式的变量来对矩阵进行计算，X必须为方阵。 X = pascal(4) p = [1,2,3] polyval(p,X) polyvalm(p,X) Y = polyvalm(p,X) * 13、多项式拟合 p = polyfit(x,y,n) 多项式曲线拟合（最小二乘法） n=1就是进行线性拟合 * x y 0 0.2 0.6 1 1.3 1.6 1.7 1.8 1.9 2.2 2.3 2.5 2.6 2.9 0 -2.5 -4 -5.7 -3.5 -2 -1 2 3.5 4 7 7.5 9.9 10.9 x y 3.1 3.4 3.8 4.1 4.4 4.7 4.8 4.9 5 5.1 5.3 11.9 13.5 13 11.9 9 6.5 4 1.5 0 -2.5 -5 * 多项式拟合时的注意事项 在进行曲线拟合时对多项式阶次的选择是任意的。 虽然高阶的多项式可以更准确地拟合数据（仅指该曲线与给定的数据之间的均方误差最小），但在进行曲线拟合时，并不需要采用太高阶的多项式，这主要基于以下原因： 1、越是高阶的多项式其数值特性越差，计算起来也越耗时； 2、随着多项式阶次的升高，拟合的曲线变的越来越不平滑，通常会出现用户不愿意看到的局部波形； 3、由于数据本身的近似性，因此在进行数据拟合时没有必要仅仅考虑使拟合的曲线无限接近数据点，而要在曲线的阶次合均方误差之间综合考虑，因为越是高阶的多项式在物理实现时越困难。 从数学原理上，n+1个数据点可以惟一定义一个n次曲线（或n阶多项式）。 * 一维数组相关函数的应用以线性拟合为例 n T(℃) R(Ω) 1 19.1 76.30 2 25.0 77.80 3 30.1 79.25 4 36.0 80.80 5 40.0 82.35 6 45.1 83.90 7 50.0 85.10 测得铜导线在温度Ti下的电阻为Ri，编写一个M脚本文件，用matlab提供的数组函数对数据进行线性拟合，求出斜率、截距和相关系数，并与origin的拟合结果进行比较。 * x=[0:0.2:2*pi]; plot(x,cos(x),’k*-’); 二、一维数组在二维绘图中的应用 1、plot(a,b,string) a, b 都是一维数组，以 a 中元素为横坐标， b 中元素为纵坐标作平面曲线。 a, b 必须具有相同长度。string 是用 单引号 括起来的字符串，用来指定图形的属性（点、线的形状和颜色） 黑色、实线离散点用星号 属性可以全部指定，也可以只指定其中某几个。 排列顺序任意 * 2、线型、标记和颜色 线型 点标记 颜色 - 实线 : 点线 -. 点划线 -- 虚线 . 点 o 小圆圈 x 叉号