不明白也,能不能详细说一说啊) ?4 [& e" H+ a& t1 E
这是我的代码:
1 P# d, Z; C c$ D# @5 K错在哪里????
, V3 G! r3 x! qvoid main()+ k& V! U! r1 U2 k. a
{n=4;" o, f- [6 {, X
cout<<"请输入基向量 bb["<' c+ h' d" H5 m7 ~, L" |
for(i=1;i<=n-1;i++)/ [1 z1 o" w; j5 o4 b
{for(j=1;j<=n-1;j++) ; g O" L6 x" a
cin>>bb[j]; T6 k8 U# x6 z9 Q. m
cout<<'\n';}
4 o2 U6 O; Y1 u; w5 F9 {7 J Cvoid computertb() //求逆矩阵tb[]. o8 s( g9 Y$ H% y& s
{ for(i=1;i<=n-1;i++)
/ s: z9 L8 e0 R/ u3 K4 j! o{for(j=n;j<=2*n-2;j++) , X7 Z* ^" z) }2 \7 @- a. K
{if (i==j-n+1) bb[j]=1;5 U% K) k0 Z/ U0 Z
else bb[j]=0;' v1 F+ Z* F6 l' ]. N' a
}}/ _: w8 E2 |4 Z5 d
for(i=1;i<=n-1;i++)
" g4 A0 ^$ p% r7 i! x" k/ P! K{if(bb!=1)
& n3 j8 r; c6 [5 s{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;! e; H$ k% k1 i; c
bb=1;}: q# f+ O3 B6 c! T, H
for(k=i+1;k$ R$ _+ `6 c3 G{if(bb[k]!=0)
; S ?' ]* f$ g{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
$ I5 Q* [6 A8 lbb[k]=0;6 A: g! v- t# a0 A8 q
}}}/ @, r2 P% O+ o* A4 F, ? _
for(i=n-1;i>=1;i--)% R; w# i5 V' _$ |
{for(k=i-1;k>=1;k--) , p. {0 Z: S6 g# X, x
if(bb[k]!=0) * \' B/ }0 L$ a
{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];; X" t k2 p% h7 ~( v( a" f# ~
bb[k]=0;
1 F- e" Q5 \, H- ]2 Y}- J: a6 r( I1 ^. `6 S" C9 l
}
; f4 P, Q" c2 O' q& y; ?3 Y- |for(i=1;i3 r5 y: }/ z. Ofor(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];5 d& t; B* K+ J: Z' Q A9 g1 E
for(i=1;i<=n-1;i++)$ V( [4 B9 O# o2 y
{for(j=1;j<= n-1;j++)
0 y3 C6 L5 x( q9 a1 w8 ]7 S1 m' Wcout<[j];
3 q8 v' K1 U# B" B7 Lcout<<'\n';
2 h! t' A6 k+ H}}+ l: A d; `" I
代码如上:. m. |6 b; r( x6 F2 Y' f5 I
输入:
) r6 u, [: V8 @+ H3 C S7 X* E" i$ e1 2 -1
8 W2 M" ]9 k- @9 q2 -2 1) D% T2 u' }+ Q7 _
1 -1 2
6 a" u D* _( ~运行结果:
% T8 o1 |; S. p A' c, s1 0 0 0.33 1 -1.3
}) E! P% l% n3 ^2 n a! x0 1 0 0.33 -1 1.7 4 L# T7 k5 J) Y
0 0 1 6e-008 -1 2% K. C- o& p# f1 J e/ R
- p) Z# s) O( j. E5 e+ R# p1 d
: v' ?% [+ ?/ n: L! A4 @
6 l* U7 G" l6 X9 n) i( y |