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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab非线性拟合模型参数，大神你在哪！
745/8返回列表首页上一页345678下一页
查看: 23955  |  回复: 73
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

月只蓝

主管区长 (职业作家)
引用回帖:
39楼: Originally posted by 梦落在草地上 at 2015-01-22 10:11:12
忽略38楼：
重启程序重新运行了下：
Parameters X1, X2, X3, Y;
Maximum = Y;
Function Y=355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2;

在Pa ...

Y本身是目标函数，限定Y在[1,100]，又求Y的最大值，那Y最大值不就是100么？
那么，不就等于是 求
Y=355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2=100
这个三元方程的解么？那么就是一个方程，三个未知数，肯定有无穷多个解。
作为验证我的说法，以下仍然按照最优化问题来求解，给出代码：
Parameters X1[1,7], X2[1,7], X3[1,7];
Maximum = Y;
Function Y=(355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2);
       (355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2)>=1;
       (355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2)<=100;

可见解显然是无穷多的，以下给出三组：

目标函数值(最大): 100
X1: 5.55127097786439
X2: 3.75469585593139
X3: 1.48118209473375

目标函数值(最大): 100
X1: 5.54702430245763
X2: 5.92326853055975
X3: 4.89868124102268

目标函数值(最大): 100
X1: 6.58142328449946
X2: 2.51431683205087
X3: 4.83067723055093
赞一下 回复此楼

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

梦落在草地上
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
41楼2015-01-22 10:20:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
41楼: Originally posted by 月只蓝 at 2015-01-22 10:20:36
Y本身是目标函数，限定Y在，又求Y的最大值，那Y最大值不就是100么？
那么，不就等于是 求
Y=355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2=100
这个 ...

我理解您意思了，三元二次方程的确有很多解，但如果就三个三元二次方程的最优呢？
Y1=(29.9+22.1*X1+11.6*X2-1.5*X1*X2+1.2*X1*X3-2.5*X1^2-0.6*X2^2-0.7*X3^2)
Y2=(-45.3+35.6*X1+24.6*X2+8.2*X3-1.1*X1*X2-4.0*X1^2-2.4*X2^2-1.3*X3^2)
Y3=(-15.3-4.0*X1+24.4*X2+22.4*X3+1.0*X1*X2+0.8*X1*X3-2.4*X2*X3-0.4*X1^2-2.4*X2^2-1.8*X3^2)
同时X1[1,7] X2[1,7] X3[1,7]

能得到一个综合方程，然后再求解最优吗？
赞一下 回复此楼
42楼2015-01-22 10:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
42楼: Originally posted by 梦落在草地上 at 2015-01-22 10:55:33
我理解您意思了，三元二次方程的确有很多解，但如果就三个三元二次方程的最优呢？
Y1=(29.9+22.1*X1+11.6*X2-1.5*X1*X2+1.2*X1*X3-2.5*X1^2-0.6*X2^2-0.7*X3^2)
Y2=(-45.3+35.6*X1+24.6*X2+8.2*X3-1.1*X1*X2-4.0 ...

这样都没有目标函数。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
43楼2015-01-22 11:15:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
送红花一朵
引用回帖:
41楼: Originally posted by 月只蓝 at 2015-01-22 10:20:36
Y本身是目标函数，限定Y在，又求Y的最大值，那Y最大值不就是100么？
那么，不就等于是 求
Y=355.45-38.35*X1-13.99*X2-85.78*X3+0.825*X1*X2+6.234*X1*X3+1.959*X2*X3+1.278*X1^2+0.272*X2^2+6.59*X3^2=100
这个 ...

您好，感激您一直给我回帖，也帮我增加了matlab语言对化学数据的非线性拟合的认识。然而，化学数据处理不是纯粹的数学计算，有很多条件限制，所以似乎目前为止，还没有真正解决我自己的问题。现想最后求助下您：
1.对于Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Variable X1,X2,X3,Y;
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;
Data;
    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
所拟合的方程前，可不可以再限制下区间X1[1,7]X2[1,7]X3[1,7]Y[0,100],使得最后直接最优化值是落在各个区间的某个极小值，且不是边界值。比如高中数学学过，极值不一定不等于最值。
赞一下 回复此楼
44楼2015-01-22 21:52:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
44楼: Originally posted by 梦落在草地上 at 2015-01-22 21:52:02
您好，感激您一直给我回帖，也帮我增加了matlab语言对化学数据的非线性拟合的认识。然而，化学数据处理不是纯粹的数学计算，有很多条件限制，所以似乎目前为止，还没有真正解决我自己的问题。现想最后求助下您：
...

主要是拟合方程时限制区间为Y[1,100]下的拟合曲线方程（X限制好像没有意义）
赞一下 回复此楼
45楼2015-01-22 21:55:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
44楼: Originally posted by 梦落在草地上 at 2015-01-22 21:52:02
您好，感激您一直给我回帖，也帮我增加了matlab语言对化学数据的非线性拟合的认识。然而，化学数据处理不是纯粹的数学计算，有很多条件限制，所以似乎目前为止，还没有真正解决我自己的问题。现想最后求助下您：
...

那我建议你先拟合出方程，然后对该方程求偏导数，令偏导数=0，求得极值点。
1stopt求得的都是最值（限制参数范围，求得的也是该范围内的最值）。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
46楼2015-01-22 22:14:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
46楼: Originally posted by 月只蓝 at 2015-01-22 22:14:31
那我建议你先拟合出方程，然后对该方程求偏导数，令偏导数=0，求得极值点。
1stopt求得的都是最值（限制参数范围，求得的也是该范围内的最值）。...

请问不以求极值为目的的情况下，拟合方程时能不能限制Y区间在0-100之间呢？
赞一下 回复此楼
47楼2015-01-23 08:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
47楼: Originally posted by 梦落在草地上 at 2015-01-23 08:46:39
请问不以求极值为目的的情况下，拟合方程时能不能限制Y区间在0-100之间呢？...

运行代码：
CODE:
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Variable X1,X2,X3,Y;
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;
Data;
    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

结果：
相关系数之平方(R^2): 0.910015554060544
决定系数(DC): 0.910015554060543
卡方系数(Chi-Square): 6.30625010587871
F统计(F-Statistic): 121.356380369046

参数        最佳估算
----------        -------------
b0                 274.297442101538
b1                 -44.818446582328
b2                 21.7155760160432
b3                 -78.7219393455877
b4                 0.820652156345272
b5                 6.79527940302928
b6                 0.732950301205965
b7                 1.89990679569777
b8                 -4.33596272234971
b9                 6.34139744840854

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

No        实测值y        计算值y
1        76.84        78.7268296
2        77.85        79.9128753
3        62.19        51.1783409
4        70.77        71.2205570
5        73.3        79.1110228
6        66.23        69.6920682
7        10.66        20.4925341
8        79.36        69.5274660
9        79.87        76.4079319
10        77.34        71.5289771
11        21.77        21.3194432
12        74.82        85.8316589
13        74.31        72.2471250
14        58.65        56.7631700

可以看到上面的三列数据，实测值y和计算值y 都在[0,100]的范围。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
48楼2015-01-23 09:16:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梦落在草地上

新虫 (小有名气)
引用回帖:
48楼: Originally posted by 月只蓝 at 2015-01-23 09:16:14
运行代码：
Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9;
Variable X1,X2,X3,Y;
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;
Data;
    1          2        4     ...

可是拟合的方程在自变量区间内，Y值却不在[0,100]之间了。因此最大值远超过100.而拟合方程后再限定区间内最大值100，毫无意义。所以在拟合出方程前就需要限定Y取值范围0-100再求最优值才行。
赞一下 回复此楼
49楼2015-01-23 09:30:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
49楼: Originally posted by 梦落在草地上 at 2015-01-23 09:30:23
可是拟合的方程在自变量区间内，Y值却不在之间了。因此最大值远超过100.而拟合方程后再限定区间内最大值100，毫无意义。所以在拟合出方程前就需要限定Y取值范围0-100再求最优值才行。...

你是说拟合出来的方程，在X1∈[1,7]，X2∈[1,7]，X3∈[1,7]的范围内，最大值大于100对吧？
以我的思路，那你先把这个拟合问题的目标函数以最小二乘的形式写出来吧，然后把 Y ≤100 加在约束条件里面，可以实现 Y值不超过100的拟合。
以我的了解，1stopt软件不能直接做约束Y值的拟合。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
50楼2015-01-23 10:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dshuanger 的主题更新
745/8返回列表首页上一页345678下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 293求调剂 +12 我爱高数高数爱� 2026-04-12 13/650 2026-04-12 21:32 by wangluhui
[考研] 计算机22408 281分，求调剂 +7 17715607211 2026-04-06 7/350 2026-04-12 00:43 by xuxiang
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 求调剂288 +6 ioodiiij 2026-04-10 8/400 2026-04-10 21:07 by zhouxiaoyu
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 考研调剂 +13 冰冰，，， 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 368化学求调剂 +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 085100建筑学 寻求跨专业调剂 一志愿南大294分 校级省级国家级奖项若干 踏实肯干 +3 1021075758 2026-04-06 4/200 2026-04-07 09:23 by 蓝云思雨
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
信息提示
关闭
请填处理意见
关闭