不明白也,能不能详细说一说啊
9 n4 t+ R0 O! d% P这是我的代码:
: T3 l/ I+ I* F* ~, a* G! W错在哪里????
) \+ J$ c/ W: l. u) G; Ovoid main()
& n, K: H3 B0 X: H{n=4;& Y7 q/ z# Y8 x- k3 U5 O0 M
cout<<"请输入基向量 bb["<8 M. a, B; e( [8 ?4 S% `; p$ a
for(i=1;i<=n-1;i++)2 i1 F) h8 _5 T5 }
{for(j=1;j<=n-1;j++) & g, m, o1 m2 s" d, N( G; n
cin>>bb[j];
; z0 I) x7 R0 A# i5 n0 m2 ]$ kcout<<'\n';} N" M2 i+ K- L9 Z- G
void computertb() //求逆矩阵tb[]- G8 y$ I5 Z) `8 g# |+ u
{ for(i=1;i<=n-1;i++) & K/ }4 L) W0 H# P; l, L% {! l
{for(j=n;j<=2*n-2;j++)
6 v& i6 e$ a- B$ I6 E. y{if (i==j-n+1) bb[j]=1;* J& \2 z7 Y; |/ \# M( h
else bb[j]=0;- U: R* t0 U' O4 C% d" A
}}
3 Y0 s4 f. O; _ `% w* \for(i=1;i<=n-1;i++)
8 }& E1 s/ e& t" N0 `. n3 d' R{if(bb!=1)
/ J% P/ V: p* h U1 r1 p$ e{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
# V, [, j- V: j% T( @* W% g! u1 Mbb=1;}% L" s& O" y, x, I0 ?# U: h
for(k=i+1;k" o0 ]# `4 e. }( [9 O8 F0 q
{if(bb[k]!=0)4 P9 u2 N+ U5 N0 `
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
) ^1 Q* W. g/ \+ ?: t. |7 hbb[k]=0;0 R6 G( m8 }: w. ~7 T
}}}+ Z$ c% m% ^8 F$ Y( l0 n8 K* f
for(i=n-1;i>=1;i--)
$ ]' ]# j. z E j! Q" G# |9 b{for(k=i-1;k>=1;k--) 8 l6 r5 U6 m3 y' ~% b' V1 W
if(bb[k]!=0)
" _$ c9 [, b* e; H& G, L{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];0 n6 H0 r2 F/ G
bb[k]=0;
" z* Y$ }3 j; @5 m, S+ V4 A}8 [- c s7 E$ O) Q) r
}, M) E' d% ~6 g" c ]9 u# H
for(i=1;i6 t# ^$ b$ v5 b, E
for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];
% g& Z" a9 F6 M P1 Rfor(i=1;i<=n-1;i++)
' a9 Y0 ]' J6 J. I e{for(j=1;j<= n-1;j++)
% \( Z: G g) l8 y9 H0 i# {cout<[j];1 n0 _; B* q$ n" U0 c9 w- e3 ^
cout<<'\n';* c' g8 G( J2 j
}}. c' D# h& ?- X* h5 d
代码如上:, `! m; b P6 [" g- P5 Q
输入: * d5 A& ?8 k/ V3 w9 {2 R4 n" m
1 2 -1
( d$ |8 `8 q$ K. j! ]2 -2 15 I# j; \5 U1 a7 m
1 -1 2
) F( o& \6 E, }: E! Y运行结果:
/ T& p; F% A5 \( @, `1 0 0 0.33 1 -1.3! h3 `- |9 x. I: w3 K7 |
0 1 0 0.33 -1 1.7 5 z0 d6 }* E7 ?0 Z R9 w2 F
0 0 1 6e-008 -1 29 J9 ?7 b$ d$ G" ?, U7 x! B
& z5 U* g, c4 e' P& q7 x 5 Q' |" R+ H# U3 Z5 R* U. ^, W: [
5 g% @9 g+ U C. H1 d/ { |