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

xiangpeng395

1. c=[1 1 0 0 0 0 0;2 J# n* D6 V4 B! K; |' F
   
2. -1 0 -1 1 0 0 0;' Q: y  Q1 e\" q+ P' I. T
   
3. 0 -1 1 0 1 1 0;
   ; _! o2 n- E! Z5 i\" ?6 s7 ?9 ?3 v
4. 0 0 0 -1 -1 0 -1;+ J\" B4 C: X2 h: ]
   
5. 0 0 0 0 0 -1 1 ]
   8 T: ?1 y8 x; K! O( Q8 N9 ^
6. for i=1:55 R- k7 a2 H\" V5 k/ F4 K4 ~
   
7. a=c;
   / e5 t9 n' w6 v5 L, h% ?, T
8. b=a(i,:)\" @! n' J  O; D( ~8 Z% L
   
9.      for j=1:4
   ) W7 G% I% K  C5 v3 V0 \+ j. k
10.         b(i,j);$ m! p6 J! r0 f. T
    
11.             if b(i,j)~=0! P% x\" v- L7 u# M
    
12.                a(:,j)=[]
    9 F8 a5 T8 X- I) k' B
13.             end) q. S\" o' ?  P5 N- ~. y3 W
    
14.      end
    3 }$ e3 b5 M' [- a, j
15. a(i,:)=[]! k* v+ J- X0 s. t3 a: Q
    
16. end

复制代码

为什么总是超出维度啊，
0 h4 Y. f& y1 y. C' d我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！

madio

1. c=[1 1 0 0 0 0 0;
   ' O1 w2 a; e9 l  K
2. -1 0 -1 1 0 0 0;& \: s' Z. g; l. W
   
3. 0 -1 1 0 1 1 0;  C4 {0 J+ Z# z. h2 U\" b\" e
   
4. 0 0 0 -1 -1 0 -1;
   , D, W5 I  `4 s. f( S% F- H( N8 x
5. 0 0 0 0 0 -1 1 ]
     ?: _; r\" [% e: P4 `
6. for i=1:5! ?6 Z! ^& _* d* A- T+ z- y
   
7. a=c;
   0 F% b6 `* C% ?/ z
8. b(i,:)=a(i,:)& ]. r! I, `7 J( O# Y# \1 g/ O
   
9.      for j=1:4/ y& q8 i\" w\" o# N\" B$ N5 R- A
   
10.         b(i,j);
    1 P' r8 E: Y- S\" j% H
11.             if b(i,j)~=0- a) F% D6 x% M  t! W
    
12.                a(:,j)=[]
    ) M, m+ U9 d5 t& T5 c
13.             end
    % a4 u2 L6 P. |0 C% S! t% f$ V
14.      end
    + i6 g* p+ [2 o\" t1 a( O2 D/ }
15. a(i,:)=[]
      V, m8 P& i% q# T+ p2 `1 Y' v
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
" u- Z- X2 Z& K, N/ }* R* P改成这样的就可以运行了

' B) }. q& r: k8 m1 H9 |不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;
   6 n5 A% ?' e0 A6 `% h
2. -1 0 -1 1 0 0 0;
   1 `% p/ b7 I; C7 x0 Z
3. 0 -1 1 0 1 1 0;2 H; o8 m& u0 \' R+ J, X7 e; }
   
4. 0 0 0 -1 -1 0 -1;
     M5 Q  z\" \* o$ s\" j' k+ }
5. 0 0 0 0 0 -1 1 ]+ o  n% ~1 V; }. p
   
6. for i=1:5# S: Q\" r( Y) K5 s$ i; {
   
7. a=c;
   4 k3 B: G\" L9 R* ~, o
8. b=[];
   # j6 R+ A4 v( @\" X- y9 C- L' X
9. del_col=[];
   ; |$ Z\" w\" r/ ?+ [8 O, Q( p
10. b(i,:)=a(i,:)
    : D5 j- ^+ P8 K! a- g0 M\" q
11.      for j=1:7
    3 @\" Q6 v& z- j  I5 A
12.         b(i,j);
    ! f+ ~( d& m6 ~5 B3 B/ L% J
13.          if b(i,j)~=03 O/ m; D& _; i6 ]  V) v5 M3 h
    
14.            del_col=[del_col,j];
    ; O; m: t, u' y2 _7 R6 X
15.          end, B  d. p7 X1 G2 U. k\" e
    
16.      end
    0 V  S! s& }% H( _
17.      a(:,del_col)=[]# a3 B) g\" b  w: `7 p- b2 N
    
18. end

复制代码

xiangpeng395

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

( y! v, ]7 u$ [& B% Q好的，谢谢了，

xiangpeng395

madio 发表于 2014-8-30 21:15
2 e$ Y# I8 s- v& H# C你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

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