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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab非线性拟合模型参数，大神你在哪！
746/8返回列表首页上一页345678下一页
查看: 23955  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 268分085602化学工程调剂 +31 月照花林。 2026-04-09 31/1550 2026-04-12 23:37 by 学zh
[考研] 2本，初试303，0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 0856专硕求调剂 希望是a区院校 +22 好好休息好不好 2026-04-09 25/1250 2026-04-12 17:09 by ajpv风雷
[考研] 286求调剂 +26 Faune 2026-04-06 26/1300 2026-04-12 10:24 by lhj2009
[考研] 化工调剂求导师收留！一志愿失利，踏实肯干，有植物提取科研经历 +20 yzyzx 2026-04-09 21/1050 2026-04-12 00:12 by 小小小小啦啦啦
[考研] 271求调剂 +20 2261744733 2026-04-11 22/1100 2026-04-11 23:14 by labixiaoqiao
[考研] 26自然地理学303分求调剂 +6 一战成硕啊啊啊� 2026-04-06 11/550 2026-04-11 21:27 by 裴宏伟
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[基金申请] 有爆料，一个青年教师卖房得400万，然后换了一个四青帽子 +9 babu2015 2026-04-08 9/450 2026-04-10 11:43 by 苏东坡二世
[考研] 青岛科技大学材料学院，环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +3 慕绝cc 2026-04-09 3/150 2026-04-09 09:57 by liuhuiying09
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 材料工程专业日语生求调剂 +9 111623 2026-04-07 9/450 2026-04-07 23:31 by 一只好果子?
[考研] 生物学363调剂求助 +7 fanzhang6666 2026-04-06 9/450 2026-04-07 17:37 by lijunpoly
[考研] 材料调剂 +17 小刘同学吖吖 2026-04-06 18/900 2026-04-07 11:41 by 诗与自由
[考研] 287分求调剂 有专利国奖一志愿哈工大085406 +6 白易辰 2026-04-06 7/350 2026-04-06 22:46 by 875465
[考研] 求助 +3 卡卡东88 2026-04-06 4/200 2026-04-06 15:28 by going home
信息提示
关闭
请填处理意见
关闭