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

xiangpeng395

1. c=[1 1 0 0 0 0 0;% U3 x# I; y$ B/ D9 l1 Q
   
2. -1 0 -1 1 0 0 0;; g6 Y9 ^% R9 ]1 x5 O
   
3. 0 -1 1 0 1 1 0;
   & Y. R* {, ?, ?; f5 }
4. 0 0 0 -1 -1 0 -1;
   2 r# F- r9 O, |$ D% U
5. 0 0 0 0 0 -1 1 ]\" F% F/ M1 g. u5 B
   
6. for i=1:5
   , J1 b1 E2 e! q( A
7. a=c;
   ; a) p5 ^3 |/ K1 g# H* B
8. b=a(i,:)
   + I$ [/ p/ j( z- Q! u- n+ X
9.      for j=1:4
   . f; Q6 @) E7 J/ Z3 @
10.         b(i,j);
    1 a+ l: M) x) M3 l
11.             if b(i,j)~=0
    5 D  z9 v2 F  {0 B  r5 ^0 b
12.                a(:,j)=[]
    + }! Q$ x- ?2 S* x! k
13.             end
    ; \+ e# J6 R6 M; g1 V
14.      end& A, e( T! l3 V9 [8 ^
    
15. a(i,:)=[]
    $ P8 I9 _$ S$ ]6 b3 M
16. end

复制代码

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

2 U' C9 m1 s$ \& |1 n% p- E# o* O

madio

1. c=[1 1 0 0 0 0 0;% N7 l9 c8 ^3 X9 P2 @# i- T
   
2. -1 0 -1 1 0 0 0;# X: T& Q3 _: w- ~
   
3. 0 -1 1 0 1 1 0;
   \" y( |; F\" V$ L
4. 0 0 0 -1 -1 0 -1;
   / l9 ^! N& f) n/ W* M6 G
5. 0 0 0 0 0 -1 1 ]9 p. v; c% f3 N
   
6. for i=1:51 o, D' T) M  ?/ i+ n
   
7. a=c;; ?4 p3 O' }. q. H; k* A
   
8. b(i,:)=a(i,:); c1 p1 Z9 I0 o
   
9.      for j=1:4
     X/ I. v, p- i# w, d6 [
10.         b(i,j);
    ! A* V0 m0 H% n( \
11.             if b(i,j)~=0
    8 g4 }$ ]# d* i7 R- G3 A: l) D, z
12.                a(:,j)=[]
    : c+ ?7 Q0 I0 D3 {' c- J8 `
13.             end7 O, U, a+ O& }8 G8 W
    
14.      end
    ' S0 y: j! C0 I; e) E7 s
15. a(i,:)=[]
    & M4 A( l+ ^3 S: s5 H) U0 E
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
4 w; k7 I8 _4 S6 q改成这样的就可以运行了

" J8 I# q7 w0 `/ |, l4 e$ Z/ z& P不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;
   $ e2 G4 i2 R* C' L/ P1 m
2. -1 0 -1 1 0 0 0;
     E5 \+ o' x! H. q) b
3. 0 -1 1 0 1 1 0;
   # l! I5 ^7 W! ^2 w5 \+ }
4. 0 0 0 -1 -1 0 -1;
     \& O: G; y3 h0 [, ~
5. 0 0 0 0 0 -1 1 ]$ J$ m, [3 |# L5 y
   
6. for i=1:5. T, T& T, g9 c! t; f* n
   
7. a=c;/ A! [6 F( V( ]- T
   
8. b=[];
   % R- J3 F9 T5 ^/ @2 n
9. del_col=[];( j* o% u% B# `. _
   
10. b(i,:)=a(i,:)
    , [2 `( n6 v2 {  Y
11.      for j=1:7
    / L3 C4 |% ^\" H: `7 n* [
12.         b(i,j);
    : O, M# p0 H. X  Y
13.          if b(i,j)~=0
    ; N; q, f/ @- C; V
14.            del_col=[del_col,j];
    - ~' m# ?$ t8 I  r; Z+ v; N
15.          end) U5 T4 `2 \6 W
    
16.      end
    ( w' m6 \* G; f
17.      a(:,del_col)=[]
    6 K* a* h; P) k1 F9 X' [* j
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
1 v0 s# M% Q5 _% q8 R, m2 L你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

$ S  o' `2 f% l好的，谢谢了，

xiangpeng395

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

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