clear;clc
%?
%?功能：使用成型濾波法仿真Clarke模型
%?作者：hhuyudi
%?參考：MIMO-OFDM?Wireless?Communications?with?MATLAB
%????????參考書中?FWGN?小節?
%
%最大多普勒頻為926Hz，仿真次數為50us，仿真次數為十萬次
fd=926;?????%?最大Doppler
ts=5e-5;????%?仿真最小時間
fs=1/ts;????%?sample?frequency
Ns=1e5;?????%?仿真次數
M=2^12;?????%?做FFT的總點數，可以使用自相關函數通過FFT計算功率譜密度，請自行完成
t=（0:Ns-1）*ts;%?仿真總時間，作圖使用
tt=（0:M-1）*ts;%?仿真總時間，作圖使用
ff=（-M/2:M/2-1）/（M*ts*fd）;%?頻率歸一化，功率譜密度作圖使用
temp=zeros（3Ns）;?%?3*Ns的矩陣
??????????????????%?temp（1，：）存儲仿真結果的自相關函數
??????????????????%?temp（2，：）存儲ones（1，Ns）的自相關函數，用于歸一化
??????????????????%?temp（3，：）存儲仿真結果其分量的互相關函數
??????????????????
%generat?channel?information
[hNfftNifftDoppler_coeff]=Clarke_model（fdfsNs）;?%?歸一化的Clarke模型，功率為1

%自相關函數，分量的互相關函數，為了結果更準確，這里50次仿真取平均值
for?ii=1:50
????[hNfftNifftDoppler_coeff]=Clarke_model（fdfsNs）;
????yy=xcorr（h）;
????yy_cs=xcorr（real（h）imag（h））;
????temp（1:）=temp（1:）+yy（Ns:length（yy））;
????temp（3:）=temp（3:）+yy_cs（Ns:length（yy_cs））;
end

%畫圖
%Clarke模型復包絡的幅度增益、幅度分布、相位分布
figure（1）;
subplot（211）
plot（（1:Ns）*ts10*log10（abs（h）））;
axis（[0?0.01?-20?10]）;
xlabel（‘時間/s‘）;ylabel（‘幅度/db‘）;title（‘信道時域特性‘）;
str=sprintf（‘channel?model?by?Clarke?with?fm=%dHzTs=%dus‘fdts*1e6）;
title（str）;
%?幅度分布
subplot（223）
[fxi]=ksdensity（abs（h））;
plot（xif）;
cs2=var（h）/2;??%?瑞利分布的方差是實部或虛部的方差，是h方差的一半
r=linspace（021000）;
fx2=r./（cs2）.*exp（-r.^2/2/（cs2））;
hold?on;plot（rfx2‘r:‘）;hold?off;
xlabel（‘幅度‘）;ylabel（‘統計次數‘）;title（‘幅度分布‘）;axis（[0?2.5?0?1.0]）;
%?相位分布
subplot（224）
[fxi]=ksdensity（angle（h））;
plot（xif）;
hold?on;plot（[-pi?pi][1/2/pi?1/2/pi]‘r:‘）;hold?off;
xlabel（‘相位/rad‘）;ylabel（‘統計次數‘）;title（‘相位分布‘）;axis（[-pi?pi?0?0.2]）;


%%自相關函數和互相關函數計算，功率譜密度請自行編程
temp（1:）=temp（1:）/50;
temp（3:）=temp（3:）/50;
%用于歸一化
yyy=xcorr（ones（1Ns））;
temp（2:）=yyy（Ns:length（yy））;
%自相關函數和分量互相關函數，取M個點，進行歸一化
for?k=1:M
????simulated_corr（k）=real（temp（1k））/temp（2k）;
????simulated_corr_cs（k）=real（temp（3k））/temp（2k）;
end
classical_corr=besselj（02*pi*fd*tt）;

%%畫圖：自相關函數和互相關函數
figure（2）;
subplot（211）
plot（ttsimulated_corr‘k-‘）;
title（‘自相關函數‘）;?grid?on;?
xlabel（‘時間差/s‘）;ylabel（‘相關系數‘）;axis（[0?0.004?-0.5?1]）;
legend（‘仿真結果‘）;
subplot（212）
plot（ttsimulated_corr_cs‘k-‘）;
title（‘互相關函數‘）;?grid?on;?
xlabel（‘時間差/s‘）;ylabel（‘相關系數‘）;axis（[0?0.004?-0.2?0.2]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-06-09?02:00??Clarke_Jakes_Zheng模型matlab仿真\
?????目錄???????????0??2018-06-09?02:01??Clarke_Jakes_Zheng模型matlab仿真\Clarke\
?????文件????????2724??2018-06-09?01:52??Clarke_Jakes_Zheng模型matlab仿真\Clarke\Clarke.m
?????文件????????1126??2018-06-09?01:45??Clarke_Jakes_Zheng模型matlab仿真\Clarke\Clarke_model.m
?????文件?????????485??2018-06-09?01:41??Clarke_Jakes_Zheng模型matlab仿真\Clarke\Doppler.m
?????目錄???????????0??2018-06-09?01:16??Clarke_Jakes_Zheng模型matlab仿真\Jakes\
?????文件????????2017??2018-06-09?01:57??Clarke_Jakes_Zheng模型matlab仿真\Jakes\Jakes.m
?????文件?????????797??2018-05-17?11:24??Clarke_Jakes_Zheng模型matlab仿真\Jakes\Jakes_model.m
?????文件?????????507??2018-06-09?01:59??Clarke_Jakes_Zheng模型matlab仿真\read_me.txt
?????目錄???????????0??2018-06-09?01:18??Clarke_Jakes_Zheng模型matlab仿真\改進Jakes\
?????文件????????1885??2018-06-09?01:58??Clarke_Jakes_Zheng模型matlab仿真\改進Jakes\Zheng.m
?????文件????????1309??2018-05-25?10:56??Clarke_Jakes_Zheng模型matlab仿真\改進Jakes\Zheng_QPSK.m
?????文件?????????510??2018-06-01?20:08??Clarke_Jakes_Zheng模型matlab仿真\改進Jakes\z_rayleigh.m