| 查看: 2072 | 回复: 11 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
好心虫友，帮忙看看这个matlab程序那里出错了 已有1人参与
|
|||
|
用matlab拟合数据，进行NRTL模型拟合：我编写了如下的程序，一直出错，麻烦有懂的虫友帮忙看看，怎么改错？？ （一）定义.m文件。 function g4=f(beta0,x1,x2,T,R) g1=beta0(1); g2=beta0(2); g3=beta0(3); g4=beta0(4); g4=x2.*x2.*((((g(1)./(R.*T)).*(exp(-g(3).*g(1)./(R.*T))).^2)./((x1+(exp(-g(3).*g(1)./(R.*T))).*x2).^2))+(((g(2)./(R.*T)).*(exp(-g(3).*g(2)./(R.*T))).^2)./((x2+(exp(-g(3).*g(2)./(R.*T))).*x1).^2))) （二）编写程序。 x1=[0.01525 0.02059 0.02685 0.03417 0.04221 0.04958]; x2=[0.98475 0.97941 0.97315 0.96583 0.95779 0.95042]; T=[298.15 303.15 308.15 313.15 318.15 323.15]; R=[8.314 8.314 8.314 8.314 8.314 8.314]; beta0=[0,0,0,0]; [g,resnorm]=lsqcurvefit(@f,beta0,x1,x2,T,R) 运行后出错： Warning: Length of lower bounds is > length(x); ignoring extra bounds. > In checkbounds at 27 In optim\private\lsqncommon at 39 In lsqcurvefit at 149 In Untitled2 at 6 Warning: Length of upper bounds is > length(x); ignoring extra bounds. > In checkbounds at 37 In optim\private\lsqncommon at 39 In lsqcurvefit at 149 In Untitled2 at 6 Exiting due to infeasibility: 4 lower bounds exceed the corresponding upper bounds. g = 0 0 0 0 resnorm = [] |
回复此楼» 猜你喜欢
今年春晚有几个节目很不错，点赞！
已经有4人回复
-
球磨粉体时遇到了大的问题，请指教！
已经有15人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有5人回复
-
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复
|
那我初始值就定义4个可以吗？我按照你说的改了程序，把初始值改为【1,1,0.3,1】不再为0，结果出来后就是Warning: Length of lower bounds is > length(x); ignoring extra bounds. > In checkbounds at 27 In optim\private\lsqncommon at 39 In lsqcurvefit at 149 In Untitled2 at 6 Warning: Length of upper bounds is > length(x); ignoring extra bounds. > In checkbounds at 37 In optim\private\lsqncommon at 39 In lsqcurvefit at 149 In Untitled2 at 6 Exiting due to infeasibility: 4 lower bounds exceed the corresponding upper bounds. g = 1.0000 1.0000 0.3000 1.0000 resnorm = [] 就是我定义的值，我也把公式做了改变的g4=exp(x2.*x2.*((((g(1)./(R.*T)).........) |
7楼2017-07-09 18:40:34
2楼2017-07-09 17:29:24
3楼2017-07-09 17:49:08
4楼2017-07-09 17:49:43