版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

greetrunners

金虫 (著名写手)

应助: 7 (幼儿园)
金币: 487.8
散金: 544
红花: 6
帖子: 1134
在线: 190.1小时
虫号: 1632555
注册: 2012-02-21
性别: GG
专业: 分离过程

[求助] 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];

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请大侠帮忙看一下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人回复

努力！努力！努力！

1楼2013-01-30 10:48:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

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没定义啊,怎么就开始用了

赞一下

回复此楼

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

2楼2013-01-30 17:17:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

greetrunners

金虫 (著名写手)

应助: 7 (幼儿园)
金币: 487.8
散金: 544
红花: 6
帖子: 1134
在线: 190.1小时
虫号: 1632555
注册: 2012-02-21
性别: GG
专业: 分离过程

引用回帖:

2楼: Originally posted by libralibra at 2013-01-30 17:17:00
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，是 ...

应该如何定义？

回复此楼

努力！努力！努力！

3楼2013-01-30 19:08:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

引用回帖:

3楼: Originally posted by greetrunners at 2013-01-30 19:08:17
应该如何定义？...

这...你自己知道吧?
起码函数的定义应该是

function ret_value = foo(input_value)
param1 = ...
param2 = ...
% then use param1 and param2 to calculate the others
...
end

你得先定义参数,然后使用,你的函数前面几行看着写的挺清楚,可是那是注释,你没有定义a2,a3....几乎所有要用的参数

赞一下

回复此楼

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

4楼2013-01-30 20:43:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖