不明白也,能不能详细说一说啊
1 _, `6 ]9 J! Q这是我的代码:3 r3 m! b N+ S) b
错在哪里????
/ `4 H; U! f2 Q! \void main()4 }( {3 y. b) W$ u2 M
{n=4;
2 ?0 g" ?6 e# kcout<<"请输入基向量 bb["<2 d, l2 n% J" K( [ ^! Dfor(i=1;i<=n-1;i++)
8 K% N9 F8 w" Q: o7 o' j3 f% d* b. a{for(j=1;j<=n-1;j++) " u2 C* |4 V' h) p0 P9 S" [1 u
cin>>bb[j];
, A( i# E2 M! Y- I! _# B: F; `cout<<'\n';}
; `0 o) V1 O8 U W" E7 bvoid computertb() //求逆矩阵tb[]
9 Y( {4 E- q) G4 d{ for(i=1;i<=n-1;i++)
. [1 F' |$ \* M! V* W5 ?! G$ }( f{for(j=n;j<=2*n-2;j++) 8 ?7 j+ G) e8 G r E$ x
{if (i==j-n+1) bb[j]=1;
: d& I Y% N! I: Jelse bb[j]=0;
0 _5 L1 x* [( R w& u4 f( I}}
: o9 O7 e0 @5 \4 Dfor(i=1;i<=n-1;i++)$ [/ @: `$ O- j4 B- M7 `
{if(bb!=1)1 S5 t5 u& S' g, g
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;1 N) o* X/ |; e' P5 w2 `
bb=1;}
# r( e, B% J$ E* y4 z( g( t/ Dfor(k=i+1;k3 R7 K- i0 W$ ]3 g{if(bb[k]!=0)% R: \, N# c0 {8 U
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];1 H+ s; X, E) E! g. F
bb[k]=0;. ~1 \; {2 U0 {0 {, q' m
}}}
$ G( m' k/ X* ]# c ]for(i=n-1;i>=1;i--)9 y% t& ]; @# a% e( k# ^: E
{for(k=i-1;k>=1;k--) 6 W+ f) _; e% ]- U5 z1 U
if(bb[k]!=0) % u4 O) i5 _' s! I
{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
9 j$ | C0 `8 @3 `- J, V" Sbb[k]=0;2 M) ~% X% V4 Z2 u
}
/ f" D' g/ w6 E }9 T3 \; a7 ]' U, k; p8 E: u5 _
for(i=1;i. {; J% T- v; Wfor(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];. y: @8 A" |% C0 t4 d
for(i=1;i<=n-1;i++)
2 S& U+ S0 m2 i/ k' F; K{for(j=1;j<= n-1;j++)
& W+ e+ Z$ l- |* q4 h' a: n" Lcout<[j];8 q5 A$ @1 b2 r9 @* f' i" ]! W+ s, ~
cout<<'\n';
, S x5 L/ A+ J/ V}}+ q0 o$ g- \; B
代码如上:4 g( j- L8 U3 ?
输入: ; `& P" ~! v+ f3 O1 ^2 c
1 2 -1
( g% m s. ^8 k1 N2 -2 1 n% [! y8 }' C f3 [) M
1 -1 2. U8 J3 h( U$ c% O. @, J
运行结果:
) `; k: g2 f6 @2 U1 0 0 0.33 1 -1.3
/ T& B4 I* L5 A j0 1 0 0.33 -1 1.7
/ K# y, D" u0 o5 p+ ~, N1 B0 0 1 6e-008 -1 2
2 L" M/ F; V% d4 B7 t* @
- c& v0 Q8 z: X/ C
7 v* X) `! \$ R8 P' j. T! d
3 S4 D+ l- F; `2 K# j5 r |