MATLAB三维曲面绘图求助！

jaj2008

如何根据下面一段代码画出如图的三维曲面：
clc;
clear all;

wn =2*500*pi; ipson=0.038;
k=4500; k_t=2400; k_r=0.9;
N=2; en=0;

u=0.01:0.01:1;
ex=acos(1-u);
ayy=0.5*(-cos(2*ex)-2*k_r*ex-k_r*sin(2*ex))-0.5*(-cos(2*en)-2*k_r*...
en-k_r*sin(2*en));

w=2000:1:10000;w=w';
s=1i*w;
G_y=wn^2./(k*(wn^2+s.^2+2*ipson*wn*s));

f=@(y) @(a) 1./((N/2*pi)*a*k_t*real(y));
ff=f(G_y);

ap=arrayfun(ff,ayy,'UniformOutput',false);
app=cell2mat(ap);
save data app;
试了很久了，就是不行。不知道是不是u取值问题，但是试了0.001:u:1,曲面画不出来，现在画出来的效果不好。如下图：

jaj2008

上面是u和ap,下面是w,
wn=2*500*pi;
ipson=0.038;
k=4500;
k_t=2400;
k_r=0.9;
N=2;
en=0;
ex=1.5708;
ayy=0.5*(-cos(2*ex)-2*k_r*ex-k_r*sin(2*ex))-0.5*(-cos(2*en)-2*k_r*en-k_r*sin(2*en));
fid=fopen('e:result2\qieshen3.txt','w');
fid1=fopen('e:result2\zhuansu3.txt','w');
for m=0:4
for w=2000:1:10000
s=i*w;
d=w/wn;
G_y=wn^2/(k*(wn^2+s^2+2*ipson*wn*s));
ap=1/((N/2*pi)*ayy*k_t*real(G_y));
fprintf(fid,'%f\r\n',ap);
n=60*w/(N*(2*m*pi+2*pi-2*atan((d^2-1)/2*ipson*d)));
fprintf(fid1,'%f\r\n',n);
end
end
fclose(fid1);
fclose(fid);