高斯过程回归GPR MATLAB

很好用的高斯过程回归Matlab代码，可直接使用，欢迎大家下载

function [ result ] = ApplyImagingCondition（ f1 f2 isrc ）

load（f1）;
wavefield1 = wavefield;
load（f2）;
wavefield2 = flipdim（wavefield 1）;

result = squeeze（sum（wavefield1 .* wavefield2 1））;

figure（）
imagesc（xzresult‘）;
axis image


% nx = length（wavefield1（1:1））;
% ny = length（wavefield1（11:））;
% 
% image1 = zeros（nx ny）;
% image2 = zeros（nx ny）;
% 
% 
% Ey11 = zeros（nxny）;
% Ey21 = zeros（nxny）;
% iEy = zeros（nxny）;
% 
% for i  = 1:nx
%     for j = 1:ny
% %         iEy（ij） = xcorr（Ey1（:ij）Ey2（end:-1:1ij）0）;
%         iEy（ij） = sum（Ey1（:ij） .* Ey2（end:-1:1ij））;
%         
%     end
% %     disp（i）;
% end
% figure（）
% Ey_fig = iEy/max（max（abs（iEy）））;
% % imagesc（iEy‘）
% xaxis = （1:nx）*dx;
% zaxis = （1:ny）*dz;
% 
% imagesc（xaxiszaxisiEy‘）;
% % imagesc（xaxiszaxisEy_fig‘）;
% xlabel（‘m‘）
% ylabel（‘m‘）
% % caxis（[-0.30.3]）
% axis image;
% colorbar;
% saveas（gcf[‘result_of_corr_‘num2str（isrc）‘.png‘]）
% 
% 
% figure（）
% % imagesc（xaxiszaxisiEy‘）;
% imagesc（xaxiszaxisEy_fig‘）;
% xlabel（‘m‘）
% ylabel（‘m‘）
% caxis（[-0.30.3]）
% axis image;
% % colorbar;
% saveas（gcf[‘result_of_corr2_‘num2str（isrc）‘.png‘]）
% 
% 
% 
% save（[‘result_corr_‘num2str（isrc）‘.mat‘]‘iEy‘）;
% 
% % nt = length（Ey（:11））;
% % for i = 1:3751
% %     figure（1）
% %     subplot（121）
% %     Ey11 = reshape（Ey1（i::）nxny）;
% %     imagesc（Ey11）
% %     subplot（122）
% %     Ey21 = reshape（Ey2（nt+1-i::）nxny）;
% %     imagesc（Ey21）
% %     pause（0.01）;                                                                                                                                                                                                                                                              
% % end




end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\
     文件          44  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\.gitignore
     目录           0  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\
     文件        3262  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\ApplyImagingCondition.m
     文件        2659  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\GetWaveformDifference.m
     文件        5510  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\RTM02.m
     文件        3583  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\RTM03.m
     文件        4008  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\RTM03_181.m
     文件        5146  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TEGenerateData.m
     文件        5133  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TEGenerateData_step.m
     文件        5143  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TERunBackward.m
     文件        5158  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TERunBackward_sig.m
     文件        5141  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TERunForward.m
     文件        5161  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TERunForward_sig.m
     文件        6630  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TE_make_model.m
     文件        6845  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TE_make_model_01.m
     文件       24758  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\TE_model2d.m
     文件         994  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\finddt.m
     文件        2196  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\finddx.m
     文件        2644  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\get_step.m
     文件         520  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\get_tapper.m
     文件         822  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\gridinterp.m
     文件         778  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\padgrid.m
     文件         232  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\ricker.m
     文件         444  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\serial.sh
     文件        1021  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\2components\tapper.m
     文件        1755  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\ApplyImagingCondition.m
     文件        1811  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\GetWaveformDifference.m
     目录           0  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\JI_GPR_modeling_codes\
     文件       19646  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\JI_GPR_modeling_codes\TE_model2d.m
     文件        4516  2016-04-09 12:29  FWI_GPR_Matlab_TE-master\JI_GPR_modeling_codes\TE_run_example.m
............此处省略22个文件信息