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

xiangpeng395

1. c=[1 1 0 0 0 0 0;
   4 E+ h8 s2 r0 E5 @
2. -1 0 -1 1 0 0 0;+ P+ i' ]' ~% I( Q! a4 p4 a8 L
   
3. 0 -1 1 0 1 1 0;
   : f3 Q. e' d, U3 R
4. 0 0 0 -1 -1 0 -1;
   ; z8 l. M- O6 s/ n
5. 0 0 0 0 0 -1 1 ]
   / l& M, m; n! }( p% F' k# w& q
6. for i=1:5
   ' p+ w1 r, r! M5 B, c
7. a=c;
   ; l% x- J$ p, |) l! A/ `
8. b=a(i,:)
   + M9 S; s: \( a3 k
9.      for j=1:4
   . B! p) G* O& Z4 w& H# C* W
10.         b(i,j);% y3 u' F1 v+ l) z' g) s; o
    
11.             if b(i,j)~=0
    $ m3 c  C% T\" @% V& v9 G
12.                a(:,j)=[]
    0 B, i( u  d# Z\" c
13.             end( e* u4 @1 {& R/ P
    
14.      end4 c  q( i4 o. [- t$ g
    
15. a(i,:)=[]' a+ v% d& J\" e/ B- q7 U5 u6 k
    
16. end

复制代码

为什么总是超出维度啊，
我想求出依次删除某一行以及与该行中非零元素的列，然后求det(a*a'),并输出结果，怎么搞？求教大神们！
- ~0 m+ N& h0 h2 B* f' |/ ]1 K- L
/ q, Y6 K! m$ _0 i

xiangpeng395

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

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

xiangpeng395

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

, e- A7 _- Q6 j, y6 v好的，谢谢了，

madio

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

1. c=[1 1 0 0 0 0 0;4 X9 `7 ?/ }1 F# w1 I0 ]: q5 q\" v
   
2. -1 0 -1 1 0 0 0;# o/ Y( T0 c) q& H$ ]
   
3. 0 -1 1 0 1 1 0;9 P4 |$ H( B- e
   
4. 0 0 0 -1 -1 0 -1;
   : k* i, T6 F& j3 C8 x* b
5. 0 0 0 0 0 -1 1 ]5 p0 D+ N3 s1 ^6 F9 O) `
   
6. for i=1:5) Z0 V' c& {/ ]- V: Q& P& R6 d3 s
   
7. a=c;+ o  J) o$ t% [$ K& {
   
8. b=[];/ h2 E8 [7 y2 M  n
   
9. del_col=[];9 c( W& T9 o, u( d
   
10. b(i,:)=a(i,:)
    ! {  y* ]/ j* q- T6 t
11.      for j=1:7/ c8 M! x3 Y( u
    
12.         b(i,j);! ?* @3 F- E, N7 o! w6 h
    
13.          if b(i,j)~=0- R6 @. J9 f9 a# u1 V; S
    
14.            del_col=[del_col,j];, Q4 o: W9 c( T5 @6 N0 W
    
15.          end% g% g% Y$ b) c6 [& |4 p
    
16.      end
    0 t/ U\" g+ s. O5 }  c$ g
17.      a(:,del_col)=[]' {: y9 z) n: k) P
    
18. end

复制代码

xiangpeng395

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

; z: N% s4 U6 D& @; N不好意思，我给j的赋值搞错了，j=1:7,我运行了下，还是维度有错啊，怎么办

madio

1. c=[1 1 0 0 0 0 0;\" P% D6 {8 ?$ |; H
   
2. -1 0 -1 1 0 0 0;
   $ K8 @; Z% X$ G8 O) o
3. 0 -1 1 0 1 1 0;- x$ u  C- W! g0 {5 l; \
   
4. 0 0 0 -1 -1 0 -1;
   & v( W\" q% v4 Z# ^% _( I
5. 0 0 0 0 0 -1 1 ]
   3 @0 c$ s1 Y% _+ [4 m: y/ r
6. for i=1:54 _$ T& V) t3 p# K
   
7. a=c;- R3 @: T# p2 k* G0 G
   
8. b(i,:)=a(i,:)! L! d  |5 \$ s+ g1 y& J\" R
   
9.      for j=1:49 }  Z( p\" P. f/ p
   
10.         b(i,j);7 }& x# x5 d3 g4 g
    
11.             if b(i,j)~=09 g) V$ ^/ G. x' v\" O2 j
    
12.                a(:,j)=[]
    ' H# y: P$ I8 ]* F2 f1 {# y, C
13.             end
    ) L, `  {! N: r\" A
14.      end
    # _3 T3 z) a/ x& C\" w
15. a(i,:)=[]- c( q% P1 g1 Q! {
    
16. end

复制代码

改成这样的就可以运行了