原创力文档- 1、本文档共43页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
* 找出每一行的最左边与最右边的图片编号：(程序名 zhipian208) clear,clc P=[3,12,23,29,50,55,58,66,92,96,119,130,142,144,179,187,189,191,193;……,5,41,90,102,103,109,114,115,118,120,124,141,147,152,155,156,186,195,208]; for i=1:11 H1=P(i,:); ZY=funzuo1you1(H1); Z1(1,i)=ZY(1,1);Y1(1,i)=ZY(1,2); Z1(2,i)=ZY(2,1);Y1(2,i)=ZY(2,2); end Z1,Y1 load gs13B2data01 Pz=[];Py=[]; for k=1:11 Pz=[Pz;X{Z1(2,k)}]; Py=[Py;X{Y1(2,k)}]; end imshow(Pz) figure(2) imshow(Py) * Z1= 5 6 16 6 6 3 9 8 3 2 3 50 62 169 39 72 15 95 126 30 8 90 Y1= 13 4 2 8 5 17 3 10 8 18 11 142 37 19 75 61 177 44 146 60 197 124 * （三）拼接每一行 1.求每一行的最短距离矩阵：(程序名 zhipian209) clear,clc P=[3,12,23,29,50,55,58,66,92,96,119,130,142,144,179,187,189,191,193;……,5,41,90,102,103,109,114,115,118,120,124,141,147,152,155,156,186,195,208]; load gs13B2data02 Z1=[5,6,16,6,6,3,9,8,3,2,3]; Y1=[13,4,2,8,5,17,3,10,8,18,11]; for k=1:11 tt=P(k,:); d=[]; for i=1:19 for j=1:19 a=YR(:,tt(i)); b=YL(:,tt(j)); d(i,j)=dist(a,b); end end * c{k}=zdfloyd(d); for i=1:19 c{k}(i,i)=500; end A=Z1(k); B=Y1(k); c{k}(B,A)=-50;c{k}=c{k} b=int2str(k); aa=[b,.txt]; fid=fopen(aa,w); fprintf(fid,%f\n,c{k}); fclose(fid); end 2. 运行lingo程序(程序名 px15204)11遍，将每行的排序结果 保存于文件bd1.txt, bd2.txt,…, bd11.txt * 3. 根据lingo输出各行排序结果文件，通过matlab编程复原各行： (程序名 zhipian210) clear,clc P=[3,12,23,29,50,55,58,66,92,96,119,130,142,144,179,187,189,191,193;……,5,41,90,102,103,109,114,115,118,120,124,141,147,152,155,156,186,195,208]; Z1=[5,6,16,6,6,3,9,8,3,2,3]; Px=[]; for k=1:11 b=int2str(k); bbd=[bd,b,.txt]; a=importdata(bbd); c=reshape(a,19,19); px=Z1(k); * for j=1:18 i=px(end); t=c(i,:); m=find(t==1,1); px=[px,m]; end Px(k,:)=px; tt=P(k,:); Pp=tt(px); Phao(k,:)=Pp;
文档评论（0）