bpsk.rar_协作检测_协作频谱感知_感知_频谱感知 两_频谱感知结果

% 2bsk 　　叠加噪声以后求出来的　循环自相关 db=0 时 位移为1600 与scd 作比较 clear all; close all; A=2; % 假定载波每比特能量为单位能量，则其振幅为sqrt(2) fc=10; % 10Hz N_sample=20; % 每个码元的抽样点数 N=200; %码元数 Ts=1; % 1 Baud/s %dt =Ts/N_sample; % 波形抽样间隔 dt =Ts/fc/N_sample; % 波形抽样间隔 t=0:dt:N*Ts-dt; Lt=length(t); % 产生二进制信源 d=sign(randn(1,N)); dd=sigexpand(d,fc*N_sample); %对序列插入fc*N-sample-1个0 %dd=sigexpand((d+1)/2,fc*N_sample); % 单 极性 NRZ波形 gt=ones(1,fc*N_sample); % 双极性 NRZ波形 d_NRZ=conv(dd,gt); %figure(1) %subplot(221); %输入 NRZ 信号波形（双极性） %plot(t,d_NRZ(1:length(t))); %axis([0 20 -2 2]); xlabel('输入双极性 NRZ 信号') %subplot(222) %[f,d_NRZf]=T2F(t,d_NRZ(1:length(t))); %plot(f,10*log10(abs(d_NRZf).^2/N*Ts)); %axis([-2 2 -50 10]); xlabel('输入双极性 NRZ 信号功率谱密度dB/Hz'); ht=sqrt(A)*cos(2*pi*fc*t); % 2psk 信号 d_2psk=d_NRZ; s_2psk=d_2psk(1:Lt).*ht; EsNodB=0 % 最小为 -2.5ｄｂ EsNo=10.^(-EsNodB); Es=1; sigma=sqrt(EsNo/2); for k=1:length(EsNodB) r=s_2psk +sigma(k)*randn(1,length(s_2psk)) %figure(2) %subplot(221) %plot(t,s_2psk); %axis([0 10 -1.2 1.2]); xlabel('2psk 信号') %subplot(222) %[f,s_2pskf]=T2F(t,s_2psk); %plot(f,10*log10(abs(s_2pskf).^2/N*Ts)); %axis([-fc-4 fc+4 -50 10]); xlabel('2psk 信号功率谱密度dB/Hz') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Mlag=2000; % 700 % [r1,lags]=xcorr(s_2psk,Mlag,'unbiased'); [r1,lags]=xcorr(r,Mlag,'unbiased'); %[f,s_2pskf]=T2F(r1,lags); %这个函数求出的效果不好。 %figure(2) % plot(f,abs(lags)); ts=0.01; df=0.8; [Corr,corr,df1]=fftseq(r1,ts,df); %s0 fs=1/ts; fy=[0:df1:df1*(length(corr)-1)]-fs/2; figure(3) plot(fy,abs(fftshift(Corr))); % 发射信号的循环自相关函数的结果 %axis([13 fc+5 0 200]); axis([0 50 0 140]) xlabel('Cycle Frequency') ylabel('Amplitude of CAF') grid on; %pause end