guassionpulse.m

关于超宽带的一系列仿真包括高斯波形各阶导函数波形频谱以及采用高斯导函数的跳时超宽带仿真程序

function cp0701_shape_factor_variation(alphamin,alphastep,N_alphavalues)
alphamin=0.414e-9;
alphastep=0.1e-9;
N_alphavalues=7;
A=1;
smp=1024;
Tmin=-4e-9;
Tmax=4e-9;
alpha=alphamin;
t=linspace(Tmin,Tmax,smp);
for i=1:N_alphavalues
    pulse=-A*exp(-2*pi*(t/alpha).^2);
    dt=(Tmax-Tmin)/smp;
    fs=1/dt
    N=smp;
    df=1/(N*dt);
    X=fft(pulse);
    X=X/N;
    E=fftshift(abs(X).^2/(df^2));
    Ess=2*E((N/2+1):N);
    figure(1);
    PT=plot(t,pulse);
    set(PT,'LineWidth',[2]);
    AX=gca;
    set(AX,'FontSize',12);
    T=title('时域');
    X=xlabel('时间[s]');
    set(X,'FontSize',14);
    Y=ylabel('幅度[V]');
    set(Y,'FontSize',14);
    alphabehaviour={'增加\alpha'};
    text(0.75e-9,-0.5,alphabehaviour,'BackgroundColor',[1 1 1]);
    axis([-2e-9 2e-9 -1.2 1.2]);
    hold on
    figure(2);
    positivefrequency=linspace(0,(fs/2),N/2);
    PF=semilogy(positivefrequency,Ess);
    set(PF,'LineWidth',[2]);
    AX=gca;
    set(AX,'FontSize',12);
    T=title('频域');
    set(T,'FontSize',14);
    X=xlabel('频率[Hz]');
    set(X,'FontSize',14);
    Y=ylabel('功率谱密度 [(V^2)*sec/Hz]');
    set(Y,'FontSize',14);
    axis([0 20e9 1e-60 1e-10]);
    text(7.5e9,1e-25,alphabehaviour,'BackgroundColor',[1 1 1]);
    hold on
    alpha=alpha+alphastep;
end