for i=1:6 %vong lap quet cac toc do bit
for j=1:73 %vong lap quet cac buoc song tu 830 -> 1550
t=1;
P_sensitivity(i,j)=A*rate(i)*10^6*6.625*10^(34)*3* 10^8/(lamda(j)*10^(-9)); %tinh P_sensivity(theo buoc song "lamda" va toc do "rate")
P_sensitivity(i,j)= 10*log10(P_sensitivity(i,j)) +30 ; %dBm
beta(j)=exp( (3.912*1000/visibility) * ((lamda(j)/550)^(-q ))*range );%suy hao tan xa chung
loss(j)=1.076*(rain^0.67)*range +10*log10(beta(j));%tinh suy hao tong(theo buoc song lamda)
for k=1:31 %vong lap quet cac gia tri cong suat dung duoc
P_receive(i,j,k)=P_transmit(k) - loss(j); %tinh P_receive theo cong suat phat va buoc song
if(P_receive(i,j,k) > P_sensitivity(i,j) +3)
SNR(i,j,k)=n*lamda(j)*10^(-9)*(10^(P_receive(i,j,k)/10)*10^(-3)/((rate(i)*10^3)*6.625*10^(-34)*3*10^8));
if (SNR(i,j,k)>0)
SNR(i,j,k)=10*log10(SNR(i,j,k));%ti so tin hieu tren nhieu (dB)
if (SNR(i,j,k)>0)
BER(i,j,k)=0.5*erfc(0.5*(SNR(i,j,k)/2)^0.5);%Bit Error Ratio
else
BER(i,j,k)=1;
end
end
if ((BER(i,j,k)<=10^(-3))&&(P_receive(i,j,k)>P_sensitivity(i,j)+3))%3 la do du tru cong suat(dBm)) %kiem tra thoa BER va so voi BERmin de tim to hop toi uu
P_transmit_optimize(i)=P_transmit(k);
lamda_optimize(i)=lamda(j);
index(i,j,k)=1;
if ((BER(i,j,k)<=10^(-3)) && t==1)
P_draw(i,j)=P_transmit(k);
t=0;
end
end
end
end
end
end
for j=1:73 %vong lap quet cac buoc song tu 830 -> 1550
t=1;
P_sensitivity(i,j)=A*rate(i)*10^6*6.625*10^(34)*3* 10^8/(lamda(j)*10^(-9)); %tinh P_sensivity(theo buoc song "lamda" va toc do "rate")
P_sensitivity(i,j)= 10*log10(P_sensitivity(i,j)) +30 ; %dBm
beta(j)=exp( (3.912*1000/visibility) * ((lamda(j)/550)^(-q ))*range );%suy hao tan xa chung
loss(j)=1.076*(rain^0.67)*range +10*log10(beta(j));%tinh suy hao tong(theo buoc song lamda)
for k=1:31 %vong lap quet cac gia tri cong suat dung duoc
P_receive(i,j,k)=P_transmit(k) - loss(j); %tinh P_receive theo cong suat phat va buoc song
if(P_receive(i,j,k) > P_sensitivity(i,j) +3)
SNR(i,j,k)=n*lamda(j)*10^(-9)*(10^(P_receive(i,j,k)/10)*10^(-3)/((rate(i)*10^3)*6.625*10^(-34)*3*10^8));
if (SNR(i,j,k)>0)
SNR(i,j,k)=10*log10(SNR(i,j,k));%ti so tin hieu tren nhieu (dB)
if (SNR(i,j,k)>0)
BER(i,j,k)=0.5*erfc(0.5*(SNR(i,j,k)/2)^0.5);%Bit Error Ratio
else
BER(i,j,k)=1;
end
end
if ((BER(i,j,k)<=10^(-3))&&(P_receive(i,j,k)>P_sensitivity(i,j)+3))%3 la do du tru cong suat(dBm)) %kiem tra thoa BER va so voi BERmin de tim to hop toi uu
P_transmit_optimize(i)=P_transmit(k);
lamda_optimize(i)=lamda(j);
index(i,j,k)=1;
if ((BER(i,j,k)<=10^(-3)) && t==1)
P_draw(i,j)=P_transmit(k);
t=0;
end
end
end
end
end
end
Leave a comment: