ascend1.m

来自「一种用于ＧＭＳＫ非线性调制的自适应均衡算法 」· M 代码 · 共 38 行

function prt=ascend1(Rs,fs,factor)

%Rs=128000;  %比特速率
%nn=5;       %每个比特的采样点数nn
%fs=nn*Rs;   %采样速率
%factor=1/4;

Bda=factor*Rs;
Ts=1/fs;
N=150;
k=4*fs/Rs;
prt=zeros(1,N);
for i=0:N-1
t=i/fs;
if t==0
   prt(1,i+1)=1;

else
   if (1-4*Bda*t*Bda*t)==0
   prt(1,i+1)=1/2*factor*sin(pi/(2*factor));
 else    
p1=sin(pi*Rs*t)/((pi*Rs*t)*(1-4*Bda*t*Bda*t));
p2=cos(pi*Bda*t);
prt(1,i+1)=p1*p2;
end
end
end


%prt1=prt(1:18);
%prt2=fliplr(prt1);
%prt=[prt2(1:17),prt1];

prt1=prt(1:21);
prt2=fliplr(prt1);
prt=[prt2(1:20),prt1];
a=1;