tic
clc;
clf;
BW=5*10^6;%bandwidth
T=51.2*10^-6;%pulse?width
fc=0;?%zero?centered?frequency
n=fix（2*BW*T+1）;%?sampling?points
sampling_interval?=?T?/?n;%sampling?interval
freqlimit?=?0.5/?sampling_interval;
alpha=BW/（T）;%chirp?rate
t=linspace（-T/2T/2n）;%time?interval
ichannel1=cos（2*pi*（（fc）*t+alpha*t.^2））;%real?part
qchannel2=sin（2*pi*（（fc）*t+alpha*t.^2））;%imaginary?part
LFMD=0.0892*（ichannel1+（1i.*qchannel2））;%Complex?base?band?time?domain?signal
LFMD（2:1:10）=0;
LFMD（504:1:511）=0;
freq?=?linspace（-freqlimit/2freqlimit/2n）;
FLFMD=fftshift（fft（LFMDn））;%Frequency?domain?of?base?band?signal
[lfmfreiflfm]=lffm1（BWnT）;%Function?calling?for?input?Lfm?signal?
figure（1）
subplot（311）
plot（freiabs（flfm））;
plot（freqabs（（FLFMD）））;
xlabel（‘Frequency?domain‘）
ylabel（‘Amplitude‘）
title（‘Input?signal?of?pulse?compression?filter?in?Frequency?domain‘）;
grid;
%plot（t（（LFMD））‘k‘）
%?plot（t（lfm）‘k‘）;
rep=fliplr（LFMD）;
%??disp（length（rep））;
j=0;
y1=zeros（110001）;%10Khzs?Sampling?domain???
for?i=2246:1:2758
????j=j+1;
????y1（i）=LFMD（j）;
end
%?????j=0;
%??????for?i=3976:1:4488
%?????j=j+1;
%?????y1（i）=LFMD（j）;
%??????end
?
%????for?i=9246:1:9758
%?????j=j+1;
%?????y1（i）=LFMD（j）;
%???end

Q=length（lfm）;%length?of?lfm?signal?h（k）
P=length（y1）;%length?of?base?band?signal?domain?X（k）?
N=P+Q-1;
N=2.^nextpow2（N）;%?FFTbuffer
x1=zeros（1N）;%creating?a??empty?16k?Buffer
%?creating?stored?replica?buffer?of?16k
for?k1=1:1:N
????if?k1<=Q
????????h1（k1）=rep（k1）;
????else
????????h1（k1）=0;
????end
????
end
hn=fftshift（fft（h1N））;?%taking?FFt?of?16k?H（k）buffer
hnc=（conj（hn））;%
%creating?the?16k?buffer?for?X（k）signal
n3=0;
for?k2=1:1:N
?????if?k2<=N-P+1
?????????x1（k2）=0;
?????else
?????????n3=n3+1;
????????x1（k2）=y1（n3）;
?????end
end

??x1=awgn（x10.1）;
?subplot（211）;
?tt=linspace（01*10^-6N）;
?plot（ttx1）;
?xlabel（‘Time?domain‘）
?ylabel（‘Amplitude‘）
title（‘Noise?along?with?the?signal?‘）;
%??grid;
?xf=（fftshift（fft（x1N）））;?%fft?of?the?16k?X（k）?buffer
%computing?the?product?of?X（k）*h（k）
zn=xf.*hnc;
%IFFT?of?the?Multiplier?output?buffer
zk=（ifft（znN））;
zk=abs（（zk））;
zk1=zk（1:1:8192）;
var1=var（zk1）;
zk2=zk（9590:1:16000）;
var2=var（zk2）;
n12=var1-var2;
t4=linspace（0NN）;
subplot（212）
plot（t4（zk））;
xlabel（‘Time?domain‘）
ylabel（‘Amplitude‘）
title（‘?Mainlobe?width?analysis?for?the?pulse?compression?filter‘）;
grid;
toc;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2478??2010-05-18?07:24??dd1.asv
?????文件????????2441??2010-05-14?00:52??dd1.m
?????文件????????2319??2010-05-13?23:57??direct.m
?????文件?????????620??2010-05-13?22:18??lffm1.m
?????文件???????42555??2010-06-20?16:51??llll.mdl
?????文件????????1883??2010-06-22?00:39??matched2.m
?????文件????????1890??2010-05-14?01:15??matched21.m
?????文件????????4741??2010-04-05?15:24??ovsav1.m
?????文件????????2500??2010-04-13?13:02??ovsav2.m
?????文件????????2050??2010-04-05?15:24??ovsav3.m
?????文件????????1836??2010-04-05?14:49??ovsav4.m
?????文件????????2150??2010-05-14?01:12??ovsav5?（2）.m
?????文件????????2033??2010-05-19?06:58??ovsav5.asv
?????文件????????2150??2010-05-14?01:12??ovsav5.m
?????文件????????2124??2010-05-14?01:16??ovsav51?（2）.m
?????文件????????2124??2010-05-14?01:16??ovsav51.m
?????文件???????22546??2010-06-20?17:27??sineee.mdl
?????文件????????1306??2014-02-12?13:15??license.txt