Ma'lum bir joyda astronomik tadqiqotlar o`tkazish imkoniyatini baqolash

 

 

 

 

 

 

 

 

 

 

 

---

62 

 

§ 3.2 “k” – koeffisientdi hisoblab  topishda  “Matlab” dasturidan  

foydalanish 

Matlab –so`zi matrsalar labaratoryasi degan ma`noni anglatadi. Bu 

dastur yordamida ko`pgina amallarni bajarish mumkin. Masalan: grafik 

chizish, hisob-kitob ishlari, matematika, fizikada,astronomiya va h.k 

ishlarni bajarish mumkin.

 

%Oy  Kun  Булутланиш - кундузи 

Булутланиш - кечаси 

load gismeteo2002.txt 

oy=gismeteo2002(:,1); 

kun=gismeteo2002(:,2); 

kunbulut=gismeteo2002(:,3); 

kechbulut=gismeteo2002(:,4); 

dt=datenum(2002,oy,kun); 

plot(dt,kunbulut) 

datetick('x') 

figure 

plot(dt,kechbulut) 

datetick('x') 

load gidmeteo_LT_behato2012 

%gives kvaqt bulut 

%datestr(kvaqt(1))      ans = 01-Jan-2002 05:00:00 

%datestr(kvaqt(end))    ans = 01-Nov-2005 02:00:00 

i=find(kvaqt

gmdt=kvaqt(i); 

gmblt=bulut(i); 

figure 

plot(gmdt,gmblt) 

datetick('x') 

%!!!!!!!!!!!!!!!!!!!!! 

%2002 uchun qorong'i payt 

%!!!!!!!!!!!!!!!!!!!!! 

load gidmeteo_LT_behato2012 %was gidmeteo_LT 

% datestr(kvaqt(1)) 

% bulut(1) 

v=kvaqt;    % was v=kvaqt(2:end); 

b=bulut;    % was b=bulut(2:end);nkr=14; %nima bu? 

db=fix(v(1));   %01-Jan-2002 

de=fix(v(end)); %01-Nov-2005 

[twssh,twssm,twsrh,twsrm]=rdtw02; 

twss=twssh+twssm./60; 

twsr=twsrh+twsrm./60; 

mon=[]; 

for k=1:365 % was twss 

    ibs=[]; 

---

63 

 

    kbs=[]; 

    thisday(k)=db+k-1; 

    kbs=find(v>twss(k)/24+thisday(k) & 

v

    if isempty(kbs) 

        sobl(k)=-1; 

    else if length(kbs)>1 

            %         sobl(k)=sum(b(kbs))./length(kbs); 

            %kbs 

            ibs=[kbs(1)-1 kbs kbs(end)+1]; 

            y=b(ibs); 

            x=(v(ibs)-fix(v(ibs(1)))-0.5)*24; %0.5d=12 soat 

ayrildi qulaylik uchun 

            xi=twss(k)-12 :0.25: twsr(k)+12;    %bu yerdayam 12 

soat ayrildi 

            yi = interp1(x,y,xi); 

            srobl(k)=sum(yi)/length(xi); 

            den(k)=fix(v(ibs(1))); 

            %2013                         

        end 

    end 

end 

figure 

plot(den,srobl,'.') 

title('Nochnaya srednyaya oblachnost') 

xlabel('Oylar') 

ylabel('o''rtacha bulutlanish ballarda') 

datetick('x',5) 

%plot(thisday,sobl) 

hold on  

plot(dt,kechbulut,'r.') 

legend('Gid','Gis') 

dt=den; 

sobl=srobl; 

% dt=v; 

% sobl=b; 

mnlength=[31 28 31 30 31 30 31 31 30 31 30 31]; 

 

for m=1:12 

    mobl=[]; 

    im=[]; 

    for y=1:1 %y=1:2005-2002+1 

        mbeg=datenum(y+2002-1,m,1,1,0,0); 

        mend=datenum(y+2002-1,m+1,1,1,0,0); 

        iii=find(dt>mbeg&dt

        im=[im iii]; 

    end 

    mobl=sobl(im); 

    i0=find(mobl<0); 

    if ~isempty(i0) 

        mobl(i0)=[]; 

    end 

    iy=find(mobl<=3);  %nkr ? 

---

64 

 

    kolyasnoch(m)=length(iy)/mnlength(m); 

end 

gidmkyanoch=kolyasnoch 

sobl=kechbulut; 

for m=1:12 

    mobl=[]; 

    im=[]; 

    for y=1:1 %y=1:2005-2002+1 

        mbeg=datenum(y+2002-1,m,1,1,0,0); 

        mend=datenum(y+2002-1,m+1,1,1,0,0); 

        iii=find(dt>mbeg&dt

        im=[im iii]; 

    end 

    mobl=sobl(im); 

    i0=find(mobl<0); 

    if ~isempty(i0) 

        mobl(i0)=[]; 

    end 

    iy=find(mobl<=3);  %nkr ? 

    kolyasnoch(m)=length(iy)/mnlength(m); 

end 

gismkyanoch=kolyasnoch; 

figure 

plot(gidmkyanoch.*100,'Linewidth',3)                 

hold on 

plot(gismkyanoch.*100,'r--','Linewidth',3) 

title(' 2002 yilda ochiq vaqt miqdori') 

xlabel('Oylar') 

ylabel('Ochiq vaqt miqdori') 

legend('UzGidroMet','GisMeteo') 

grid on 

gismkyanoch*100' 

gidmkyanoch*100' 

k=gidmkyanoch./gismkyanoch 

save k2002 k gismkyanoch gidmkyanoch 

 

 

 

 

 

 

 

 

 

 

---

65 

 

Bu  ma`lumotlardan  foydalanib  gismeteo  ma`lumotlarini  yuqoridagi 

pragramma asosida hisoblab chiqdik va quydagicha ma`lumotlarga ega bo`ldik.  

UzGidMet (ko'k) 

GisMet (qizil) 

k

-koeffisient 

 

 

 

0.19 

0.39 

0.50 

0.07 

0.25 

0.29 

0.32 

0.55 

0.59 

0.33 

0.33 

1.00 

0.52 

0.61 

0.84 

0.87 

0.80 

1.08 

0.87 

0.87 

1.00 

0.94 

0.87 

1.07 

0.97 

0.97 

1.00 

0.77 

0.81 

0.96 

0.47 

0.50 

0.93 

0.19 

0.32 

0.60 

                                        7-jadval 

Bu jadvalda UzGidMet va GisMet ma`lumotlar asosida 

k-

koeffisient aniqlandi. Bu 

koeffisient  asosida  ixtiyoriy  yil  uchun  qo`llab  ko`rish  mumkin  bo`ladi.  Biz 

malumotlarni solishtirib quydagi natijani oldik 

 

34-rasm. 2002 yil uchun Gismeteo va Uzgidromet ma`lumotlarini solishtirish 

---

66 

 

Tahlil  natijasi  shuni  ko’rsatdiki,  aprel-noyabr  oylarida  bu  ikki  manbadan 

olingan  ma`lumotlar  bir-biriga  juda  yaqin  qiymatlarga  ega  bo’ladi.  Qish  oylarida 

esa  Gismeteo  ma`lumotlari  asosida  baholangan  ochiq  tungi  vaqt  miqdori  kattaroq 

qiymatlarga ega. Bu nomutanosiblik k-koeffisient orqali bartaraf etililishi mumkin. 

Bu  degani,  endilikda  ihtiyoriy  yil  yoki  vaqt  oralig`i  uchun  Gismeteodan  olingan 

ma’lumotlarni  ushbu  ishda  olingan  koeffisientlar  orqali  haqiqatga  yaqin 

qiymatlarga  o’tkazish  mumkin  va  ular  Samarqand  va  boshqa  astronomik  puntlar 

astroiqlimining  asosiy  parametri  –  ochiq  vaqt  miqdorini  katta  statistik  vazn  bilan 

o’lchash imkonini yaratadi.