2016-12-31 32 views
0

私はL1> 0、L2、10非線形方程式を持ってMATLABを使用して非線形代数方程式の解を見つけるには?

L1 + 1.3*3*(P1^0.3) + 2*P1 = 12 
L2 + 1.2*5*(P2^0.2) + 3*P2 = 20 
L3 + 1.15*6*(P3^0.15) + 5*P3 = 28 
L4 - L1*0.9*0.4*(X1^-0.1) = 0 
L4 - L2*0.8*0.5*(X2^-0.2) = 0 
L4 - L3*0.7*0.6*(X3^-0.3) = 0 
P1 - 0.4*(X1^0.9) = 0 
P2 - 0.5*(X2^0.8) = 0 
P3 - 0.6*(X3^0.7) = 0 
X1 + X2 + X3 = 10 

私は解決策のための任意の初期推測を持っていない、しかし、すべての変数は、本質的に非負すなわちあり、

> 0、L3> 0、L4> 0、P1> 0、P2> 0、P3> 0、X1> 0、X2> 0、X3> 0

私はこれらの方程式を以下のコマンド、

clear 
clc 

syms L1 L2 L3 L4 P1 P2 P3 X1 X2 X3 

sol=solve([ L1 + 1.3*3*(P1^0.3) + 2*P1 == 12, L2 + 1.2*5*(P2^0.2) + 3*P2 == 20, L3 + 1.15*6*(P3^0.15) + 5*P3 == 28, ... 
L4 - L1*0.9*0.4*(X1^-0.1) == 0, L4 - L2*0.8*0.5*(X2^-0.2) == 0, L4 - L3*0.7*0.6*(X3^-0.3) == 0, ... 
P1 - 0.4*(X1^0.9) == 0, P2 - 0.5*(X2^0.8) == 0, P3 - 0.6*(X3^0.7) ==0, X1+X2 +X3 == 10, ... 
L1>0, L2>0, L3>0, L4>0, P1>0, P2>0, P3>0, X1>0, X2>0, X3>0], [L1, L2, L3, L4, P1, P2, P3, X1, X2, X3]); 

しかし、これらの非線形方程式を解く方法

Warning: 32 equations in 11 variables. 
> In C:\Program Files\MATLAB\R2013a\toolbox\symbolic\symbolic\symengine.p>symengine at 56 
    In mupadengine.mupadengine>mupadengine.evalin at 97 
    In mupadengine.mupadengine>mupadengine.feval at 150 
    In solve at 170 
Warning: Explicit solution could not be found. 
> In solve at 179 
>> 

としてエラーを示して?

更新:DIFFは実際には数字( ,10,20、または30)の違いです。私はここからDIFFを削除しました。

+0

[this](https://nl.mathworks.com/matlabcentral/answers/60852-warning-explicit-solution-could-not-be-found)の質問をご覧ください。あなたの問題を再検討する必要があるかもしれません。 –

答えて

0

あなたには10個の式がありますが、変数は10個です。

手助けができます。

P1 = 0.4(x^(0.9)) 
P2 = 0.5(x^(0.8)) 
P3 = 0.6(x^(0.7)) 

L4 = 0.36L1*X1^(-0.1) = 0.4L2*X2^(-0.2) = 0.42L3*X3^(-0.3) 

は、我々は4つの変数に4次方程式を得る: 線形方程式を削減

L1 + 3.9*((0.4^(0.3))*X1^(0.27)) + 0.8*X1^(0.9) - 12 = 0 
0.9*(X1^(-0.1))*(X2^(0.2))*L1 + 6*((0.5^0.2)*(x2^(0.16))) + 1.5*X2^(0.8) - 20 = 0 
(6/7)*(X1^(-0.1))*(X3^(0.3))*L1 + 6.9*((0.6^0.15)*(X3^(0.105))) + 3*(x3^0.7) - 28 = 0 
x1 + X2 + X3 = 10 

と私がしようと

sol=solve([X1+X2+X3 == 10,L1 + 3.9*((0.4^(0.3))*X1^(0.27)) + 0.8*X1^(0.9) - 12 == 0, 0.9*(X1^(-0.1))*(X2^(0.2))*L1 + 6*((0.5^0.2)*(X2^(0.16))) + 1.5*X2^(0.8) - 20 == 0,(6/7)*(X1^(-0.1))*(X3^(0.3))*L1 + 6.9*((0.6^0.15)*(X3^(0.105))) + 3*(X3^0.7) - 28 == 0, X1 > 0, X2 > 0, X3 > 0, L1 > 0]); 

を解決しようそれを解決する、bu私は記憶がない。 スマートな方法が欲しい。

+0

これは半分の答えです。少なくとも1つの方法を使用して完全なソリューションを提示してください。 –

+0

@petr質問を更新しました。あなたは今助けてください。 – dSb

関連する問題