va=get(gcbo,'Value'),
if va==1,%%時域波形
    nnn=1:N;
    figure(16)
    set(16,'Position',[10,50,500,200])
    plot(nnn,s16);
    axis([0,10*4*L,-4,4]);		
    grid;
    
    
elseif va==2,%%頻域波形
    S=t2f(s16);    %S是s的傅氏變換
    P=S.*conj(S)/T;    %功率譜
    aa=30+10*log10(P+eps);    %加eps以避免除以零
    figure(16)
     set(16,'Position',[10,50,500,200])
    plot(f,aa,'b');
    grid
   axis([-4,+4,-60,60])
    xlabel('f  (MHz)')
    ylabel('Ps(f)   (dBm/MHz)')
    
    
else  %%%眼圖
    figure(16)
     set(16,'Position',[10,50,500,200]) 
    tt=[0:dt:Na*L*dt];
    for jj=1:Na*L:N-Na*L
        plot(tt,s16(jj:jj+Na*L));
        hold on
    end
    grid
end