不明白也,能不能详细说一说啊$ K+ j7 q: Y# [
这是我的代码:
r% a! `2 S; e8 G错在哪里????
% u" U+ Q8 {. k" q* Q* M _1 wvoid main()% ~& l+ ], D1 q% ?( a& f
{n=4;) Y/ v5 i# _$ \! r/ p
cout<<"请输入基向量 bb["<' }- P2 I7 a: X# s$ X# \
for(i=1;i<=n-1;i++)
! [- i' L+ S7 H) e$ `: W: r{for(j=1;j<=n-1;j++) , Y1 S7 f6 R. W
cin>>bb[j];
" d$ B4 d+ P' d* Y, W) Ocout<<'\n';}
; i$ q$ C! @ Z) f7 r; Cvoid computertb() //求逆矩阵tb[]: W* B( i/ _9 r. i9 F1 R
{ for(i=1;i<=n-1;i++) ) a- k0 D" R. O
{for(j=n;j<=2*n-2;j++) * \* E- w# p( j) A4 m
{if (i==j-n+1) bb[j]=1;
j% A8 a0 p* Oelse bb[j]=0;+ s$ p% g9 o, z$ r2 r1 y
}}% |/ i# }% C; T- P @5 G |
for(i=1;i<=n-1;i++)
/ l5 J8 |2 D ?- H# Y{if(bb!=1)2 \* f8 d" \: S8 C- S
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;, r3 H( D3 y0 e+ K6 R
bb=1;}/ \9 z2 ~( N5 G. `! \
for(k=i+1;k. D- e0 T# B) G* W( O% n{if(bb[k]!=0)7 e6 C0 s% o0 f% ]: H- f
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];1 G5 d5 q# x6 @% E a: d; v$ ^
bb[k]=0;5 y1 @, M6 [5 W3 U8 t) C
}}}
) L' t4 O, W. Q4 k) Jfor(i=n-1;i>=1;i--)
/ I: q6 L' u0 M* V. K3 q! f7 G$ y" u) o{for(k=i-1;k>=1;k--)
8 e. z: E5 w1 I/ _/ |0 ?1 oif(bb[k]!=0) 3 q8 u+ u; }, B* g* U/ n S5 i
{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
, Q; e# o# `5 Gbb[k]=0;
0 X6 ^' h8 |) B0 D2 c}
6 D: G6 @9 J. U7 X }
& m2 {( S+ d0 { T: h5 h/ t6 Lfor(i=1;i1 J4 B$ X5 z- W) F4 ?for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];& I% i2 l# ]1 O! R+ v- L9 e
for(i=1;i<=n-1;i++)' v, K" Q( N% O4 e7 g5 e
{for(j=1;j<= n-1;j++) $ |3 l# z( [' {6 R6 T
cout<[j];$ r6 ?( q# w$ o7 h
cout<<'\n';0 l, o! C, ~! g' H1 \. ^/ j p2 n, f
}}
3 `: T6 P; c; I/ p代码如上:( H S4 j/ [3 c b1 J1 }: _* N* I
输入: ; n6 k# E- R+ m# p4 N: w
1 2 -1
- K) S& X% [3 a) n2 -2 1; t$ W+ p" i1 J$ L( h
1 -1 20 }: Q9 U4 a9 N% v& \" p
运行结果:0 |$ @4 N$ |. K Y( t
1 0 0 0.33 1 -1.3
* g) g) C* g; u$ j/ B/ p$ e0 1 0 0.33 -1 1.7
6 E( C: R5 x5 j8 \9 b p5 a) z" Q0 0 1 6e-008 -1 2: m, I7 {+ ]5 F. H4 b
. M5 u8 H, j- N0 t+ x
1 u* U' u2 ^: U& V+ \5 y- U8 `- I9 J' \) \0 G8 K) S; \) i; |
|