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

24小时热门版块排行榜    

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

梦落在草地上

新虫 (小有名气)
引用回帖:
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 的主题更新
746/8返回列表首页上一页345678下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 BMTdnp4OJO9f 2026-11-06 4/200 2026-11-07 18:52 by Xc2FxHsHkBNc
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 BMTdnp4OJO9f 2026-11-06 5/250 2026-11-07 17:56 by Xc2FxHsHkBNc
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 BMTdnp4OJO9f 2026-11-06 5/250 2026-11-07 17:45 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 17:01 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[基金申请] 各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
关闭
请填处理意见
关闭