不明白也,能不能详细说一说啊
8 g! O8 n/ g' v( S9 g这是我的代码:4 u2 r0 p- G; A8 j( x! X
错在哪里????
7 i; M7 @5 o$ H% p. e+ uvoid main()
) k u1 ^/ `$ U+ P5 r; N7 d) d3 c{n=4;
! f {: @7 S7 _5 ucout<<"请输入基向量 bb["<( F; m, c) Q! F' o& K' Q* W; n
for(i=1;i<=n-1;i++)
: N: q/ x9 p! q: @( q{for(j=1;j<=n-1;j++)
% S" o% w' D4 A* Zcin>>bb[j];
! Q& m- ~$ B- _$ k8 U6 m; dcout<<'\n';}5 ]- w3 F: ]" x8 }" x2 g
void computertb() //求逆矩阵tb[]
: w3 T/ O( h5 M+ ^{ for(i=1;i<=n-1;i++)
' I: C# I2 A& A$ L# }{for(j=n;j<=2*n-2;j++) 4 ]; ^) l+ E* n* M2 w" r
{if (i==j-n+1) bb[j]=1;0 T9 d/ b3 t% s5 Z
else bb[j]=0;
- q/ z# c# d8 Z( {' f2 S8 y}}. ~: _7 Y' f$ U: n5 K D" x$ ~
for(i=1;i<=n-1;i++)
& z! I: [; g0 z+ _{if(bb!=1)
/ V8 V p; g) A{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
" C( T, L2 P4 nbb=1;}
& I2 L6 ^2 F1 efor(k=i+1;k$ O) m* _$ z K2 V5 |
{if(bb[k]!=0)
* k' i( U) q0 X9 E" W2 @1 m% z{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];5 L$ `, ~- F/ Z
bb[k]=0;; w6 s& g8 R: P5 \, N4 B
}}}
% S Y- |4 L6 t' B; \for(i=n-1;i>=1;i--)
3 {9 U# |" z& C3 z{for(k=i-1;k>=1;k--)
5 ?) R, t: Z% t' F7 b5 Pif(bb[k]!=0)
[9 p) k* N6 \3 N( a{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];4 Z- _3 C/ z2 a$ n
bb[k]=0;
9 @3 U) A6 F b Q4 ` y- ^2 K# w}
4 `3 n4 p$ P, T# D }
. |# ~* _ }+ F% S6 g' L3 W9 dfor(i=1;i. O; Y+ G) r: ?# j- k1 s: zfor(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];
* G% ?% {: v; l: b! vfor(i=1;i<=n-1;i++)5 W8 P. _3 a* k! E8 f5 A
{for(j=1;j<= n-1;j++) 7 N E7 u- J; u; L9 B
cout<[j];4 A7 {4 V; p$ [- n) h) r9 n
cout<<'\n';! d! L. G9 y* o9 s2 P1 t
}}
v* p, {) f2 v \/ F代码如上:# k- f1 D. Z5 T/ \
输入:
' x b5 V# e& O$ Z& ^; E+ |1 2 -1' O/ \7 k. k, L2 D
2 -2 1( c% R4 B0 Q) M9 D' X4 _
1 -1 2
. s/ _4 Q% }$ l1 O9 Y& s( H1 n运行结果:3 E4 S5 G# u, B) _; h' l& S8 ^
1 0 0 0.33 1 -1.3# }8 H9 a( A6 ?8 O
0 1 0 0.33 -1 1.7 " L; w/ O) }. c/ j
0 0 1 6e-008 -1 2, P4 q& k a$ @- Z4 {2 I* q
. A) g& m' z% k' x3 l1 U/ |) }2 t
- P, ]3 L$ F8 f& ?5 G! K* M# Y
0 @2 F7 G" a% U9 x8 r3 S5 A* C, t
|