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

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
A(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
1 S+ q5 f  e; ?8 {4 o+ l12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
) t. w7 e1 w/ |5 @# o% W可以用D(4:end,5:end)=1很明显的看出来：
D =

3 c/ J' G* n% [7 O1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
- m2 s4 Z1 O, _& ~8 {5 o8 x0 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
; i' a& E7 W6 Z; Z) K( x7 A) r6 d0 0 0 4 1 1 1 1 1 1 1 1 1 1 1
+ A* _  B# h8 ?2 t9 A/ y8 |+ f! x0 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
# D9 s6 T8 ~( X) @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
9 `9 m5 @) x. h& z! C0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
, d- Z2 m; Q: o% `) f0 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
& d3 S$ _; [1 W+ a8 O而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
+ R! e- C8 t# Y* o% y- d: U' [D(4:end,=[],D(:,5:end)=[]
我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
D =
  P7 @( d' E* J9 a! u5 z; Z* S
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
1 1 1 1 1 1 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 H* N' N5 I! C6 c, v% q1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
. E8 L3 n3 t' i) t7 ]% 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
7 X, G( {) [% C5 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
- k& {+ K( N% i8 L7 h: z" Y2 j13、插入行或列，要注意反向依次调用数据的时候，步长为-1
8 F4 q0 c, z$ U# a! v( pa=[1:15]
B=diag(a)
B(:,7)=[]
B(:,end+1:-1:7)=B(:,end:-1:6)
B(:,7)=0,B(7,7)=7
1 y' u, N5 h/ `# m  X