理論內(nèi)容基本在上一堂課講清楚了，那么現(xiàn)在就該來(lái)看仿真程序了！通過(guò)程序來(lái)進(jìn)一步的消化位定時(shí)算法！

每一次寫仿真程序都是一次完善的過(guò)程，仿真程序?qū)懙脑俸?，它也?huì)出現(xiàn)或多或少的問(wèn)題，同時(shí)還要必須不斷地向?qū)嶋H的情況去靠攏，這才是寫程序最終的目的仿真，這是基礎(chǔ)，最終能夠和實(shí)際產(chǎn)品應(yīng)用才是目的，所以我們?cè)趯W(xué)仿真的時(shí)候，一定要清楚程序和實(shí)際產(chǎn)品最終還有多大的距離。因此，每次寫程序都是一個(gè)長(zhǎng)時(shí)間的不斷修訂的過(guò)程！大家如果走上算法研究的道路，務(wù)必記?。翰粩嗟耐晟品抡姝h(huán)境是一個(gè)漫長(zhǎng)的過(guò)程，也是不斷深化自己理解的過(guò)程！

先看看網(wǎng)絡(luò)上的一段程序！

close all;clear all;

%%************** 程序主體 **************%%

N=20000; % 符號(hào)數(shù)

K=4; % 每個(gè)符號(hào)采4個(gè)樣點(diǎn)

Ns=K*N; % 總的采樣點(diǎn)數(shù)

w=[0.5,zeros(1,N-1)];

%環(huán)路濾波器輸出寄存器，初值設(shè)為0.5

n=[0.7 zeros(1,Ns-1)];

%NCO寄存器，初值設(shè)為0.9

n_temp=[n(1),zeros(1,Ns-1)];

u=[0.6,zeros(1,2*N-1)];

%NCO輸出的定時(shí)分?jǐn)?shù)間隔寄存器，初值設(shè)為0.6

yI=zeros(1,2*N);

%I路內(nèi)插后的輸出數(shù)據(jù)

yQ=zeros(1,2*N);

%Q路內(nèi)插后的輸出數(shù)據(jù)

time_error=zeros(1,N);

%Gardner提取的時(shí)鐘誤差寄存器

i=1;

%用來(lái)表示Ts的時(shí)間序號(hào),指示n,n_temp,nco,

k=1;

%用來(lái)表示Ti時(shí)間序號(hào),指示u,yI,yQ

ms=1;

%用來(lái)指示T的時(shí)間序號(hào),用來(lái)指示a,b以及w

strobe=zeros(1,Ns);

% 環(huán)路濾波器系數(shù)

c1=5.41*10^(-3);

c2=3.82*10^(-6);

%%% %%%PSK基帶數(shù)據(jù)%%% %%%

bitstream=(randi(2,N,1)-1);

psk2=pskmod(bitstream,2);

xI=zeros(1,Ns);

xQ=zeros(1,Ns);

xI(18*N) = real(psk2);

%8倍插值，后面2倍抽取實(shí)現(xiàn)輸入采樣信號(hào)是輸出插值符號(hào)速率的4倍

%即4個(gè)采樣點(diǎn)確定一個(gè)插值

xQ(18*N) = imag(psk2);

% 截短后的根升余弦匹配濾波器

h1 = rcosfir(0.8,[-8,8],4,1,'sqrt');

aI1 = conv(xI,h1);

bQ1 = conv(xQ,h1);

L = length(aI1);

%仿真輸入數(shù)據(jù)

aI = [aI1(22L),0,0]; %2倍抽取

bQ = [bQ1(22L),0,0];

%%%%%%% Gardner算法

ns=length(aI)-2;

while(i

n_temp(i+1) = n(i)-w(ms);

if(n_temp(i+1)>0)

n(i+1)=n_temp(i+1);

else

n(i+1)=mod(n_temp(i+1),1);

% 內(nèi)插濾波器

FI1 = 0.5*aI(i+2)-0.5*aI(i+1)-0.5*aI(i)+0.5*aI(i-1);

FI2 = 1.5*aI(i+1)-0.5*aI(i+2)-0.5*aI(i)-0.5*aI(i-1);

FI3 = aI(i);

yI(k) = (FI1*u(k)+FI2)*u(k)+FI3;

FQ1 = 0.5*bQ(i+2)-0.5*bQ(i+1)-0.5*bQ(i)+0.5*bQ(i-1);

FQ2 = 1.5*bQ(i+1)-0.5*bQ(i+2)-0.5*bQ(i)-0.5*bQ(i-1);

FQ3 = bQ(i);

yQ(k) = (FQ1*u(k)+FQ2)*u(k)+FQ3;

strobe(k) = mod(k,2);

% 時(shí)鐘誤差提取模塊

if(strobe(k)==0)

% 每個(gè)符號(hào)計(jì)算一次時(shí)鐘誤差

if(k>2)

time_error(ms) = yI(k-1)*(yI(k)-yI(k-2))+yQ(k-1)*(yQ(k)-yQ(k-2));

else

time_error(ms) = (yI(k-1)*yI(k)+yQ(k-1)*yQ(k));

end

%環(huán)路濾波器,每個(gè)數(shù)據(jù)符號(hào)計(jì)算一次環(huán)路濾波器輸出

if(ms>1)

w(ms+1)=w(ms)+c1*(time_error(ms)-time_error(ms-1))+c2*time_error(ms);

else

w(ms+1)=w(ms)+c1*time_error(ms)+c2*time_error(ms);

end

ms=ms+1;

end

k=k+1;

u(k)=n(i)/w(ms);

end

i=i+1;

編輯：jq