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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助：以下非线性拟合方程可以用1stopt求解么？用什么方法
北京石油化工学院2026年研究生招生接收调剂公告
171/2返回列表12下一页
查看: 1670  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hanmei98232

木虫 (小有名气)

[求助] 求助：以下非线性拟合方程可以用1stopt求解么？用什么方法

以下非线性拟合方程可以用1stopt求解么？用什么方法
0.07272703=-(m*0.17809/(n*0.17809^0.5+1))*8.314*298.15*(ln(0.17809)-0.784*(1.7321*0.17809^0.5/(1+2.0785*0.17809^0.5)+1.6667*ln(1+2.0785*0.17809^0.5))+2.6667*0.17809*(c+d*(1+(-1-3.464*0.17809^0.5+6*0.17809)*e0.17809p(-3.464*0.17809^0.5))/(12*0.17809))+2.8289*0.17809^2*e)+0.07214;
          0.07319177=-(m*0.36061/(n*0.36061^0.5+1))*8.314*298.15*(ln(0.36061)-0.784*(1.7321*0.36061^0.5/(1+2.0785*0.36061^0.5)+1.6667*ln(1+2.0785*0.36061^0.5))+2.6667*0.36061*(c+d*(1+(-1-3.464*0.36061^0.5+6*0.36061)*e0.36061p(-3.464*0.36061^0.5))/(12*0.36061))+2.8289*0.36061^2*e)+0.07214;
          0.07391741=-(m*0.53672/(n*0.53672^0.5+1))*8.314*298.15*(ln(0.53672)-0.784*(1.7321*0.53672^0.5/(1+2.0785*0.53672^0.5)+1.6667*ln(1+2.0785*0.53672^0.5))+2.6667*0.53672*(c+d*(1+(-1-3.464*0.53672^0.5+6*0.53672)*e0.53672p(-3.464*0.53672^0.5))/(12*0.53672))+2.8289*0.53672^2*e)+0.07214;
          0.07451667=-(m*0.72049/(n*0.72049^0.5+1))*8.314*298.15*(ln(0.72049)-0.784*(1.7321*0.72049^0.5/(1+2.0785*0.72049^0.5)+1.6667*ln(1+2.0785*0.72049^0.5))+2.6667*0.72049*(c+d*(1+(-1-3.464*0.72049^0.5+6*0.72049)*e0.72049p(-3.464*0.72049^0.5))/(12*0.72049))+2.8289*0.72049^2*e)+0.07214;
          0.07521378=-(m*0.92046/(n*0.92046^0.5+1))*8.314*298.15*(ln(0.92046)-0.784*(1.7321*0.92046^0.5/(1+2.0785*0.92046^0.5)+1.6667*ln(1+2.0785*0.92046^0.5))+2.6667*0.92046*(c+d*(1+(-1-3.464*0.92046^0.5+6*0.92046)*e0.92046p(-3.464*0.92046^0.5))/(12*0.92046))+2.8289*0.92046^2*e)+0.07214;
          0.07661342=(m*1.3003/(n*1.3003^0.5+1))*8.314*298.15*(ln(1.3003)-0.784*(1.7321*1.3003^0.5/(1+2.0785*1.3003^0.5)+1.6667*ln(1+2.0785*1.3003^0.5))+2.6667*1.3003*(c+d*(1+(-1-3.464*1.3003^0.5+6*1.3003)*e1.3003p(-3.464*1.3003^0.5))/(12*1.3003))+2.8289*1.3003^2*e)+0.07214;
          0.07871152=(m*1.67471/(n*1.67471^0.5+1))*8.314*298.15*(ln(1.67471)-0.784*(1.7321*1.67471^0.5/(1+2.0785*1.67471^0.5)+1.6667*ln(1+2.0785*1.67471^0.5))+2.6667*1.67471*(c+d*(1+(-1-3.464*1.67471^0.5+6*1.67471)*e1.67471p(-3.464*1.67471^0.5))/(12*1.67471))+2.8289*1.67471^2*e)+0.07214;
          0.07934748=(m*1.86291/(n*1.86291^0.5+1))*8.314*298.15*(ln(1.86291)-0.784*(1.7321*1.86291^0.5/(1+2.0785*1.86291^0.5)+1.6667*ln(1+2.0785*1.86291^0.5))+2.6667*1.86291*(c+d*(1+(-1-3.464*1.86291^0.5+6*1.86291)*e1.86291p(-3.464*1.86291^0.5))/(12*1.86291))+2.8289*1.86291^2*e)+0.07214;
          0.08119555=(m*2.42891/(n*2.42891^0.5+1))*8.314*298.15*(ln(2.42891)-0.784*(1.7321*2.42891^0.5/(1+2.0785*2.42891^0.5)+1.6667*ln(1+2.0785*2.42891^0.5))+2.6667*2.42891*(c+d*(1+(-1-3.464*2.42891^0.5+6*2.42891)*e2.42891p(-3.464*2.42891^0.5))/(12*2.42891))+2.8289*2.42891^2*e)+0.07214;
          0.08388612=(m*2.94245/(n*2.94245^0.5+1))*8.314*298.15*(ln(2.94245)-0.784*(1.7321*2.94245^0.5/(1+2.0785*2.94245^0.5)+1.6667*ln(1+2.0785*2.94245^0.5))+2.6667*2.94245*(c+d*(1+(-1-3.464*2.94245^0.5+6*2.94245)*e2.94245p(-3.464*2.94245^0.5))/(12*2.94245))+2.8289*2.94245^2*e)+0.07214;
          0.08544611=(m*3.49978/(n*3.49978^0.5+1))*8.314*298.15*(ln(3.49978)-0.784*(1.7321*3.49978^0.5/(1+2.0785*3.49978^0.5)+1.6667*ln(1+2.0785*3.49978^0.5))+2.6667*3.49978*(c+d*(1+(-1-3.464*3.49978^0.5+6*3.49978)*e3.49978p(-3.464*3.49978^0.5))/(12*3.49978))+2.8289*3.49978^2*e)+0.07214;
          0.08798993=(m*4.05405/(n*4.05405^0.5+1))*8.314*298.15*(ln(4.05405)-0.784*(1.7321*4.05405^0.5/(1+2.0785*4.05405^0.5)+1.6667*ln(1+2.0785*4.05405^0.5))+2.6667*4.05405*(c+d*(1+(-1-3.464*4.05405^0.5+6*4.05405)*e4.05405p(-3.464*4.05405^0.5))/(12*4.05405))+2.8289*4.05405^2*e)+0.07214;
即x<1
function y=-(m*X/(n*X^0.5+1))*8.314*298.15*(ln(x)-0.784*(1.7321*x^0.5/(1+2.0785*x^0.5)+1.6667*ln(1+2.0785*x^0.5))+2.6667*x*(c+d*(1+(-1-3.464*x^0.5+6*x)*exp(-3.464*x^0.5))/(12*x))+2.8289*x^2*e)+0.07214;
x>1
y=(m*X/(n*X^0.5+1))*8.314*298.15*(ln(x)-0.784*(1.7321*x^0.5/(1+2.0785*x^0.5)+1.6667*ln(1+2.0785*x^0.5))+2.6667*x*(c+d*(1+(-1-3.464*x^0.5+6*x)*exp(-3.464*x^0.5))/(12*x))+2.8289*x^2*e)+0.07214;
data;
0.17809         0.07272703
0.36061        0.07319177
0.53672        0.07391741
0.72049           0.07451667
0.92046        0.07521378
1.3003         0.07661342
  1.67471       0.07871152
1.86291        0.07934748
  2.42891       0.08119555
2.94245        0.08388612
3.49978        0.08544611
4.05405        0.08798993
如果1stopt解不了用matlab应该怎么解？
求助求助！！！！非常感谢
Parameters c,d,e;
Variable x1,x2,y1,y2 ;
function y1=(m*X1/(n*X1^0.5+1))*8.314*298.15*(ln(x1)-0.784*(1.7321*x1^0.5/(1+2.0785*x1^0.5)+1.6667*ln(1+2.0785*x1^0.5))+2.6667*x1*(c+d*(1+(-1-3.464*x1^0.5+6*x1)*exp(-3.464*x1^0.5))/(12*x1))+2.8289*x1^2*e)+0.07214;
         y2=-(m*X2/(n*X2^0.5+1))*8.314*298.15*(ln(x2)-0.784*(1.7321*x2^0.5/(1+2.0785*x2^0.5)+1.6667*ln(1+2.0785*x2^0.5))+2.6667*x2*(c+d*(1+(-1-3.464*x2^0.5+6*x2)*exp(-3.464*x2^0.5))/(12*x2))+2.8289*x2^2*e)+0.07214;
Data;
//x1 y1 x2 y2
0.17809         0.07272703                    1.3003         0.07661342
0.36061        0.07319177                    1.67471       0.07871152
0.53672        0.07391741                    1.86291        0.07934748
0.72049           0.07451667                    2.42891       0.08119555
0.92046        0.07521378                    2.94245        0.08388612
                                                              3.49978        0.08544611
                                                              4.05405        0.08798993
求教 这个样编写对么？为什么出来的图和y对不上呢？

[ Last edited by hanmei98232 on 2014-3-15 at 13:46 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-03-15 12:55:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
1stopt软件拟合、解方程的标准代码参见附件。
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 1stOpt使用手册.pdf
    • 2014-03-15 14:17:42, 1.67 M
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-03-15 14:18:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanmei98232

木虫 (小有名气)
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-03-15 14:18:20
1stopt软件拟合、解方程的标准代码参见附件。

使用手册我有 也看了 程序发上去了 但是拟合不出来啊 不知道应该怎么办
赞一下 回复此楼
3楼2014-03-15 14:43:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
hanmei98232: 金币+5, ★有帮助, 虽然没有解答 还是谢谢你 2014-03-19 09:49:10
引用回帖:
3楼: Originally posted by hanmei98232 at 2014-03-15 14:43:14
使用手册我有 也看了 程序发上去了 但是拟合不出来啊 不知道应该怎么办...

共享参数拟合么？你找本版的实习版主，请他帮你做
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2014-03-15 14:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
3楼: Originally posted by hanmei98232 at 2014-03-15 14:43:14
使用手册我有 也看了 程序发上去了 但是拟合不出来啊 不知道应该怎么办...

高版本的1stopt可以做共享参数拟合。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2014-03-15 14:47:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
0.07272703=-(m*0.17809/(n*0.17809^0.5+1))*8.314*298.15*(ln(0.17809)-0.784*(1.7321*0.17809^0.5/(1+2.0785*0.17809^0.5)+1.6667*ln(1+2.0785*0.17809^0.5))+2.6667*0.17809*(c+d*(1+(-1-3.464*0.17809^0.5+6*0.17809)*e0.17809p(-3.464*0.17809^0.5))/(12*0.17809))+2.8289*0.17809^2*e)+0.07214;
          0.07319177=-(m*0.36061/(n*0.36061^0.5+1))*8.314*298.15*(ln(0.36061)-0.784*(1.7321*0.36061^0.5/(1+2.0785*0.36061^0.5)+1.6667*ln(1+2.0785*0.36061^0.5))+2.6667*0.36061*(c+d*(1+(-1-3.464*0.36061^0.5+6*0.36061)*e0.36061p(-3.464*0.36061^0.5))/(12*0.36061))+2.8289*0.36061^2*e)+0.07214;
          0.07391741=-(m*0.53672/(n*0.53672^0.5+1))*8.314*298.15*(ln(0.53672)-0.784*(1.7321*0.53672^0.5/(1+2.0785*0.53672^0.5)+1.6667*ln(1+2.0785*0.53672^0.5))+2.6667*0.53672*(c+d*(1+(-1-3.464*0.53672^0.5+6*0.53672)*e0.53672p(-3.464*0.53672^0.5))/(12*0.53672))+2.8289*0.53672^2*e)+0.07214;
          0.07451667=-(m*0.72049/(n*0.72049^0.5+1))*8.314*298.15*(ln(0.72049)-0.784*(1.7321*0.72049^0.5/(1+2.0785*0.72049^0.5)+1.6667*ln(1+2.0785*0.72049^0.5))+2.6667*0.72049*(c+d*(1+(-1-3.464*0.72049^0.5+6*0.72049)*e0.72049p(-3.464*0.72049^0.5))/(12*0.72049))+2.8289*0.72049^2*e)+0.07214;
          0.07521378=-(m*0.92046/(n*0.92046^0.5+1))*8.314*298.15*(ln(0.92046)-0.784*(1.7321*0.92046^0.5/(1+2.0785*0.92046^0.5)+1.6667*ln(1+2.0785*0.92046^0.5))+2.6667*0.92046*(c+d*(1+(-1-3.464*0.92046^0.5+6*0.92046)*e0.92046p(-3.464*0.92046^0.5))/(12*0.92046))+2.8289*0.92046^2*e)+0.07214;
          0.07661342=(m*1.3003/(n*1.3003^0.5+1))*8.314*298.15*(ln(1.3003)-0.784*(1.7321*1.3003^0.5/(1+2.0785*1.3003^0.5)+1.6667*ln(1+2.0785*1.3003^0.5))+2.6667*1.3003*(c+d*(1+(-1-3.464*1.3003^0.5+6*1.3003)*e1.3003p(-3.464*1.3003^0.5))/(12*1.3003))+2.8289*1.3003^2*e)+0.07214;
          0.07871152=(m*1.67471/(n*1.67471^0.5+1))*8.314*298.15*(ln(1.67471)-0.784*(1.7321*1.67471^0.5/(1+2.0785*1.67471^0.5)+1.6667*ln(1+2.0785*1.67471^0.5))+2.6667*1.67471*(c+d*(1+(-1-3.464*1.67471^0.5+6*1.67471)*e1.67471p(-3.464*1.67471^0.5))/(12*1.67471))+2.8289*1.67471^2*e)+0.07214;
          0.07934748=(m*1.86291/(n*1.86291^0.5+1))*8.314*298.15*(ln(1.86291)-0.784*(1.7321*1.86291^0.5/(1+2.0785*1.86291^0.5)+1.6667*ln(1+2.0785*1.86291^0.5))+2.6667*1.86291*(c+d*(1+(-1-3.464*1.86291^0.5+6*1.86291)*e1.86291p(-3.464*1.86291^0.5))/(12*1.86291))+2.8289*1.86291^2*e)+0.07214;
          0.08119555=(m*2.42891/(n*2.42891^0.5+1))*8.314*298.15*(ln(2.42891)-0.784*(1.7321*2.42891^0.5/(1+2.0785*2.42891^0.5)+1.6667*ln(1+2.0785*2.42891^0.5))+2.6667*2.42891*(c+d*(1+(-1-3.464*2.42891^0.5+6*2.42891)*e2.42891p(-3.464*2.42891^0.5))/(12*2.42891))+2.8289*2.42891^2*e)+0.07214;
          0.08388612=(m*2.94245/(n*2.94245^0.5+1))*8.314*298.15*(ln(2.94245)-0.784*(1.7321*2.94245^0.5/(1+2.0785*2.94245^0.5)+1.6667*ln(1+2.0785*2.94245^0.5))+2.6667*2.94245*(c+d*(1+(-1-3.464*2.94245^0.5+6*2.94245)*e2.94245p(-3.464*2.94245^0.5))/(12*2.94245))+2.8289*2.94245^2*e)+0.07214;
          0.08544611=(m*3.49978/(n*3.49978^0.5+1))*8.314*298.15*(ln(3.49978)-0.784*(1.7321*3.49978^0.5/(1+2.0785*3.49978^0.5)+1.6667*ln(1+2.0785*3.49978^0.5))+2.6667*3.49978*(c+d*(1+(-1-3.464*3.49978^0.5+6*3.49978)*e3.49978p(-3.464*3.49978^0.5))/(12*3.49978))+2.8289*3.49978^2*e)+0.07214;
          0.08798993=(m*4.05405/(n*4.05405^0.5+1))*8.314*298.15*(ln(4.05405)-0.784*(1.7321*4.05405^0.5/(1+2.0785*4.05405^0.5)+1.6667*ln(1+2.0785*4.05405^0.5))+2.6667*4.05405*(c+d*(1+(-1-3.464*4.05405^0.5+6*4.05405)*e4.05405p(-3.464*4.05405^0.5))/(12*4.05405))+2.8289*4.05405^2*e)+0.07214;
以上放到function中，qiuyicicanshu
将data 分成两段，针对不同的拟合函数，分别做一次拟合。
将一个程序分成三个小程序处理，低版本1stop就可以胜任了。
赞一下 回复此楼
6楼2014-03-15 15:48:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
6#  qiuyicicanshu = 求解一次参数。
赞一下 回复此楼
7楼2014-03-15 15:49:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanmei98232

木虫 (小有名气)
引用回帖:
6楼: Originally posted by feixiaolin at 2014-03-15 15:48:21
0.07272703=-(m*0.17809/(n*0.17809^0.5+1))*8.314*298.15*(ln(0.17809)-0.784*(1.7321*0.17809^0.5/(1+2.0785*0.17809^0.5)+1.6667*ln(1+2.0785*0.17809^0.5))+2.6667*0.17809*(c+d*(1+(-1-3.464*0.17809^0.5+6*0. ...

我用的低版本1stopt 每次输入都提示错误 不知道为什么
赞一下 回复此楼
8楼2014-03-15 16:55:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanmei98232

木虫 (小有名气)
引用回帖:
5楼: Originally posted by 月只蓝 at 2014-03-15 14:47:20
高版本的1stopt可以做共享参数拟合。...

没有高版本1stopt想请教您一下如果用matlab需要什么程序去拟合？
赞一下 回复此楼
9楼2014-03-15 16:56:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanmei98232

木虫 (小有名气)
引用回帖:
4楼: Originally posted by 月只蓝 at 2014-03-15 14:45:56
共享参数拟合么？你找本版的实习版主，请他帮你做...

谢谢你
回复此楼
10楼2014-03-15 16:58:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hanmei98232 的主题更新
171/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕311分求调剂 +7 梁富贵险中求 2026-04-04 7/350 2026-04-05 11:11 by imissbao
[考研] 求调剂 +3 电气小神童 2026-04-04 3/150 2026-04-05 10:17 by barlinike
[考研] 285求调剂 +11 哦呦呼o 2026-04-04 11/550 2026-04-05 08:15 by 544594351
[考研] 349求调剂 +11 zwjjjjjj 2026-03-31 11/550 2026-04-04 19:52 by 蓝云思雨
[考研] 一志愿085404，总分291，四级已过，求调剂 +5 阿俊阿俊阿俊 2026-04-04 7/350 2026-04-04 13:23 by 莲菜就是藕吧
[考研] 本9一志愿2 0854低分专硕286求调剂 +9 芒种111 2026-04-04 9/450 2026-04-04 11:01 by tangruihua
[考研] 材料调剂 +11 吴棂颖！ 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 本科985，专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 085600专硕材料与化工348分求调剂 +10 上学啦！ 2026-04-01 11/550 2026-04-03 14:13 by 百灵童888
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[基金申请] 请问共同通讯和共同一作的认可度问题 10+4 psa1234 2026-04-01 10/500 2026-04-03 11:08 by Kittylucky
[考研] 考研调剂 +3 李木子0120 2026-04-02 5/250 2026-04-02 21:45 by dongzh2009
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 一志愿安徽大学计算机科学与技术学硕，331分求调剂 +5 蒋昌鹏qtj 2026-04-01 5/250 2026-04-02 08:10 by fxue1114
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
信息提示
关闭
请填处理意见
关闭