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

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
A(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
& E+ O" l( k9 V% J& s. G( x12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
可以用D(4:end,5:end)=1很明显的看出来：
2 b5 g1 U8 \6 T4 u" p; P' jD =
# a8 T6 z" w* X0 y! A
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
- Y" B2 c) s5 G3 m5 i0 2 0 0 0 0 0 0 0 0 0 0 0 0 0
% ]; t+ O/ c" i3 N; A: ^  w' H1 ^4 f0 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
+ P0 `, z" f7 r: P, t- |. |- c! M0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
& ]" O- r  Y- L/ w3 i7 N0 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
) g6 T6 q2 n; \: B. i6 ~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
1 D2 S- K! ?, a- z, `! u0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
% p  M+ x& Q' W) c6 c- 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
8 k5 S# H2 n" F$ K0 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
而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
D(4:end,=[],D(:,5:end)=[]
我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
7 `# ?, s- i3 O2 r& ~4 T& ID =

$ X& k3 ~+ v8 ]' b1 0 0 0 1 1 1 1 1 1 1 1 1 1 1
& V9 Q6 ?) k) N1 K4 G; ^) v: V# w0 2 0 0 1 1 1 1 1 1 1 1 1 1 1
, W  U( D: S; ]. i0 0 3 0 1 1 1 1 1 1 1 1 1 1 1
# P# o# L& u. G1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
' G& ?8 Q4 G& w  C, L$ w/ y1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
  h, w$ @( d7 ~, Y6 m1 F1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
. y/ m( e0 b/ Q9 r1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
" S0 N; q/ O1 N$ t$ z* H1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
6 P8 A; e( i2 k+ 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
5 m* B# u9 L2 ?1 }1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
( Z! ]+ n# l& z8 X$ `1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
9 V2 S. i: r( q' A1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
# x( G: Z% P5 r# x6 ~; n+ E13、插入行或列，要注意反向依次调用数据的时候，步长为-1
) `7 W0 ~) Z0 z- i# o& O, Va=[1:15]
B=diag(a)
) g2 a; b; b/ e! o' w) U0 R9 z6 ZB(:,7)=[]
1 r5 g! f1 `& o! y/ i: xB(:,end+1:-1:7)=B(:,end:-1:6)
- L2 ~" }9 E& k/ tB(:,7)=0,B(7,7)=7