24小时热门版块排行榜

>论坛更新日志 (2949)
>虫友互识 (186)
>论文投稿 (180)
>文献求助 (162)
>基金申请 (152)
>导师招生 (115)
>考研 (89)
>休闲灌水 (81)
>学术会议 (74)
>硕博家园 (74)
>博后之家 (65)
>论文道贺祈福 (51)
>教师之家 (48)
>公派出国 (43)
>考博 (38)
>科研资料 (37)

返回列表

【悬赏金币】回答本帖问题，作者Evoly_z将赠送您 88 个金币

wlfc

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 77
帖子: 19
在线: 12小时
虫号: 26920802
注册: 2021-08-01

【答案】应助回帖

用OpenLu求解，Lu脚本代码：
!!!using["luopt","math"]; //使用命名空间
f(x,y,dy, params::k1,k2)=
{
dy=k1*(0.701*exp(-3.211*x)+0.299*exp(-0.067*x))-k2*y,
0 //必须返回0
};
目标函数(_k1,_k2 : i,s,tyz : tyArray,tA,max,k1,k2)=
{
k1=_k1, k2=_k2, //传递优化变量
//最后一个参数50表示gsl_ode函数在计算时,最多循环计算50次,这样可以提高速度
tyz=gsl_ode[@f,nil,0.0,tA,ra1(0), 1e-6, 1e-6, 2, 1e-6,50],
i=0, s=0, while{++i<max, s=s+[tyz(i,1)-tyArray(i,1)]^2},
s
};
main(::tyArray,tA,max)=
{
tyArray=matrix{ //存放实验数据xi,yi
      "0       0
0.33    0.061043
1       0.03675
2       0.05932
4       0.095993
6       0.072057
10       0.05085
15       0.04678
24       0.047673
30       0.034973
48       0.030375
72       0"
},
len[tyArray,0,&max], tA=tyArray(all:0), //用len函数取矩阵的行数，tA取矩阵的列
Opt1[@目标函数] //Opt1函数全局优化
};

结果：
0.1325616120497096       0.358508535326814       3.242596823546582e-003

赞一下

回复此楼

11楼2021-08-09 21:27:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 Evoly_z 的主题更新

返回列表

不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)

普通表情龙兔虎猫

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] Filecode 又变了，巨变 +3	WH3796 2026-08-12	4/200	2026-08-13 14:13 by 小木虫6752397
[基金申请] 不应该看fileCode +6	且听虎啸 2026-08-12	8/400	2026-08-13 12:07 by 小木虫6752397
[基金申请] 关于Filecode分析方法 +9	majunge000 2026-08-10	11/550	2026-08-13 10:40 by ZH远山
[基金申请] 重要来源：本周末出结果 +9	瞬息宇宙 2026-08-12	9/450	2026-08-13 10:12 by Vivilian
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +9	家与远方 2026-08-10	14/700	2026-08-13 09:41 by 臭臭不臭01
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3	Tide man 2026-08-12	4/200	2026-08-13 08:35 by ZJTJZ
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22	Lanmanbaby 2026-08-09	36/1800	2026-08-12 22:50 by sdfapple719
[基金申请] filecode +13	documentary 2026-08-10	15/750	2026-08-12 15:56 by 云上清扬
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4	工厂打螺丝 2026-08-08	4/200	2026-08-12 09:18 by hansi2025
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +7	布布和一二 2026-08-10	10/500	2026-08-12 09:12 by 布布和一二
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8	Tide man 2026-08-10	9/450	2026-08-11 20:39 by beefly
[基金申请] 帮忙看看fileCode +7	wwncly 2026-08-10	13/650	2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果，有咨询渠道？？？ +3	Tide man 2026-08-11	3/150	2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题，想知道的进来 +17	且听虎啸 2026-08-07	24/1200	2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4	Vierhys 2026-08-10	8/400	2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11	布布和一二 2026-08-08	22/1100	2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5	布布和一二 2026-08-07	5/250	2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode +4	布布和一二 2026-08-07	7/350	2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了，今天 +6	archvillain 2026-08-06	10/500	2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6	布布和一二 2026-08-07	22/1100	2026-08-07 14:45 by 且听虎啸