不明白也,能不能详细说一说啊
" W: y6 f: l2 I这是我的代码:
# y9 ]/ }' R, h错在哪里????
* G0 ^# R1 y" T7 G6 J. H; zvoid main(); `' b7 h6 l) d
{n=4;# i" T8 V" u( j* s
cout<<"请输入基向量 bb["<5 z g: F+ f3 M4 H" G
for(i=1;i<=n-1;i++) g0 Y m$ Q5 E
{for(j=1;j<=n-1;j++) , O5 a- A( m @& U: j' T& _ c8 n
cin>>bb[j];! k8 }" w' A, k7 C4 U0 @
cout<<'\n';}
$ Z/ ^2 K; Z$ n9 m" ]5 `2 Ivoid computertb() //求逆矩阵tb[]
/ u, G$ B3 _1 b; M{ for(i=1;i<=n-1;i++)
7 I/ R; i; |1 p! ?3 B4 H{for(j=n;j<=2*n-2;j++) 3 g' i( B; ~$ e, W
{if (i==j-n+1) bb[j]=1;# h! p$ t; z- Y2 v
else bb[j]=0; W _' p$ k. D9 i4 @# V( Y
}}
$ t& ~" z7 e& Wfor(i=1;i<=n-1;i++)
% x0 e% ?( M& e. v{if(bb!=1)
4 s7 W$ v) ~) u" T, @{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
! m, z8 C. ^0 j6 Z! P+ R/ jbb=1;}
5 n$ `, X) S: X& y1 K: ^for(k=i+1;k \4 w2 q# t x# c
{if(bb[k]!=0)" Y' Q$ o! ], w& f& Z% N! g ]' W
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
. U: P. K# D4 Tbb[k]=0;
( V0 f7 P5 I' E& R y: [# }}}}
) @1 R3 Q7 M! ?. x) }for(i=n-1;i>=1;i--)
* R w2 @. ]) K# E4 [+ F. s/ t{for(k=i-1;k>=1;k--) 8 E/ q4 q0 L3 w$ d* A" u4 a8 Z+ x6 f
if(bb[k]!=0)
; e G7 z7 Z$ S8 _+ D+ m2 Q{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];, r/ M( m" U9 I0 Q
bb[k]=0;
) l0 n* i% A. J! C7 C}
0 C M2 [ R5 s. x" o) s1 q1 D0 j }
, V9 W! U2 j, K- H1 [2 _7 D) e' Tfor(i=1;i( h' {2 Z8 K" e7 r5 D, ]8 d
for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];& B& S) n8 ~! ?, q
for(i=1;i<=n-1;i++)
3 Q6 n1 D3 J$ i+ z' m8 H0 x2 m+ e{for(j=1;j<= n-1;j++) 4 U7 V- b/ N: k
cout<[j];
5 E% H5 f: X/ L# o/ [- q& Hcout<<'\n';" c/ H! J0 L. L I; {8 r
}}1 f3 E1 a% H: Q# W3 d* q
代码如上:
+ g$ \) v, K$ J2 C6 U. ^输入: V2 T2 i2 o) g1 S: \, q
1 2 -1
, l* U, M6 w [2 -2 1
4 B! H) d8 R/ y6 x1 -1 21 F$ Z( ?) V0 M
运行结果: Z ~- E9 [: @/ h
1 0 0 0.33 1 -1.32 e4 U. U+ Q; R% K' } \3 T5 V* O
0 1 0 0.33 -1 1.7 z0 E, a2 l: Z, s7 p& z4 e
0 0 1 6e-008 -1 2
3 y0 E0 V. i5 Z! `
; M: N" L* R0 U+ ]
6 w7 D' ]3 Z3 V4 \$ H- I+ |4 m* @/ |; {0 w' m! F2 G4 t
|