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

limopcb

请教大神，
; l0 j# C5 [/ X
3 J% Y# j; i$ |# M) k: Q求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！

. A' x# a8 V! P- M2 zclc;
clear all;
- }0 y9 ]: \: ^2 A/ m, Fclose all;
syms v1 v2 v4 P3ref
syms i12 i13 i14 i24 i34

r = 0.01;   
5 g% n$ U, J' M, ^) e) ^8 f% TR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;

P1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
' ], k  Y8 [1 {" o' {& tv3 = v3ref;
& K" \! L9 Q" l9 ~
8 O  r4 }( A3 l# l# oeq1 = P1ref/v1 - (i13 + i12 + i14);
& i7 b5 W4 h9 w( qeq2 = P2ref/v2 - (-i12 + i24);
eq3 = P3ref/v3 - (i13 - i34);
eq4 = P4ref/v4 - (i34 + i14 + i24);
& [9 G5 W3 X) `eq5 = (v1 - v3)/R13 - i13;
eq6 = (v1 - v4)/R14 - i14;
eq7 = (v1 - v2)/R12 - i12;
& q2 Z, m$ v( ~$ f, i& x: ?6 Zeq8 = (v2 - v4)/R24 - i24;
  R: o2 J9 d1 J$ ]4 |) i, geq9 = (v3 - v4)/R34 - i34;

8 L, ]) x: @# m; [s = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
v1: 101005.230297198
i13: 804.184237758687
i12: 314.30767957579
0 U/ c5 ^# R* R6 Z% d7 S4 ti14: 614.09162794411
v2: 100769.499537517
i24: 661.634997470918
p3ref: 67678254.7078767
) Z7 I2 K8 o5 X2 A# o6 A/ W/ si34: 127.40169067992
& v6 E. E. O7 W3 Zv4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
  p2 U, k2 m6 ?% J: l0 gsyms i12 i13 i14 i24 i34

r = 0.01;   
. `- d4 h) |/ p6 u8 ~9 [R12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
5 B3 r: K) N& {' H
/ k. |: [/ v' ]% P9 g- yP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
6 [+ }  q" Q! @- ^- rv3 = v3ref;

eq1 = P1ref/v1 - (i13 + i12 + i14)==0;
eq2 = P2ref/v2 - (-i12 + i24)==0;
eq3 = P3ref/v3 - (i13 - i34)==0;
5 X( x1 j! E+ z3 Q4 [% qeq4 = P4ref/v4 - (i34 + i14 + i24)==0;
: L- b' A1 U2 F1 oeq5 = (v1 - v3)/R13 - i13==0;
' {: h- Z* N: l/ meq6 = (v1 - v4)/R14 - i14==0;
eq7 = (v1 - v2)/R12 - i12==0;
/ |4 t# d4 O4 Deq8 = (v2 - v4)/R24 - i24==0;
( a6 A3 l+ Q3 s. m& w: @eq9 = (v3 - v4)/R34 - i34==0;
eqns=[eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9]
' X8 I/ p/ W* @9 Mvars=[v1,v2,v4,P3ref,i12,i13,i14,i24,i34]

s = vpasolve(eqns,vars)