2017年4月19日足协杯第2轮江西联盛vs北京北控直播

%% clear; close all; clc; %参数设置 data_length=5000; fc= 120000;%频偏120000Hz fb=10*10^6;%符号速率10MHz fs=40*10^6;%采样速率40MHz %% %产生码、串并转换、产生I、Q两路 rand_data=(randn(1,data_length)>0.5)*2-1; %seriel to parallel同时单极性码转为双极性码 for i=1:data_length p=floor((i-1)/2)*2; a(i)=rand_data(p+1); b(i)=rand_data(p+2); end; %% %调制 Delta_Freq=fc; Simulation_Length=floor((data_length*1/fb)*fs); for i=1:Simulation_Length-1 temp=floor((i*1/fs)*fb); I_Data(i)=a(temp+1); Q_Data(i)=b(temp+1); end Simulation_Length=Simulation_Length-1; Signal_Source=I_Data + j*Q_Data; Carrier=exp(j*(2*pi*Delta_Freq/fs*(1:Simulation_Length))); %Imod(i)=I(i)*cos(2*pi*fc*t(i)); %Qmod(i)=Q(i)*sin(2*pi*fc*t(i)); %y=Imod-j*Qmod; %QPSK_Out=awgn(y,SNR); %I_Data=I; %Q_Data=Q; %调制处理 Signal_Channel_1=Signal_Source.*Carrier; %Signal_Channel=Signal_Source.*Carrier; SNR=1000; Signal_Channel=awgn(Signal_Channel_1,SNR); %Signal_Channel=Signal_Channel.'; %% %costas环载波恢复 %频偏，Hz close all; Time_Sample=1/fs; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %以下为锁相环处理过程 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %参数清零 Signal_PLL=zeros(Simulation_Length,1); NCO_Phase = zeros(Simulation_Length,1); Discriminator_Out=zeros(Simulation_Length,1); Freq_Control=zeros(Simulation_Length,1); PLL_Phase_Part=zeros(Simulation_Length,1); PLL_Freq_Part=zeros(Simulation_Length,1); %环路处理 C1=0.22013; C2=0.0024722; for i=2:Simulation_Length Signal_PLL(i)=Signal_Channel(i)*exp(-j*mod(NCO_Phase(i-1),2*pi));%%? I_PLL(i)=real(Signal_PLL(i)); Q_PLL(i)=imag(Signal_PLL(i)); Discriminator_Out(i)=(sign(I_PLL(i))*Q_PLL(i)-sign(Q_PLL(i))*I_PLL(i))... /(sqrt(2)*abs(Signal_PLL(i))); PLL_Phase_Part(i)=Discriminator_Out(i)*C1; Freq_Control(i)=PLL_Phase_Part(i)+PLL_Freq_Part(i-1); PLL_Freq_Part(i)=Discriminator_Out(i)*C2+PLL_Freq_Part(i-1); NCO_Phase(i)=NCO_Phase(i-1)+Freq_Control(i); end %画图显示结果 figure subplot(2,2,1) plot(-PLL_Freq_Part(2:Simulation_Length)*fs/(2*pi)); grid on; title('锁相环频率响应曲线'); % axis([1 Simulation_Length -100 100]); subplot(2,2,2) plot(PLL_Phase_Part(2:Simulation_Length)*180/pi); title('锁相环相位响应曲线'); axis([1 Simulation_Length -2 2]); grid on; %设定显示范围 Show_D=1; %起始位置 Show_U=Simulation_Length; %终止位置 % Show_D=300; %起始位置 % Show_U=900; %终止位置 Show_Length=Show_U-Show_D; subplot(2,2,3) plot(Signal_Channel(Show_D:Show_U),'*'); title('进入锁相环的数据星座图'); axis([-2 2 -2 2]); grid on; hold on; subplot(2,2,3) plot(Signal_PLL(Show_D:Show_U),'r*'); grid on; subplot(2,2,4) plot(Signal_PLL(Show_D:Show_U),'r*'); title('锁相环锁定及稳定后的数据星座图'); axis([-2 2 -2 2]); grid on; %% figure %设定显示范围 Show_D=Simulation_Length-200; %起始位置 Show_U=Simulation_Length-1; %终止位置 Show_Length=Show_U-Show_D; subplot(2,2,1) plot(I_Data(Show_D:Show_U)); grid on; title('I路信息数据'); axis([1 Show_Length -2 2]); subplot(2,2,2) plot(Q_Data(Show_D:Show_U)); grid on; title('Q路信息数据'); axis([1 Show_Length -2 2]); subplot(2,2,3) plot(I_PLL(Show_D:Show_U)); grid on; title('锁相环输出I路信息数据'); axis([1 Show_Length -2 2]); subplot(2,2,4) plot(Q_PLL(Show_D:Show_U)); grid on; title('锁相环输出Q路信息数据'); axis([1 Show_Length -2 2]); figure(3); plot(NCO_Phase(Show_D:Show_U))