1stopt拟合参数比较多,希望有5.0以上版本的大神帮助运算一下,程序见附件,非常感谢! 返回小木虫查看更多
自己发完贴都不带再看一眼的啊?附件在哪儿?
感觉楼主的公式有问题:exp(-(48900/8.314*x2)*(x2/436.45-1)+....,其中“-(48900/8.314*x2)*(x2/436.45-1)”项,对应x2系列最小(278.15)和最大值(323.15)分别为:593368.271762467,493399.16621225,这两个数值还能进行指数计算吗?试下exp(493399.16621225),早就计算溢出了。
自己发完贴都不带再看一眼的啊?附件在哪儿?
非常抱歉,当时发完贴以为附件上传了,没有仔细看。我这是个2个因变量一个自变量的非线性拟合,函数是个隐函数。关系式是y=exp(-(48900/8.314*x2)*(x2/436.45-1)+((a4*(1-y)*x1*exp(-a7*a4)+a5*(1-y)*(1-x1)*exp(-a8*a5))*((1-y)*x1*exp(-a7*a4)+(1-y)*(1-x1)*exp(-a8*a5)))/(y+(1-y)*x1*exp(-a7*a4)+(1-y)*(1-x1)*exp(-a8*a5))^2+((1-y)*x1*exp(-a7*a1)*(a1*(1-y)*x1+a1*(1-y)*(1-x1)*exp(-a9*a6)-a6*(1-y)*(1-x1)*exp(-a9*a6)))/(y*exp(-a7*a1)+(1-y)*x1+(1-y)*(1-x1)*exp(-a9*a6))^2+((1-y)*(1-x1)*exp(-a8*a2)*(a2*(1-y)*(1-x1)+a2*(1-y)*x1*exp(-a9*a3)-a3*(1-y)*x1*exp(-a9*a3)))/(y*exp(-a8*a2)+(1-y)*x1*exp(-a9*a3)+(1-y)*(1-x1))^2) x1,x2是自变量,y是因变量。a1-a9是参数,其中a7,a8,a9这三个参数取值是0-1,在论坛上看matlab对于隐函数的拟合好像不是太友好,然后我自己也用低版本的1stopt试了试,并不能出结果。不清楚高版本的1stopt能不能运行出来,求大佬帮忙指点一下
,
感觉楼主的公式有问题:exp(-(48900/8.314*x2)*(x2/436.45-1)+....,其中“-(48900/8.314*x2)*(x2/436.45-1)”项,对应x2系列最小(278.15)和最大值(323.15)分别为:593368.271762467,493399.16621225,这两个数值还能进行指数计算吗?试下exp(493399.16621225),早就计算溢出了。