每个节点的生成树数目程序

xiangpeng395

1. c=[1 1 0 0 0 0 0;- i( h: }$ V: [  d) G
   
2. -1 0 -1 1 0 0 0;  Y8 y- y* u2 T4 ~/ A
   
3. 0 -1 1 0 1 1 0;
   3 @/ o! R5 p, x
4. 0 0 0 -1 -1 0 -1;
   ) i& H- m$ o! t/ X$ ]0 i$ \2 r
5. 0 0 0 0 0 -1 1 ]
   0 \! y( |$ g2 h& o3 K- r$ v1 d
6. for i=1:5
   * G' {' D/ l\" m$ f7 p$ z9 d2 Y
7. a=c;
   2 D% s( P8 D  M( r: }
8. b=a(i,:)8 E3 J+ E4 N2 D1 y* F% G
   
9.      for j=1:46 \& l! [$ ]) P- C  m- ?; G) X
   
10.         b(i,j);\" I; a1 p- h1 r
    
11.             if b(i,j)~=0
    ! g\" J# S  j1 U8 {9 d( J
12.                a(:,j)=[]
    / U! l4 q, T* F- R6 O
13.             end' B. l2 w  @\" M* J/ A9 i- ^
    
14.      end2 X& D/ F; k* Z9 U0 h# P
    
15. a(i,:)=[]* e+ j  x( A. e
    
16. end

复制代码

为什么总是超出维度啊，
我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
5 v5 T) b2 y9 |( M; p" C
+ O8 J& T, d$ H; X& i% W

madio

1. c=[1 1 0 0 0 0 0;3 `$ P3 M$ s# v\" ~# ~
   
2. -1 0 -1 1 0 0 0;
   . V8 d, ~- ?7 L
3. 0 -1 1 0 1 1 0;  P\" a$ ]! n0 g3 p
   
4. 0 0 0 -1 -1 0 -1;
   * G7 u  M8 T  Q4 G4 l
5. 0 0 0 0 0 -1 1 ]
   2 {$ k; x3 ]3 @( h$ K. C
6. for i=1:5
   3 o; B7 z% O) {) G- z7 Q
7. a=c;
     r& l8 A# ~& h
8. b(i,:)=a(i,:)/ O: S( g8 j2 w6 u; n
   
9.      for j=1:4; V  \$ C1 l* k( `  P
   
10.         b(i,j);\" Q; i3 n, \4 A* C
    
11.             if b(i,j)~=09 s* ]: X- `9 i$ W
    
12.                a(:,j)=[]  N7 K3 N+ x: e* P1 @' @3 H! k
    
13.             end; O+ U. J6 M+ _* e# U3 A& v
    
14.      end& R% t, Z+ M$ J: g) C: G+ c
    
15. a(i,:)=[]% A0 v) I. g* i! b
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
3 r% |0 _- Q" G0 q2 f" Q改成这样的就可以运行了

不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

1. c=[1 1 0 0 0 0 0;  V; x, r, D! c
   
2. -1 0 -1 1 0 0 0;
   & B4 K\" i5 V2 u; O  A$ U; x+ q& E
3. 0 -1 1 0 1 1 0;. G3 f' B) C! |4 v3 S- B
   
4. 0 0 0 -1 -1 0 -1;. m1 s2 T. k$ H) w9 q: h) S$ z
   
5. 0 0 0 0 0 -1 1 ]/ D0 I5 j1 f/ \5 J- E2 Q
   
6. for i=1:5& ]5 G2 l\" R$ x  s2 ^
   
7. a=c;
   4 O/ I' }+ r8 t# i
8. b=[];
   + i- s0 Y  e6 i; [8 l4 I
9. del_col=[];- C) n2 i6 v! `; x  Q\" i% z. G
   
10. b(i,:)=a(i,:)4 S, W2 n; S* |. A9 V  \/ P. l
    
11.      for j=1:7
    4 y' a  k. S8 f4 i7 B+ T# Z
12.         b(i,j);8 Y% y# Y4 R: Z- |% z
    
13.          if b(i,j)~=0
    + J  ^( ~0 r! Z' Y
14.            del_col=[del_col,j];
    * X4 S2 g: n0 [* P, q7 x0 d. l$ L0 D/ w
15.          end& p& j* d% A8 h5 C4 B& o9 ~8 L
    
16.      end
    - x/ @4 R5 U9 O1 I, T7 F; ~- O
17.      a(:,del_col)=[]% \: G+ J$ S) ~
    
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
4 {. J. `: ^5 M' }6 n' F) X! _你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

  V: s( e4 e" ?" D& z3 L可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！