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

xiangpeng395

1. c=[1 1 0 0 0 0 0;/ Z) D% O4 q) U) B\" a
   
2. -1 0 -1 1 0 0 0;
   ) j9 ?. z+ I- u& X. S& o
3. 0 -1 1 0 1 1 0;+ I) f, O4 ]: @6 ]5 ]
   
4. 0 0 0 -1 -1 0 -1;8 ~: _& Q. f) U: i; P3 E
   
5. 0 0 0 0 0 -1 1 ]
   & H. B) K* P& H' G* @
6. for i=1:5
   ! e) U& R& l, s0 S- }9 ^) y
7. a=c;
   / f$ E) Z\" y\" p7 u
8. b=a(i,:)3 j. w, n8 W6 I7 S6 \' @
   
9.      for j=1:4
   8 G& y7 K: ^\" {
10.         b(i,j);7 `* \8 }' Z7 D/ `; a7 a3 _4 z
    
11.             if b(i,j)~=0
    ' f; d  F- [5 i+ t: t; J# ?
12.                a(:,j)=[]
    \" {$ P* j1 J0 _- C' _7 P  F
13.             end% ^# s( J6 h- m8 H8 q3 F, V7 L0 `
    
14.      end! A& w( [: `) i7 b3 e1 E$ c
    
15. a(i,:)=[]
    + O' e+ H5 [$ w6 m
16. end

复制代码

为什么总是超出维度啊，
  c5 o! I4 e; E6 d% t; v( @$ d; p我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！

; P4 \9 ^$ }6 e" }) {, ~" E

madio

1. c=[1 1 0 0 0 0 0;
   ; o% y4 K$ i. \* Q9 z& r, L: [* k
2. -1 0 -1 1 0 0 0;
   7 b  H! H; ~6 d# [
3. 0 -1 1 0 1 1 0;
   & M. l! v5 E& b9 M
4. 0 0 0 -1 -1 0 -1;
   ! ^- {' K: o3 v\" Y$ h
5. 0 0 0 0 0 -1 1 ]1 z5 W* b1 `+ V) i$ \+ Q
   
6. for i=1:5
   & h  X; R/ W0 v2 }- p' ^
7. a=c;
   ' V) g8 z4 g. B# _0 I4 [
8. b(i,:)=a(i,:)
   , i% r- Q3 `% D) t( i2 B4 u
9.      for j=1:4
   - K: N( a1 }9 C1 _
10.         b(i,j);
    : g6 |& M. V& f2 A, I# ~
11.             if b(i,j)~=06 m% G& I4 u& v. Y% S
    
12.                a(:,j)=[]) N# X\" ?! j. R* i$ {  y! ~
    
13.             end
    7 I; d) l. X2 |. w
14.      end
    7 z4 Y+ W, u7 [5 d3 b' V
15. a(i,:)=[]' Y3 {2 N- D7 |+ d' `
    
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

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

# j; C/ }% W+ z2 v  |& P不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;$ [4 r; Y( |. x% a9 J/ r' r. B' x& E
   
2. -1 0 -1 1 0 0 0;
   ) f\" q$ Y# a: V2 Z4 D: o
3. 0 -1 1 0 1 1 0;3 R1 N, o% }) l2 k) ^3 q/ m
   
4. 0 0 0 -1 -1 0 -1;% G! d1 d7 y\" ?/ D
   
5. 0 0 0 0 0 -1 1 ]
   2 Z' X8 s+ U! U
6. for i=1:5
   ( e# I& {8 f9 o0 Z* p- C- @
7. a=c;' K5 M# a* P* z5 G% g
   
8. b=[];2 w! f$ C2 k\" x0 ?1 f3 O! i9 |: w
   
9. del_col=[];& ?6 |8 Y5 c( B- W
   
10. b(i,:)=a(i,:)* r$ d2 g; q/ ^  s2 L! ^
    
11.      for j=1:7
    5 P& W+ @; p6 z8 O  M7 i
12.         b(i,j);, ^; H& T1 z' l+ }8 C
    
13.          if b(i,j)~=08 Z2 Z9 q! q' _+ U
    
14.            del_col=[del_col,j];: b) o( o\" R0 X+ k$ c# c, `9 |
    
15.          end
    ( D* e- ^/ ~7 j2 C& M
16.      end
    5 V6 V$ q0 h$ M
17.      a(:,del_col)=[]% y9 @( j\" Y4 i0 r2 A; t. T. i- Q
    
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
0 n9 w7 F, c) v* R3 x* B$ K' o你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

好的，谢谢了，

xiangpeng395

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

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