matlab中适合表示场强的colormap

我在用matlab画三维场强图时，找不到合适的colormap。最接近的colormap是matlab自带的jet，即随着函数值的增大颜色从蓝过渡到红。但由于场强图主要表征场的幅值特性，颜色区分只与函数值的绝对值相关，直接用jet函数就会导致场强最小位置的颜色为中间过渡色，而非蓝色。因此我利用colormap颜色RGB值渐变过渡的关系，编写了一个适用于场强表示的colormap。
程序效果见下图，本程序参考jet的色彩过渡关系，为了编写方便提升了色彩的纯度，采用blue[0,0,1]/cyan[0,1,1]/ yellow[1,1,0]/red[1,0,0]的纯色rgb值。
科研初学者，程序没啥含量，就是自己的一点小理解。

%% 电场colormap设定
n_step=9;
color_row=n_step*6+7;
color_map=zeros(color_row,3); 
color_r = 0:1/(n_step+1):1; 
color_g = 0:1/(n_step+1):1;
color_b = 0:1/(n_step+1):1;
% red to yellow
color_map(1:1+(n_step+1),1)=1;
color_map(1:1+(n_step+1),2)=color_g;
a=1+(n_step+1);
% yellow to cyan
color_map(a:a+(n_step+1),1)=fliplr(color_r);
color_map(a:a+(n_step+1),2)=1;
color_map(a:a+(n_step+1),3)=color_b;
a=a+(n_step+1);
% cyan to blue
color_map(a:a+(n_step+1),2)=fliplr(color_b);
color_map(a:a+(n_step+1),3)=1;
a=a+(n_step+1);
% blue to cyan
color_map(a:a+(n_step+1),2)=color_b;
color_map(a:a+(n_step+1),3)=1;
a=a+(n_step+1);
% cyan to yellow
color_map(a:a+(n_step+1),1)=color_r;
color_map(a:a+(n_step+1),2)=1;
color_map(a:a+(n_step+1),3)=fliplr(color_b);
a=a+(n_step+1);
% yellow to red
color_map(a:a+(n_step+1),1)=1;
color_map(a:a+(n_step+1),2)=fliplr(color_g);
colormap(color_map);