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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   ; s% s8 K4 r\" F1 E6 ~( c5 f. P+ A
2. -1 0 -1 1 0 0 0;
   9 Z2 K\" o2 B# B
3. 0 -1 1 0 1 1 0;: s! s\" w) Z/ P4 B! u) b- J
   
4. 0 0 0 -1 -1 0 -1;$ J6 O& T7 Y# m\" q' s9 i+ {; T
   
5. 0 0 0 0 0 -1 1 ]5 ~8 R  h2 {0 ]! o: g; N3 H) E$ ^
   
6. for i=1:5
   ! D# H7 L7 X/ i2 K; N3 t( P
7. a=c;2 p, V! A& l! `7 w
   
8. b=a(i,:), e* \* s1 G* Q* f% ~# ^
   
9.      for j=1:4
   ' \  y\" l5 J; i9 j. i7 M8 z* s
10.         b(i,j);
    : h0 F5 j8 E9 p) E1 j* m- a* O
11.             if b(i,j)~=0$ e. y, |2 J% D. M# k
    
12.                a(:,j)=[]# K3 |: v  J+ p; _8 L
    
13.             end  G* c& w: \) l4 U7 W
    
14.      end\" O) O$ d1 S: D  B
    
15. a(i,:)=[]
    ' V. a0 N$ c) l9 ?2 w
16. end

复制代码

为什么总是超出维度啊，
  G" r. a) `2 \7 H/ A3 z- M$ E1 U" Z8 k; q我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！

1 j8 K& J7 Z9 `* ?4 i3 [* @8 k

madio

1. c=[1 1 0 0 0 0 0;, d& B\" C' a6 F: m9 b. l
   
2. -1 0 -1 1 0 0 0;
   1 x* l  D, Z0 q  c4 k: ^3 d
3. 0 -1 1 0 1 1 0;
   & r7 A  j6 @8 X! i0 u# ~
4. 0 0 0 -1 -1 0 -1;7 j) D% |8 x+ n. ^; \9 e: E
   
5. 0 0 0 0 0 -1 1 ]2 w6 p; q0 r, @$ m/ ]
   
6. for i=1:5; i4 E; \8 f3 A% O) \) t) ]
   
7. a=c;
   . U! e5 Z$ f. u# o/ S
8. b(i,:)=a(i,:)
   : h, c- J$ F3 }5 I: v\" C
9.      for j=1:45 l& Y! m$ |3 T& J6 {6 i
   
10.         b(i,j);
    6 C2 Q% C# {* F3 x6 x
11.             if b(i,j)~=0+ F$ L* C  c  {
    
12.                a(:,j)=[]- S% g$ Y& \: U' I# l7 n5 d& t
    
13.             end  ~! h: ^+ I. b3 l
    
14.      end
    , B2 `6 ^' ]7 B+ o
15. a(i,:)=[]
    5 k) u1 R9 a. _2 P\" i, e! i
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
/ m8 a6 n) e2 q, T) t4 q$ c改成这样的就可以运行了

/ K. L4 l+ B  L9 m  x7 K/ Q1 X不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;
   ' Z1 D5 e9 T5 w4 L8 h
2. -1 0 -1 1 0 0 0;
   8 u- p: j8 U' ~' f2 k# K  O# e! e
3. 0 -1 1 0 1 1 0;
   ( C1 f5 B4 F* O. {% n( g
4. 0 0 0 -1 -1 0 -1;( {0 h. k6 }9 Z: g/ \
   
5. 0 0 0 0 0 -1 1 ]
   9 }6 }* n+ c! z4 _
6. for i=1:5
   - C/ ]. P0 H/ H
7. a=c;/ h. ?7 t; _% R/ a8 w) M# Q( F
   
8. b=[];, i, }4 [; m- H' }. Q; y
   
9. del_col=[];  j1 @( K) j' a\" g) S/ M0 O& ^  c
   
10. b(i,:)=a(i,:)& C: _* {# O1 k- K! {% p0 w
    
11.      for j=1:77 y( A2 V# h: }4 B
    
12.         b(i,j);% w1 O. y4 h/ a' e
    
13.          if b(i,j)~=0# J  V6 u# J/ M. i
    
14.            del_col=[del_col,j];
    1 V5 K% V. B5 `
15.          end
    # ^0 o9 |& H1 K4 D7 d  x) P
16.      end, `% ]7 q* q2 Z& }' A+ i4 v& y! v+ O
    
17.      a(:,del_col)=[]6 f; N9 V) N* P, e+ t% e1 h, X
    
18. end

复制代码

xiangpeng395

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

5 t0 D5 o: a# s3 w* @好的，谢谢了，

xiangpeng395

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

可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！