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

百年孤独

11、如果A(r,c)出现在等号的左边，并且用(r,c)声明的一个或者多个元素并不存在，就根据需要将A扩展，并将扩展位置上的元素置为0，这样就使得A(r,c)均指向已知元素，但是，如果数组A出现在等号右边，那么A(r,c)所指向的所有元素都必须存在，否则就返回一个出错信息。
简而言之，A(r,c)出现在等号左边时，若该值存在，则输出，若该值不存在，则赋值为0；
A(r,c)出现在等号右边时，若该值存在，则赋值，若该值不存在，则报错；
! u: n; |$ m' G4 F8 b- j6 b12、D(4:end,5:end)=[]会报错，因为这个区域不是拐角形状，而是这个矩阵的右下角；
% c/ c! l3 X( `. U/ e' Y& g可以用D(4:end,5:end)=1很明显的看出来：
4 r5 ^7 l9 y4 v! e: ^D =

1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
8 Y+ `) N+ i: t& Z3 \" j0 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
" Y6 M8 X8 R2 U. U+ x  v" O0 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
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
6 x6 T+ @8 ^% n% E7 t: d0 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
$ F0 L, ]- f- ?# e8 d0 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
而如果想要只保留矩阵的左上角的部分，正确的删除命令如下：
- t6 p/ U) |2 ?D(4:end,=[],D(:,5:end)=[]
我们可以用D(4:end,=1,D(:,5:end)=1将删去的部分显而易见的标示出来：
D =

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
* g1 Z' H* q* e6 c! q3 R8 ?1 1 1 1 1 1 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 i$ R7 ?$ A1 U3 o! N# ^1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
& j4 ^7 s1 b7 \1 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
- l8 V; H& l% z" T+ V2 c# w  H, f1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
7 W; Z- b1 @: a- ^1 1 1 1 1 1 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 [+ Q2 Y9 O% s& W' l0 ^; w1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
( A! t) F; ^* O0 F  I( h1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
% i# n7 K$ p2 Y, Z$ X3 E4 \" A, D13、插入行或列，要注意反向依次调用数据的时候，步长为-1
a=[1:15]
B=diag(a)
& j, b5 n) L" a1 L7 P+ o$ }B(:,7)=[]
$ d% `0 P4 `! a6 FB(:,end+1:-1:7)=B(:,end:-1:6)
7 L4 U# `3 u* e6 i+ ?! A% I1 CB(:,7)=0,B(7,7)=7