不明白也,能不能详细说一说啊5 |! I% d1 F; B2 _. D( V
这是我的代码:
- b/ c6 d* w% F7 Q3 x! y L$ E错在哪里????
. w: U8 K* a" ^) I( Rvoid main(); a$ V# k8 X4 l
{n=4;3 v! N$ i! p0 u; u! B
cout<<"请输入基向量 bb["<; [1 @1 e* b3 B# B& G' }+ ~* Pfor(i=1;i<=n-1;i++)( K+ o3 H: z: U% u$ M# S
{for(j=1;j<=n-1;j++) , l/ x2 f4 ~8 y& h V
cin>>bb[j];/ A; P( q) ?5 @3 V& H5 p9 b* ]
cout<<'\n';}4 S! [+ x6 D8 \: P
void computertb() //求逆矩阵tb[]/ c3 s4 [6 G' Y' I
{ for(i=1;i<=n-1;i++) " D; f6 N! V- ?3 Z" [, t
{for(j=n;j<=2*n-2;j++)
6 k: l; {% O8 F7 d{if (i==j-n+1) bb[j]=1;. t& g; X9 i4 h% U8 a1 e8 w) I
else bb[j]=0;
# ~) D$ o/ t) {' h/ {$ {2 {8 A}}
; R* i- ?9 G$ c6 _for(i=1;i<=n-1;i++)* a* H. k9 m5 J$ e B
{if(bb!=1)
1 d: E( o' J, D) U8 D{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
% f9 `. o/ _% U: C- J" Cbb=1;}
6 i6 [: B h. O( V+ xfor(k=i+1;k X" p: z4 U! R
{if(bb[k]!=0)
4 C4 V. p. u; \0 [{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
; i2 ?- h- K+ j2 sbb[k]=0;6 o' i$ @% l# [* m( e3 a
}}}+ k1 ]6 Q7 s7 E9 M/ ~9 r
for(i=n-1;i>=1;i--)
* i7 f/ [0 u+ E{for(k=i-1;k>=1;k--)
4 M1 B. p4 P1 N) U3 Uif(bb[k]!=0) - t" L8 }. M7 K2 J( S
{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
# y" y, Z4 x7 y1 ^; wbb[k]=0;
1 W! A( ]+ G% } }+ V# {}$ J r ^+ {1 N" Z/ r; X5 R1 |) i
}
6 J$ t# d0 A; c: Afor(i=1;i0 J( o* G' k- ]7 H$ s& W, A5 A
for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];
+ S" _/ R4 ?. Yfor(i=1;i<=n-1;i++)
/ q6 Y# R3 P1 Y1 Y1 g& ^0 V{for(j=1;j<= n-1;j++)
$ ?- l7 x) k/ C, gcout<[j];
8 A7 t$ a7 v7 U; x b+ Ucout<<'\n';" C7 S( G" ?5 v/ A. S! y
}}3 }8 m' `; k- G7 ` m* C( _ T' [3 z% R
代码如上:
! f3 \9 {" ]# g4 o' P* ]输入: 3 T5 ]5 _/ F1 `
1 2 -1, g1 Y+ q3 @$ ]; H5 o* j
2 -2 1
- F& u; Y& r+ Z3 p3 |& B1 -1 2
/ b5 ^$ b! _& R8 m- G运行结果:* v& y+ ]0 K( u. o% ?6 v: }! Q# a' J
1 0 0 0.33 1 -1.30 N0 @$ `+ L( x, o
0 1 0 0.33 -1 1.7 ; b5 S2 Z1 q3 B5 E
0 0 1 6e-008 -1 2
$ h; B2 H) i z) f$ K$ O( @- W% W7 z, a4 Y
* a$ h7 F; }$ S( U4 b! W1 _
% F' h. M7 |& F, i! v; T9 H |