1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

% 脈沖雷達仿真

clc;closeall;clearall;

%%

Fs = 125e6;%數據采樣率125M SPS

beta= 20e6;%接收中頻帶寬，20MHz

tao = 10e-6;%脈沖時長10uS

%%

c = 3e8;%光速，單位米

l = 500;%目標距離，單位米

t_d = l/c;%回波延遲，單位秒

a_t = 0.088;%發(fā)射功率系數

a_r = 1e-6;%信道衰減系數

a_n = 5e-6;%噪聲增益

%%

N = tao * Fs;%信號采樣點數

t = 0 : tao/(N-1) : tao;%時間序列

%發(fā)射信號

s_t = a_t*cos(pi*(beta/tao)* t.*t);

%噪聲信號

s_n = a_n*a_t*rand(1,N);

%接收信號

s_r0 = a_r*a_t*cos(pi*(beta/tao)* (t + t_d).*(t + t_d));%衰減部分

s_r = s_r0 + s_n;%實際波形

%相干接收

rr = xcorr(s_t,s_r);

%檢測峰值并標記

[M I] =max(rr);

%計算延遲

delta_t = (I - N) / Fs;

%換算距離

l_ra = delta_t * c;

%計算誤差

err = (l - l_ra)/l * 100;

str =sprintf('實際距離%d米,測算距離%d米,誤差%f%%',l,l_ra,err);

disp(str);

%發(fā)射信號內積

p_t =dot(s_t,s_t);

%接收信號內積

p_r =dot(s_r0,s_r0);

%噪聲內積

p_n =dot(s_n,s_n);

%信噪比

snr = db(p_r/p_n,'power');

str =sprintf('發(fā)射功率 %e，接收功率 %e，噪聲功率 %e,信噪比 %fdB',p_t,p_r,p_n,snr);

disp(str);

%%

figure;

subplot(311);

plot(t,s_t);title('發(fā)射波形');

subplot(312);

plot(t,s_r);title('接收波形');

subplot(313);

plot(rr);title('相干接收');

text(I,M,'leftarrow 回波位置','Color','red');