[发明专利]一种利用模拟分子动力学方法确定后装源驻留时间的方法有效

专利信息
申请号: 201410714374.X 申请日: 2014-11-28
公开(公告)号: CN104383639A 公开(公告)日: 2015-03-04
发明(设计)人: 勾成俊;吴章文;侯氢 申请(专利权)人: 四川大学
主分类号: A61N5/00 分类号: A61N5/00
代理公司: 成都顶峰专利事务所(普通合伙) 51224 代理人: 李崧岩
地址: 610000 四*** 国省代码: 四川;51
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种利用模拟分子动力学方法确定后装源驻留时间的方法,用于解决现有技术因临床情况的复杂性,而导致驻留位较多、剂量参考点不能与驻留位一一对应等问题。该方法包括以下步骤:(一)根据驻留点和剂量节制点的位置关系,采用下式算出每个驻留点在剂量节制点处的剂量:(二)根据剂量节制点处的计算剂量与处方剂量的偏差和相邻驻留点的时间差构成评价函数:(三)根据上式对时间tk求偏导数,得到式(3):(四)将上式近似为一个分子动力学方程,得到分子间的相互作用力fk(五)采用下式进行迭代,迭代过程中,通过不断改变τ的值,得出不同的驻留时间。
搜索关键词: 一种 利用 模拟 分子 动力学 方法 确定 后装源 驻留 时间
【主权项】:
一种利用模拟分子动力学方法确定后装源驻留时间的方法,其特征在于,包括以下步骤;(I)根据驻留点和剂量节制点的位置关系,采用式(1)算出每个驻留点在剂量节制点处的剂量:<mrow><msub><mi>d</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>(II)根据剂量节制点处的计算剂量与处方剂量的偏差和相邻驻留点的时间差构成评价函数:<mrow><mi>O</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>d</mi><mi>j</mi><mi>i</mi></msubsup><msub><mi>t</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>d</mi><mi>j</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><munder><mi>&Sigma;</mi><mi>i</mi></munder><mi>&alpha;</mi><msup><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>t</mi><mi>i</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>(III)根据式(2)对时间tk求偏导数,得到式(3):<mrow><mfrac><mrow><mo>&PartialD;</mo><mi>O</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msub><mi>t</mi><mi>k</mi></msub></mrow></mfrac><mo>=</mo><mn>2</mn><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>d</mi><mi>j</mi><mi>i</mi></msubsup><msub><mi>t</mi><mi>i</mi></msub><msubsup><mi>d</mi><mi>j</mi><mi>k</mi></msubsup><mo>-</mo><mn>2</mn><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><msubsup><mi>d</mi><mi>j</mi><mn>0</mn></msubsup><msubsup><mi>d</mi><mi>j</mi><mi>k</mi></msubsup><mo>+</mo><mn>2</mn><mi>&alpha;T</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>(IV)将式(3)类似于一个分子动力学方程,分子间的相互作用力fk表示如下:<mrow><mo>-</mo><msub><mi>f</mi><mi>k</mi></msub><mo>=</mo><mfrac><mrow><mo>&PartialD;</mo><mi>O</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><msub><mi>t</mi><mi>k</mi></msub></mrow></mfrac><mo>=</mo><mn>2</mn><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>d</mi><mi>j</mi><mi>i</mi></msubsup><msub><mi>t</mi><mi>i</mi></msub><msubsup><mi>d</mi><mi>j</mi><mi>k</mi></msubsup><mo>-</mo><mn>2</mn><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><msubsup><mi>d</mi><mi>j</mi><mn>0</mn></msubsup><msubsup><mi>d</mi><mi>j</mi><mi>k</mi></msubsup><mo>+</mo><mn>2</mn><mi>&alpha;T</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>(V)采用下式(5)和(6)进行迭代,迭代过程中,通过不断改变τ的值,得出不同的驻留时间:<mrow><msub><mi>v</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>+</mo><mi>&delta;&tau;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>v</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>m</mi><mi>k</mi></msub></mrow></mfrac><mi>&delta;&tau;</mi><mo>[</mo><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>+</mo><mi>&delta;&tau;</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>t</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>+</mo><mi>&delta;&tau;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>t</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><mi>&delta;&tau;</mi><msub><mi>v</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msub><mi>m</mi><mi>k</mi></msub></mrow></mfrac><mi>&delta;</mi><msup><mi>&tau;</mi><mn>2</mn></msup><msub><mi>f</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>(VI)当O(t)小于设定值时或者迭代次数达到设定值时退出迭代过程,得出最终的驻留时间解;其中,M是剂量节制点的数目,N是驻留点的数目,ωj是第j个剂量节制点的权重,ti是第i个驻留点的驻留时间,是第i驻留点对第j个剂量节制点的剂量率,是第j个剂量节制点的节制剂量,α表示驻留时间平滑因子,vk(τ)是第k个原子在τ时刻的速度,δτ是时间增量。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于四川大学,未经四川大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201410714374.X/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top