%頻率選擇性rayleigh信道測試
num=5;
%假設功率延遲譜服從負指數分布~exp(-t/trms),trms=(1/4)*cp時長；
%t在0~cp時長上均勻分布
%若cp時長為16e-6s，可以取5徑延遲如下
delay=[0 2e-6 4e-6 8e-6 12e-6];
trms=4e-6;
var_pow=10*log10(exp(-delay/trms));

fd=132;
t_interval=1e-6;counter=200000;count_begin=0;
t_shift=floor(delay/t_interval);%歸一化各徑延時
%theta_shift=2*pi*fc*delay;

chann_l=100000;%信道采樣點數，若一個調制符號采樣一個信道點，則信道采樣點數等于輸入信號中的調制符號個數
selec_ray_chan=zeros(num,chann_l);%初始化頻率選擇性信道，徑數＝num
pow_per_channel=10.^(var_pow/10);%各徑功率線性化，從dB轉變成線性
total_pow_allchan=sum(pow_per_channel);%各徑功率之和
%以下for循環產生相互獨立的num條rayleigh信道
for k=1:num
    atts=sqrt(pow_per_channel(k));
    selec_ray_chan(k,:)=atts*rayleighnew(chann_l,t_interval,fd,count_begin+k*counter)/sqrt(total_pow_allchan);
end
k=1:chann_l;
%plot(k,abs(selec_ray_chan(1,:)),'r',k,abs(selec_ray_chan(2,:)),'g',k,abs(selec_ray_chan(3,:)),'b',k,abs(selec_ray_chan(4,:)),'k',k,abs(selec_ray_chan(5,:)),'m');
figure(1)
plot(k,abs(selec_ray_chan(1,:)));
figure(2)
plot(k,abs(selec_ray_chan(2,:)));
figure(3)
plot(k,abs(selec_ray_chan(3,:)));
figure(4)
plot(k,abs(selec_ray_chan(4,:)));
figure(5)
plot(k,abs(selec_ray_chan(5,:)));