计算药学第2章.ppt

第2章 化学中常用的计算方法 矩阵及其基本运算 线性方程组和回归分析 高次方程的求解 插值和拟合 一、 矩阵及其基本运算 矩阵的基本概念 矩阵的基本运算 MATLAB 程序： 矩阵的建立 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] a=[] b=zeros(2,3) c=ones(2,3) d=eye(2,3) e=eye(3,3) MATLAB 程序： 矩阵元素的操作 a=[1 2 3;4 5 6;7 8 9] a1=a(2,:) a2=a(:,2) a3=a(:) a4=a(1:2,2:3) a5=a(2:-1:1,:) a6=a(:,3:-1:2) a7=a;a7(1:2,:)=[] a8=a;a8(:,1)=[] a9=[a a2] a10=[a;a1] 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) 矩阵的基本运算 二、 线性方程组和回归分析 矩 阵 解 法 正规方程组 （m=n，方阵） AX=B A-1AX=A-1B （ A-1A=In ） InX=A-1B （ InX=X） X=A-1B 应努力建立条件数小的方程组，避免因解病态方程组造成的误差。 由于方程组的条件数取决于系数矩阵，根据研究体系的特征，选择适当的实验点，是避免产生病态方程组的关键。 如计算分光光度法中当各组分光谱完全相同，将得到无解的奇异矩阵；但假如虽然有差别，可差别很小，则条件数必然很大，则将得到病态方程组。 分光光度法中波长的选择十分重要。 一元线性回归及有关计算 一元线性回归 - 二变量间x和y的线性关系 线性相关系数的求算 三、高次方程的求解 迭代法 弦截法 牛顿-雷扶生法 四、插值和拟合 对表2-2 的数据，如仅用末两列数据，则只能用线性插值得到正确的结果， 上式中x应为滴定剂体积V, y为对应的电位二次微商值(△2E/△V2),这里要求结果为0, 计算滴定终点的matlab 命令为  Ve=interp1(x,y,0,‘method’,) 0 可为数组， ‘method’为linear（缺省），nearest(最近邻点插值), cubic (三次插值) polyfit 可广泛用于各种拟合 help polyfit POLYFIT Polynomial curve fitting. POLYFIT(x,y,n) finds the coefficients of a polynomial p(x) of degree n that fits the data, p(x(i)) ~= y(i), in a least-squares sense. [p,S] = POLYFIT(x,y,n) returns the polynomial coefficients p and a matrix S for use with POLYVAL to produce error estimates on predictions. If the errors in the data, y, are independent normal with constant variance, POLYVAL will produce error bounds which contain at least 50% of the predictions. See also POLY, POLYVAL, ROOTS. 设初始剂量为D0,每次注射剂量为D,注射间隔时间为τ, 给药方案为[D0,D,τ]。 血药浓度 应保持在10～25 ug/mL之间, 则 D0 =vc2 = 375 mg, D=v(c2-c1)=225 mg 可制定给药方案首次注射375mg，其余每次注射225mg，注射的间隔时间为4小时。 上机习题 第一部分 MATLAB 基本操作 熟习菜单各部分功能 基本命令练习 demo intro help who whos mkdir dir cd save load clear clf 第二部分： 教材28页 习题1