非相参积累.zip_信号积累_相参中频信号_雷达_雷达信号_非相参

%中频非相参积累 简版 clc clear all close all k=100; %Hz 设置基本频率单位 A=4; % 脉冲串幅度 fs=40*k; %采样率 fc=10*k; %载波频率 N=10; %脉冲个数 fai=pi*rand(1,N); %载波随机相位 只要是在零到pi内依然会积累 fr=10; tr=1/fr; %脉冲频率 duty=5; %占空比 正比负 tao=(1/fr)*(duty/100);%脉宽 t=1/fs:1/fs:N/fr; %十个周期 ps=A*square(2*pi*fr*t,duty); %脉冲串 for i=1:length(ps) if ps(i)==-A; ps(i)=0; end end NTr=length(ps)/N; %一个周期 %产生射频信号 t1=0:1/fs:(NTr-1)/fs; c=zeros(1,length(ps)); for i=1:N a=cos(2*pi*fc*t1+fai(i)); %载波 c(NTr*(i-1)+1:NTr*(i))=a; end t=1/fs:1/fs:N/fr; %十个周期 %调制 m=c.*ps; figure(1) n=t*fr; plot(n,ps,n,m) %产生高斯白噪声 ns=wgn(1,length(m),0); figure(2) plot(n,ns) xlabel('t');ylabel('noise(t)');title('高斯白噪声'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %经过信道 figure(3) mr=m+ns; plot(n,mr) xlabel('t');ylabel('mr(t)');title('叠加噪声后的调制信号'); %非相参积累 mrn=zeros(1,NTr); for i=1:NTr for j=0:N-1 mrn(i)=mr(i+j*NTr)+mrn(i); end end %计算积累前信噪比 Nav=sum(ns.^2)/length(ns); %噪声功率 Ns=Nav*NTr; Ss=sum(mr(1:NTr).^2)-Ns; %单次信号功率 SNRs=Ss/Ns; %积累后 Sn=sum(mrn.^2)-sum(ns.^2);%N次信号功率 Nn=N*Ns; SNRn=Sn/Nn; %解调 t1=0:1/fs:(length(mr)-1)/fs; Is=mr.*cos(2*pi*fc*t1); Qs=mr.*(-sin(2*pi*fc*t1)); t2=0:1/fs:(length(mrn)-1)/fs; In=mrn.*cos(2*pi*fc*t2); Qn=mrn.*(-sin(2*pi*fc*t2)); %滤波 Fs = fs; % Sampling Frequency Fpass = 4*fr; % Passband Frequency Fstop = 20*fr; % Stopband Frequency Dpass = 0.057501127785; % Passband Ripple Dstop = 0.0001; % Stopband Attenuation flag = 'scale'; % Sampling Flag [n,Wn,BETA,TYPE] = kaiserord([Fpass Fstop]/(Fs/2), [1 0], [Dstop Dpass]); b = fir1(n, Wn, TYPE, kaiser(n+1, BETA), flag); Hd = dfilt.dffir(b); %psc=filter(Hd,psc); Is=filter(Hd,Is);Qs=filter(Hd,Qs); In=filter(Hd,In);Qn=filter(Hd,Qn); psc=sqrt(Is.^2+Qs.^2); figure(4) n=0:length(psc)-1; plot(n/fr,psc) %接收解调信号 xlabel('t');ylabel('psc(t)');title('信号对比图'); hold on plot(n/fr,ps) %发射信号 hold on pscn=sqrt(In.^2+Qn.^2); n2=0:length(pscn)-1; plot(n2/fr,pscn) %积累信号 Ss Sn SNRs SNRn