24小时热门版块排行榜    

查看: 1052  |  回复: 8

wcx1192002

新虫 (初入文坛)

[求助] 模型参数求解

小弟是学药剂的,在做药物释放方面的模型拟合。选用模型如下所示。
想问问怎么求解参数D啊?其中时间t为自变量,Mt/Mo为因变量。L为已知参数。已经测得在不同t下的Mt/Mo的数值了,用matlab如何实现D的求解?请各位大侠帮忙一下,谢谢!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoqing8569

铁杆木虫 (著名写手)

奥林匹亚光学院院长

【答案】应助回帖


wcx1192002(金币+25): 2011-05-04 10:06:13
sunyang1988(金币+1): 谢谢交流 2011-05-05 09:57:27
1 - (2.94 D t)/L^2
2楼2011-05-04 08:30:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wcx1192002

新虫 (初入文坛)

引用回帖:
Originally posted by xiaoqing8569 at 2011-05-04 08:30:19:
1 - (2.94 D t)/L^2

请问一下详细的推导步骤啊,呵呵······
3楼2011-05-04 10:07:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoqing8569

铁杆木虫 (著名写手)

奥林匹亚光学院院长

【答案】应助回帖


sunyang1988(金币+1): 谢谢回帖交流 2011-05-05 09:57:35
引用回帖:
Originally posted by wcx1192002 at 2011-05-04 10:07:49:
请问一下详细的推导步骤啊,呵呵······

我用mathematica计算的。
  
1 - 8 D t /L^2 NSum[1/(2 n + 1)^2 Exp[-(2 n + 1)^2], {n, 0, Infinity}]
4楼2011-05-04 10:32:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wcx1192002

新虫 (初入文坛)

引用回帖:
Originally posted by xiaoqing8569 at 2011-05-04 10:32:45:
我用mathematica计算的。
  
1 - 8 D t /L^2 NSum[1/(2 n + 1)^2 Exp[-(2 n + 1)^2], {n, 0, Infinity}]

请问一下,t和L可以从指数函数里面直接提到外面么?呵呵···
5楼2011-05-04 10:46:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoqing8569

铁杆木虫 (著名写手)

奥林匹亚光学院院长

【答案】应助回帖

★ ★
sunyang1988(金币+2): 谢谢热心帮助 2011-05-05 09:57:43
1 - 8/Pi^2 Exp[-D Pi^2 t/L^2] NSum[
   1/(2 n + 1)^2 Exp[-(2 n + 1)^2 ], {n, 0, Infinity}]
1- 0.298203 E^(-((D Pi^2 t)/L^2))
6楼2011-05-04 12:34:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wcx1192002

新虫 (初入文坛)

引用回帖:
Originally posted by xiaoqing8569 at 2011-05-04 12:34:29:
1 - 8/Pi^2 Exp[-D Pi^2 t/L^2] NSum[
   1/(2 n + 1)^2 Exp[-(2 n + 1)^2 ], {n, 0, Infinity}]
1- 0.298203 E^(-((D Pi^2 t)/L^2))

是不是应该是
1-8/Pi^2Exp[-D Pi^2t/L^2]^{NSum[1/(2n+1)^2Exp[-(2n+1)^2],{n,0,Infinity}}?exp(a*b)应该是[exp(a)]^b才对啊···
7楼2011-05-04 12:53:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoqing8569

铁杆木虫 (著名写手)

奥林匹亚光学院院长

引用回帖:
Originally posted by wcx1192002 at 2011-05-04 12:53:26:
是不是应该是
1-8/Pi^2Exp[-D Pi^2t/L^2]^{NSum[1/(2n+1)^2Exp[-(2n+1)^2],{n,0,Infinity}}?exp(a*b)应该是[exp(a)]^b才对啊···

这样也不可以,能不能提供一点D、L的限定条件?
8楼2011-05-04 16:08:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wcx1192002

新虫 (初入文坛)

引用回帖:
Originally posted by xiaoqing8569 at 2011-05-04 16:08:16:
这样也不可以,能不能提供一点D、L的限定条件?

D大概在1e-8到-9之间,L是0.005。关键是那个级数求和比较麻烦,谢谢!
9楼2011-05-04 16:27:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wcx1192002 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +9 majunge000 2026-08-11 11/550 2026-08-15 23:20 by botar
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见