数学建模十类经典算法(6)

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
* t2 ^  E+ p% |- a* {) G  O简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
! K5 I9 s8 Z6 e6 ^* hA(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
可以用D(4:end,5:end)=1很明显的看出来：
8 P6 M: [/ C5 f* rD =
1 G2 D# h4 I$ v: ]
; Y1 @$ m" w! W* e1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
  l# L2 ~; l9 D, V1 R2 W2 u0 2 0 0 0 0 0 0 0 0 0 0 0 0 0
) r. y; d+ d& X6 S7 [# G& p/ _0 0 3 0 0 0 0 0 0 0 0 0 0 0 0
! h  [+ L* N+ j3 f# O0 0 0 4 1 1 1 1 1 1 1 1 1 1 1
+ o9 J7 Q# [4 Z: g8 y# D0 M0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
' @: K( K8 y& F! C/ x) Z3 P0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
$ [9 e2 m7 P/ v2 Z0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
& |4 w( i5 N0 _3 s% u9 K0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
% o; s2 ^' S1 `' r1 F0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
  ^% f1 ]; c7 D0 H+ A* J# g0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
D(4:end,=[],D(:,5:end)=[]
+ v* ?; Y- B' i# ^, p( l% w. I我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
D =
# E4 E4 @, ^& \. ^' `
' r! u' k- ]: `, ?  n" f% |1 0 0 0 1 1 1 1 1 1 1 1 1 1 1
6 X, H, b4 C4 }7 x5 y! d. G- Z/ d0 2 0 0 1 1 1 1 1 1 1 1 1 1 1
4 M1 S- ]: d8 j' j" a# h0 0 3 0 1 1 1 1 1 1 1 1 1 1 1
2 h. h# {8 j+ N7 i# n8 I3 A, k0 f- ^1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
! a! E# s  T6 l+ K- p  x# Z1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
# q& n; D' |, X& F9 h1 |1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
" _2 E) _. J6 W4 l7 f" J5 O& J1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
# _% V, U. C# T  n& Y4 N8 F1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
) T0 M* p+ k3 I0 w. G# k/ F1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
5 T0 i2 B8 o1 q6 s5 ?; [1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
13、插入行或列，要注意反向依次调用数据的时候，步长为-1
$ U8 _' Y( l) C/ Z$ @8 p7 B* La=[1:15]
B=diag(a)
B(:,7)=[]
3 ^# _/ S0 l' S* n' NB(:,end+1:-1:7)=B(:,end:-1:6)
" L& R# G: P3 |: m/ v2 n4 g9 K  x8 }B(:,7)=0,B(7,7)=7
9 r9 d3 t# b0 m5 j( `2 ~$ ?