24小时热门版块排行榜    

查看: 692  |  回复: 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

新虫 (著名写手)

求大神指点,我只是增加定量

发自小木虫Android客户端
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 的主题更新
信息提示
请填处理意见