|
原帖由 qpalz 于 2008-5-31 23:11 发表
a =
[ 1/r1+1/r3, -1/r1, -1/r3]
[ -1/r1+g, 1/r1+1/r2, 0]
[ -1/r3-g, 0, 1/r3+1/r4]
b =
[ 0, 0, is]
x=a\b';
得
x =
1/(-r4*g*r1+r4+r2*g*r3+ ...
下面的结果是lz想要的吗?
-
- syms r1 r2 r3 r4 g
- IS = sym('IS','real');
- a =[1/r1+1/r3,-1/r1,-1/r3;-1/r1+g,1/r1+1/r2,0;-1/r3-g,0,1/r3+1/r4];
- b =[ 0, 0, IS];
- x=a\b'
复制代码
x =
-IS*r4*(r2+r1)/(r4*g*r1-r4-r2*g*r3-r2-r3-r1)
r4*IS*r2*(-1+g*r1)/(r4*g*r1-r4-r2*g*r3-r2-r3-r1)
-(r2*g*r3+r2+r3+r1)*IS*r4/(r4*g*r1-r4-r2*g*r3-r2-r3-r1) |
|