matlab的costas环载波恢复,环路滤波器系数有推导.doc

《matlab的costas环载波恢复,环路滤波器系数有推导.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、fs=50e6;%采样频率ts=1/fs; num=2e6;  %数据长度SNR=-15;real_fc=;%实信号频率data=sin(2*pi*real_fc*(0:num-1)*ts+pi/4)+sqrt(10^(SNR/10))*randn(1,num);%科斯塔斯环的输入信号fc=;%本地频率n=fs/10000;%累积时间为0.1msnn=[0:n-1];nf=floor(length(data)/n);%将输入数据分成1ms的多个数据块wfc=2*pi*fc;  %本地信号phi_prv=0;temp=0;frame=0;carrie

2、r_phase=0;phase=0;%环路滤波器的参数c1=1203.9*1.5;c2=16.22*1.5; %数据率为4kbps,1bit数据采样点数为12.5K，每5K个点累加一次，相当于4个数据累加清零一次，%w=0.01*4K,T=1/4K,k0=2*pi*T*fs/2^32,c1=2*0.707*w*T/k0,c2=(wT)^2/k0,forframe=1:nf %产生本地的sin和cos函数expcol=exp(j*(wfc*ts*nn+phase));sine=imag(expcol);  cosine=real(expcol);x=

3、data((1:n)+((frame-1)*n));%将数据转换到基带x_sine=x.*sine;x_cosine=x.*cosine;Q=sum(x_sine);    %经过滤波器I=sum(x_cosine);phase_discri(frame)=atan(Q/I);  %得到锁相环的输入%锁相环dfrq=c1*phase_discri(frame)+temp;%经过环路滤波器temp=temp+c2*phase_discri(frame);wfc=wfc-dfrq*2*pi;  %改变本地频率dfrq_frame(frame)=wfc;

4、 phase=wfc*ts*n+phase;  %得到不同块的相位dphase_frame(frame)=phase; endplot(dfrq_frame/(2*pi));holdonplot([1:length(dfrq_frame)],real_fc,'r');legend('锁相环跟踪','实际的载波频率');gridmean_freq=mean(dfrq_frame/2/pi)p=abs(real_fc-mean_freq)/real_fc;