[求助]matlab解非线性方程组出错

mathtype · 发表于 2006-2-27 22:48

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

用matlab解非线性方程组出错了，出错信息如下: ?? In an assignment A(I) = B, the number of elements in B and I must be the same. Error in ==> thermal>NonlinEqs at 24 f(4)=x(3).^2.*x(4).^2-K1.*x(1).*x(2)./(P./xtotal).^2; Error in ==> fsolve at 180 fuser = feval(funfcn{3},x,varargin{:}); Error in ==> thermal at 14 x1 =fsolve(@NonlinEqs,x0,[],K1,K2,K3,n10,n20) 程序如下 function xFsolve % xFsolve.m % 求解热力学平衡计算
clear all clc
n10=0.2;n20=1-n10; T=500:100:1500; %体系温度，K K1=exp(7.9305.*log(T)-5.6347e-03.*T+5.8622e-07.*T.^2-2.7848e+04./T-19.03390284); K2=exp(0.248617.*log(T)-1.4693e-03.*T+2.28e-07.*T.^2-4.1539e+04./T+4.213066); K3=exp(6.833053.*log(T)-3.7606e-03.*T+3.48e-07.*T.^2-7.276675e+03./T-33.8969); x0 = [n10 n20 0 0 0 0]'; x1 =fsolve(@NonlinEqs,x0,[],K1,K2,K3,n10,n20)
 % ------------------------------------------------------------------ function f = NonlinEqs(x,K1,K2,K3,n10,n20) P=101325; xtotal=x(1)+x(2)+x(3)+x(4)+x(5)+x(6); f(1) = x(1)+x(2)+x(3)+x(6)-1; f(2) = 2*x(1)+x(4)+x(5)-2*n10; f(3)=2.*x(2)+x(3)+x(5)-2*n20; f(4)=x(3).^2.*x(4).^2-K1.*x(1).*x(2)./(P./xtotal).^2; f(5)=x(3).*x(5)-K2.*x(2).*x(4); f(6)=x(6).*x(4).^2-K3.*x(1);

happy · 发表于 2006-2-28 08:44

<DIV class=quote>以下是引用mathtype在2006-2-27 22:48:38的发言： 
用matlab解非线性方程组出错了，出错信息如下: ?? In an assignment A(I) = B, the number of elements in B and I must be the same. Error in ==> thermal>NonlinEqs at 24 f(4)=x(3).^2.*x(4).^2-K1.*x(1).*x(2)./(P./xtotal).^2; Error in ==> fsolve at 180 fuser = feval(funfcn{3},x,varargin{:}); Error in ==> thermal at 14 x1 =fsolve(@NonlinEqs,x0,[],K1,K2,K3,n10,n20) 程序如下 function xFsolve % xFsolve.m % 求解热力学平衡计算
clear all clc
n10=0.2;n20=1-n10; T=500:100:1500; %体系温度，K K1=exp(7.9305.*log(T)-5.6347e-03.*T+5.8622e-07.*T.^2-2.7848e+04./T-19.03390284); K2=exp(0.248617.*log(T)-1.4693e-03.*T+2.28e-07.*T.^2-4.1539e+04./T+4.213066); K3=exp(6.833053.*log(T)-3.7606e-03.*T+3.48e-07.*T.^2-7.276675e+03./T-33.8969); x0 = [n10 n20 0 0 0 0]'; x1 =fsolve(@NonlinEqs,x0,[],K1,K2,K3,n10,n20)
 % ------------------------------------------------------------------ function f = NonlinEqs(x,K1,K2,K3,n10,n20) P=101325; xtotal=x(1)+x(2)+x(3)+x(4)+x(5)+x(6); f(1) = x(1)+x(2)+x(3)+x(6)-1; f(2) = 2*x(1)+x(4)+x(5)-2*n10; f(3)=2.*x(2)+x(3)+x(5)-2*n20; f(4)=x(3).^2.*x(4).^2-K1.*x(1).*x(2)./(P./xtotal).^2; f(5)=x(3).*x(5)-K2.*x(2).*x(4); f(6)=x(6).*x(4).^2-K3.*x(1); </DIV>
 x(3).^2.*x(4).^2 是一个数 K1.*x(1).*x(2)./(P./xtotal).^2; 是向量，两者无法相加

realyyy · 发表于 2006-5-15 14:44

高！

账号		自动登录	找回密码
密码			我要加入

[求助]matlab解非线性方程组出错

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（mathtype）[求助]matlab解非线性方程组出错