MIMO雷达OFDM Chirp波形设计

产生OFDM Chirp波形，包括时域表示和频域表示，用于mimo雷达，且附有参考文献，里面有波形的时频域表示，以便进一步研究

% mimo sar ofdm chirp waveform diversity design with random matrix modulation 
% Wen Qin Wang
clear all
clc
%%  波形参数  %%
deltaf = 50e6;%50Mhz
Tb = 1.0e-6;  %1us
k = 50e12;    %调频率
M = 8;        %列，8个chip
N = 8;        %行，8个subcarrier
Tp = 8*Tb;    %脉冲时宽
B = 8*deltaf;    %脉冲带宽 400e6
fs = 400e6; 
n = round（Tp*fs）;%采样点个数（信号点数）
t = linspace（（Tp-0）/nTpn）;
j = sqrt（-1）;
f18 = 7*deltaf;
f26 = 5*deltaf;
f31 = 0*deltaf;
f43 = 2*deltaf;
f55 = 4*deltaf;
f64 = 3*deltaf;
f72 = 1*deltaf;
f87 = 6*deltaf;

%% 定义每个码片内信号的幅度   （截取码片对应的数据）
a = ones（1length（t））;
% a1
for i = 1:length（t）/8
    a1（i） = a（i）;
end
for i = length（t）/8+1:length（t）
    a1（i） = 0;
end

%a2
for i = 1:length（t）/8
    a2（i） = 0;
end
for i = length（t）/8+1:2*length（t）/8
    a2（i） = a（i）;
end
for i = 2*length（t）/8+1:length（t）
    a2（i） = 0;
end

%a3
for i = 1:2*length（t）/8
    a3（i） = 0;
end
for i = 2*length（t）/8+1:3*length（t）/8
    a3（i） = a（i）;
end
for i = 3*length（t）/8+1:length（t）
    a3（i） = 0;
end

%a4
for i = 1:3*length（t）/8
    a4（i） = 0;
end
for i = 3*length（t）/8+1:4*length（t）/8
    a4（i） = a（i）;
end
for i = 4*length（t）/8+1:length（t）
    a4（i） = 0;
end

%a5
for i = 1:4*length（t）/8
    a5（i） = 0;
end
for i = 4*length（t）/8+1:5*length（t）/8
    a5（i） = a（i）;
end
for i = 5*length（t）/8+1:length（t）
    a5（i） = 0;
end

%a6
for i = 1:5*length（t）/8
    a6（i） = 0;
end
for i = 5*length（t）/8+1:6*length（t）/8
    a6（i） = a（i）;
end
for i = 6*length（t）/8+1:length（t）
    a6（i） = 0;
end

%a7
for i = 1:6*length（t）/8
    a7（i） = 0;
end
for i = 6*length（t）/8+1:7*length（t）/8
    a7（i） = a（i）;
end
for i = 7*length（t）/8+1:length（t）
    a7（i） = 0;
end

%a8
for i = 1:7*length（t）/8
    a8（i） = 0;
end
for i = 7*length（t）/8+1:8*length（t）/8
    a8（i） = a（i）;
end

%% 重新定义t1 t2 t3 t4 t5 t6 t7 t8   %%【正确】
t1 = t;t2 = t-Tb;t3 = t-2*Tb;t4 = t-3*Tb;t5 = t-4*Tb;t6 = t-5*Tb;t7 = t-6*Tb;t8 = t-7*Tb;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%    
phi18 = a1.*exp（j*2*pi*（f18*t1+0.5*（k）*t1.^2））;
phi26 = a2.*exp（j*2*pi*（f26*t2+0.5*（k）*t2.^2））;
phi31 = a3.*exp（j*2*pi*（f31*t3+0.5*（k）*t3.^2））;
phi43 = a4.*exp（j*2*pi*（f43*t4+0.5*（k）*t4.^2））;
phi55 = a5.*exp（j*2*pi*（f55*t5+0.5*（k）*t5.^2））;
phi64 = a6.*exp（j*2*pi*（f64*t6+0.5*（k）*t6.^2））;
phi72 = a7.*exp（j*2*pi*（f72*t7+0.5*（k）*t7.^2））;
phi87 = a8.*exp（j*2*pi*（f87*t8+0.5*（k）*t8.^2））;

%%  时域表示  %%
s1 = phi18+phi26+phi31+phi43+phi55+phi64+phi72+phi87;
s1 = s1/max（abs（s1））;
figure
plot（ts1）
set（gca‘ytick‘-1:0.5:1）
set（gca‘xtick‘0:2.0e-6:8.0e-6）
xlabel（‘time[second]‘）
ylabel（‘amplitude‘）
title（‘waveform a1‘）

%%  频域表示  %%
fs1 = fftshift（fft（s1））;
fs1 = abs（fs1）;
fs1 = fs1/max（fs1）;
fs1 = 20*log10（fs1）;
fx = linspace（（B-0）/nBn）;
figure
plot（fxfs1）
set（gca‘ytick‘-50:10:0）
set（gca‘xtick‘0:100e6:400e6）
axis（[0 400e6 -50 0]）
xlabel（‘frequency[Hz]‘）
ylabel（‘normalized modulus[dB]‘）
title（‘waveform a1‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    2666302  2017-12-20 18:26  mimo radar OFDM Chirp waveform design\MIMO SAR OFDM Chirp Waveform Diversity Design With Random Matrix Modulation.pdf

     文件       3036  2017-12-20 12:12  mimo radar OFDM Chirp waveform design\ofdmchirpa1.m

     文件       3017  2017-12-20 12:13  mimo radar OFDM Chirp waveform design\ofdmchirpa2.m

     目录          0  2017-12-27 12:20  mimo radar OFDM Chirp waveform design

----------- ---------  ---------- -----  ----

              2672355                    4