%張潤(rùn)寧論文中多普勒補(bǔ)償問題?巴克碼（加噪聲）
%改變fd的值可以改變多普勒頻率的大小
%改變SNR的值可以改變信噪比的大小
clear?all;
T=0.2*10.^-6;?????????????????????????%子脈沖寬度T
B=1/T;????????????????????????????????%子脈沖帶寬B=1/T
fd=0.05*B??????????????????????????????%多普勒頻率fd????

SNR=20?????????????????????????????????%信噪比SNR?
Pnoise=0.01;??????????????????????????%噪聲功率Pnoise
Psignal=Pnoise*10.^（SNR/10）;??????????%信號(hào)功率Psignal
%Psignal=1;
%Pnoise=0

sequence=[1?1?1?1?1?-1?-1?1?1?-1?1?-1?1];

SeqLength=length（sequence）;??????????%序列長(zhǎng)度P=SeqLength?

%%%%%%%%%%%%%%%%%
%加多普勒信息和噪聲
Vector=1:SeqLength;
%SignalR=BPCSeq.*exp（1i*2*pi*fd*T*Vector）;????????????????????????????????????????????????????%回波信號(hào)加多普勒信息
rand（‘state‘sum（100*clock））;
SignalR=sqrt（Psignal）*sequence.*exp（1i*2*pi*fd*T*Vector）+sqrt（Pnoise）*randn（1SeqLength）;????????????%SignalR_N為加噪聲和多普勒信息的回波信號(hào)
%%%%%%%%%%%%%%%%%

s1=SignalR.*SignalR;

?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%用兩次內(nèi)插的方法來實(shí)現(xiàn)分頻處理
?u=s1（1:（SeqLength+1）/2）;
?for?k=0:（SeqLength-1）/2
?????%k
?????Even=2*k;
?????Odd=2*k+1;
?????v（Even+1）=u（Even/2+1）;
?????if?（Odd?????v（Odd+1）=（u（（Odd-1）/2+1）+u（（Odd+1）/2+1））/2;
?????end
?end
?s2=v;
?%s2為分頻后的結(jié)果
?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
?s3=SignalR.*conj（s2）;

matchfilter=sequence（SeqLength:-1:1）;??????????	%全景匹配濾波器

MFOut3=conv（s3matchfilter）;
[PSL_AISL_A]=PSLandISL（MFOut3）
mm_A=max（abs（MFOut3））

MFOut1=conv（SignalRmatchfilter）;
[PSL_1ISL_1]=PSLandISL（MFOut1）
mm_1=max（abs（MFOut1））

figure（1）
?subplot（223）
?plot（20*log10（0.01+abs（MFOut1）/max（abs（MFOut1））））;
?title（‘C（補(bǔ)償前）S/N=6dBfd/B=0.0‘）
?y=length（MFOut1）;
?axis（[0y-405]）;
?grid?on
?subplot（224）
?plot（20*log10（0.01+abs（MFOut3/max（abs（MFOut3）））））;
?title（‘D（補(bǔ)償后）S/N=6dBfd/B=0.0‘）
?y=length（MFOut3）;
?axis（[0y-405]）;
?grid?on
?%subplot（223）
?%plot（s）
?%subplot（224）
?%plot（s3）

%匹配濾波器的輸出基本是和論文中的圖一致的（P26）
%但是補(bǔ)償后的脈壓輸出有很大的不同
%問題估計(jì)存在于插值實(shí)現(xiàn)分頻的計(jì)算中
%經(jīng)過檢查，不是插值的問題
%完成，開始錯(cuò)在公式s3=s.*conj（s2）寫成s3=s2.*conj（s2）了
%如何在信號(hào)中加噪聲？？x=s+n??解決了見公式[s=s+Pnoise*rand（1SeqLength）]
%2002/3/25增加了計(jì)算SNR損失的函數(shù)PSLandISLandMSSLandLSNRandLPG

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-12-27?07:58??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\DopplerCom\
?????文件????????2349??2002-09-25?11:49??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\DopplerCom\DoplleComBarker.m
?????文件????????2501??2003-01-22?15:38??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\DopplerCom\DoplleComMSeries.m
?????文件????????3349??2003-01-02?11:44??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\DopplerCom\MultiSignal.m
?????文件?????????969??2002-08-26?17:08??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\DopplerCom\PSLandISL.m
?????文件?????????997??2002-09-24?14:32??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\DopplerCom\PSLandISLandMSSL.m
?????文件????????1437??2002-09-24?14:33??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\DopplerCom\PSLandISLandMSSLandLSNRandLPG.m
?????目錄???????????0??2015-12-27?07:58??相位編碼脈沖壓縮雷達(dá)的多普勒補(bǔ)償算法\