sqrc_time.m

本程序仿真一个从卫星到区域数据中心的下行链路（由一组6个载波构成）

function[imp,impw]=sqrc_time(R,beta,nsamples,nsymbols)
beta=beta*(R/2);
a=R+(2*beta);b=R-(2*beta);
length1=fix(nsamples*nsymbols/2);
ts=(1/R)/nsamples;
time=[-length1*ts:ts:(length1-1)*ts]+0.00000013;
n=length(time);w=hanning(n);
term1=a*pi*time;term2=b*pi*time;
cos1=cos(a*pi*time);sin1=sin(b*pi*time);
denominator=(pi*sqrt(R))*((1-((8*beta*time).^2)));
for k=1:n
numerator=(8*beta*(cos1(k)))+(sin1(k)/time(k));
if(denominator(k)==0)
imp(k)=1;
else
 imp(k)=numerator/denominator(k);
end
 impw(k)=imp(k)*w(k);
end
yy=sum(impw.^2);
impw=impw/(yy^0.5);
impw=impw/sqrt(nsamples);
yy=sum(imp.^2);
imp=imp/(yy^0.5);
imp=imp/sqrt(nsamples);