不明白也,能不能详细说一说啊7 R) v" z" |: ?+ B9 G
这是我的代码:- C) b2 b: H4 c0 O+ u5 F: d
错在哪里????
, S9 o& l: V0 `void main()
1 C% |& N* ~. g% \{n=4;( A- o8 c' w, z6 Y% `7 o3 n
cout<<"请输入基向量 bb["<- Q# F* H0 ~) d8 F$ T1 {for(i=1;i<=n-1;i++), F- r4 c) e# m0 f; m+ f: {
{for(j=1;j<=n-1;j++)
* v# [! S" c, }+ l8 e6 `9 ecin>>bb[j];
& e9 @+ n+ Y( D0 C) l, Lcout<<'\n';}
$ m m2 K* u, S( @) s; M( N2 b) Pvoid computertb() //求逆矩阵tb[]
' @ o2 ^& L$ U{ for(i=1;i<=n-1;i++) $ N! T' Q5 Y2 q' V6 b# l- h: S, I+ {) \
{for(j=n;j<=2*n-2;j++) : C3 w) K7 U1 W. e
{if (i==j-n+1) bb[j]=1;
, R/ o' y3 ~. K6 _: ielse bb[j]=0;
% l9 Y2 [3 ^: q; |2 F, N}}: x9 M& ?1 K+ q7 _
for(i=1;i<=n-1;i++)2 q, _- f3 P6 \0 e7 r) e
{if(bb!=1)0 j0 u* ]5 v: \
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;! M6 Q8 _# |% I7 w; W9 ~# l. H/ T& q
bb=1;}
7 C7 S! \- I4 X+ ]# X! cfor(k=i+1;k) q' i) M" F* G+ B" A8 j4 W
{if(bb[k]!=0)0 {: R% _1 L" ^/ f( ^9 {- L
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
. k- {# ]4 m8 g8 H. F8 f/ Kbb[k]=0;
, _7 V$ Z, z/ a2 H: a* Q" o7 E5 O}}}
5 f7 m/ U% y c- c5 B/ D4 u' w. `for(i=n-1;i>=1;i--)
M0 U$ J. X! }5 X{for(k=i-1;k>=1;k--) ( d3 A/ e+ n6 k& v0 r* e
if(bb[k]!=0)
, g9 Q7 b6 e- G{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];/ _* x" r2 r3 b4 d0 m) n
bb[k]=0;: y7 P# i. |; D
}
3 u7 r8 Z5 g6 J }
' t5 e8 F' q9 h8 a; g& j3 j2 efor(i=1;i( g0 C9 z1 m5 F! }! @
for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];
. z( U( P' p2 D, vfor(i=1;i<=n-1;i++)
( R R6 _7 {- ?1 u1 X1 e8 [{for(j=1;j<= n-1;j++)
% \9 i9 ?' F1 v) _/ T) O/ Scout<[j];0 ^( ?+ ]3 {; V- L
cout<<'\n';
6 U! N- p6 c L; m0 P/ L! I' Z4 R}}
' C( T% Z4 g0 |4 Z代码如上:
0 X1 d! g( g9 y4 x输入:
4 A; D& Y( v/ c0 d9 `: d1 2 -1+ e) F! d' c, H7 p. ^- E
2 -2 1
& g$ f/ b1 w2 Y2 p: s ^0 _# p. K1 -1 2
4 O! X8 N8 E, S运行结果: Q5 b) h' T1 C2 Y5 F& ]& Z
1 0 0 0.33 1 -1.3
$ p; w" S- Z1 a( F: m0 1 0 0.33 -1 1.7
1 ?! {* K) I3 j' q1 ~' W0 0 1 6e-008 -1 2; k( X% R5 ]- Z- u [, b8 F
5 V/ X+ @4 W: ? l( q" j
7 y u& m1 h' m
$ M7 _( q! m' x1 P& N) f) L |