24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者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能看出啥？ +8	要乐观耀哥 2026-08-10	23/1150	2026-08-13 00:35 by 要乐观耀哥
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22	Lanmanbaby 2026-08-09	36/1800	2026-08-12 22:50 by sdfapple719
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +10	Tide man 2026-08-11	10/500	2026-08-12 20:59 by 陈晨晨陈啊
[基金申请] 长年满屏的广告，版主太不责任了。基金也等的急 +4	gltch 2026-08-12	4/200	2026-08-12 16:08 by eulota
[基金申请] filecode +13	documentary 2026-08-10	15/750	2026-08-12 15:56 by 云上清扬
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11	yufeiwaner 2026-08-09	13/650	2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果？ +3	yuleib84 2026-08-11	3/150	2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8	Tide man 2026-08-10	9/450	2026-08-11 20:39 by beefly
[基金申请] 什么时候出结果，有咨询渠道？？？ +3	Tide man 2026-08-11	3/150	2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3	布布和一二 2026-08-11	4/200	2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26	布布和一二 2026-08-06	49/2450	2026-08-11 13:18 by WH3796
[基金申请] 确定了，国自然21号放榜 +6	布布和一二 2026-08-10	7/350	2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17	且听虎啸 2026-08-07	24/1200	2026-08-10 18:35 by zhangduo2008
[基金申请] 2026国自然放榜时间 +9	布布和一二 2026-08-08	9/450	2026-08-10 11:22 by xxxx2020
[基金申请] filecode与中标关系的预测 +5	布布和一二 2026-08-07	5/250	2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10	Tide man 2026-08-08	18/900	2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5	yang182083 2026-08-06	7/350	2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7	a089 2026-08-07	7/350	2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4	布布和一二 2026-08-07	7/350	2026-08-07 22:55 by zhanghaozhu
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5	Tide man 2026-08-06	7/350	2026-08-07 13:11 by 医学老男孩