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

end