24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hanmei98232

木虫 (小有名气)

应助: 1 (幼儿园)
金币: 2315.9
红花: 2
帖子: 285
在线: 95.7小时
虫号: 1357946
注册: 2011-07-31
性别: MM
专业: 催化化学

[求助] 求助：以下非线性拟合方程可以用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 ]

回复此楼