应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 基础数学 » 【求助】用matlab中的高斯牛顿法或单纯形法进行模型嵌合
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 759  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hht666

金虫 (正式写手)

[交流] 【求助】用matlab中的高斯牛顿法或单纯形法进行模型嵌合

各位高手, 我最近碰到一问题,需要分析人体内的药物含量随时间变化,需要用上述方法进行模型嵌合, 只知道matlab 软件可以计算，想自己来时间不够，只能求助各位高手。
t        X
0.5        29701602
1        28952373
2        27578765
3        26202931
4        24616875
5        23031635
6        21218045
7        19869353
8        18343096
9        16833890
10        15504818
11        14335500
12        13110579
13        12048990
14        10987643
模型公式如下：X=A×e∧(-a×t)+B× e∧(-b×t)+C×e∧(-c×t)+D× e∧(-d×t)
  A，a, B,b,C,c, D,d 需满足如下条件。
A+B+C+D=29792400， A×a+B×b+C×c+D×d=0,
a,b,,c,d 均大于零，且 b>a,b>c,b>d;
请高手解决，30金币敬请笑纳。

[ Last edited by 点点风 on 2009-3-21 at 12:44 ]
回复此楼

» 猜你喜欢
1楼 2009-03-20 15:20:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hht666

金虫 (正式写手)
我用残差法,假设在T 足够大时,前面三项都趋近于零,就可以得出那把D 和  d固定下来.D=EXP(17.449)=37844643, d=0.0882.请再算算情况怎么样.

[ Last edited by hht666 on 2009-3-21 at 13:54 ]
赞一下 回复此楼
7楼2009-03-21 13:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

fspdlh

金虫 (正式写手)
★ ★ ★ ★ ★
hht666(金币+5,VIP+0):多谢啊.有一个条件不满足,请帮忙再算算 3-21 12:02
以下是我用非线性拟合做的，结果与几个量的初值关系比较大，你能不能说一下各值大约在什么范围内啊，或者数据再多些，因为待定参数太多，相对来说数据有点少
abcd =

    0.0631    0.9828    0.0593    0.0650


ABCD =

  1.0e+009 *

    9.3546   -0.0019   -3.1108   -6.2121


r =

    0.9999



abcd =

    0.0949  184.3219   37.3425    0.2386


ABCD =

  1.0e+007 *

    4.3306    0.0086   -0.0452   -1.3148


r =

    0.9999

[ Last edited by fspdlh on 2009-3-20 at 19:53 ]
赞一下(6人) 回复此楼
2楼2009-03-20 19:50:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)
★ ★
Doctorcbw(金币+2,VIP+0):谢谢参与 12-8 15:07
这是程序
datafit.m中：
function [abcd ABCD r ]=datafit(beta)
load data.txt
t=data(:,1);
X=data(:,2);
beta=nlinfit(t,X,'fun',beta);
a1=beta(1);
delta_b=beta(2);
c1=beta(3);
d1=beta(4);
A=beta(5);
B=beta(6);
a=a1^2;
c=c1^2;
d=d1^2;
b=max(max(a,c),d)+delta_b^2;
C=-(-B*d+A*a+B*b+29792400*d-A*d)/(c-d);
D=(-c*B+A*a+29792400*c-c*A+B*b)/(c-d);
X1=fun(beta,t);
r=corrcoef(X,X1);
r=r(1,2);
abcd=[a b c d];
ABCD=[A B C D];

fun.m中：
function X=fun(beta,t)
a1=beta(1);
delta_b=beta(2);
c1=beta(3);
d1=beta(4);
A=beta(5);
B=beta(6);
a=a1^2;
c=c1^2;
d=d1^2;
b=max(max(a,c),d)+delta_b^2;
C=-(-B*d+A*a+B*b+29792400*d-A*d)/(c-d);
D=(-c*B+A*a+29792400*c-c*A+B*b)/(c-d);
X=A.*exp(-a.*t)+B.*exp(-b.*t)+C.*exp(-c.*t)+D.*exp(-d.*t);

data.txt中：
0.5        29701602
1        28952373
2        27578765
3        26202931
4        24616875
5        23031635
6        21218045
7        19869353
8        18343096
9        16833890
10        15504818
11        14335500
12        13110579
13        12048990
14        10987643
赞一下(6人) 回复此楼
3楼2009-03-20 19:58:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)
★
点点风(金币+1,VIP+0):感谢参与讨论，加分鼓励！ 3-21 12:45
我编的函数的参数为
sqrt(a) sqrt(b-max(a,c,d)) sqrt(c) sqrt(d) A B 六个
由这六个参数来确定a b c d A B C D

[ Last edited by fspdlh on 2009-3-20 at 20:02 ]
赞一下(10人) 回复此楼
4楼2009-03-20 19:59:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-01 21:48 by chyhaha
[考研] 070300一志愿211，312分求调剂院校 +14 小黄鸭宝 2026-03-30 14/700 2026-04-01 20:19 by 赖春艳
[考研] 调剂 +3 好好读书。 2026-04-01 3/150 2026-04-01 17:06 by zhouyuwinner
[考研] 0703求调剂 +4 zizimo 2026-03-31 4/200 2026-04-01 16:04 by yanflower7133
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +14 丹青奶盖 2026-03-26 16/800 2026-04-01 15:58 by oooqiao
[考研] 0856，材料与化工321分求调剂 +13 大馋小子 2026-03-27 14/700 2026-04-01 15:58 by yanflower7133
[考研] 314求调剂 +6 1xiaojun23 2026-03-31 6/300 2026-04-01 11:28 by 今天星期八LQ
[考研] 材料与化工调剂一志愿大连海事085600，349 +9 吃的不少 2026-03-30 9/450 2026-04-01 11:24 by wangjy2002
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk￥
[考研] 289求调剂 +7 BrightLL 2026-03-29 7/350 2026-03-31 22:05 by 544594351
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 329求调剂，一志愿西北工业大学，材料工程（085601） +6 小小机灵虫 2026-03-29 12/600 2026-03-31 16:58 by 记事本2026
[考研] 286求调剂 +6 Faune 2026-03-30 6/300 2026-03-31 14:37 by jp9609
[考研] 266分，求材料冶金能源化工等调剂 +8 哇呼哼呼哼 2026-03-27 10/500 2026-03-31 13:35 by Huaxue_Wang
[考研] 一志愿大连理工大学，机械工程学硕，341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
[考研] 085600 286分 材料求调剂 +11 麻辣鱿鱼 2026-03-27 12/600 2026-03-30 19:33 by Wang200018
[考研] 材料专硕 085600求调剂 +7 BBQ233 2026-03-30 7/350 2026-03-30 17:44 by oooqiao
[考研] 327求调剂 +6 汲亦昊 2026-03-29 6/300 2026-03-29 13:40 by peike
[考研] 085602 化工专硕 338分 求调剂 +12 路痴小琪 2026-03-27 12/600 2026-03-28 15:41 by L135790
信息提示
关闭
请填处理意见
关闭