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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 1stopt 语句修改
242/3返回列表上一页123下一页
查看: 1616  |  回复: 23
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xxz903

金虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;
赞一下 回复此楼
11楼2014-12-08 09:15:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
引用回帖:
11楼: Originally posted by xxz903 at 2014-12-08 09:15:10
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;

您好。我的是高版本的，不能运行，请问a值如何体现？不用体现？
赞一下 回复此楼
12楼2014-12-08 10:50:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
引用回帖:
11楼: Originally posted by xxz903 at 2014-12-08 09:15:10
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;

请指导，我运行结果都是0.   a是初值，y是某一点（t=2）函数值，咋么写程序嫩狗狗体现出来这些。请高手指点
赞一下 回复此楼
13楼2014-12-08 10:56:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
chengxl: 金币+2 2014-12-09 08:33:43
引用回帖:
11楼: Originally posted by xxz903 at 2014-12-08 09:15:10
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;

Parameters c, k;
Variable t, y;
Function  y=c*(exp(k*t)-1)+1;
data;
0  1
1  3
2  4
3  5
// ----------------------------------------
"Type your title here"

迭代数: 66
计算用时(时:分:秒:微秒): 00:00:00:94
优化算法: 粒子群算法PSO
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0895640117970001
残差平方和(SSE): 0.0320868488366927
相关系数(R): 0.998194643831721
相关系数之平方(R^2): 0.996392546974336
决定系数(DC): 0.996332931561521
卡方系数(Chi-Square): 0.00435149314601671
F统计(F-Statistic): 552.407773509935

参数        最佳估算
----------        -------------
c                  -5.45183726787717
k                  -0.428004961121383
====== 结果输出 =====
No        实测值y        计算值y
1        1        1.0000000
2        3        2.8982851
3        4        4.1356029
4        5        4.9420968
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

chengxl
14楼2014-12-08 15:08:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
送红花一朵
引用回帖:
14楼: Originally posted by feixiaolin at 2014-12-08 15:08:18
Parameters c, k;
Variable t, y;
Function  y=c*(exp(k*t)-1)+1;
data;
0  1
1  3
2  4
3  5
// ----------------------------------------
"Type your title here"

迭代数: 66
计算用时( ...

y'=k*y
已知不同初始条件下（t=0）初始值及其对应的t=10时函数值，如何用1stopt 语句拟合k值。下列语句运行不了，请高手指导，感激不尽。假定函数不能积分得到解析解。
Parameters k;
Variables a,y;
Function y=int(k*y,t=0,10)+a;
RowData;
0 1 2 3;
1 3 4 5;
赞一下 回复此楼
15楼2014-12-08 19:40:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
15楼: Originally posted by chengxl at 2014-12-08 19:40:52
y'=k*y
已知不同初始条件下（t=0）初始值及其对应的t=10时函数值，如何用1stopt 语句拟合k值。下列语句运行不了，请高手指导，感激不尽。假定函数不能积分得到解析解。
Parameters k;
Variables a,y;
Function ...

我的版本也很低，跑不了。
赞一下 回复此楼
16楼2014-12-08 20:35:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chengxl

木虫 (小有名气)
引用回帖:
16楼: Originally posted by feixiaolin at 2014-12-08 20:35:10
我的版本也很低，跑不了。...

谢谢了，麻烦您了
回复此楼
17楼2014-12-09 08:33:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雕刻你的拾光

新虫 (初入文坛)
SharedModel;
Parameters X0, Xm, um,a,b,P0,Yx,Yp,m,S0;
Variables t, X, P,S;
Function
           X= X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))
           P=P0+a*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)+(b*Xm/um)*ln((Xm-X0+X0*exp(um*t))/Xm)
           S=S0-(1/Yx+a/Yp)*(X0*Xm*exp(um*t)/(Xm-X0+X0*exp(um*t))-X0)-Xm/um*(b/Yp+m)*ln((Xm-X0+X0*exp(um*t))/Xm);
Data;
t    X          P        S
0   0.05        0        79
4   0.16        0        78
8   0.6515      0        76
12  1.5426      1.5      74
16  2.4047      4.9      69
20  3.1066      8.7      60
24  3.6693      12.4     50
28  3.8165      16.1     43
32  3.8874      17.2     40
36  3.9618      17.9     38
三个函数表达式分开可以算出来，可是合在一起就说表达式有误，分开算出的数据是错的，应该怎么弄，求大神指点
赞一下 回复此楼
18楼2014-12-16 22:30:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

CODE:
Parameters k;
Variable a,y;
Function y=int(k*y,y=0,10)+a;
RowData;
0,1,2,3;
1,3,4,5;

用微分方程的形式拟合应该也可以实现。
1stopt 语句修改
附图1.png
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

chengxl
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
19楼2014-12-17 09:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

6.0运行结果（稳定唯一）：

均方差(RMSE):0.483374374706425
残差平方和(SSE):7.00952358368483
相关系数(R): 0.981724799896839
相关系数之平方(R^2): 0.963783582732488
确定系数(DC): 0.937297250322943
F统计(F-Statistic): 0

参数                  最佳估算
--------------------        -------------
x0        0.0334076309586811
xm        4.64024768139639
um        0.241057508591346
a        4.1046342230345
b        -0.00415265747668047
p0        -0.287687677048546
yx        0.127392818019361
yp        5.42816890517109
m        0.0260506659407243
s0        78.5148748503121

====== 结果输出 ======

文件: 数据文件-1
No        目标 x        计算 x        目标 p        计算 p        目标 s        计算 s
1        0.05        0.0334076309586811        0        -0.287687677048546        79        78.5148748503121
2        0.16        0.0866085836080747        0        -0.0702457278129587        78        78.0513786216106
3        0.6515        0.220476302701316        0        0.476847080403129        76        76.8848015882738
4        1.5426        0.536863715228011        1.5        1.76956429951415        74        74.1258474300948
5        2.4047        1.18549279397836        4.9        4.41819545471953        69        68.4600563772828
6        3.1066        2.19800926890448        8.7        8.54648023793028        60        59.5776135163638
7        3.6693        3.25942090148948        12.4        12.8576041989121        50        50.165648205155
8        3.8165        3.99496251671989        16.1        15.8159219548877        43        43.4653609261927
9        3.8874        4.3710537088627        17.2        17.2897532001951        40        39.8032194191773
10        3.9618        4.53378908930007        17.9        17.8835669783185        38        37.9511990644952
赞一下 回复此楼
20楼2014-12-17 09:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chengxl 的主题更新
242/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 18:22 by skyjudy
[基金申请] 国自然面上复盘~欢迎讨论 +3 晴天加油 2026-08-26 4/200 2026-08-26 17:59 by Tide man
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
信息提示
关闭
请填处理意见
关闭