不明白也,能不能详细说一说啊
* F. a$ n, ^$ X3 }# g这是我的代码:* O$ b) i8 o6 v6 y0 {7 g# e
错在哪里????0 \( |* l% H- [" |, |. V: C* D
void main()
, q4 @; y$ w0 j: t# v& o{n=4;: l8 g* q( g/ i% b* N
cout<<"请输入基向量 bb["<3 I% s2 a5 s& \* {; o$ B* E7 U2 {for(i=1;i<=n-1;i++)
7 ?0 @7 ]6 q3 `( }- B{for(j=1;j<=n-1;j++) 1 r" m3 C7 q) C" ?2 u/ }7 l. ~- t9 S4 k
cin>>bb[j];
6 S5 D: Z, J2 C" |1 o3 j( Lcout<<'\n';}4 a; W) r( _. }5 l
void computertb() //求逆矩阵tb[]; t/ \' A0 D8 H7 m v ^
{ for(i=1;i<=n-1;i++) 6 D: G# ~+ H( ?. [
{for(j=n;j<=2*n-2;j++)
) E2 b* Q c6 ?* r1 H/ ^{if (i==j-n+1) bb[j]=1;8 ]( z: D& U5 l$ c R y' V$ f. {
else bb[j]=0;
/ `; y. e* E# `5 m}}
% q; e* Y, N" ]: S. D( e/ a4 T) ofor(i=1;i<=n-1;i++)& `+ @( \+ E8 x" \6 M
{if(bb!=1)) ]. [7 r" V) B
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
6 S0 \. F# j. r4 R/ b9 K) ubb=1;}
- f4 x) X6 D! lfor(k=i+1;k4 \5 v. B& h# b5 s{if(bb[k]!=0), m7 r P2 n; |, o+ R% h
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];# y! \& Q4 L: w
bb[k]=0;# E& \2 B5 `( i: `, R& R" A
}}}- Z( Q' }2 u$ P
for(i=n-1;i>=1;i--)
8 F5 { _3 M2 l2 g+ A{for(k=i-1;k>=1;k--) 2 s. `/ J; ?% |/ B t" L9 r! L/ [
if(bb[k]!=0) 9 l% i) ?3 K3 x4 T5 w6 s( L' H7 T
{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];: U" n& a" A% a5 ~. a/ Q3 K& R/ T
bb[k]=0;3 e Y7 S. L" a- r
}4 G6 T- f- w+ ?' N
}
! U3 }/ O( m+ e( M! P2 vfor(i=1;i0 q) i, h9 \2 F" M) V1 U u& dfor(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];9 e/ D4 k1 W" I
for(i=1;i<=n-1;i++)
* S; X- u: \, I4 L( D* R' Z{for(j=1;j<= n-1;j++) 1 c; J* f# V/ l7 {, t
cout<[j];
: O; q4 G) w$ F* H1 u' c% Y% F+ K- Qcout<<'\n';: i/ Z' ~/ y+ H0 G
}}" R4 B' O; g' ]( x2 `, z0 e
代码如上:
& c4 J: p# P2 b x, d输入:
2 o$ _6 J6 z% D, ~* s# x1 2 -1
* u s6 M$ n: h( \2 -2 1
, @6 W) N3 |! ?- {3 B8 W N4 p1 -1 2
* g6 Q7 S/ R/ z6 q/ [6 U运行结果:
/ \2 X/ S0 q% J W# G& H1 0 0 0.33 1 -1.3
V; U6 {4 G& r% z0 1 0 0.33 -1 1.7
' d$ y; B: v8 v. Y/ Z% T0 0 1 6e-008 -1 2
- ^- L. k, e( O4 @& ~0 \
" q+ i1 ^2 ]0 Z4 q+ \ 3 l1 [ e- @4 }7 C7 R* m2 t
: b' r! ?! N, J# P, N6 C% E
|