makeimpulse.m

matlab的源程序

function [h] = MakeImpulse(Source,Symbole,NbTrajet,Lambda,Tau,Te,k)

if (strcmp(Source.Modulation,'Lineaire')==1)
    if (Source.Debit==Te)
        h = 0;
        for (iTrajet=1:NbTrajet)
            if (k+2-iTrajet>0)
                h = h+Lambda(iTrajet)*Symbole(k+2-iTrajet);
            end
        end
    else
        alpha = 10; % Nombre de p閞iode dans la racine de Nyquist.
        h = 0;
        iTsmin = max(-alpha+floor(k*Te/Source.Debit),0);
        iTsmax = min(length(Symbole)-1,alpha+floor(k*Te/Source.Debit));
      
        for (iTs=iTsmin:iTsmax)
            Temp = 0;
            for (iTrajet=1:NbTrajet)
                Temp =  Temp+Lambda(iTrajet)*GetNyquist(Source.RollOff,Source.Debit,k*Te-iTs*Source.Debit-Tau(iTrajet));   
            end
	    
                Temp = Temp * Symbole(1,iTs+1);
            h = h + Temp;
        end
    end
    
elseif (strcmp(Source.Modulation,'CPM')==1)

    h = 0;
    for (iTrajet=1:NbTrajet)
        n = floor((k*Te-Tau(iTrajet))/Source.Debit);
        if (n>=Source.LongueurFiltre)
            Termes_Precedents = sum(Symbole(1:n-Source.LongueurFiltre+1));
        else
            Termes_Precedents = 0;
        end
        Symb = Symbole(n+1:-1:max(1,n-Source.LongueurFiltre+2));
        t = k*Te-Tau(iTrajet)-(n-(0:length(Symb)-1))*Source.Debit;
        if (strcmp(Source.FiltreMiseEnForme,'LRC')==1)
            Fonc = (t)/(Source.LongueurFiltre*Source.Debit)-1/(2*pi)*sin(2*pi*t/(Source.LongueurFiltre*Source.Debit))/(Source.LongueurFiltre*Source.Debit);
        else
            Fonc = (t)/(Source.LongueurFiltre*Source.Debit);
        end
        Terme_Courants = sum(Symb.*Fonc);
        h = h + Lambda(iTrajet)*exp(i*pi*Source.IndiceModulation*(Termes_Precedents+Terme_Courants));
    end
end