求解如下方程组，MATLAB算不出来怎么办？

limopcb

$ R" o# y- \+ L1 y
/ p% U8 p1 a+ U& @" n0 z请教大神，

- Q/ c5 b6 h$ [" i" S/ I. A0 g/ J求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！

2 V6 c" D2 q8 y  X- rclc;
1 [+ t' B4 n  Y/ F. Cclear all;
- [; B9 q) n0 M* t1 V6 e* hclose all;
  w2 T' X; Y. s$ B' usyms v1 v2 v4 P3ref
& m( @% p( d# ^syms i12 i13 i14 i24 i34

r = 0.01;   
R12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
' H: m0 p- C# O; w
P1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
v3 = v3ref;
) G' t4 V7 m, m3 F4 ]
) ^6 A/ p% G9 j; {5 Q1 neq1 = P1ref/v1 - (i13 + i12 + i14);
  @6 H' S8 X4 [1 X9 O+ yeq2 = P2ref/v2 - (-i12 + i24);
eq3 = P3ref/v3 - (i13 - i34);
) C  Z4 K2 x; v5 O+ D0 ~eq4 = P4ref/v4 - (i34 + i14 + i24);
eq5 = (v1 - v3)/R13 - i13;
eq6 = (v1 - v4)/R14 - i14;
* W; Z" }+ @/ Q6 Peq7 = (v1 - v2)/R12 - i12;
eq8 = (v2 - v4)/R24 - i24;
eq9 = (v3 - v4)/R34 - i34;
: v6 e6 o0 J0 c- {
s = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
v1: 101005.230297198
# D! _; K5 M& {1 \: {) Q5 Di13: 804.184237758687
i12: 314.30767957579
i14: 614.09162794411
9 ~$ S& p- o! V1 ~v2: 100769.499537517
% p- _' G6 R# z8 e6 k& hi24: 661.634997470918
0 k+ W+ Q+ n6 P! Vp3ref: 67678254.7078767
i34: 127.40169067992
v4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
syms i12 i13 i14 i24 i34
9 `; l' o" ?) k$ w5 j7 ]# j
# r+ E! V. d5 \9 ^( c8 y8 B! nr = 0.01;   
/ `  Z8 S6 T+ b0 r: w4 lR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;

, \! N8 f- G$ NP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
v3 = v3ref;
% S% U6 N# P5 s
eq1 = P1ref/v1 - (i13 + i12 + i14)==0;
; c) n% _# m( L% z6 T- @eq2 = P2ref/v2 - (-i12 + i24)==0;
3 ?. X7 r: H6 ~eq3 = P3ref/v3 - (i13 - i34)==0;
eq4 = P4ref/v4 - (i34 + i14 + i24)==0;
& z  {  Y" ^9 V# Z) Yeq5 = (v1 - v3)/R13 - i13==0;
. ?, c8 u1 g  r2 z* K: k& q" deq6 = (v1 - v4)/R14 - i14==0;
; y4 H& H5 }* s  c0 `eq7 = (v1 - v2)/R12 - i12==0;
eq8 = (v2 - v4)/R24 - i24==0;
. w( d% a3 z( aeq9 = (v3 - v4)/R34 - i34==0;
eqns=[eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9]
vars=[v1,v2,v4,P3ref,i12,i13,i14,i24,i34]
  C5 I7 ~9 j; ?; i* x
s = vpasolve(eqns,vars)