计算药物分析(药学与生物信息学)第二章_图文.pptVIP

* （2）矩阵-矩阵运算 [1] 元素对元素的运算，同数组-数组运算。 3、矩阵的运算 （1）标量-矩阵运算 同标量-数组运算。 [2]矩阵运算： 矩阵加法：A+B 矩阵乘法：A*B 方阵的行列式：det（A） 方阵的逆：inv（A） MATLAB 程序：矩阵元素的操作 a=[1 2 3 4 5 6] b=[1 2 1 2 1 2] c1=a+a c2=a*b c=[2 7 3;3 9 4;1 5 3] c3=det(c) c4=inv(c) 矩阵的基本运算 转置: AT或A’ (矩阵的行与列互换) 对称矩阵（A=AT） 求逆: A，B均为方阵，如AB=BA，则A是可 逆的，记为B=A-1 矩 阵 解 法 正规方程组 （m=n，方阵） AX=B A-1AX=A-1B （ A-1A=In ） InX=A-1B （ InX=In ） X=A-1B 例 解： 超定方程组 （ ） AX=B (ATA) X=ATB (ATA)-1 (ATA) X= (ATA)-1 ATB InX =(ATA)-1 ATB X= (ATA)-1 ATB 矩 阵 解 法 应用示例：（光谱分析中的多组分测定 ） A=ECL A=A1+A2+…+An 应努力建立条件数小的方程组，避免因解病态方程组造成的误差。 由于方程组的条件数取决于系数矩阵，根据研究体系的特征，选择适当的实验点，是避免产生病态方程组的关键。 如计算分光光度法中当各组分光谱完全相同，将得到无解的奇异矩阵；但假如虽然有差别，可差别很小，则条件数必然很大，则将得到病态方程组。 分光光度法中波长的选择十分重要。 二、 线性方程组和回归分析 克莱姆公式 高斯消去法 矩阵解法 矩阵解法 系数矩阵-A 未知数矩阵-X 常数矩阵-B 矩阵形式：AX=B X=inv(A)*B=A\B 一元线性回归及有关计算 一元线性回归 - 二变量间x和y的线性关系 y x * * * * * 线性相关系数的求算 总 变 差 平 方 和 0 S=Q(残差平方和)+U(回归差平方和) 三、高次方程的求解 迭代法 三、高次方程的近似求解 对方程f(x) = 0 求近似解，使f(x*)≈0 设初值x0，按一定规则生成新值x1 依次计算生成数列: x0,x1,x2,x3……xn lim xn =x* │ xn-xn-1│ ε 弦截法 基本原理 迭代通式 收敛指标 牛顿-雷扶生法的基本思想 设 是f(x)=0的一个近似根，把f(x)在 处作泰勒展开 若取前两项来近似代替f(x)(称为f(x)的线性化)，则得近似的线性方程 设 ，令其解为 ，得 这称为f(x)=0的牛顿迭代格式。 牛顿-雷扶生法 切线逼近法 特点：一个初始值；收敛速度快，求根方便 x0=1, x1 = 1- (1-2-5)/(3*1-2) =7 x2=4.7655 x3=3.3487 x4=2.5316 x5=2.1739 x6=2.09788 x7 =2.094552 x8=2.094552 注意点： 1、 收敛标准： 2、 初始值： 3、 可能有多个解 设函数 在区间 上有定义，且已知在点 上的 函数 ， 成立，就称 为 的插值函数，点 称为插 值节点，包含节点的区间 称为插值区间，求插值函数 若存在一简单 使 的方法称为插值法。 四、插值和拟合 若 是次数不超过 的代数多项式， 其中 为实数，就称 为插值多项式， 若 为分段的多项式，就称为分段插值. 若 为三角多项式，就称为三角插值. 即 相应的插值法称为多项式插值. 从几何上看，插值法就是就曲线 ，使其通过 给定的 个点 ，并用它近似已知曲线 . 见图 对表2-2 的数据，如仅用末两列数据，则只能用线性插值得到正确的结果。 上式中x应为滴定剂体积V, y为对应的电位二次微商值(△2E/△V2)，这里要求结果为