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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
     a; `4 T* W$ D6 F+ e
2. -1 0 -1 1 0 0 0;; w( O$ B$ S/ a( C$ w# d  [, ~
   
3. 0 -1 1 0 1 1 0;
   / a: B5 \& _* Q' C( M; }
4. 0 0 0 -1 -1 0 -1;+ b( ]% X  a% Y$ |
   
5. 0 0 0 0 0 -1 1 ]* P: X; U& ]+ e- B
   
6. for i=1:5
   ! s) h4 L: \% p2 t4 m
7. a=c;# L: {. v5 E5 e& `
   
8. b=a(i,:)
   2 W2 x' [+ G1 ]# Q: Y* b
9.      for j=1:4) k- R9 u6 V8 B* D9 f# @8 g
   
10.         b(i,j);0 _) x3 t9 }2 [$ ~+ L. g# ^\" p
    
11.             if b(i,j)~=0, H. \4 I; Y8 s. e' g' }: H
    
12.                a(:,j)=[]4 T. E4 R4 M4 M% l0 e1 R, H0 \  K: `: _
    
13.             end; D: [! [! Y6 ^+ |# q4 M2 ^1 q\" t3 W
    
14.      end
    3 ?8 `$ ]) }7 E8 y# ]1 h2 }
15. a(i,:)=[]
    ! v3 c0 g: ]4 d2 F3 \- J
16. end

复制代码

为什么总是超出维度啊，
我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
2 J& Y7 f+ R% s5 Y3 F
$ G: K( X1 W2 ]) z* B

madio

1. c=[1 1 0 0 0 0 0;/ A( l4 T2 u3 k: R' I
   
2. -1 0 -1 1 0 0 0;
   3 w- u# d9 i2 I4 t; {
3. 0 -1 1 0 1 1 0;
   ' g7 H  n9 Q- [2 e+ v& C
4. 0 0 0 -1 -1 0 -1;: C3 I% H! _  g
   
5. 0 0 0 0 0 -1 1 ]  `. E5 `- B  D7 e& b8 G
   
6. for i=1:5
   * a0 r, m/ g4 ?* v
7. a=c;. ?' C7 [& U. [( ]\" E
   
8. b(i,:)=a(i,:); [% Y% y0 s) W8 ?
   
9.      for j=1:4( j+ K' v5 G( r# T2 u
   
10.         b(i,j);
    # Q( N$ V. }) P: x: s7 c2 p7 b
11.             if b(i,j)~=0
    ' Z* \0 p' N  K
12.                a(:,j)=[]' D, t; c, U  [% Y% {! A0 V
    
13.             end$ _5 _3 |; q& m. C
    
14.      end8 }( c$ v2 z: Q! s
    
15. a(i,:)=[]
    ' R9 _) V\" c/ j' c1 Y$ Y
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
改成这样的就可以运行了

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

madio

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

1. c=[1 1 0 0 0 0 0;
   5 S* @* q' L; _; _0 _2 {& J
2. -1 0 -1 1 0 0 0;7 n$ w7 ]9 m8 b4 ?  X4 v
   
3. 0 -1 1 0 1 1 0;
   % X# B/ }0 U1 [& R$ W
4. 0 0 0 -1 -1 0 -1;
   & l* x/ [+ O5 f  r) E\" B\" ]2 v
5. 0 0 0 0 0 -1 1 ]
   2 r) L4 D, \4 H/ o- A) n( E% K
6. for i=1:5; [; s7 a% q. Y  z# Z- n5 D
   
7. a=c;
   8 n5 J' l6 k0 M, I2 j* B: C
8. b=[];
   1 M5 {9 i. b: Q( L# F
9. del_col=[];4 k8 l9 P8 }: @8 R7 W' [- l  x
   
10. b(i,:)=a(i,:)
    * f3 @1 l/ G7 ^; g
11.      for j=1:7& Z* m  K$ b; T4 I/ O6 }
    
12.         b(i,j);
    3 K7 @; r. Z7 q% P& J
13.          if b(i,j)~=0\" F% B# E) n! t- r. v
    
14.            del_col=[del_col,j];  T* R4 b# S/ z( e8 z
    
15.          end
    1 B$ P2 M/ w/ l4 D
16.      end
    ' E+ N; W$ I1 E  l  ~6 i$ t
17.      a(:,del_col)=[]
    ; t( o+ T* o- R' ~
18. end

复制代码

xiangpeng395

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

好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
- f6 S, E6 I% g: h# A, u你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

1 y. o! F/ C) Z4 `可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！