ks_twosample.m

program ini untuk membuat matrik dinamis

%============== two-sample Kolmogorov-Smirnov Test==========
clear all; clc; format long
load  redaman_disdroSST2 diskum5  % load curah hujan disdro link B-T
load kolmogorov_R_raingauge RA RB RD % load curah hujan raungauge

A=diskum5;      %sampel 1
A=nonzeros(A);
B=RD;           %sampel 2
B=nonzeros(B); 
m=mean(A);s=std(A);v=var(A);a=sort(A);
mb=mean(B);sb=std(B);vb=var(B);b=sort(B);
G=log(A);Anor=sort(G);% data diurutkan
Gb=log(B);Bnor=sort(Gb);


% A2
% A3
% A4
% A5
% A6

na=length(Anor); %jmlh data sample
nb=length(Bnor);

for i=1:length(Anor);
    F1(i)=i./na;
end

for j=1:length(Bnor);
    F2(j)=j./nb;
end

%%%%%%%%%%%%%%%%%%%%%BuaT RanGe%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


   
% for i=1:length(Anor)
%    u(i)=Anor(i)
%    Fu(i)=normcdf(u(i)); %pembangkitan fungsi lognormal untuk teori
% end
% 
% for i=1:length(Anor)
% Dplus(i)=i./na-Fu(i);
% Dmin(i)=Fu(i)-(i-1)./na;
% end
% 
% Da=[Dplus' Dmin'];
% D=max(Da); %Nilai D-KS
% 
% %MenenTukan Nilai D taBeL Berdasarkan Derajat Kepercayaan
% % d1=1.07/sqrt(na);% 0.8 derajat kepercayaan
% % d1=1.14/sqrt(na);%0.85
% % d1=1.22/sqrt(na);%0.9
% d1=1.36/sqrt(na);%0.95
% % d1=1.63/sqrt(na);%0.99
% 
nna=0:1/(na-1):1;
nnb=0:1/(nb-1):1;
% teori11=Fu+d1; %Batas Atas
% teori12=Fu-d1; %Batas BawaH

figure(1)%gambar CDF
semilogy(a,1-nna','k--',b,1-nnb','r-.','linewidth',2)
% legend('Teoritis SST[1]','Teoritis SST[2]','Empiris SST link 4km  ')
% xlabel('Redaman (dB)')
% % xlabel('Curah Hujan')
% ylabel('Prob.[Redaman < absis]')
% % ylabel('Prob.[Curah Hujan < absis]')
% title({'Kolmogorov-Smirnov fitting lognormaL ';'Link B-T'})
grid on;
% 
% figure(2)% gambar CCDF
% semilogy(a,1-teori12','k--',a,1-teori11','k-.',a,1-nna,'k-','linewidth',2)
% legend('Teoritis SST[1]','Teoritis SST[2]','Empiris SST link 4km ')
% xlabel('Redaman (dB)')
% % xlabel('Curah Hujan')
% ylabel('Prob.[Redaman > absis]')
% % ylabel('Prob.[Curah Hujan > absis]')
% title({'Kolmogorov-Smirnov fitting lognormaL ';'Link B-T'})
% grid on;