24小时热门版块排行榜

返回列表

liyincumt

金虫 (职业作家)

小学生

应助: 4 (幼儿园)
金币: 22263.2
散金: 730
红花: 5
帖子: 3223
在线: 664.1小时
虫号: 379584
注册: 2007-05-22
性别: GG
专业: 理论和计算化学

[交流] 用Matlab求解主客体相互作用的稳定常数已有2人参与

请教有哪位大虾做过主客体分子的稳定常数的计算？
目前市场上虽然有现成的软件hyperquad，但是价格昂贵得惊人。小弟虽不才，却想用matlab实现稳定常数的计算，但一直苦于对模型的线性化处理。
实验条件：
UV-VIS紫外可见光光谱仪。
其中纯主体分子在紫外区有吸收，摩尔吸光系数可测得；客体分子无紫外吸收；当主客体结合后，吸光度会增加，即络合物的摩尔吸光系数升高。公式表达如下：
H+G->HG
稳定常数计算
K=HG/[(H0-HG)*(G0-HG)]
其中H0表示初始主体分子浓度，G0表示初始客体分子浓度。
络合物HG浓度计算如下：
HG=0.5*(H0+G0+1/K-((H0+G0+1/K)^2-4*H0*G0)^0.5)
(HG还有另一个解，与实际不符，故舍弃)
εh表示主体分子摩尔吸光率；εhg表示络合物分子摩尔吸光率；l表示光程长度，A表示总吸光度。

目标：建立总吸光度A（因变量）与初始主体分子浓度H0（自变量1）和初始客体分子浓度G0（自变量2）之间关系，其中已知参数是主体分子吸光率εh与光程l，待估参数是络合物吸光率εhg与稳定常数K。
图片显示的就是该模型。
残差U=Sum[实际观测值Aobs-理论计算值A=f(H0,G0,εh,εhg,l,K)]^2
采用Newton-Gauss进行非线性最小二乘进行迭代，使得残差U最小，并且计算出K值与εhg值。

问题：如要采用matlab计算的话，得线性化该模型，从而建立起Y=XB的关系。但是模型中有一个根式，不知该如何处理？
欢迎相关课题的童鞋们一起进行讨论。

附件是一篇关于紫外、荧光得到稳定常数的文献。

[ Last edited by liyincumt on 2011-8-4 at 19:47 ]

回复此楼