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

%% （target  1）    Get_Optimal_Position_Band    
% 1_5b
X_best=[P_1_best(1);P_2_best(1);P_3_best(1);P_4_best(1);P_5_best(1)];
Y_best=[P_1_best(2);P_2_best(2);P_3_best(2);P_4_best(2);P_5_best(2)];
Z_best=[P_1_best(3);P_2_best(3);P_3_best(3);P_4_best(3);P_5_best(3)];
Delta_k_best=[Max_1;Max_2;Max_3;Max_4;Max_5];
figure(N);
N=N+1;
subplot(1,3,1);
plot3(X_best,Y_best,Z_best,'*','markersize',12);
xlabel('x/b'), ylabel('y/b'), zlabel('z/b');
grid on;
hold on;
% 稀疏
subplot(1,3,2);
[yi,zi] = meshgrid(-0.375:precision:0.375, -0.375:precision:0.375);
xi = griddata(Y_best,Z_best,X_best,yi,zi,'cubic');
surf(xi,yi,zi);
xlabel('x/b'), ylabel('y/b'), zlabel('z/b')
% 密
subplot(1,3,3);
[x_i,y_i] = meshgrid(1.4:precision:15.1, -0.375:precision:0.375);
z_i = griddata(X_best,Y_best,Z_best,x_i,y_i,'cubic');
surf(x_i,y_i,z_i);
xlabel('x/b'), ylabel('y/b'), zlabel('z/b')
[m,n] = size(x_i);
Optimal_Position_Band_x=reshape(x_i,m*n,1);
Optimal_Position_Band_y=reshape(y_i,m*n,1);
Optimal_Position_Band_z=reshape(z_i,m*n,1);
Optimal_Position_Band_P=[Optimal_Position_Band_x,Optimal_Position_Band_y,Optimal_Position_Band_z];
[Optimal_Position_Band_P] = IsNan_Delte_row(Optimal_Position_Band_P,3);

% v=
% F = scatteredInterpolant(dddd,v);
% [m,n] = size(x_i);
% dddd=reshape(x_i,m*n,1);
% % % 
% figure(N);
% N=N+1;
% subplot(1,2,1);
% [yi,zi] = meshgrid(-0.375:precision:0.375, -0.375:precision:0.375);
% Delta_k_best_i_7 = griddata(Y_best,Z_best,Delta_k_best,yi,zi,'cubic');
% surf(Delta_k_best_i_7,yi,zi);
% xlabel('Delta_k_best_i_6'), ylabel('y/b'), zlabel('z/b')
% % 
% subplot(1,2,2);
% [x_i,y_i] = meshgrid(1.4:precision:15.1, -0.375:precision:0.375);
% Delta_k_best_i_7m = griddata(X_best,Y_best,Delta_k_best,x_i,y_i,'cubic');
% surf(x_i,y_i,Delta_k_best_i_7m);
% xlabel('x/b'), ylabel('y/b'), zlabel('Delta_k_best_i_6')