不明白也,能不能详细说一说啊1 Z5 o" P9 X" Y0 I6 L
这是我的代码:
! o" h9 V% m, L- h错在哪里????
3 K+ c, t+ `9 c: Qvoid main()
: R9 P7 M' V2 o. X( u7 T! W( I, f; C{n=4;
+ ~/ e( L' F! Q, q- Z& ~cout<<"请输入基向量 bb["<2 W# J9 U6 e( D$ Zfor(i=1;i<=n-1;i++)
/ _3 g* ^4 O0 S4 y8 h5 E{for(j=1;j<=n-1;j++) 7 \$ H8 X) y) q A- E6 v# a
cin>>bb[j];" Q0 g" Q0 [, a( F$ p8 D
cout<<'\n';}% v7 R' ~( g9 e
void computertb() //求逆矩阵tb[]
* B1 d9 f6 J: Q' U1 g3 R, {2 D{ for(i=1;i<=n-1;i++)
; F. N$ G6 E6 l" a. [+ p{for(j=n;j<=2*n-2;j++) ( X$ o% [9 s; ^! W0 l5 }8 g
{if (i==j-n+1) bb[j]=1;+ \( G1 `$ ?. G9 K6 X4 Q% L
else bb[j]=0;) v/ W, g' m+ I% R3 D
}} }0 Q( O1 o# P8 e L" ^" B) z
for(i=1;i<=n-1;i++)
# }2 X9 o" v# z# R{if(bb!=1)6 D$ M* ^. I3 u6 t! P
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
% W+ @4 I" N7 M# ?) [4 Fbb=1;}
" Z4 Z$ N6 r% E8 v6 h& kfor(k=i+1;k/ K9 ?# z1 M g{if(bb[k]!=0)
# C& N+ |4 } c$ q& R9 r( B{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
* M1 \0 c$ M, X$ Xbb[k]=0;9 W# {; N$ T" s
}}}9 n& ^% a4 N9 [' a" F( m) R
for(i=n-1;i>=1;i--)
/ m" N7 L2 }) N0 Y{for(k=i-1;k>=1;k--) 4 W. V7 g2 N* M1 s: [
if(bb[k]!=0)
1 F" s2 D s6 C5 C& E: U{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
9 h) f. e2 s, g. ] `bb[k]=0;
1 |: a& b2 z S( ~. n# [$ b}
: w8 S0 r, w" ?' Q) n' n" l }
5 j) F* y9 Z! xfor(i=1;i! W) r! D. q" {for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];+ F# _& M; K5 N' ~0 j) R
for(i=1;i<=n-1;i++)
+ N4 R! Z4 T) v7 x0 p{for(j=1;j<= n-1;j++)
; g2 V* I; K, Q- h5 Y4 icout<[j];. E% m8 ^' t( a& I( C+ Z0 j
cout<<'\n';
* U( D8 V4 E& r: _; F}}7 r5 z7 I% q$ |- K) d
代码如上:
9 h9 ?8 b0 g# g输入:
# r6 q6 J" K4 n' J$ n1 2 -13 p& b6 }% Z* z8 j
2 -2 1
1 [6 h- f: n) g& p) i8 a6 G1 c1 -1 2+ g2 ?, U( J$ m% i$ w7 `8 f
运行结果:
( p% y1 k' H, x y$ Z; F; V1 0 0 0.33 1 -1.3
7 S5 e3 w% T# F# O0 1 0 0.33 -1 1.7 ) R* h4 S* n% b3 ^1 l" E
0 0 1 6e-008 -1 2- D( [' }4 T; C3 U3 p i5 G$ {* ~
" b) _+ f/ _8 B% J0 H9 `, y# y7 D 5 x7 h# E! _4 b# z2 l
7 j0 R7 u+ f7 i4 g1 i* r2 O |