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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   4 ]) e\" s7 h. q9 A3 v
2. -1 0 -1 1 0 0 0;
   1 H7 u: |% W  T$ @8 Z! z\" n( ?\" T. n
3. 0 -1 1 0 1 1 0;$ J\" _1 {6 F, W; X/ F5 h7 }* }
   
4. 0 0 0 -1 -1 0 -1;
   / S0 v3 m# N' @) i2 ^* X) g+ r
5. 0 0 0 0 0 -1 1 ]$ F  W% X9 l, u1 G0 i8 Q
   
6. for i=1:5
   7 a; n* l\" A# ]/ g* Z* v- j
7. a=c;\" m3 E; w\" o! U; x
   
8. b=a(i,:), _( c1 x& o& I. K
   
9.      for j=1:47 Q# n1 X0 V  z% p9 d2 V
   
10.         b(i,j);' |- S# ^- E, g6 ~
    
11.             if b(i,j)~=0& t  L& T2 M; l! A  ~5 n0 |# H
    
12.                a(:,j)=[]
    3 D( z6 m( Z9 c  k9 c/ K5 l  A4 _
13.             end
    # D+ l8 L! t/ }9 r\" j0 R' P
14.      end
    ! N6 f( P& y1 f( K
15. a(i,:)=[]
    - w* I! Y3 L- T: q
16. end

复制代码

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

madio

1. c=[1 1 0 0 0 0 0;! S4 Y- b4 O3 ~- x\" [
   
2. -1 0 -1 1 0 0 0;7 ^! X9 K+ e6 _! O. ~& `' v7 O
   
3. 0 -1 1 0 1 1 0;
   ' R* o, |2 L9 j9 O\" f, M5 Y
4. 0 0 0 -1 -1 0 -1;7 K  i9 @# ]8 |8 M
   
5. 0 0 0 0 0 -1 1 ]
     G1 A: q+ a) Y' D$ _9 `, c
6. for i=1:5
   9 D  v5 E4 H+ k; r
7. a=c;
   / p! y4 z: J$ b2 C
8. b(i,:)=a(i,:)8 p6 f% u2 w) ?\" Z3 S5 D7 y6 d
   
9.      for j=1:4
     q, ^! t5 |! o6 R0 R4 R
10.         b(i,j);, E  B$ J: v3 s9 I7 ]/ ~\" Y1 p
    
11.             if b(i,j)~=0
    6 G$ X3 u: J' t; [; r
12.                a(:,j)=[]2 \; n- {! ~) m  c! X$ u
    
13.             end! p% R$ }- L& ?/ n' R
    
14.      end1 M( C+ Z) z5 r) [& l
    
15. a(i,:)=[]
    $ ]9 H/ W% e. u6 V9 g7 Q
16. end

复制代码

改成这样的就可以运行了

xiangpeng395

madio 发表于 2014-8-30 20:26
3 C- M( N, {/ R1 G* q; C$ ]- {改成这样的就可以运行了

. h6 T+ [9 ~; O) S% T5 t不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

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

1. c=[1 1 0 0 0 0 0;
   3 p' z3 s* D/ L- U! s* D  t
2. -1 0 -1 1 0 0 0;1 `1 X4 r% G) d
   
3. 0 -1 1 0 1 1 0;) V3 ~; v; P- f& F& q6 q; j
   
4. 0 0 0 -1 -1 0 -1;+ e: Z# Y4 I  A( x/ e2 B7 _, ~7 p
   
5. 0 0 0 0 0 -1 1 ]
   & [# x9 u+ [: i' S5 @$ i
6. for i=1:5
   2 T; p  @5 h* o' y
7. a=c;% e# y/ l; g2 t
   
8. b=[];
   . _4 L\" J) F# Z
9. del_col=[];! _\" ]# t- ^& U& I
   
10. b(i,:)=a(i,:)
    3 a8 [4 K+ c, {. s
11.      for j=1:7
    : V/ n+ ]# a* B4 |0 |
12.         b(i,j);) Q. M5 G  h2 K2 [: H  Z) r
    
13.          if b(i,j)~=07 h; D% O/ k% \% |& z) w: P
    
14.            del_col=[del_col,j];
    ' z' @3 f+ Y7 O8 Z+ n  u
15.          end/ m. v- x9 S5 |+ N# g# f
    
16.      end
    $ M; {) A5 L2 p$ a0 X
17.      a(:,del_col)=[]
    . Q  K( q; Y/ h1 y8 u
18. end

复制代码

xiangpeng395

madio 发表于 2014-8-30 21:15
" i$ _+ O& l( U. ~你这个删除代码写的有问题，不能这样实现，我给你改了一下，你看看结果对不对。

2 }6 t) f/ d7 c( M# r好的，谢谢了，

xiangpeng395

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

  e, ]* ?9 P. i* K* ^1 y0 L0 i8 y; z可以了，就是没有删除对应行，我加了一句，现在可以了，十分感谢！！！