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

limopcb

% P1 H& v2 h5 r# i7 ?* {4 n
! d$ w2 M( u3 R. z9 z% W! ~请教大神，
. S+ K6 W3 R* d' I* A: O  I
求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！

$ q  W9 `/ L. o' c5 `3 kclc;
* M2 _4 {# W2 Aclear all;
0 s# i. J. a. d/ j8 a4 nclose all;
. w: `! L0 w- m# r5 Z; d0 Bsyms v1 v2 v4 P3ref
syms i12 i13 i14 i24 i34

r = 0.01;   
R12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;

P1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
v3 = v3ref;

eq1 = P1ref/v1 - (i13 + i12 + i14);
eq2 = P2ref/v2 - (-i12 + i24);
/ P2 E( [" U+ l0 }% leq3 = P3ref/v3 - (i13 - i34);
, H; i) o6 e0 m* ~eq4 = P4ref/v4 - (i34 + i14 + i24);
; t4 \! [0 h1 ~3 i0 weq5 = (v1 - v3)/R13 - i13;
eq6 = (v1 - v4)/R14 - i14;
! g( j" ^5 j5 D$ ]5 M" `eq7 = (v1 - v2)/R12 - i12;
eq8 = (v2 - v4)/R24 - i24;
eq9 = (v3 - v4)/R34 - i34;
$ f, N! _0 q5 }/ x5 w
* s2 @6 v% c' T% Gs = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
v1: 101005.230297198
7 O/ J3 P( o5 ii13: 804.184237758687
i12: 314.30767957579
% c7 E8 [' ?9 l% C3 Bi14: 614.09162794411
3 O1 {5 x, n% g' k: _- M/ iv2: 100769.499537517
. n; R/ U3 c2 @/ |0 E; ni24: 661.634997470918
" @, _" f8 ?% D  F/ @* bp3ref: 67678254.7078767
i34: 127.40169067992
v4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
3 k6 V) Q4 p8 i- D/ Csyms i12 i13 i14 i24 i34
1 J/ ^2 q- G4 {* l
r = 0.01;   
  q' Y) ~+ I# P4 [4 J; C/ aR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;

1 g0 _: I- ^. V6 t. Q. sP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
- ]; m* Z9 ?- W4 I  R3 e6 ~/ X9 tv3 = v3ref;
4 U8 d. X& S* g2 w5 k' w
8 l- Z8 @2 N; H$ q. e: Eeq1 = P1ref/v1 - (i13 + i12 + i14)==0;
eq2 = P2ref/v2 - (-i12 + i24)==0;
) n: Q2 z- D' `4 U& c2 I4 weq3 = P3ref/v3 - (i13 - i34)==0;
/ G3 i& i# P8 c6 N7 d1 ]eq4 = P4ref/v4 - (i34 + i14 + i24)==0;
% @) K1 `6 u5 ?2 ^  ]9 H4 Teq5 = (v1 - v3)/R13 - i13==0;
eq6 = (v1 - v4)/R14 - i14==0;
eq7 = (v1 - v2)/R12 - i12==0;
' P3 O# I! ]; x0 teq8 = (v2 - v4)/R24 - i24==0;
eq9 = (v3 - v4)/R34 - i34==0;
eqns=[eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9]
& i2 I) [, v& Svars=[v1,v2,v4,P3ref,i12,i13,i14,i24,i34]
( v8 c' x3 \$ i; A) n
s = vpasolve(eqns,vars)