不明白也,能不能详细说一说啊
8 P& B M9 w2 s6 |7 C这是我的代码:! H* V8 c5 z. f& {3 W! ~
错在哪里????) E2 s6 K1 \ A; ?3 k- u
void main()4 {0 _6 O. Z1 d4 [
{n=4;$ R$ I( l5 t1 o1 ] @) [
cout<<"请输入基向量 bb["<, {2 p; B5 B8 I8 t/ n g; Yfor(i=1;i<=n-1;i++)' k4 ]* A* W- M3 j( @
{for(j=1;j<=n-1;j++) $ d+ Z O# S' Q2 M6 p
cin>>bb[j];
" R) i! d! ^) ~ C) E. L# K3 |cout<<'\n';}
% V, v2 M" u; ivoid computertb() //求逆矩阵tb[]
3 R/ P# r& y2 D{ for(i=1;i<=n-1;i++)
' s/ e {& ], w{for(j=n;j<=2*n-2;j++) - g+ `- Z9 G2 K/ ~
{if (i==j-n+1) bb[j]=1;
R! h2 ^: K# ]) ^; }7 q) telse bb[j]=0;
; r# t7 i5 o3 k4 k& _}}5 f0 C o2 u) n O/ z$ G5 \
for(i=1;i<=n-1;i++)$ _3 ~* n) v: n# H' E; e2 i
{if(bb!=1)3 O |9 \4 a' w0 F0 i0 ~
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
( c: Q2 J& v6 Nbb=1;}5 G4 `/ R# n R% V |
for(k=i+1;k. g! K6 ^: S9 ~- u' y" K9 Q k{if(bb[k]!=0)" ?* ~9 f3 T& C0 q( l
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
4 b8 }, j/ k: r, n( Z- Jbb[k]=0;
9 t d0 J1 H4 Z& g( D}}}4 r4 x( r5 j- i6 \" F6 s% D
for(i=n-1;i>=1;i--). s1 M7 I5 H+ y% N6 {
{for(k=i-1;k>=1;k--)
/ S! \& C3 w1 `, Z# c6 L+ {3 c; uif(bb[k]!=0)
5 n5 q0 G' Y+ K. e; }3 ^' x{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
" T. [5 j% w2 v/ B1 t5 p! Sbb[k]=0;. x0 u# \2 B/ W( _
}
5 b) Z- q, h5 U. U* E3 J }
9 p' j( @! j N# k3 g) M' qfor(i=1;i8 ?5 Y3 \! M/ o. k& ^for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];
+ [7 o5 e' a5 y" J4 D* }for(i=1;i<=n-1;i++)
) z! J! o0 Z/ L; a{for(j=1;j<= n-1;j++)
; z5 F" S' T8 L7 h/ Jcout<[j];
" p7 p. [8 Z5 p% V; ^cout<<'\n';5 z4 I! S9 L) b8 A2 j
}}
; w2 _" J, n2 p代码如上:
# I. {2 V" v$ N4 I输入: , @+ c5 o8 G# b& H
1 2 -1
+ V7 t# e) I* g6 ]2 -2 1
% i; r1 c2 U: J3 a4 _1 -1 25 M$ Y" `6 B1 o$ T" D! N2 F- M
运行结果:
$ d: C" k8 J% R# X R1 0 0 0.33 1 -1.3
5 j* F) m4 S9 G0 1 0 0.33 -1 1.7
: C2 u4 s5 R! W- B0 0 1 6e-008 -1 24 u4 ~- k. l- I' ~; C0 u6 _: e' I) j
4 Z8 l4 ?+ C; c/ h% ?) G6 f , L, W, l7 J7 ? r8 B$ Y
& Q( m$ r0 e* g$ \0 A/ y |