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

limopcb

请教大神，

" A0 R& t4 z" ^) v/ F) t求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！

9 I: Y  r( F0 n! }, cclc;
clear all;
close all;
" Y3 X/ Y/ E  Ksyms v1 v2 v4 P3ref
2 H9 U# N9 ~7 r6 v9 ?syms i12 i13 i14 i24 i34
2 a* F0 C8 k9 x/ ~8 m/ |
/ r( o, T& }9 ?& d5 ]# O2 ^. x, t7 Ur = 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;
$ c& u* e; s/ B9 L
5 E; }4 v3 s9 g5 ]eq1 = P1ref/v1 - (i13 + i12 + i14);
eq2 = P2ref/v2 - (-i12 + i24);
& k3 r* P8 z  O4 O1 ?3 N! Heq3 = P3ref/v3 - (i13 - i34);
eq4 = P4ref/v4 - (i34 + i14 + i24);
eq5 = (v1 - v3)/R13 - i13;
, H4 ?/ D& @9 R8 M. y  Eeq6 = (v1 - v4)/R14 - i14;
$ ], f8 t7 I* w/ A- S* W$ Neq7 = (v1 - v2)/R12 - i12;
eq8 = (v2 - v4)/R24 - i24;
eq9 = (v3 - v4)/R34 - i34;
; c" \  l: j9 J, q2 s, }& E9 }. g
s = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
v1: 101005.230297198
i13: 804.184237758687
i12: 314.30767957579
i14: 614.09162794411
0 |* T; y, D. }; Kv2: 100769.499537517
i24: 661.634997470918
/ g5 Q* }' Y# k$ ?" [2 y2 Gp3ref: 67678254.7078767
i34: 127.40169067992
v4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
; P  i6 L7 ?8 z/ F; N" Qsyms i12 i13 i14 i24 i34
2 i+ l6 g+ t1 K' z) J
5 }: X3 `9 K% I' Jr = 0.01;   
R12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
  u- r! H* ^: \$ W$ L
P1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
v3 = v3ref;
& Y; N9 h' E8 d. S1 q
7 |: Y" a3 H$ }! W! H) ?* Oeq1 = P1ref/v1 - (i13 + i12 + i14)==0;
! `0 Y! W) r4 r( Teq2 = P2ref/v2 - (-i12 + i24)==0;
eq3 = P3ref/v3 - (i13 - i34)==0;
eq4 = P4ref/v4 - (i34 + i14 + i24)==0;
eq5 = (v1 - v3)/R13 - i13==0;
eq6 = (v1 - v4)/R14 - i14==0;
! g1 H- N4 a1 ~9 O% Q/ Keq7 = (v1 - v2)/R12 - i12==0;
: R4 r  F0 Q1 Req8 = (v2 - v4)/R24 - i24==0;
$ Y( T6 x* ?* M8 H) _- Oeq9 = (v3 - v4)/R34 - i34==0;
8 b+ Z! V+ [6 J' e3 y2 Ueqns=[eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9]
vars=[v1,v2,v4,P3ref,i12,i13,i14,i24,i34]
& l/ k- W+ p3 B! `5 z6 ~
s = vpasolve(eqns,vars)