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

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
A(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
可以用D(4:end,5:end)=1很明显的看出来：
D =
3 N( M& o$ {% i
: q+ k& a7 ^. W# `. q$ ?; U1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 2 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 3 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 4 1 1 1 1 1 1 1 1 1 1 1
8 T5 {' v, q- q0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
8 O# n3 Q+ B% U$ j: g" @" _; y. A0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
+ z$ j# ?) `- V. C/ _& {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
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
$ F, o* i- ^/ {. G5 i! w0 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
% C) R4 f, p" `+ X' y0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
5 Q8 z4 Y. G& d2 a, ?5 o4 gD(4:end,=[],D(:,5:end)=[]
我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
D =

1 0 0 0 1 1 1 1 1 1 1 1 1 1 1
0 2 0 0 1 1 1 1 1 1 1 1 1 1 1
0 0 3 0 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 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
. J. g8 B7 X& W! ^# Y# X1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
) ?; ?6 x, d5 v1 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
: h; v7 V9 C1 G/ L) Z) d- 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
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
4 L! `2 C+ G8 ~1 g. y1 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
13、插入行或列，要注意反向依次调用数据的时候，步长为-1
) }# p. j" J! W# L# Q6 qa=[1:15]
$ r5 j2 P& [9 E" h6 [B=diag(a)
B(:,7)=[]
3 J; r( n% ~0 y% q: aB(:,end+1:-1:7)=B(:,end:-1:6)
B(:,7)=0,B(7,7)=7