不明白也,能不能详细说一说啊1 P+ X1 S8 c) v# C/ M
这是我的代码:
& }: g8 U: [! W; A错在哪里????* ? C! X6 Y8 [) w; N8 M, V- v3 F7 S; M
void main()2 s8 Y( U) A; {$ ]* C }6 ~2 n
{n=4;
, v$ ?( K; V5 m; {$ icout<<"请输入基向量 bb["<: _8 j. |# s+ N2 S* t, ^' V
for(i=1;i<=n-1;i++)
7 J( p6 J* {, u8 Z! m{for(j=1;j<=n-1;j++) ! |* A+ S* a$ b$ y
cin>>bb[j];& z4 Z7 f7 \% k: d0 b
cout<<'\n';}
' D1 d3 j @$ q$ ~! I4 B Gvoid computertb() //求逆矩阵tb[]2 J L' ~8 B* s) @& ]5 C
{ for(i=1;i<=n-1;i++) . L, V: v7 G/ V4 ]
{for(j=n;j<=2*n-2;j++) $ }3 l2 {, B( ~4 v. i
{if (i==j-n+1) bb[j]=1;
3 [( {! y" L$ Melse bb[j]=0;' G2 Y4 S- s. y0 s' B9 S
}}
% S" u6 e) A6 U+ |for(i=1;i<=n-1;i++)
( J9 W" q4 R! @{if(bb!=1)2 \. P* n! T0 o4 r2 |9 j
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;/ f) p7 N5 R7 K; H
bb=1;}' Z' j7 ]" N3 y7 f7 o6 ?( W0 w( C# d$ h
for(k=i+1;k5 w# T, C8 f4 N5 W{if(bb[k]!=0)
; X) C L+ ~. A{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];: x8 o5 a* T; b( Z+ t
bb[k]=0;1 ]% m" \/ m& e( q6 N$ \5 n
}}}( \) o6 V) J0 Q
for(i=n-1;i>=1;i--)
7 u# Q8 i% O* V' e6 U{for(k=i-1;k>=1;k--)
, o! ~2 J u: j0 `if(bb[k]!=0)
2 a4 F3 q1 r# Q1 _# r" J{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
9 j/ d- D, g7 obb[k]=0;
# u' i* U L( K# A* e}
5 P/ ~5 J1 r& D1 n% I% j7 E: g }& s. |# c2 F7 E
for(i=1;i! |6 o f+ ]: W! U5 E
for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];& d6 f. V* G" ?% ^! ?$ u
for(i=1;i<=n-1;i++) y2 L( l0 D6 L: r- w/ k; J: n
{for(j=1;j<= n-1;j++) ! E1 ]6 M& X9 T- s1 X+ L% }2 b
cout<[j];% F% b5 M/ U' J
cout<<'\n';$ A4 D4 k8 I+ q0 ^
}}. M ?% T5 y% _- U! A0 |# b4 c3 q
代码如上:! E8 H- T: }9 H4 X# M
输入:
4 B+ `" {4 i* p" X+ k1 2 -15 l: `7 n# i9 `3 X: x4 E5 P, D
2 -2 1! T7 M* K4 A) h3 J2 J c4 s! d
1 -1 2
2 f& R1 P E1 s. D( l, u1 `* K9 I运行结果:
/ n' X% ~4 A/ |. C9 S8 j1 0 0 0.33 1 -1.3
' k: _/ @$ [$ \) [, G0 1 0 0.33 -1 1.7 9 i' f, ?# m) N8 W
0 0 1 6e-008 -1 2
% E$ [ g' J1 D; y
0 M8 _! L$ {6 j4 k% s' R
" ]- y6 \0 z6 A0 M9 z3 ~3 z- Y9 N7 }, p0 y& }7 \- G
|