%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%多徑電力線信道模型已經(jīng)建立了電力線信道的頻域響應。但是頻域到時間域。
%%下面是建立5徑、15徑電力線信道模型?
%Lohani?Satyendra?Nath.Communication?Network?Analysis?in?Smart?Grid?[D].
%???Linnaeus?University.2012（電力線典型信道）.。
%mqm?tested?it?on?Nov.32012
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%衰減參數(shù)?k=1a0=0a1=7.8*10^（-10）
%?多徑參數(shù)：
%?%N=15徑仿真信道參數(shù)
%?%衰減參數(shù)?k=1a0=0a1
%?%多徑參數(shù)：
%?N?為最大路徑數(shù)目;?
%?gi為路徑i?的多徑衰減系數(shù)
%?di為路徑長度
%?τi?為信號通過路徑i?的延%?時;?
%??a0、a1?衰減系數(shù)
%??k?是衰減因子的指?數(shù)，取0.5-1；
%?現(xiàn)仿真4?徑信道，令真空中的傳播速度p?=3×10^8?m/s?[季建玲2011在低壓電力線信道]，
%[張惠敏2008在中壓電力線信道]，傳播速度vp?=1.5×10^8?m/s
%?介電常數(shù)εr?=?3.8，信號的頻率范圍從0.01MHz?～?20?MHz，

%%%%%%%%%%%%%%%%低壓電力線信道相關參數(shù)
clear?all;
close?all;
clc
k=1;????a0=0;???a1=1.5*10^（-9）;
g1=[0.6-0.080.08-0.080.15];%5徑信道模型參數(shù)gi為衰減系數(shù)
d1=[100130160190?300];%5徑信道模型參數(shù)路徑長度di
g2=[0.0290.0430.103-0.058-0.045-0.0400.038-0.0380.071?-0.035?...
0.065?-0.0550.042-0.0590.049];%15徑
d2=[9010211314314820026032241149056774096011301250];%15徑
%設定頻帶寬度10kHz-20MHz
bandwidth=20e+6;%20MHz?
samplefrequency=5*1e+7;%采樣頻率50MHz
step_frequency=2048;%頻率步長，Hz
N_subcarrier=2048;
f=10000:step_frequency:bandwidth;
h1=zeros（1length（f））;
h2=zeros（1length（f））;
%----------------------
vp=3*10^8;%低壓電力線信道
for?i=1:1:5
????h1=h1+g1（i）*exp（-（a0+a1*f.^k*d1（i））-j*2*pi*f*2*d1（i）/vp）;
end;
for?i=1:1:15
????h2=h2+g2（i）*exp（-（a0+a1*f.^k*d2（i））-j*2*pi*f*2*d2（i）/vp）;
end;

%-----低壓電力線信道畫圖-------------
figure;
subplot（321）;
plot（f/100000020*log10（abs（h1）））;
axis（[0?bandwidth/1e+6?-50?0]）;
xlabel（‘低壓電力線信道，頻率（MHz）‘）;
ylabel（‘|H（f）|（dB）‘）;
title（‘5徑電力線信道傳輸函數(shù)幅度‘）;
subplot（322）;
plot（f/100000020*log10（abs（h2）））;
axis（[0?bandwidth/1e+6?-50?0]）;
xlabel（‘低壓電力線信道，頻率（MHz）‘）;
ylabel（‘|H（f）|（dB）‘）;
title（‘15徑電力線信道傳輸函數(shù)幅度‘）;
subplot（323）;
plot（f/100000angle（h1）/pi）;
axis（[0?bandwidth/1e+6?-1.5?1.5]）;
xlabel（‘頻率（MHz）‘）;
ylabel（‘相位‘）;
title（‘5徑電力線信道傳輸函數(shù)相位‘）;
subplot（324）;
plot（f/100000angle（h2）/pi）;
axis（[0?bandwidth/1e+6?-1.5?1.5]）;
xlabel（‘頻率（MHz）‘）;
ylabel（‘相位‘）;
title（‘15徑電力線信道傳輸函數(shù)相位‘）;

%?d3=d1*1/vp*1e+6;%us
%?d4=d2*1/vp*1e+6;%us
h1_time=ifft（h1）;
subplot（325）;
t1=1/samplefrequency*[1:length（h1_time）];%us
h1_main=abs（h1_time）;
Length_show=128/2;
plot（1e+6*t1（1:Length_show）h1_main（1:Length_show））;
%?axis（[0?10?-1.5?1.5]）;
xlabel（‘時間（us）‘）;
ylabel（‘幅度‘）;
title（‘5徑電力線信道沖擊響應響應‘）;
h2_time=ifft（h2）;
h2_main=abs（h2_time）;
subplot（326）;
plot（1e+6*t1（1:Length_show）h2_main（1:Length_show））;
%?axis（[0?10?-0.1?0.1]）;
xlabel（‘時間（us）‘）;
ylabel（‘幅度‘）;
title（‘15徑電力線信道沖擊響應響應‘）;
%
%%%%%%%%%%%%%%%%中壓電力線信道相關參數(shù)
%----------------------
vp=1.5*10^8;%中壓電力線信道
for?i=1:1:4
????h1=h1+g1（i）*exp（-（a0+a1*f.^k*d1（i））-j*2*pi*f*2*d1（i）/vp）;
end;
for?i=1:1:15
????h2=h2+g2（i）*exp（-（a0+a1*f.^k*d2（i））-j*2*pi*f*2*d2（i）/vp）;
end;

%-----中壓電力線信道畫圖-------------
figure;
subplot（321）;
plot（f/100000020*log10（abs（h1）））;
axis（[0?bandwidth/1e+6?-5