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

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
A(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
$ ~  Z( T4 @: V' |5 b12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
- T0 B! x; e- T# f  ]1 t1 d可以用D(4:end,5:end)=1很明显的看出来：
D =

0 o8 a4 W( q; X, {, V1 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
' P+ s4 C) E1 ^7 V0 0 3 0 0 0 0 0 0 0 0 0 0 0 0
; z+ X( v/ h  I! r0 0 0 4 1 1 1 1 1 1 1 1 1 1 1
9 @( q: t3 w6 f3 A/ A" U0 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' V' b+ m9 g1 \, b3 q9 t0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
8 E$ I6 }( }# G$ z1 y, R: B0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
2 G; N$ `. p9 k# y0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
$ D6 r9 L8 i6 t0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
$ Z% I: \7 n. S7 \) |0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
& n) u/ {& G7 I  W% j+ F0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
- S5 X' U/ ~9 g+ }% C2 h* E! b: _& M' w0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
3 q  M8 B7 J* R# w1 _' o0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
) D( Y- E* S, k$ U( G; h# O而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
D(4:end,=[],D(:,5:end)=[]
; X7 }) V* t7 i2 R. t8 \我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
+ a0 [9 X2 x" P# ND =
+ H) @/ A) m7 {7 @+ P0 Q, K
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
$ ~. {- w* b# f- U1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
) ~; `& R4 q% Y6 ~+ D2 |6 q1 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
+ \' V1 Y* x& b3 U6 N8 H0 \1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
/ S, c2 _, v/ [2 ]1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
7 b7 E- v$ ]- ^, p4 U1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
6 W' d  _# T* ~9 i' x  r3 K: \8 ^1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
9 s/ U* X0 O; F* b! ]# a1 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
13、插入行或列，要注意反向依次调用数据的时候，步长为-1
a=[1:15]
& V& j2 y% Y- r8 B# p9 IB=diag(a)
) N. ]& N3 K6 ]5 j, eB(:,7)=[]
6 p( U: [$ L0 Z) R0 n* nB(:,end+1:-1:7)=B(:,end:-1:6)
. t1 p4 V6 M. cB(:,7)=0,B(7,7)=7

( \3 N9 X2 E% C0 \