clc;
clear?all;
close?all;
fs=2e6;%采樣率
fc=1e10;%載波
R1=50000;R2=30000;v1=10;v2=3;c=3e8;
Tm=500e-6;%調制周期
T0=5e-6;%脈寬
period_counter=64;%表示經歷的周期的數目
t=0:1/fs:period_counter*Tm;
tm=0:1/fs:Tm;
counter_tm=size（tm）;%一個周期采樣點的個數
t0=1/fs:1/fs:T0;%一個脈寬采樣點
counter_t0=size（t0）;%一個脈寬采樣點數
%發射信號
send=0;
for?i=1:period_counter
????send=[send?ones（1T0*fs）?zeros（1（Tm-T0）*fs）];%
end;
figure（1）;
plot（send）;
axis（[0?1e4?-0.5?2.0]）
title（‘發射信號‘）;xlabel（‘采樣點數‘）;ylabel（‘幅度‘）;
%回波信號
%有遲延
receive=0;
t_delay1=2*R1/c;
t_delay2=2*R2/c;
delay_array1=zeros（1fix（t_delay1*fs））;%fix:去尾
delay_array2=zeros（1fix（t_delay2*fs））;
receive1=[delay_array1?send];
receive2=[delay_array2?send?zeros（1fix（t_delay1*fs）-fix（t_delay2*fs））];
receive=0.3*receive1+receive2;
figure（2）;
plot（receive）;
axis（[0?1e4?-2.0?2.0]）
title（‘固定目標反射的回波‘）;xlabel（‘采樣點數‘）;ylabel（‘幅度‘）;
%有多普勒調制
fd1=200;fd2=800;
tt=0:1/fs:period_counter*Tm+t_delay1;
receive_fd1=receive1.*cos（2*pi*fd1*tt）;
receive_fd2=receive2.*cos（2*pi*fd2*tt）;
receive_fd=0.3*receive_fd1+receive_fd2;
figure（3）;
plot（receive_fd）;
axis（[0?5e4?-1.5?1.5]）
title（‘兩個移動目標反射的回波‘）;xlabel（‘采樣點數‘）;ylabel（‘幅度‘）;

%匹配濾波
cankao=[ones（1T0*fs）?zeros（1（Tm-T0）*fs）];
cankao=fliplr（cankao）;
MF=conv（cankaoreceive_fd）;
figure（4）;
plot（MF）;
axis（[0?5e4?-32?32]）
title（‘匹配濾波結果‘）;xlabel（‘采樣點數‘）;ylabel（‘幅度‘）;

%數據重排
zz=1000-mod（length（MF）1000）;
MF1=[MF?zeros（1zz）];
for?m=1:1:（length（MF）/1000）
????for?n=1:1:1000
????????speed（mn）=MF1（11000*（m-1）+n）;
????end
end
real_speed=real（speed）;
figure（5）;
meshc（real_speed）
title（‘數據重排三維圖‘）;xlabel（‘距離門‘）;ylabel（‘速度門‘）;
speed_out=fft（speed）;
figure（6）;
meshc（abs（speed_out））
title（‘FFT處理‘）;xlabel（‘距離門‘）;ylabel（‘速度門‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1882??2010-01-28?19:44??pd_duopule.m

-----------?---------??----------?-----??----

?????????????????1882????????????????????1