24小时热门版块排行榜    

CyRhmU.jpeg
查看: 23379  |  回复: 73

梦落在草地上

新虫 (小有名气)

引用回帖:
50楼: Originally posted by 月只蓝 at 2015-01-23 10:33:44
你是说拟合出来的方程,在X1∈,X2∈,X3∈的范围内,最大值大于100对吧?
以我的思路,那你先把这个拟合问题的目标函数以最小二乘的形式写出来吧,然后把 Y ≤100 加在约束条件里面,可以实现 Y值不超过100的拟合 ...

对的。
最小二乘的形式∑(Y --Y平)^2吧?
matlab可以实现吗?
51楼2015-01-23 11:02:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
51楼: Originally posted by 梦落在草地上 at 2015-01-23 11:02:10
对的。
最小二乘的形式∑(Y --Y平)^2吧?
matlab可以实现吗?...

最小二乘的形式∑ (Y公式 --Y实验值)^2
这样就是最优化问题了,MATLAB和1stopt都可以。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
52楼2015-01-23 11:04:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)

引用回帖:
52楼: Originally posted by 月只蓝 at 2015-01-23 11:04:09
最小二乘的形式∑ (Y公式 --Y实验值)^2
这样就是最优化问题了,MATLAB和1stopt都可以。...

麻烦您示范下怎么写,可以吗?以这组数据为例:1          2        4       76.84
    1          4        7        77.85
    2          6        3        62.19
    2          1        7        70.77
    3          3        3        73.30
    3          5        6        66.23
    4          7        2        10.66
    4          1        6        79.36
    5          3        2        79.87
    5          5        5        77.34
    6          7         1        21.77
    6          2        5        74.82
    7          4        1        74.31
    7          6        4        58.65
53楼2015-01-23 11:21:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
53楼: Originally posted by 梦落在草地上 at 2015-01-23 11:21:43
麻烦您示范下怎么写,可以吗?以这组数据为例:1          2        4       76.84
    1          4        7        77.85
    2          6        3        62.19
    2          1        7        70.77
    3          3        3        73.30
    3          5        6        66.23
    4          7        2        10.66
    4          1        6        79.36
...

第一项:
(b0+b1*1+b2*2+b3*4+b4*1*2+b5*1*4+b6*2*4+b7*1^2+b8*2^2+b9*4^2 - 76.84)^2 + ......
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
54楼2015-01-23 11:31:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)

引用回帖:
54楼: Originally posted by 月只蓝 at 2015-01-23 11:31:14
第一项:
(b0+b1*1+b2*2+b3*4+b4*1*2+b5*1*4+b6*2*4+b7*1^2+b8*2^2+b9*4^2 - 76.84)^2 + .........

Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9,X1[1,7],X2[1,7],X3[1,7];
Function Y=b0+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+b8*X2^2+b9*X3^2;
Minimum Z
Z=(b0+b1*1+b2*2+b3*4+b4*1*2+b5*1*4+b6*2*4+b7*1^2+b8*2^2+b9*4^2-76.84)^2
+(b0+b1*1+b2*4+b3*7+b4*1*4+b5*1*7+b6*4*7+b7*1^4+b8*4^2+b9*7^2-77.85)^2
+(b0+b1*2+b2*6+b3*3+b4*2*6+b5*2*3+b6*6*3+b7*2^4+b8*6^2+b9*3^2-62.19)^2
+(b0+b1*2+b2*1+b3*7+b4*2*1+b5*2*7+b6*1*7+b7*2^4+b8*1^2+b9*7^2-70.77)^2+(b0+b1*3+b2*3+b3*3+b4*3*3+b5*3*3+b6*3*3+b7*3^3+b8*3^2+b9*3^2-73.30)^2
+(b0+b1*3+b2*5+b3*6+b4*3*5+b5*3*6+b6*5*6+b7*3^4+b8*5^2+b9*6^2-66.23)^2
+(b0+b1*4+b2*7+b3*2+b4*4*7+b5*4*2+b6*7*2+b7*4^4+b8*7^2+b9*2^2-10.66)^2
+(b0+b1*4+b2*1+b3*6+b4*4*1+b5*4*6+b6*1*6+b7*4^4+b8*1^2+b9*6^2-79.36)^2
+(b0+b1*5+b2*3+b3*2+b4*5*3+b5*5*2+b6*3*2+b7*5^4+b8*3^2+b9*2^2-79.87)^2
+(b0+b1*5+b2*5+b3*5+b4*5*5+b5*5*5+b6*5*5+b7*5^4+b8*5^2+b9*5^2-77.34)^2
+(b0+b1*6+b2*7+b3*1+b4*6*7+b5*6*1+b6*7*1+b7*6^4+b8*7^2+b9*1^2-21.77)^2
+(b0+b1*6+b2*2+b3*5+b4*6*2+b5*6*5+b6*2*5+b7*6^4+b8*2^2+b9*5^2-74.82)^2
+(b0+b1*7+b2*4+b3*2+b4*7*4+b5*7*1+b6*4*1+b7*7^4+b8*4^2+b9*1^2-74.31)^2
+(b0+b1*7+b2*6+b3*4+b4*7*6+b5*7*4+b6*6*4+b7*7^4+b8*6^2+b9*4^2-58.65)^2
Data;
...
然后就出现定义Z有误了...
能不能麻烦你给出可执行的命令啊?谢谢
matlab非线性拟合模型参数,大神你在哪!
16.png

55楼2015-01-23 15:31:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
55楼: Originally posted by 梦落在草地上 at 2015-01-23 15:31:21
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9,X1,X2,X3;
Function Y=b0+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+b8*X2^2+b9*X3^2;
Minimum Z
Z=(b0+b1*1+b2*2+b3*4+b4*1*2+b5*1*4+b6*2*4+b7*1^2+b ...

现在已经不做拟合了,不必给出关键词Data和具体数据。按以下即可构造目标函数了。1stopt代码的写法和格式,请详见 1stopt的帮助手册,按手册来就不会错的。

Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Minimum= Z;
Function Z=(b0+b1*1+b2*2+b3*4+b4*1*2+b5*1*4+b6*2*4+b7*1^2+b8*2^2+b9*4^2-76.84)^2
+(b0+b1*1+b2*4+b3*7+b4*1*4+b5*1*7+b6*4*7+b7*1^4+b8*4^2+b9*7^2-77.85)^2
+(b0+b1*2+b2*6+b3*3+b4*2*6+b5*2*3+b6*6*3+b7*2^4+b8*6^2+b9*3^2-62.19)^2
+(b0+b1*2+b2*1+b3*7+b4*2*1+b5*2*7+b6*1*7+b7*2^4+b8*1^2+b9*7^2-70.77)^2
+(b0+b1*3+b2*3+b3*3+b4*3*3+b5*3*3+b6*3*3+b7*3^3+b8*3^2+b9*3^2-73.30)^2
+(b0+b1*3+b2*5+b3*6+b4*3*5+b5*3*6+b6*5*6+b7*3^4+b8*5^2+b9*6^2-66.23)^2
+(b0+b1*4+b2*7+b3*2+b4*4*7+b5*4*2+b6*7*2+b7*4^4+b8*7^2+b9*2^2-10.66)^2
+(b0+b1*4+b2*1+b3*6+b4*4*1+b5*4*6+b6*1*6+b7*4^4+b8*1^2+b9*6^2-79.36)^2
+(b0+b1*5+b2*3+b3*2+b4*5*3+b5*5*2+b6*3*2+b7*5^4+b8*3^2+b9*2^2-79.87)^2
+(b0+b1*5+b2*5+b3*5+b4*5*5+b5*5*5+b6*5*5+b7*5^4+b8*5^2+b9*5^2-77.34)^2
+(b0+b1*6+b2*7+b3*1+b4*6*7+b5*6*1+b6*7*1+b7*6^4+b8*7^2+b9*1^2-21.77)^2
+(b0+b1*6+b2*2+b3*5+b4*6*2+b5*6*5+b6*2*5+b7*6^4+b8*2^2+b9*5^2-74.82)^2
+(b0+b1*7+b2*4+b3*2+b4*7*4+b5*7*1+b6*4*1+b7*7^4+b8*4^2+b9*1^2-74.31)^2
+(b0+b1*7+b2*6+b3*4+b4*7*6+b5*7*4+b6*6*4+b7*7^4+b8*6^2+b9*4^2-58.65)^2;
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
56楼2015-01-23 15:46:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)

引用回帖:
56楼: Originally posted by 月只蓝 at 2015-01-23 15:46:04
现在已经不做拟合了,不必给出关键词Data和具体数据。按以下即可构造目标函数了。1stopt代码的写法和格式,请详见 1stopt的帮助手册,按手册来就不会错的。

Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Minimu ...

这样做的目的,得到了目标函数值(最小): 535.046842913012和b0-b9的值。然后再带入公式,用minimum Y求最优就可以了吗?这种方法和之前的有什么意义上的不同?才导致了结果变了啊
57楼2015-01-23 16:05:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)

引用回帖:
56楼: Originally posted by 月只蓝 at 2015-01-23 15:46:04
现在已经不做拟合了,不必给出关键词Data和具体数据。按以下即可构造目标函数了。1stopt代码的写法和格式,请详见 1stopt的帮助手册,按手册来就不会错的。

Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Minimu ...

这样构造的目标函数,对于其R、F-statistic、P-Value就没有
58楼2015-01-23 16:36:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

引用回帖:
57楼: Originally posted by 梦落在草地上 at 2015-01-23 16:05:03
这样做的目的,得到了目标函数值(最小): 535.046842913012和b0-b9的值。然后再带入公式,用minimum Y求最优就可以了吗?这种方法和之前的有什么意义上的不同?才导致了结果变了啊...

上述代码只是把目标函数写出来,约束条件没有写进去,把约束条件加进去,得到的 b0-b9才能满足你在49楼说的条件。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
59楼2015-01-23 16:41:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)

引用回帖:
59楼: Originally posted by 月只蓝 at 2015-01-23 16:41:38
上述代码只是把目标函数写出来,约束条件没有写进去,把约束条件加进去,得到的 b0-b9才能满足你在49楼说的条件。...

这里得到目标函数,同时也得到了b0-b9了,难道不能带入我们的拟合方程后再在X1[1,7]X2[1,7]X3[1,7]Y[0,100]求最优?还有拟合方程的R2、F值该怎么求?
60楼2015-01-24 17:37:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dshuanger 的主题更新
信息提示
请填处理意见