%?mimo?sar?ofdm?chirp?waveform?diversity?design?with?random?matrix?modulation?
%?Wen?Qin?Wang
clear?all
clc
%%??波形參數(shù)??%%
deltaf?=?50e6;%50Mhz
Tb?=?1.0e-6;??%1us
k?=?50e12;????%調(diào)頻率
M?=?8;????????%列，8個chip
N?=?8;????????%行，8個subcarrier
Tp?=?8*Tb;????%脈沖時寬
B?=?8*deltaf;????%脈沖帶寬?400e6
fs?=?400e6;?
n?=?round（Tp*fs）;%采樣點(diǎn)個數(shù)（信號點(diǎn)數(shù)）
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;

%%?定義每個碼片內(nèi)信號的幅度???（截取碼片對應(yīng)的數(shù)據(jù)）
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（‘a(chǎn)mplitude‘）
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