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

limopcb

/ @$ u% C, Z9 ?5 p3 j
请教大神，
4 B/ z. T3 l6 f% B. y
求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！

" Y' f! @2 r( R0 |  X0 gclc;
clear all;
8 A/ D+ x- V% S1 S  p7 M( Q6 k0 b0 L7 @) kclose all;
# c% ]3 R8 U# {. y8 `syms v1 v2 v4 P3ref
! L4 p- U: M8 M, psyms i12 i13 i14 i24 i34

# n' m, J$ r4 e( Gr = 0.01;   
( t5 ?7 \# q7 t2 B" vR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;

  p/ S. J# F! Y6 ^) j9 ~# FP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
0 T, I1 `/ I* W" T  [$ S! Ov3 = v3ref;

# t% z, m0 T2 d" `: [; Q' S# Ueq1 = P1ref/v1 - (i13 + i12 + i14);
: Z2 x( U1 Y% Xeq2 = P2ref/v2 - (-i12 + i24);
eq3 = P3ref/v3 - (i13 - i34);
eq4 = P4ref/v4 - (i34 + i14 + i24);
eq5 = (v1 - v3)/R13 - i13;
eq6 = (v1 - v4)/R14 - i14;
( w- p2 M9 c: q8 y; meq7 = (v1 - v2)/R12 - i12;
" v, E8 V) X0 s1 \' }eq8 = (v2 - v4)/R24 - i24;
eq9 = (v3 - v4)/R34 - i34;
( p. h% O+ F2 G0 d
s = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
. ]/ _+ n5 `5 b4 a. Y% e- A* D1 wv1: 101005.230297198
i13: 804.184237758687
0 B! H2 G5 ^6 L" ki12: 314.30767957579
, _: f, T- z3 w- t$ ki14: 614.09162794411
v2: 100769.499537517
& _! D; f0 t: w1 i4 z; h+ bi24: 661.634997470918
( w% V3 L5 f) c! |6 q2 h/ g" lp3ref: 67678254.7078767
, w% e( g, q0 Z4 mi34: 127.40169067992
5 A& B. @: T) \v4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
7 C5 ~: M( ~% n+ [* i$ s% \syms i12 i13 i14 i24 i34
, B! [; o) m# t2 d7 P; y7 Y! B
  a% H9 B5 Q5 G2 ^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;
7 ]" a+ b7 O/ g& b3 k, K
2 K$ S8 {$ F+ k; }$ heq1 = P1ref/v1 - (i13 + i12 + i14)==0;
eq2 = P2ref/v2 - (-i12 + i24)==0;
+ [6 C# N/ G+ S4 `! \eq3 = P3ref/v3 - (i13 - i34)==0;
6 J0 _7 R- R2 G! Deq4 = P4ref/v4 - (i34 + i14 + i24)==0;
+ `& a+ o  X, [4 Neq5 = (v1 - v3)/R13 - i13==0;
1 @2 h' f8 i  N7 x! Leq6 = (v1 - v4)/R14 - i14==0;
( v. u& \+ R; Y$ Teq7 = (v1 - v2)/R12 - i12==0;
6 U/ ?  u; [% @+ E5 @4 o  Ceq8 = (v2 - v4)/R24 - i24==0;
5 u; D2 H% w6 Y* r: w! J! seq9 = (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]
' i# Q. ?2 J0 c  t% T! }
7 L) f8 ?7 C/ ^7 ms = vpasolve(eqns,vars)