不明白也,能不能详细说一说啊7 g1 e9 x6 v' L8 s/ q9 {" D
这是我的代码:3 K: q8 ?8 n; ]$ b J T+ P, L* X1 a4 a% z
错在哪里????
+ Z+ S! D3 |% {void main()
2 d0 ?) U7 V) b& s2 h+ Z+ w1 b{n=4;0 S6 T# |( V5 H2 b3 |
cout<<"请输入基向量 bb["<8 x* v; V) }) f& Q/ r" `. |
for(i=1;i<=n-1;i++)5 I( C# @2 c& L% e* M4 Z" N
{for(j=1;j<=n-1;j++)
7 l2 @9 `8 k+ W- ?7 Jcin>>bb[j];
% z, m [0 g$ N3 M3 e" F: J8 L9 t5 [cout<<'\n';}# u% M. b! ~2 l6 P! j
void computertb() //求逆矩阵tb[]! v* A5 ~' x8 c
{ for(i=1;i<=n-1;i++)
- Z, r" S# `5 z6 u, F* M( K{for(j=n;j<=2*n-2;j++)
, ?8 h! e$ s5 l$ M _5 |& B{if (i==j-n+1) bb[j]=1;
5 X4 X3 Y, b1 Q }7 D& o! Nelse bb[j]=0;* X8 K8 G! q) J" c% Y+ y/ H
}}( P0 }* @& ?6 d6 q7 H: Q
for(i=1;i<=n-1;i++)
. F. Z( j2 Q4 r% T{if(bb!=1)
0 Z" r4 T( ?. b% O# {5 {{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;' V7 f/ O" T: z1 Y
bb=1;}
$ f0 D; ~/ f( R' p6 afor(k=i+1;k7 x: o3 N0 u+ Y) x{if(bb[k]!=0)4 R7 w& M$ a% }. A8 ?
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
9 z+ w8 W! w8 {/ h' t# @8 ^& Vbb[k]=0;
; j' L. `! d% J; D& p. b) g3 J: x}}}9 P3 K4 x* i5 q% A) e
for(i=n-1;i>=1;i--)- W1 r n& M( W
{for(k=i-1;k>=1;k--)
/ D9 }" U/ [( C; C0 mif(bb[k]!=0)
& ?! `4 G {2 B) `+ S{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];+ l: |4 K3 x. }6 U: X" h$ |! f% k, d
bb[k]=0;
; }: G! w- K* \6 Y* e6 T}
9 v9 m' X& U0 K: o- Y6 b6 m8 d" m }
4 u0 t4 v$ h7 I1 qfor(i=1;i+ u3 B% A+ z$ G( O* p
for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];: x) b6 I" _% n! O E3 _0 r
for(i=1;i<=n-1;i++)
) C& u0 {5 T7 ?/ K, ^9 M8 C- L{for(j=1;j<= n-1;j++) 2 `3 |) Q! A6 }
cout<[j];+ V {8 J5 |, U- K
cout<<'\n';
# J0 z/ S$ X5 Q X" V6 `}}2 r2 H, g, f$ B/ ]
代码如上:
# Q# P# L& x5 u- ~0 H0 ?输入:
( w/ Q. U$ [2 t5 X1 2 -1
! Q+ ]3 N3 C5 i+ e! v2 -2 1
, @" Q& b! z# h) `1 -1 2
) k2 _! D/ A! M运行结果:
$ K; M3 e0 x3 P; l$ b8 l9 ]8 @1 0 0 0.33 1 -1.3
. B5 x1 n2 {, T! X8 B0 1 0 0.33 -1 1.7
) ?6 ]+ ?" x' B4 Q0 0 1 6e-008 -1 2
. E. K/ t. R4 p$ X5 v
% Y v. V* o% ~* L8 N( `% t ! N. I) }+ z/ S. y9 d o! j* q
6 j4 W+ n( u+ E2 J) b |