24小时热门版块排行榜    

查看: 915  |  回复: 9

何fox

新虫 (著名写手)

[求助] 1stopt代码问题

为什么拟合会出现k0,k1,k2.......
Title "Type your title here";
Parameters x,y,z;
Consistant x0= 0.218737114,y0= 0.710746186,z0=0.029325119 ,k1=1612,k2=300.976,k3=89.59,k4=400.54,k5=6.421328887,k6=234.8123341,k7=239.4058681,k8=218.096192,k9=210.1498107,k10=216.7968469;
FullLoopModel;
Function
x1=k1*y*x;
x2=k1*x*k2*y^2;
x3=k1*x*k2^2*y^3;
x4=k1*x*k2^3*y^4;
x5=k1*x*k2^4*y^5;
x6=k1*x*k2^5*y^6;
x7=k1*x*k2^7*y^7;
x8=k1*x*k2^8*y^8;
x9=k1*x*k2^9*y^9;
x10=k1*x*k2^10*y^10;

z1=k3*z*y;
z2=k3*z*k4*y^2;
z3=k3*z*k4*k5*y^3;
z4=k3*z*k4*k5^2*y^4;
z5=k3*z*k4*k5^3*y^5;
z6=k3*z*k4*k5^4*y^6;
z7=k3*z*k4*k5^5*y^7;
z8=k3*z*k4*k5^6*y^8;
z9=k3*z*k4*k5^7*y^9;
z10=k4*z*k4*k5^8*y^10;

a1=k5*x1*x/z;
a2=k6*a1*y;
a3=k7*a2*y;
a4=k8*a3*y;
a5=k9*a4*y;
a6=k10*a5*y;

x0=x+x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+2*(a1+a2+a3+a4+a5+a6);
z0=z+z1+z2+z3+z4+z5+z6+z7+z8+z9+z10+a1;
y0=y+x1+2*x2+3*x3+4*x4+5*x5+6*x6+7*x7+8*x8+9*x9+10*x10+z1+2*z2+3*z3+4*z4+5*z5+6*z6+7*z7+8*z8+9*z9+10*z10+a1+2*a2+3*a3+4*a4+5*a5+6*a6;
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何fox

新虫 (著名写手)

2楼2016-06-16 12:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

Consistant 应该是Constant

» 本帖已获得的红花(最新10朵)

一团和气
3楼2016-06-16 13:03:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何fox

新虫 (著名写手)

4楼2016-06-16 13:04:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何fox

新虫 (著名写手)

送红花一朵
引用回帖:
3楼: Originally posted by lipenggg at 2016-06-16 13:03:18
Consistant 应该是Constant

可以跑出来麽?请问

发自小木虫Android客户端
5楼2016-06-16 13:04:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

1.5破解版,没有FullLoopModel命令
去掉FullLoopModel命令,可以出来结果
目标函数值: 0.00785754731254706
x: 0.000580335405842577
y: 0.00207990585706349
z: 2.4593461250917E-5
x1: 0.00023321284886932
x2: 0.00121957126836834
x3: 0.000757766950931213
x4: 0.000502273137502681
x5: 0.000296444398768115
x6: 0.000186185503599189
x7: 0.0352459162929691
x8: 0.0220573557922947
x9: 0.0138136116582057
x10: 0.00863076558610986
z1: 1.1545585510466E-6
z2: -7.65791554488516E-6
z3: -2.53234979813311E-5
z4: -1.40485723010301E-5
z5: -2.90752697472699E-5
z6: -1.34208609761919E-5
z7: -0.000153422533363949
z8: -9.08045576677666E-5
z9: -5.5591901216861E-5
z10: -0.00565533245570862
a1: 0.0353461470995159
a2: 0.0172637642073741
a3: 0.00859910527174406
a4: 0.00390610655639767
a5: 0.00171468647787587
a6: 0.000776859739219732

» 本帖已获得的红花(最新10朵)

一团和气
6楼2016-06-16 13:07:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何fox

新虫 (著名写手)

送红花一朵
引用回帖:
6楼: Originally posted by lipenggg at 2016-06-16 13:07:43
1.5破解版,没有FullLoopModel命令
去掉FullLoopModel命令,可以出来结果
目标函数值: 0.00785754731254706
x: 0.000580335405842577
y: 0.00207990585706349
z: 2.4593461250917E-5
x1: 0.00023321284886932 ...

请问你用的什么方法?类牛顿法?我这个答案主要没负数。。。

发自小木虫Android客户端
7楼2016-06-16 13:09:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

标准简面体爬山法 + 通用全局优化法
一团和气
8楼2016-06-16 13:10:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lipenggg

专家顾问 (正式写手)

【答案】应助回帖

不要负数可以加取值范围

» 本帖已获得的红花(最新10朵)

一团和气
9楼2016-06-16 13:15:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何fox

新虫 (著名写手)

送红花一朵
引用回帖:
9楼: Originally posted by lipenggg at 2016-06-16 13:15:23
不要负数可以加取值范围

谢了,非常感谢

发自小木虫Android客户端
10楼2016-06-16 13:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 何fox 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 5/250 2026-08-13 16:10 by majunge000
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 6/300 2026-08-13 15:34 by 高飞的梦想
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 26/1300 2026-08-13 12:01 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然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
[基金申请] 国自然结果 +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有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见