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

limopcb

" S3 h- m+ y3 b5 q. Q请教大神，

求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！
) z. A; B  N, B" d9 d  a
clc;
clear all;
# q2 r7 p/ q5 L9 P( S2 uclose all;
- A# `, z3 Z1 d9 @syms v1 v2 v4 P3ref
syms i12 i13 i14 i24 i34

r = 0.01;   
& N" k8 v$ K, f- \R12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
( ?9 B" }# u2 o* f6 `# u
- n0 ?5 n7 l# q. }- ]/ C. v) V% JP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
v3 = v3ref;

eq1 = P1ref/v1 - (i13 + i12 + i14);
eq2 = P2ref/v2 - (-i12 + i24);
eq3 = P3ref/v3 - (i13 - i34);
eq4 = P4ref/v4 - (i34 + i14 + i24);
" G; J% n9 u; F5 Xeq5 = (v1 - v3)/R13 - i13;
; p7 P# S' l- k7 E, ceq6 = (v1 - v4)/R14 - i14;
eq7 = (v1 - v2)/R12 - i12;
' R; r# P4 _5 A/ u, Meq8 = (v2 - v4)/R24 - i24;
eq9 = (v3 - v4)/R34 - i34;
5 h7 k+ R4 P1 Y; [2 V; i  P( C
; ]* Q( K5 V% V2 I: Es = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
* \" y5 e5 d7 _- Z5 F8 m4 Dv1: 101005.230297198
& k3 s/ M6 f" _3 ~/ i& z+ N  Xi13: 804.184237758687
i12: 314.30767957579
" A" G( Q% K! f- N+ `i14: 614.09162794411
: M+ }- f) w3 n, Tv2: 100769.499537517
i24: 661.634997470918
p3ref: 67678254.7078767
i34: 127.40169067992
; g, X: Z$ Y7 Wv4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
; K, J1 [3 s5 P- i( l1 Dsyms i12 i13 i14 i24 i34
! V% H6 b9 U7 J0 f9 g
r = 0.01;   
/ ~! J5 ^1 T& VR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
8 s( H( {4 T, a  k; Q7 u
; f! a" l1 _, w/ jP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
' h0 k' y; U. ~, L) }v3 = v3ref;
6 x$ a5 L0 z" ]# C) ^! {
" v% h; h$ S) l3 e, W3 neq1 = P1ref/v1 - (i13 + i12 + i14)==0;
; y" K9 q- n$ X4 r6 oeq2 = P2ref/v2 - (-i12 + i24)==0;
5 I3 J1 ?$ z% q7 n. g7 Req3 = P3ref/v3 - (i13 - i34)==0;
eq4 = P4ref/v4 - (i34 + i14 + i24)==0;
eq5 = (v1 - v3)/R13 - i13==0;
eq6 = (v1 - v4)/R14 - i14==0;
. D5 l/ {: K9 D  P* ?& Eeq7 = (v1 - v2)/R12 - i12==0;
eq8 = (v2 - v4)/R24 - i24==0;
eq9 = (v3 - v4)/R34 - i34==0;
eqns=[eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9]
6 S* C. S( Q; _  {5 W1 g, n2 {( tvars=[v1,v2,v4,P3ref,i12,i13,i14,i24,i34]

s = vpasolve(eqns,vars)