陈雄达-数学软件介绍-本2015-数模讲座101.pdf

数学软件介绍-MATLAB

同济大学陈雄达

2015年9月5日

提纲

•2013B碎纸片拼接复原

•赛题的基本分析

•MATLAB图形处理

•MATLAB数据及文件处理

•MATLAB的一些组合计算方法

2013B碎纸片拼接复原

•破碎文件的拼接

•司法证物复原

•历史文献修复

•军事情报获取

•人工VS智能拼接

•问题1：中文+纵切

•问题2：英文+纵切

•问题3：中文+纵横切（19X11）

•问题4：英文+纵横切（19X11）

•问题5：英文双面+纵横切（19X11）

2013B碎纸片拼接复原

1.总体思路

分行+行内排序+行间排序

全局生长算法

邻近生长算法

其它思路

读取图片

文件编号001.bmp到209.bmp

fork=119*11,

[A{k},map]=imread([num2str(k,'%03d')'.bmp‘]);

end

每个A{k}大小为180x72

输出图片

假设文件013.BMP最后被定位在第2行第3列，()

kij

fori1:19,

forj1:11,

kINDEX(i,j);

B(180*(i-1)+1:180*i,72*(j-1)+1:72*j)A{k};

end

end;image(B);colormap(map);axisequal;

单个纸片可以提取的信息

1.行距

1.判别是否空行ifall(A{k}(i,:)0),b(k,i)….end

2.两纸片空行对应ifall(b(k1)b(k2)),….end

2.每行灰度平均值

1.mean(A{k}(i,:))

3.边缘

1.A{k}(1,:)A{k}(end,:)A{k}(:,1)A{k}(:,end)

4.边缘走势

1.diff(A{k}(:,[1,2]))

第一步：分行

1.行距信息

1.英文精度差（为什么？）

2.之后算法是否都需要考虑中英文的区别？

2.聚类算法

1.每一行的像素个数

1.sum(A{k}(i,:)~0)

2.聚类

1.余弦定理cR(I,:)’*R(J,:)/norm(R(I,:))/norm(R(J,:))

2.相关系数

3.欧氏距离倒数

第一步：分行

•规划算法

•假设每一组最左边一块可以识别出来(i)，相似度为C(i,j),则

maxCX

•MATLAB命令ijij