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

xiangpeng395

1. c=[1 1 0 0 0 0 0;# T8 r, k& c5 y2 A6 v1 j( h
   
2. -1 0 -1 1 0 0 0;& Q. u* j: c; o# y( g
   
3. 0 -1 1 0 1 1 0;
   0 j- X! c\" o, n) D, S/ y  E& [
4. 0 0 0 -1 -1 0 -1;
   & _6 d- c4 V# m2 l6 I# m
5. 0 0 0 0 0 -1 1 ]; ], i\" p9 k8 i, k) L+ u! ?
   
6. for i=1:5. C, |: j- J! E\" V9 M2 y: ?
   
7. a=c;! r- t/ ~7 O. A  F! g, v
   
8. b=a(i,:)
   ( V2 s/ f5 u: k4 o) N
9.      for j=1:40 Y4 E( E5 q9 M$ Y7 {+ F, g
   
10.         b(i,j);  v, `: P( b7 I& N8 B  r7 X# }
    
11.             if b(i,j)~=0
    5 h# w- o6 t' Q4 ^) v2 Y
12.                a(:,j)=[]
    0 N% y4 g; }# W7 O: o- m; S
13.             end! C- f, r0 w( o9 @
    
14.      end
    % [# M3 V/ V9 B8 v! h
15. a(i,:)=[]
    0 G7 x, {7 J4 M. C' I5 t
16. end

复制代码

为什么总是超出维度啊，
7 ?. b; U+ p" @7 X& Y; i% ]! ^我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
3 B" W: O3 W: J! G8 K' u0 `
3 A2 \6 e: b+ l  ]( }9 b

xiangpeng395

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

" k1 k- w+ d- r& e" f' T( s+ [可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！

xiangpeng395

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

好的，谢谢了，

madio

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

1. c=[1 1 0 0 0 0 0;
   ) N! w. p2 r/ c; Z/ O
2. -1 0 -1 1 0 0 0;, o5 Z3 n1 ]' ~. D% M
   
3. 0 -1 1 0 1 1 0;
   / \( {: u  Z' W3 n+ d. G5 R9 F
4. 0 0 0 -1 -1 0 -1;$ M& o- b% Z0 ]0 [+ s) Q
   
5. 0 0 0 0 0 -1 1 ]
   # r; O5 |- |* B
6. for i=1:5. l& R0 u9 W7 z2 I3 y
   
7. a=c;
   1 I5 s% J; r4 n4 G
8. b=[];
   0 j- b- ^0 V# Q/ Y/ P; k
9. del_col=[];
     @' i6 [$ [$ X3 \4 [# \( q
10. b(i,:)=a(i,:)
    ; E8 |% o% k2 D& ~
11.      for j=1:7% x4 {3 t3 U# Y% G1 m2 i
    
12.         b(i,j);
    2 U9 k- ?$ I! ]4 C; H\" x8 i
13.          if b(i,j)~=0/ K4 ?& P6 j) w+ z0 ]( v
    
14.            del_col=[del_col,j];8 `2 }: J! Q# ~# F* C
    
15.          end% J$ z/ O$ H5 c  q\" N- u
    
16.      end
    1 s$ Z/ X. U; }( v2 Y
17.      a(:,del_col)=[]
    , c# a6 R' Q1 a/ g0 F5 V
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 20:26
/ I1 G) r8 ?( ]/ r改成这样的就可以运行了

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

madio

1. c=[1 1 0 0 0 0 0;
     W9 B) }4 a6 b. V' d( K
2. -1 0 -1 1 0 0 0;
   ( u7 u: `- u5 i5 K\" l0 ]
3. 0 -1 1 0 1 1 0;  C1 b( ~) F1 a: i\" N' S
   
4. 0 0 0 -1 -1 0 -1;
     g0 I/ p( c5 }* c$ s+ y: `! _( K
5. 0 0 0 0 0 -1 1 ]. Q\" r+ e% k8 d# x0 r; l3 u
   
6. for i=1:5- ]% P, Z4 ^1 z5 a% Y9 f7 R$ B\" T1 L
   
7. a=c;0 x$ y0 o( h: w, Z' c4 }* s( x
   
8. b(i,:)=a(i,:)
   - l' Z- C+ ^; L3 n  f\" P5 R8 }( K
9.      for j=1:4# q& K( h% n2 L, G9 s) a) A
   
10.         b(i,j);
    - h% @2 u7 f+ h
11.             if b(i,j)~=09 j+ X5 B& L# ?6 o! i
    
12.                a(:,j)=[]
    ( `1 b+ w! O- w8 }% i2 i8 G
13.             end
    / n2 O: s* Q\" y3 H
14.      end
    7 `  P( j! ~- ]1 T6 b6 _6 n. g
15. a(i,:)=[]- _* v! G- a# {3 |+ k
    
16. end

复制代码

改成这样的就可以运行了