matlab图片拼接

clear
clc
a1=imread（‘F:\第6题\附件1\0001.bmp‘‘bmp‘）;
a2=imread（‘F:\第6题\附件1\0002.bmp‘‘bmp‘）;
a3=imread（‘F:\第6题\附件1\0003.bmp‘‘bmp‘）;
a4=imread（‘F:\第6题\附件1\0004.bmp‘‘bmp‘）;
a5=imread（‘F:\第6题\附件1\0005.bmp‘‘bmp‘）;
a6=imread（‘F:\第6题\附件1\0006.bmp‘‘bmp‘）;
a7=imread（‘F:\第6题\附件1\0007.bmp‘‘bmp‘）;
a8=imread（‘F:\第6题\附件1\0008.bmp‘‘bmp‘）;
a9=imread（‘F:\第6题\附件1\0009.bmp‘‘bmp‘）;
a10=imread（‘F:\第6题\附件1\0010.bmp‘‘bmp‘）;
a11=imread（‘F:\第6题\附件1\0011.bmp‘‘bmp‘）;
a12=imread（‘F:\第6题\附件1\0012.bmp‘‘bmp‘）;
a13=imread（‘F:\第6题\附件1\0013.bmp‘‘bmp‘）;
a14=imread（‘F:\第6题\附件1\0014.bmp‘‘bmp‘）;
a15=imread（‘F:\第6题\附件1\0015.bmp‘‘bmp‘）;
[b1c1]=size（a1）;
A1=a1（1:）;B1=a1（b1:）;
[b2c2]=size（a2）;
A2=a2（1:）;B2=a2（b2:）;
[b3c3]=size（a3）;
A3=a3（1:）;B3=a3（b3:）;
[b4c4]=size（a4）;
A4=a4（1:）;B4=a4（b4:）;
[b5c5]=size（a5）;
A5=a5（1:）;B5=a5（b5:）;
[b6c6]=size（a6）;
A6=a6（1:）;B6=a6（b6:）;
[b7c7]=size（a7）;
A7=a7（1:）;B7=a7（b7:）;
[b8c8]=size（a8）;
A8=a8（1:）;B8=a8（b8:）;
[b9c9]=size（a9）;
A9=a9（1:）;B9=a1（b9:）;
[b10c10]=size（a10）;
A10=a10（1:）;B10=a10（b10:）;
[b11c11]=size（a11）;
A11=a11（1:）;B11=a11（b11:）;
[b12c12]=size（a12）;
A12=a12（1:）;B12=a12（b12:）;
[b13c13]=size（a13）;
A13=a13（1:）;B13=a13（b13:）;
[b14c14]=size（a14）;
A14=a14（1:）;B14=a14（b14:）;
[b15c15]=size（a15）;
A15=a15（1:）;B15=a15（b15:）;
A=[A1;A2;A3;A4;A5;A6;A7;A8;A9;A10;A11;A12;A13;A14;A15];
B=[B1;B2;B3;B4;B5;B6;B7;B8;B9;B10;B11;B12;B13;B14;B15];
[g1h1]=size（A）;
for i=1:g1
    for j=1:g1
        s1（ij）=abs（sum（B（i:））-sum（A（j:）））;
    end
end
s1;
ss1=s1;
for i=1:g1
    ss1（ii）=inf;
end
ss1;
[m1n1]=find（ss1==min（min（ss1）））;
d=[m1‘;n1‘]‘;
[MN]=size（d）;
d1=d（1:）
for i=2:M
    if d（12）==d（i1）
        d2=d（i:）;
    end
end
d2
dd=[d1 d2（2）]
A（d1（2）:）=-inf;
A（d2（2）:）=-inf;
B（d1（1）:）=inf;
B（d2（1）:）=inf;
[g2h2]=size（A）;
for i=1:g2
    for j=1:g2
        s2（ij）=abs（sum（B（i:））-sum（A（j:）））;
    end
end
s2;
ss2=s2;
for i=1:g2
    ss2（ii）=inf;
end
ss2;
[m2n2]=find（ss2==min（min（ss2）））;
d3=[m2 n2]
A（d3（2）:）=-inf;
B（d3（1）:）=inf;
[g3h3]=size（A）;
for i=1:g3
    for j=1:g3
        s3（ij）=abs（sum（B（i:））-sum（A（j:）））;
    end
end
s3;
ss3=s3;
for i=1:g3
    ss3（ii）=inf;
end
ss3;
[m3n3]=find（ss3==min（min（ss3）））;
d4=[m3 n3]
A（d4（2）:）=-inf;
B（d4（1）:）=inf;
[g4h4]=size（A）;
for i=1:g4
    for j=1:g4
        s4（ij）=abs（sum（B（i:））-sum（A（j:）））;
    end
end
s4;
ss4=s4;
for i=1:g4
    ss4（ii）=inf;
end
ss4;
[m4n4]=find（ss4==min（min（ss4）））;
d5=[m4 n4]
A（d5（2）:）=-inf;
B（d5（1）:）=inf;
[g5h5]=size（A）;
for i=1:g5
    for j=1:g5
        s5（ij）=