| 查看: 1121 | 回复: 2 | |||
liyincumt金虫 (职业作家)
小学生
|
[交流]
用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 ] |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有189人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 速求结合常数方法
已经有11人回复
dinghm2012
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 101.4
- 散金: 47
- 帖子: 82
- 在线: 57.6小时
- 虫号: 2141353
- 注册: 2012-11-22
- 性别: MM
- 专业: 有机分子功能材料化学
2楼2014-01-20 09:58:25
戴钢盔的猪头
木虫 (知名作家)
- 应助: 200 (高中生)
- 金币: 5083.9
- 散金: 1310
- 红花: 18
- 帖子: 6001
- 在线: 319.9小时
- 虫号: 2854040
- 注册: 2013-12-07
- 性别: GG
- 专业: 波谱分析与成像分析
3楼2014-01-20 17:21:21
投票:
10