| 查看: 1238 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
greetrunners金虫 (著名写手)
|
[求助]
MATLAB, 代码错误，请大侠帮忙看看。
|
|
|
[t12 t21 t13 t31 t23 t32]=lsqnonlin(@lsqfun3,1,1,1,1,1,1) Undefined function or variable 'a2'. Error in lsqfun3 (line 8) k1=(a2.^2*t21*exp(-0.3*t21)*exp(-0.3*t21)+a3.^2*t31*exp(-0.3*t31)*exp(-0.3*t31)+a2.*a3*exp(-0.3*t21)*exp(-0.3*t31)*(t21+t31))./((a1+a2.*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2); Error in lsqnonlin (line 197) initVals.F = feval(funfcn{3},xCurrent,varargin{:}); Caused by: Failure in initial user-supplied objective function evaluation. LSQNONLIN cannot continue. 为了减少输出的参数，编写以下函数，出现了上述问题。 function obj=lsqfun3(t12,t21,t13,t31,t23,t32) %t12,t21,t13,t31,t23,t32，是所需要关联的NRTL方程参数； %a1,a2,a3,是组分1,2,3，液相平衡数据； %y1,y2,y3,是组分1,2,3，气相平衡数据； %ps1，ps2，ps3，是组分1,2,3的饱和蒸汽压； %p==101.35; p=101.35; k1=(a2.^2*t21*exp(-0.3*t21)*exp(-0.3*t21)+a3.^2*t31*exp(-0.3*t31)*exp(-0.3*t31)+a2.*a3*exp(-0.3*t21)*exp(-0.3*t31)*(t21+t31))./((a1+a2.*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2); k2=(a2.^2*t12*exp(-0.3*t12)+a2.*a3*exp(-0.3*t12)*exp(-0.3*t32)*(t12-t32))./((a1*exp(-0.3*t12)+a2+a3*exp(-0.3*t32)).^2); k3=(a3.^2*t13*exp(-0.3*t13)+a2.*a3*exp(-0.3*t13)*exp(-0.3*t23)*(t13-t23))./((a1*exp(-0.3*t13)+a2*exp(-0.3*t23)+a3).^2); k4=(a1.^2*t21*exp(-0.3*t21)+a1.*a3*exp(-0.3*t21)*exp(-0.3*t31)*(t21-t31))./((a1+a2*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2); k5=(a1.^2*t12*exp(-0.3*t12)^2+a3.^2*t32*exp(-0.3*t32)^2+a1.*a3*exp(-0.3*t12)*exp(-0.3*t32)*(t12+t32))./((a1*exp(-0.3*t12)+a2+a3*exp(-0.3*t32)).^2); k6=(a3.^2*t23*exp(-0.3*t23)+a1.*a3*exp(-0.3*t13)*exp(-0.3*t23)*(t23-t13))./((a1*exp(-0.3*t13)+a2*exp(-0.3*t23)+a3).^2); k7=(a1.^2*t31*exp(-0.3*t31)+a1.*a2*exp(-0.3*t21)*exp(-0.3*t31)*(t31-t21))./((a1+a2*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2); k8=(a2.^2*t32*exp(-0.3*t32)^2+a1.*a2*exp(-0.3*t12)*exp(-0.3*t32)*(t32-t12))./((a1*exp(-0.3*t12)+a2+a3*exp(-0.3*t32)).^2); k9=(a1.^2*t13*exp(-0.3*t13)^2+a2.^2*t23*exp(-0.3*t23)^2+a1.*a2*exp(-0.3*t13)*exp(-0.3*t23)*(t13+t23))./((a1*exp(-0.3*t13)+a2*exp(-0.3*t23)+a3).^2); gama1=exp(k1+k2+k3); gama2=exp(k4+k5+k6); gama3=exp(k7+k8+k9); y1c=ps1.*gama1.*a1./p; y2c=ps2*gama2*a2./p; y3c=ps3*gama3*a3./p; obj=sum(sum((y1c-y1).^2+(y2c-y2).^2+(y3c-y3).^2)); >> y1=[0.8907 0.8924 0.733 0.7731 0.698 0.7778 0.7622 0.7796 0.7934]; y2=[0.1082 0.1067 0.6626 0.5623 0.732 0.4326 0.3978 0.3766 0.3473]; y3=[0 0 0 0 0 0 0 0 0]; a1=[0.8918 0.8933 0.3374 0.4377 0.268 0.5674 0.6022 0.6234 0.6527]; a2=[0.10723 0.10649 0.62613 0.53391 0.70928 0.42518 0.39054 0.37007 0.34364]; a3=[0.00097 0.00021 0.03647 0.02839 0.02272 0.00742 0.00726 0.00653 0.00366]; |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
今天维护系统维护 祝所有人 高中
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有98人回复
-
2027广东省杰青
已经有8人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有5人回复
-
投稿咨询
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有6人回复
-
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请大侠帮忙看一下matlab利用lsqcurvefit曲线拟合问题
已经有5人回复
- 请大侠帮忙看一下我这个matlab小程序结果为什么不返回矩阵
已经有7人回复
- MATLAB关于lsqcurvefit程序几个错误，求大神指导一下
已经有6人回复
- matlab中关于非线性拟合，用nlinfit.程序总出错。谢谢各位大侠
已经有3人回复
- 用matlab 的 fminu 函数做非线性最优化 出现错误， 求高手帮忙
已经有4人回复
- 帮忙给看看这个matlab优化函数 问题
已经有8人回复
- 帮忙修改Matlab中出现的错误
已经有4人回复
- MATLAB运行错误
已经有5人回复
- 采用matlab求解微分方程遇到问题，请帮忙分析一下错误-重谢
已经有3人回复
- 下载了MATLAB2012可是安装不了。大侠帮帮忙谢谢
已经有13人回复
- 怎么用matlab画极坐标三维图？各位大侠帮帮忙啊！
已经有5人回复
- 求助matlab一程序运行的问题，帮忙看看哪里不对
已经有4人回复
- 求大侠帮忙解决一个B类错误！！谢谢
已经有5人回复
- 最优化计算求助！（1stopt或者Matlab）
已经有16人回复
- matlab 问题求助，是否计算结果错误，如何改正？
已经有3人回复
- matlab程序错误
已经有5人回复
- 用matlab求解方程出问题，请帮忙看看
已经有3人回复
- 帮忙看看这个平面图用matlab怎么画
已经有19人回复
- 求助大侠帮忙，有关matlab模拟数据
已经有29人回复
- matlab程序运行总是出现下面的错误，请各位高手指点一下，谢谢了！急 急急
已经有5人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB错误
已经有6人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
greetrunners
金虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 487.8
- 散金: 544
- 红花: 6
- 帖子: 1134
- 在线: 190.1小时
- 虫号: 1632555
- 注册: 2012-02-21
- 性别: GG
- 专业: 分离过程
3楼2013-01-30 19:08:17
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
greetrunners: 金币+10, ★★★很有帮助, 谢谢 2013-01-31 09:39:29
xzhdty: 金币+1, 谢谢骠骑将军 2013-01-31 11:08:38
感谢参与，应助指数 +1
greetrunners: 金币+10, ★★★很有帮助, 谢谢 2013-01-31 09:39:29
xzhdty: 金币+1, 谢谢骠骑将军 2013-01-31 11:08:38
|
function obj=lsqfun3(t12,t21,t13,t31,t23,t32) %t12,t21,t13,t31,t23,t32，是所需要关联的NRTL方程参数； %a1,a2,a3,是组分1,2,3，液相平衡数据； %y1,y2,y3,是组分1,2,3，气相平衡数据； %ps1，ps2，ps3，是组分1,2,3的饱和蒸汽压； %p==101.35; p=101.35; k1=(a2.^2*t21*exp(-0.3*t21)*exp(-0.3*t21)+a3.^2*t31*exp(-0.3*t31)*exp(-0.3*t31)+a2.*a3*exp(-0.3*t21) =============== 函数中a2,a3没定义啊,怎么就开始用了 |
2楼2013-01-30 17:17:00
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
4楼2013-01-30 20:43:55
greetrunners
金虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 487.8
- 散金: 544
- 红花: 6
- 帖子: 1134
- 在线: 190.1小时
- 虫号: 1632555
- 注册: 2012-02-21
- 性别: GG
- 专业: 分离过程
5楼2013-01-31 12:14:29
