matlab画三维伯德图，bode图

做了一个简单的小例子，有什么疑惑可以留言。

clc
clear
close all
%% 初始化
R_Sta=1;    %R的起始，结束范围以及间隔
R_end=10;
Nap_R=1;
% 数组初始化
num=[R_Sta 1];
den=[5e-4 1];
w={1e-2,1e4};   %确定频率范围，单位：rad/s
[mag,phase,wout]=bode(num,den,w);   % 先要运行一次bode函数以确定输出量的长度

R=R_Sta:Nap_R:R_end;
Num_R=length(R);
Num_w=length(wout);
Z_mag = ones(Num_R,Num_w);
Z_phase = ones(Num_R,Num_w);
[X,Y]=meshgrid(wout,R);

R=R_Sta;
for i=1:Num_R
    num=[R 1]
    den=[5e-3 1]
    [mag,phase,wout]=bode(num,den,w);
    Z_mag(i,:)=20*log10(mag);
    Z_phase(i,:)=phase;
    R=R+Nap_R;
end

% figure(1)
% bode(num,den);
figure(2)
C_mag=Z_mag;
surf(X,Y,Z_mag,C_mag)
set(gca,'XScale','log')
colorbar    
shading interp   % 去掉网格
xlabel('rad/s');ylabel('R');zlabel('mag')
title('Magnitude')
figure(3)
C_phase=Z_phase;
surf(X,Y,Z_phase,C_phase)
set(gca,'XScale','log')
colorbar    
shading interp   % 去掉网格
xlabel('rad/s');ylabel('R');zlabel('Phase')
title('Phase')

输出图形如下：