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

limopcb

- E2 k( O5 e. H) u
请教大神，

求解如下方程组，九个方程九个未知数，MATLAB算不出来，很奇怪，按道理计算出来应该没问题的。感谢！

3 p# n( _7 }" ]. j; vclc;
clear all;
/ E6 t3 o% b3 s- _0 N; g; Mclose all;
% j. i" E# E0 b/ x, y7 k5 E% Tsyms v1 v2 v4 P3ref
6 Y- j4 B8 g. A" }& }syms i12 i13 i14 i24 i34
. S. S& z) {0 w
' G# t6 Q! m. C" c, @1 e9 |' z) jr = 0.01;   
% y* d  r# R% oR12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;

# L$ Q) Q7 Z! \3 Q/ qP1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
v3 = v3ref;

eq1 = P1ref/v1 - (i13 + i12 + i14);
eq2 = P2ref/v2 - (-i12 + i24);
1 ?* c, D* ?* Keq3 = P3ref/v3 - (i13 - i34);
eq4 = P4ref/v4 - (i34 + i14 + i24);
eq5 = (v1 - v3)/R13 - i13;
eq6 = (v1 - v4)/R14 - i14;
: w4 I- R, Q, u* C5 N+ c. weq7 = (v1 - v2)/R12 - i12;
eq8 = (v2 - v4)/R24 - i24;
eq9 = (v3 - v4)/R34 - i34;

: o/ s7 k" n7 T( |4 @2 l  t, Ms = solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9);

House

数值解试试：
v1: 101005.230297198
  a6 Z' w1 C) `. G7 D8 C* wi13: 804.184237758687
- ]/ O' x+ Y9 c! M" U' G3 ji12: 314.30767957579
. W8 \: Q6 J* G: u" N, Hi14: 614.09162794411
v2: 100769.499537517
i24: 661.634997470918
p3ref: 67678254.7078767
. l, |9 P3 ^# }2 k/ m9 W; Pi34: 127.40169067992
v4: 99777.0470413101

artic

syms v1 v2 v4 P3ref
5 B  ?# b' s. ^, t; e1 U5 G' psyms i12 i13 i14 i24 i34
0 S% ]! z( ?  W4 O! W# I/ d
* Z; D1 {( Q$ g' S5 Kr = 0.01;   
R12 = 75*r; R13 = 125*r; R14 = 200*r; R24 = 150*r; R34 = 175*r;
6 K( m# z8 q+ O! Q
, I! h+ B$ L0 f3 q4 k. i5 @P1ref = 175e6; P2ref = 35e6; P4ref = 140e6; v3ref = 100e3;
9 ~+ F, O# S. t/ D8 Q. Hv3 = v3ref;
& T9 B: |( J; g
7 j; f- o3 i* R4 W) keq1 = P1ref/v1 - (i13 + i12 + i14)==0;
8 r, u0 A2 v0 |( G8 ceq2 = P2ref/v2 - (-i12 + i24)==0;
# m- {$ I+ Q: z# o1 h$ z  heq3 = P3ref/v3 - (i13 - i34)==0;
eq4 = P4ref/v4 - (i34 + i14 + i24)==0;
5 p: U' q! B" _; {5 F6 u, Weq5 = (v1 - v3)/R13 - i13==0;
eq6 = (v1 - v4)/R14 - i14==0;
5 _: ~4 A/ h/ _5 U* }# g- |' g8 n) heq7 = (v1 - v2)/R12 - i12==0;
- V8 T0 l% K. p7 ^" D: Teq8 = (v2 - v4)/R24 - i24==0;
( c6 l6 {4 e! q' ]4 }eq9 = (v3 - v4)/R34 - i34==0;
eqns=[eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9]
8 @" v1 b3 M2 j3 ]  m* Z) Cvars=[v1,v2,v4,P3ref,i12,i13,i14,i24,i34]
" _9 `' G3 z2 l9 J  a
+ n3 C* a( w! _" t2 e3 Rs = vpasolve(eqns,vars)