Drukuj

Zaloguj się lub zarejestruj.
1 Godzina 1 Dzień 1 Tydzień 1 Miesiąc Zawsze
Zaloguj się podając nazwę użytkownika, hasło i długość sesji

[greg]

Zrobiłem wreszcie ten matematyczny model psychometryczny. Arkusz jest do pobrania z serwera forum. Jest co prawda w formacie xls (Excel), ale i tak został przez mnie zrobiony w Open Office i ten darmowy pakiet biurowy z powodzeniem go odczyta. Prawdę mówiąc spodziewałem się, że wyjdzie trochę lepiej. O ile współczynnik kierunkowy prostej dało się dość gładko przybliżyć za pomocą wielomianu drugiego stopnia to problem był z wyrazem wolnym. Tu dopiero pomógł wielomian piątego stopnia, a i efekt dopasowania był gorszy (r²=0.99 dla współczynnika kierunkowego i r²=0.91 dla wyrazu wolnego). Model został sporządzony dla zakresu temperatur od 0 C do 12 C i w tym właśnie zakresie jest użyteczny.

W czterech pierwszych komórkach jest kalkulator różnicy temperatury pomiędzy termometrem suchym i mokrym.

Komórki poniżej zawierają wartości tabeli wygenerowane przygotowanym modelem. Zrobiłem to po to byś mógł ocenić jego jakość. Wartości liczbowe w kolorze czarnym to temperatura w stopniach Celsjusza oraz wilgotność względna. Wartości na zielono to są wartości oryginalne z Twojej tabeli. Wartości na żółto są wyliczone moim modelem. Czyli właśnie zielone z żółtym trzeba porównywać.

Marcinie, oceń proszę czy z taką dokładnością to się będzie nadawało do dalszego wykorzystania w naszym przyszłym „kalkulatorze przymrozków”.

Dla bardzo dociekliwych. Formuła arkusza wyliczająca różnicę pomiędzy termometrem suchym a mokrym na podstawie wilgotności względnej i temperatury termometru suchego ma postać:

=(B2-(3,773*10^-4*B1^5-1,166*10^-2*B1^4+1,211*10^-1*B1^3-4,46*10^-1*B1^2+1,966*10^-1*B1+9,883*10^1))/(-2,805*10^-2*B1^2+9,369*10^-1*B1-18,29)

gdzie B2 to wilgotność względna, a B1 to temperatura w stopniach Celsjusza.

[Marcin D]

Greg jestem pod wrażeniem tego co zrobiłeś.
Co mi tu oceniać? Sam napisałeś, że dopasowanie nie jest idealne. Jeżeli nie jesteś wstanie lepiej tego zrobić to musi nam to wystarczyć.
Sprawdziłem kilka parametrów i na pierwszy rzut oka wydaje się w miarę ok.

Model został sporządzony dla zakresu temperatur od 0 C do 12 C i w tym właśnie zakresie jest użyteczny.

Dlaczego tylko do 12C?

Należy chyba zaznaczyć, że zabawę z wzorem traktujemy to jako ciekawostkę. Na pewno wiadomo już dziś, że wzór nie najlepiej sprawdza się w temperaturach poniżej 0C. Także określenie na jego podstawie czy będzie -10 czy -15C na pewno może być kłopotliwe.

Jestem natomiast niemal pewny, że wzór mógłby służyć także w okresie zimowym tylko należałoby jakoś zmodyfikować współczynnik C w mocno ujemnych temperaturach.

Nie mniej jednak jest to dosyć ciekawe bynajmniej dla mnie. Nie dawno zrobiłem sobie prostą formułę w excelu i po prostu podstawiam sobie wartość wilgotnego termometru manualnie. Nie ukrywam, że jest to dosyć męczące bo za każdym razem trzeba się grzebać w tablicach. Super by było gdybyś tak zrobił, że wystarczy wpisać tylko temp i wilgotność.

Błąd raczej nie przekracza 1-2C. Jednak w zimę już tak dobrze nie będzie chyba.

[greg]

Marcinie, formuła jest użyteczna od 0 C do 12 C, bo dla takiego właśnie zakresu została sporządzona. Nic nie stoi na przeszkodzie by zrobić policzyć formułę dla innych zakresów. Rachunkowo nie jest to trudne, lecz niestety trochę pracochłonne. Trzeba policzyć równania wszystkich prostych (różnica temperatury suchy/mokry jako funkcja wilgotności) a następnie znaleźć w miarę dokładną zależność pomiędzy współczynnikami równań a temperaturą termometru suchego).

Obawiam się, że taka empiryczna formuła musi być podzielona na zakresy, bo dalsze zwiększenie zakresu będzie skutkować obniżeniem dokładności i/lub jeszcze większą komplikacją wzoru.

Jakie zakresy temperatur uznałbyś za potrzebne?

[Marcin D]

Zakres?
od -20 do +35C hehe.

To może lepiej zainteresujmy się wx sim? Mniej czasu może zajmie konfiguracja programu jak obliczanie tych równań.

http://www.wxsim.com/

Z tego co wiem to Biedny Miś intensywnie już działa w tym kierunku. Program jest trochę drogi, ale na pewno bardzo ciekawy. Prognozy będą generowana na podstawie danych ze stacji meteo Biednego Misia co pozwoli na dokładniejsze prognozy, które będą uwzględniać jego lokalny klimat.

 

[greg]

Marcinie, program kosztuje. To ja wolę sam coś zrobić, a za zaoszczędzone pieniądze jakąś roślinę kupić.

Może inaczej. Idzie jesień i zima. Niech Marek spróbuje zabawy z programem a my spróbujmy zabawy z wzorami. To, że marudzę, że to żmudne to tym się nie przejmuj. W sumie to lubię robić takie rzeczy tylko nie zawsze mam tyle czasu by móc spokojnie nad tym usiąść.

Proponowałbym zatem zrobić coś takiego. Tej jesieni i zimy potraktujmy to jako narzędzie do obliczania możliwości wystąpienia przymrozków oraz temperatury minimalnej (bo do tego się to chyba sprowadza). Zatem do tego celu wystarczyłoby dorobić dół tej Twojej tabeli. Czyli na początek mielibyśmy zakres od 12 C do -17 C. Zobaczymy jak się to sprawdzi. Jak będzie dobrze to w przyszłym roku możemy i górę tabeli dorobić.

Masz na pewno sporo mądrych książek na temat meteorologii. Daj mi tylko algorytmy na prognozy czego tam chcesz i możemy nasze narzędzie dowolnie rozbudowywać. Jakieś wykresy, obrazki słoneczek, chmurek itp. Idzie zima, wieczory coraz dłuższe, na ogród się nie wychodzi. Będzie czas na zabawę. No chyba, że wolisz kupić program z USA zamiast samemu zostać współtwórcą programu.

[Marcin D]

Program i tak kupie a we wzory swoją drogą można się bawić. Natomiast zakres od -17 do +12C trochę mi nie pasuje. Ja bym zrobił zakres od -10 do +20 chociaż. Dlaczego? Dlatego, że jest to program głównie do przewidywania  przymrozków, który potrzebuje dane z 12UTC. (czyli 13oo czasu zimowego i 14oo czasu letniego. Przymrozki najbardziej niespodziewane są wtedy gdy dzień jest ciepły (np 20C). Nawet późna jesienią często zdarza się, że o 12UTC jest w okolicach 20C natomiast podczas bardzo srogiej zimy rzadko zdarza się by było o 12UTC -17C.

[andres]

Przeniosłem z tematu pogodowego:

Przymrozki często są niespodziewane, ale chyba warto obserwować pogodę i brać pod uwagę różne czynniki:
- pora wystąpienia (pierwsze przymrozki u mnie nie wcześniej niż ok. połowy października, ale czasami dopiero w listopadzie)
- zapowiedzi meteo (nadciąga wyż)
- czyste niebo
- jeżeli przymrozek ma być solidny to już wieczorem czuć, że tak będzie
- położenie danego miejsca (osłona drzew, wystawa, itp).

Na wiosnę podobnie, ale wielką niewiadomą jest to, czy np. po ostatnich przymrozkach w marcu (u mnie ok. 24 - 28 marca) nie pojawią się przymrozki w maju (tak jak w 2011 r.) (w 2009 i 2010 r. nie było u mnie przymrozków po 24 - 28 marca).

Rośliny nieodporne na najlżejszy mróz i tak trzeba wnieść zanim dojdzie do zera, bo one i tak nie lubią chłodu,
a np. bananowce i canny zostaną wykopane po przymrozkach).

Zawsze chcę, żeby poważne przymrozki nastąpiły jak najpóźniej (12 X w 2010, 28 X w 2009, 17 XI w 2008).
W ubiegłym roku pospieszyłem się i za wcześnie wykopałem canny: po pierwszych przymrozkach, które trochę je uszkodziły, było kilka względnie ciepłych dni.

[greg]

e = 6.112 10^[7.5 t / (237.7 + t)] 0.01 f

Wyżej podany wzór należy rozumieć oczywiście jako:  6.112 * 10^{[7.5 t / (237.7 + t)]} * 0.01 * f?

[Marcin D]

Zaraz podstawie jakieś wartości i zobaczę, ale raczej tak. Ptaszek to potęga a nawias kwadratowy oznacza co ma być w potędze. To jakiś komputerowy zapis bo miałem go w pomocy programu bioklima.

http://www.igipz.pan.pl/geoekoklimat/blaz/BioKlima.htm

---

Jak podstawisz 7C i 100% to 10hpa powinno wyjść.
Wyszło mi 10.02368 wiec potwierdzam prawidłowość wzoru.

---

Scaliłem posty.
greg

[greg]

Dzięki za informację. Jutro wieczorem postaram się dołożyć ten trzeci algorytm.