数学建模十类经典算法(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很明显的看出来：
+ @! x8 f7 N$ j7 R; \6 [, w) V+ r( ^D =

  g" w0 V# ^3 m& f: B1 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
0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
8 n; x+ v2 n9 d9 D! u9 F. _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
/ M. R) D5 k( H0 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
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
而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
7 U+ G/ j7 T' k; F; E2 D+ V- }D(4:end,=[],D(:,5:end)=[]
我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
4 j' `1 D" S* H) ~8 B6 YD =
; K3 A9 ^$ A; m0 c0 k1 \6 j
8 R0 X0 g9 }' K5 d5 K- D0 I1 0 0 0 1 1 1 1 1 1 1 1 1 1 1
% ^+ b, z* j1 K4 d0 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
" v$ K' t$ j3 {+ Q2 k" T1 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
0 F. k# n! e/ \* v  N4 V, i1 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
, y6 a) k0 _+ I4 O1 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
9 h+ w; ]6 |. K# g* N9 T! h  }1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
7 n! v0 k  `) p/ @5 H& B2 r7 Y1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
13、插入行或列，要注意反向依次调用数据的时候，步长为-1
a=[1:15]
7 F* c" x) @( {, v; xB=diag(a)
B(:,7)=[]
1 u, R2 `# f2 |0 p) Q" o1 ~B(:,end+1:-1:7)=B(:,end:-1:6)
B(:,7)=0,B(7,7)=7
) `. l6 d7 X: {! |
/ H. a- s/ T/ U