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

limopcb

; I' x* n6 n; Y0 D4 v
0 y" p( ]" L: E; @3 H请教大神，

3 _. s; c- R6 Z/ x求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！

clc;
clear all;
close all;
4 f8 A1 T+ h- hsyms v1 v2 v4 P3ref
' j" c' L. m! u9 D2 |' Gsyms i12 i13 i14 i24 i34
" U5 I/ Z8 J3 ]# j1 J" f
r = 0.01;   
0 C9 Y, R6 K* x: M4 n9 I9 Q/ tR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
% e1 i1 O3 \: }1 @' A. F
P1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
/ z* _2 i, D$ f  h9 _v3 = v3ref;
2 t8 n& W0 ?+ r( l/ H: R  R
eq1 = P1ref/v1 - (i13 + i12 + i14);
eq2 = P2ref/v2 - (-i12 + i24);
eq3 = P3ref/v3 - (i13 - i34);
eq4 = P4ref/v4 - (i34 + i14 + i24);
$ C, m8 M% Q! _9 V1 h5 I' yeq5 = (v1 - v3)/R13 - i13;
7 T  O' p! ~4 h4 q+ k3 G! R! Neq6 = (v1 - v4)/R14 - i14;
( j: B& d2 j! l  m8 Geq7 = (v1 - v2)/R12 - i12;
eq8 = (v2 - v4)/R24 - i24;
  `, @5 Y/ M8 ~eq9 = (v3 - v4)/R34 - i34;

+ p5 p. r% X1 |) r8 G4 S3 B2 Cs = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
v1: 101005.230297198
/ i2 ?% V9 l# t+ r$ ], m: S4 g/ S- G, fi13: 804.184237758687
i12: 314.30767957579
i14: 614.09162794411
4 c" K9 z! Q# `; o6 |v2: 100769.499537517
! S: \+ o" l6 ?) b4 ^0 ?i24: 661.634997470918
p3ref: 67678254.7078767
i34: 127.40169067992
: y8 i( I5 v6 _3 n9 c; _v4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
0 |8 @% H/ l/ h# t; N9 asyms i12 i13 i14 i24 i34
1 c/ {' z9 I8 \
' u) B. I$ O' M6 F; R6 C" Ur = 0.01;   
/ N- T9 J8 H2 y# m. iR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;

) C: n1 Z  ]6 g( `2 g) [P1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
v3 = v3ref;

. S- c" h3 S* X( `& o* Seq1 = P1ref/v1 - (i13 + i12 + i14)==0;
eq2 = P2ref/v2 - (-i12 + i24)==0;
/ p) M. D( w4 F& `) x  z# K1 d: zeq3 = P3ref/v3 - (i13 - i34)==0;
* H7 b' o, H/ A* e+ Beq4 = P4ref/v4 - (i34 + i14 + i24)==0;
7 I# V) f" w. `% N& g* J+ Weq5 = (v1 - v3)/R13 - i13==0;
eq6 = (v1 - v4)/R14 - i14==0;
eq7 = (v1 - v2)/R12 - i12==0;
' e8 c3 R9 @& v' Zeq8 = (v2 - v4)/R24 - i24==0;
3 E( M* h( f1 @6 yeq9 = (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]

/ n( V3 `" t/ M1 G; Is = vpasolve(eqns,vars)