有大神能教我怎么画出下面这个三维图的么？

wdz

有大神能教我怎么画出下面这个三维图的么？
命令如下：
clear;
clc;
%用Clough-Penzien模型生成第一条功率谱函数
%8度设防，采用第I0类场地，设计地震分组第三组,多遇
wg=20.94;                        
wf=0.2*wg;
kesig=0.64;
kesif=kesig;
S0=4.99;%还没进行修正
w0=0;
wn=40;
NN=500;
dw=(wn-w0)/NN;
for i=1:NN
    w(i)=w0+(i-1)/NN*wn;%论文P93，三角级数第四种方法，ak根据已知的sw求得
    sw1(i)=(wg^4+4*kesig^2*wg^2*w(i)^2)/((wg^2-w(i)^2)^2+4*kesig^2*wg^2*w(i)^2)*w(i)^4/((wf^2-w(i)^2).^2+4*kesif^2*wf^2*w(i)^2)*S0;
    A(i)=sqrt(4*sw1(i)*dw);%论文P93，三角级数第四种方法，ak根据已知的sw求得
end;
phi=2*pi*rand(1,500);%随机相位角
t1=0.5;
t2=5.5;
c=0.45;
t=zeros(1,669);%1-1500的0矩阵，一行？相当于6.69秒
y1=t;
for i=1:669
    t(i)=0+(i-1)*0.01;
    if t(i)<t1
       g(i)=(t(i)/t1)^2;
    elseif t(i)>=t1 & t(i)<t2
       g(i)=1;
    elseif t(i)>=t2
       g(i)=exp(-c*(t(i)-t2));
    end;
    for j=1:NN
        as(j)=A(j)*cos(w(j)*t(i)-phi(j));%三角级数法
    end;
    a(i)=sum(as);
    y1(i)=a(i)*g(i);%  公式5-27
end;
请指出画图命令。。。。

wdz

频率轴其实就是功率谱密度函数，时间轴就是时间包络函数，不知道怎么画成这样的三维图，求大神帮助