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

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
A(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
% q9 S0 f  m8 Q" f( T12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
4 R( F' e- Z3 t& s可以用D(4:end,5:end)=1很明显的看出来：
' s. D" c1 C, L7 }% zD =

1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 Z. n# D6 `9 c( [/ o" j; O9 @% W0 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
3 h" K& n7 g( X+ k; i' R( j2 B0 0 0 4 1 1 1 1 1 1 1 1 1 1 1
! L: u( X) o' P% R% k4 `0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
/ ^) u: ]3 h2 ?+ I' p; H- f2 H0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
: e, c- h3 `& T. K& J/ M) ]4 ?( 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
) B( o7 q' u+ u  d5 p7 K  P2 c0 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
3 m. ]" ~, S! 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
# g0 G# ?8 J% ^) i9 U0 0 0 0 1 1 1 1 1 1 1 1 1 1 1
* L, W% z% o1 _: {6 S$ Y6 g0 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/ O% M; Y# @" I* y& a5 ?而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
3 E0 I/ e( {2 c* wD(4:end,=[],D(:,5:end)=[]
我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
! j, n& _3 N# ~* C, m4 \5 _( VD =

1 0 0 0 1 1 1 1 1 1 1 1 1 1 1
4 L/ h; {3 A( l' |# K/ o" q" z0 2 0 0 1 1 1 1 1 1 1 1 1 1 1
; u- f1 x1 Q2 u5 Z! e5 [9 E4 K: }) W0 0 3 0 1 1 1 1 1 1 1 1 1 1 1
3 ]9 \; M6 F* v1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
$ q4 [  K) t2 ]! l$ v5 p( ]1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
) V3 E& N& `. a4 i2 ?- 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
/ |5 o% s, |$ g7 B. e( G1 e& b1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
: o" v# C8 b! L  S5 i1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
13、插入行或列，要注意反向依次调用数据的时候，步长为-1
a=[1:15]
2 O) b# x, G9 ^" K, Y$ e" R1 }B=diag(a)
3 }5 `7 @* l& _! `B(:,7)=[]
4 H, L/ l/ W  fB(:,end+1:-1:7)=B(:,end:-1:6)
, ^5 g# d: L, j" RB(:,7)=0,B(7,7)=7

: E) I( P0 p( S0 k$ c3 b