2018年研究生数模竞赛D题

function p=Codidal_Map（）
warning off;
load location; 
i=1;
    fid1=figure（1）; 
    load shg1;
    cx=[99:0.5:122]; cy=[2:0.5:25];
    x1=location（:1）‘; y1=location（:2）‘; z1=shg（:2）‘; 
    cz1=griddata（x1y1z1cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz120‘r‘）; hold on;
   
    x2=location（:1）‘; y2=location（:2）‘; z2=shg（:3）‘; 
    cz2=griddata（x2y2z2cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz220‘b*‘）; hold on;

    m_grid（‘linest‘‘none‘‘box‘‘on‘）; hold on;
    m_proj（‘miller‘‘lon‘[99122]‘lat‘[225]）;
    m_gshhs_f（‘patch‘[0.30.30.3]）; hold on;
    title（[‘第1个分潮同潮图‘]）;
    saveas（fid1[‘F:建模20180915MyProjectsMyProjects结果图‘‘第‘num2str（i）‘个分潮的同潮图‘]‘fig‘）;

i=i+1;
    fid2=figure（2）; 
    load shg2;
%     m_grid（‘linest‘‘none‘‘box‘‘on‘）; hold on;
%     m_proj（‘miller‘‘lon‘[99122]‘lat‘[225]）;
%     m_gshhs_f（‘patch‘[0.30.30.3]）; hold on;
    
    cx=[99:0.5:122]; cy=[2:0.5:25];
    x1=location（:1）‘; y1=location（:2）‘; z1=shg（:2）‘; 
    cz1=griddata（x1y1z1cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz120‘r‘）; hold on;
   
    x2=location（:1）‘; y2=location（:2）‘; z2=shg（:3）‘; 
    cz2=griddata（x2y2z2cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz220‘b*‘）; hold on;
    title（[‘第2个分潮同潮图‘]）;
    saveas（fid2[‘F:建模20180915MyProjectsMyProjects结果图‘‘第‘num2str（i）‘个分潮的同潮图‘]‘fig‘）;

i=i+1;
    fid3=figure（3）; 
    load shg3;
%     m_grid（‘linest‘‘none‘‘box‘‘on‘）; hold on;
%     m_proj（‘miller‘‘lon‘[99122]‘lat‘[225]）;
%     m_gshhs_f（‘patch‘[0.30.30.3]）; hold on;
    
    cx=[99:0.5:122]; cy=[2:0.5:25];
    x1=location（:1）‘; y1=location（:2）‘; z1=shg（:2）‘; 
    cz1=griddata（x1y1z1cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz120‘r‘）; hold on;
   
    x2=location（:1）‘; y2=location（:2）‘; z2=shg（:3）‘; 
    cz2=griddata（x2y2z2cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz220‘b*‘）; hold on;
    title（[‘第3个分潮同潮图‘]）;
    saveas（fid3[‘F:建模20180915MyProjectsMyProjects结果图‘‘第‘num2str（i）‘个分潮的同潮图‘]‘fig‘）;


i=i+1;
    fid4=figure（4）; 
    load shg4;
%     m_grid（‘linest‘‘none‘‘box‘‘on‘）; hold on;
%     m_proj（‘miller‘‘lon‘[99122]‘lat‘[225]）;
%     m_gshhs_f（‘patch‘[0.30.30.3]）; hold on;
    
    cx=[99:0.5:122]; cy=[2:0.5:25];
    x1=location（:1）‘; y1=location（:2）‘; z1=shg（:2）‘; 
    cz1=griddata（x1y1z1cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz120‘r‘）; hold on;
   
    x2=location（:1）‘; y2=location（:2）‘; z2=shg（:3）‘; 
    cz2=griddata（x2y2z2cxcy‘‘cubic‘）;   % 注意一定要有一个是转置
    contour（cxcycz220‘b*‘）; hold on;
    title（[‘第4个分潮同潮图‘]）;
    saveas（fid4[‘F:建模20180915MyProjectsMyProjects结果图‘‘第‘num2str（i）‘个分潮的同潮图‘]‘fig‘）;

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-09-21 10:33  D题
     文件      619089  2018-09-15 08:41  D题2018年中国研究生数学建模竞赛D题.docx
     文件     2725901  2018-09-21 10:17  D题D18102800333.pdf
     目录           0  2018-09-19 10:11  D题MyProjects
     文件        2982  2018-09-19 10:08  D题MyProjectsCodidal_Map.m
     文件     1228218  2018-06-03 16:41  D题MyProjectsdepth.dat
     文件       51612  2018-09-19 06:31  D题MyProjectsENAG_M2.mat
     文件      209920  2018-09-19 06:36  D题MyProjectsENAG_M2.xls
     文件        2042  2018-09-18 21:44  D题MyProjectsEvaluation_East_China.m
     文件        2006  2018-09-18 18:05  D题MyProjectsEvaluation_World_seas.m
     文件         146  2018-09-18 13:05  D题MyProjectsfun1.m
     文件        3379  2018-09-18 15:03  D题MyProjectsHarmonic_Analysis.m
     文件        3683  2018-09-19 07:19  D题MyProjectsInitilization.m
     文件         227  2018-09-16 17:05  D题MyProjectsLeap_Year.m
     文件       43220  2018-09-19 05:44  D题MyProjectslocation.mat
     文件       53789  2018-09-19 05:43  D题MyProjectslocation0.mat
     文件       77977  2018-09-18 13:40  D题MyProjectsshg1.mat
     文件       78926  2018-09-18 14:05  D题MyProjectsshg2.mat
     文件       79154  2018-09-18 13:32  D题MyProjectsshg3.mat
     文件       78140  2018-09-18 14:09  D题MyProjectsshg4.mat
     文件        1842  2018-09-19 06:14  D题MyProjectsTide_Separate.m
     文件    61639767  2018-09-04 18:26  D题MyProjectsTP.dat
     文件    23409877  2018-09-19 05:44  D题MyProjects p.mat
     文件        1031  2018-09-19 09:59  D题MyProjectsTrack_Drawing.m
     文件         620  2018-09-19 10:13  D题MyProjects代码及数据说明.txt
     目录           0  2018-09-19 10:11  D题MyProjects结果图
     目录           0  2018-09-19 10:11  D题MyProjects结果图同潮图
     文件      230087  2018-09-19 10:03  D题MyProjects结果图同潮图K1同潮图.fig
     文件      156782  2018-09-19 10:01  D题MyProjects结果图同潮图M2同潮图.fig
     文件      204224  2018-09-19 10:04  D题MyProjects结果图同潮图O1同潮图.fig
     文件      170219  2018-09-19 10:04  D题MyProjects结果图同潮图S2同潮图.fig
............此处省略28个文件信息