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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab非线性拟合模型参数，大神你在哪！
746/8返回列表首页上一页345678下一页
查看: 24575  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭