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

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
A(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
% j0 F3 h: c) s  X( i7 ~2 s12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
! b# W7 P% N1 c2 a可以用D(4:end,5:end)=1很明显的看出来：
) u& e, I' U( E; xD =

1 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
" d3 i% b( ^# J) Z# d6 X0 0 0 4 1 1 1 1 1 1 1 1 1 1 1
: t- m2 r9 r( x" M* G+ I0 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
/ V# d& |8 x, d5 V: x+ ~. a* ^. E0 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
% L/ D) |3 u# y0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
- Y: P) D0 U# @  j2 g$ j0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
8 M! E# B& b; z% N: _' e0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
% E5 F- ?5 j6 K4 C6 Z+ [D(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
: u! M8 U: K/ x/ x& T; N' f! z0 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
3 e  ?2 ~% Q5 G$ X1 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 \  B7 C  \3 e, L4 [8 V  \1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
& T6 e7 f6 t, D5 U# S1 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& S5 X/ z3 @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
. `, k7 s7 }5 {* H- 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
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
: E$ E; X. R& `3 k: ?13、插入行或列，要注意反向依次调用数据的时候，步长为-1
$ f: F, {; j7 h( [9 `( K7 I. Ma=[1:15]
% g. Z& \) L3 B. a6 z6 k! T  oB=diag(a)
B(:,7)=[]
! b( v1 T% m5 t2 }$ z5 Y2 o1 S; YB(:,end+1:-1:7)=B(:,end:-1:6)
B(:,7)=0,B(7,7)=7
9 L: C/ m, u# m2 L, C8 l& y