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

limopcb

& n+ F, m: R9 z/ _
请教大神，
; Y/ N* Y: Q4 E9 i- X
求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！
. T7 s0 U  Z. P- O) i: w
clc;
clear all;
3 K4 j- m7 e+ C% Q/ k' j9 }close all;
syms v1 v2 v4 P3ref
0 r% {% s: G; z, qsyms i12 i13 i14 i24 i34

* V! X+ _- k: n& a! Fr = 0.01;   
R12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
5 h  ~# B& Q" _/ E+ p( o) |
# N& C- `8 v2 oP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
7 j  L. Y& m) K* E2 Lv3 = v3ref;

eq1 = P1ref/v1 - (i13 + i12 + i14);
eq2 = P2ref/v2 - (-i12 + i24);
eq3 = P3ref/v3 - (i13 - i34);
eq4 = P4ref/v4 - (i34 + i14 + i24);
eq5 = (v1 - v3)/R13 - i13;
, W, t2 e% y0 C7 a# ?7 ?  geq6 = (v1 - v4)/R14 - i14;
9 F8 N! U4 d. a. q. j6 z7 s( Weq7 = (v1 - v2)/R12 - i12;
6 y( E/ q2 p8 |7 peq8 = (v2 - v4)/R24 - i24;
eq9 = (v3 - v4)/R34 - i34;

  o) w+ x% j; k! j8 z$ b! `4 as = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
( v) h# G( Y3 W! i7 @9 {v1: 101005.230297198
i13: 804.184237758687
i12: 314.30767957579
$ e3 W* \# j6 Q% S2 `i14: 614.09162794411
v2: 100769.499537517
i24: 661.634997470918
; N, Z" D$ l9 R9 a5 ~p3ref: 67678254.7078767
i34: 127.40169067992
v4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
9 m$ `; e; L& u$ a( S* e! p: }syms i12 i13 i14 i24 i34

( b  B2 N4 r8 _. E! v) k( v+ v* or = 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;

2 O& _" m, }8 H* Z* L2 ^# b/ Y' beq1 = P1ref/v1 - (i13 + i12 + i14)==0;
eq2 = P2ref/v2 - (-i12 + i24)==0;
0 d( g- z: f7 t6 O9 xeq3 = P3ref/v3 - (i13 - i34)==0;
eq4 = P4ref/v4 - (i34 + i14 + i24)==0;
eq5 = (v1 - v3)/R13 - i13==0;
! B- O% U/ P8 h# ^6 K  Xeq6 = (v1 - v4)/R14 - i14==0;
eq7 = (v1 - v2)/R12 - i12==0;
eq8 = (v2 - v4)/R24 - i24==0;
eq9 = (v3 - v4)/R34 - i34==0;
0 r9 Z" e% |. b" _9 xeqns=[eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9]
% x7 G" \* z: Q% z$ V8 G3 cvars=[v1,v2,v4,P3ref,i12,i13,i14,i24,i34]

  A! X, w3 r8 Us = vpasolve(eqns,vars)