不明白也,能不能详细说一说啊# k( c }2 z; ~! r1 h5 b' \
这是我的代码:
) s9 w- i0 L6 O错在哪里????
* I5 t7 r' j! B6 b7 p. s' Zvoid main()
; T2 e# t% I2 e- P6 M{n=4;' S0 j5 W9 I W, `. c! `2 J
cout<<"请输入基向量 bb["<1 ]0 T* f. d$ w9 O* F. o+ c/ ^1 n
for(i=1;i<=n-1;i++). O: s2 ~- a$ g; L' W0 P
{for(j=1;j<=n-1;j++) , E6 J" p* u4 R* ?$ T! ^
cin>>bb[j];
9 W) a5 v" ^# ?! v" h2 Z. c7 Xcout<<'\n';}
9 w5 Y3 V. ^+ y+ ]0 i( s+ b2 Gvoid computertb() //求逆矩阵tb[]
- A6 N* E6 K W$ g3 N9 W; Y+ o{ for(i=1;i<=n-1;i++)
# @ T9 z. s9 _+ b* H6 R{for(j=n;j<=2*n-2;j++)
4 e/ e( T: }* r/ o7 |3 y{if (i==j-n+1) bb[j]=1;/ m8 D$ l/ u8 j) ?3 s
else bb[j]=0;% P- O; y( v* \7 I8 P
}}
8 c% K3 Z0 P. l! m; R, ~2 }for(i=1;i<=n-1;i++)
' l1 _( A# s$ |! F{if(bb!=1). R h$ F/ @* d( ]. u0 n
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;6 E( {: f1 P- b+ G0 ]
bb=1;}
" |. i* ?( n kfor(k=i+1;k# y# W' T# s3 Z0 [{if(bb[k]!=0)+ | Z3 d: t0 X; l4 J, H1 T/ b
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
, K/ G& b: K! N. K. ~! z* O' Kbb[k]=0;
( X I( H: P6 q# ?}}}
" Q! s$ i6 N4 p Efor(i=n-1;i>=1;i--)
; |- c+ `2 d0 w, ~+ r& }{for(k=i-1;k>=1;k--) 9 w# U8 T: D( o% p D, O
if(bb[k]!=0)
. g; b) K+ S$ f4 i& F( l) \{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];9 ]5 [0 K3 i4 A! I, L8 J. q% d
bb[k]=0;, p8 Q( y% i, {- Y8 V! P
}
7 x) _( T& w& q! L }
+ @$ H4 ]+ v5 X$ Mfor(i=1;i; m1 q! d7 j& g& P4 l" F
for(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];1 I8 c! D: l# I8 T3 y+ h4 s3 \
for(i=1;i<=n-1;i++)
/ C1 A, ?! o8 |: E! t4 {{for(j=1;j<= n-1;j++) / z7 _ @. x! g0 L
cout<[j];1 T& a7 }& n3 }& R" d; \
cout<<'\n';5 A5 D$ ~+ X9 ~- ]
}}
$ u6 p1 @& N8 B1 A/ m代码如上:& U: v; c' C3 W# {
输入: & L) K! @8 G" {5 P4 R: L9 b0 e
1 2 -1) {) s( N, ^/ F a
2 -2 1
: U1 [0 {+ h% T! B8 n8 \ Y* v1 -1 2
9 N6 Z6 ]9 |9 X3 V2 n. \' E. B运行结果:
1 D; p2 K0 ~( F, h& z1 0 0 0.33 1 -1.3
6 @: f" }$ L0 P8 y7 y0 1 0 0.33 -1 1.7
1 g( Q8 N" q) F, T0 |' S0 q0 0 1 6e-008 -1 24 d, @9 o; u+ ?% h9 x9 m
/ ^$ ]* Z7 `6 R8 o4 b' X' k: L 5 d- `1 ?( E! @* K9 O
% h7 n" [8 j( @# P/ L |