使用matlab对一组数据的处理以及研究思路路线（4）--Get_Map_the_xxx_of_The_Facet.m\Printing.m

%% 得到初步 
N=1;
[Delta_k_1,N]=Printing(y_1_5b,z_1_5b,Delta_K_1_5b,N,'Delta_K_1.5b',precision);   
[Delta_k_2,N]=Printing(y_3b  ,z_3b  ,Delta_K_3b  ,N,'Delta_K_3b'  ,precision);  
[Delta_k_3,N]=Printing(y_5b  ,z_5b  ,Delta_K_5b  ,N,'Delta_K_5b'  ,precision);  
[Delta_k_4,N]=Printing(y_10b ,z_10b ,Delta_K_10b ,N,'Delta_K_10b' ,precision);  
[Delta_k_5,N]=Printing(y_15b ,z_15b ,Delta_K_15b ,N,'Delta_K_15b' ,precision); 

function [Delta_k_i,N]= Printing(Y,Z,Delta_k,N,b,precision)
% N为第几幅图片 b为z轴标识  precision ：为精度
%PRINTING 此处显示有关此函数的摘要
%   此处显示详细说明
%% 1_5b
figure (N)
subplot(1,2,1)
plot3(Y,Z,Delta_k,'.','markersize',8)
xlabel('y/b'), ylabel('z/b'), zlabel(b)
grid on
hold on
%%
subplot(1,2,2)
[yi,zi] = meshgrid( -0.375:precision:0.375, -0.375:precision:0.375);
Delta_k_i = griddata(Y,Z,Delta_k, yi, zi,'cubic');
surf(yi,zi,Delta_k_i);
xlabel('y/b'), ylabel('z/b'), zlabel(b)
N=N+1;
end

% %% 1_5b
% figure (1)
% subplot(1,2,1)
% plot3(y_1_5b,z_1_5b,Delta_K_1_5b,'.','markersize',12)
% xlabel('y'), ylabel('z'), zlabel('Delta_K_1.5b')
% grid on
% hold on
% %%
% subplot(1,2,2)
% [yi,zi] = meshgrid( -1:0.01:1, -1:0.01:1);
% Delta_K_1_5b_i = griddata(y_1_5b,z_1_5b,Delta_K_1_5b, yi, zi);
% surf(yi,zi,Delta_K_1_5b_i);
% xlabel('y'), ylabel('z'), zlabel('Delta_K_1.5b')