Regrese ve VB - 1. díl
 x   TIP: Přetáhni ikonu na hlavní panel pro připnutí webu
Reklama
Reklama

Regrese ve VB - 1. dílRegrese ve VB - 1. díl

 

Regrese ve VB - 1. díl

Google       Google       28. 2. 2007       13 974×

Podíváme se na tuto zajímavou metodu používanou nejen v matematice. Regresní analýza slouží k odhadování výsledků z jiných známých hodnot.

Reklama
Reklama

Regrese je metoda, která je, řekl bych, docela zajímavá a hlavně prakticky použitelná. Regresí je několik druhů – lineární, exponenciální, logaritmická, mocninná a další. Předmětem toho dvoudílného miniseriálu budou první tři zmíněné druhy. A rovnou se pustíme do lineární regrese.

Lineární regrese:

Lineární regrese nám slouží pro výpočty hodnot Y ze zadané hodnoty X a ze „vzorových“ dat. Zadáme-li tedy data (alespoň dvě dvojce) ve tvaru (x,y) (4,3 a 5,1) a budeme-li potom chtít zjistit hodnotu y pro x = 6, získáme číslo −1. To vyplývá i z grafu:

Pokud ovšem zadáme hodnot více a ony nejsou úplně lineární (spojnice bodů nevytvoří přímku), potom se celý graf, řekněme, nějak zprůměruje (zjednodušeně řečeno). Zadáme-li tedy hodnoty 1,2; 2,4; 3,8, nebude průběhem zelená křivka, ale červená přímka.

Nyní ale k programování. Prvně si vytvoříme regresi pomocí, řekněme, statického kódu, který bude y vypočítávat pouze podle dvou vzorových dvojic dat. No a nebylo by úplně od věci uvést vzorečky, podle kterých to všechno počítá.

No a teď tedy kód:

    xa = 1 'vzorky x
    xb = 2
    
    ya = 3 'vzorky y
    yb = 4
    
    n = 2 'počet vzorků
    
    'výpočet regresního koeficientu B
    B = ((n * ((xa * ya) + (xb * yb))) - ((xa + xb) * (ya + yb))) / ((n * ((xa * xa) + (xb * xb)) - ((xa + xb) * (xa + xb))))
    
    'výpočet regresního koeficientu A
    A = ((ya + yb) - (B * (xa + xb))) / (n)
    
    x = 0 'hodnota, kterou chceme dopočítat
    
    y = A + (B * x) 'závěrečný vzorec pro výpočet y

Pokud chceme používat tři parametry, potom je zde upravený kód, který je schopný toliko parametrů pojmout.

    xa = 1 'vzorky x
    xb = 2
    xc = 3
    
    ya = 2 'vzorky y
    yb = 2.5
    yc = 3
    
    n = 3 'počet vzorků
    
    'výpočet regresního koeficientu B
    B = ((n * ((xa * ya) + (xb * yb) + (xc * yc))) - ((xa + xb + xc) * (ya + yb + yc))) / ((n * ((xa * xa) + (xb * xb) + (xc * xc)) - ((xa + xb + xc) * (xa + xb + xc))))
    
    'výpočet regresního koeficientu A
    A = ((ya + yb + yc) - (B * (xa + xb + xc))) / n
    
    x = 0 'hodnota, kterou chceme dopočítat
    
    y = A + (B * x) 'závěrečný vzorec pro výpočet y

Ale nyní přistupme k variabilnějšímu programu. Nejenže bude umožňovat vložit n dat, ale bude umět vykreslit i graf (i když ne v tomto článku – na to si budete muset ob článek počkat).

    Dim data(256, 1) As Single
    Dim sumxy As Single
    Dim sumx As Single
    Dim sumy As Single
    Dim sumx2 As Single
    Dim n As Byte
    Dim x As Single
    Dim y As Single
    Dim A As Single
    Dim B As Single
    
    n = InputBox("Zadej počet vzorů", "Požadavek na data", "2")
    temp = n - 1
    
    For i = 0 To temp
    x = InputBox("Zadej hodnotu x dat pro n = " & i + 1, "Požadavek na data", "0")
    y = InputBox("Zadej hodnotu y dat pro n = " & i + 1, "Požadavek na data", "0")
    data(i, 0) = x
    data(i, 1) = y
    Next i
    
    sumxy = 0
    sumx = 0
    sumy = 0
    sumx2 = 0
    
    For i = 0 To temp
    sumxy = sumxy + data(i, 0) * data(i, 1)
    sumx = sumx + data(i, 0)
    sumy = sumy + data(i, 1)
    sumx2 = sumx2 + (data(i, 0) * data(i, 0))
    Next i
    
    B = ((n * sumxy) - (sumx * sumy)) / ((n * sumx2) - (sumx * sumx))
    Print B
    A = (sumy - (B * sumx)) / n
    Print A
    
    Min = InputBox("Zadej interval grafu hodnotu min: ", "Požadavek na data", "0")
    Max = InputBox("Zadej interval grafu hodnotu max: ", "Požadavek na data", "0")
    
    For x = Min To Max
        y = A + (B * x)
        Print "x = " & x & ", y = " & y
    Next x

Tedy to by byla základní část bez grafu. K tomu grafu zatím jen tak naznačím. Funkce: Line − (offsetx + (x * zoom),offsety − (y * zoom)) to všechno jistí. Pomocí proměnných offsetx a offsety si zvolíte střed a pomocí zoomu potom měřítko.

V cyklu For … Next nezapomínejte na možnost použití Step. Určitě by vám to při práci hodně pomohlo a možná nebudete potřebovat tolik zoomu. Pozor ale, Step mění měřítko pouze na jedné ose.

Není úplně od věci si při grafu zároveň tisknout na obrazovku kótování, je to sice celkem dost kódu, ale stojí to za to.

No v příštím díle se podíváme na exponenciální regresi.

×Odeslání článku na tvůj Kindle

Zadej svůj Kindle e-mail a my ti pošleme článek na tvůj Kindle.
Musíš mít povolený příjem obsahu do svého Kindle z naší e-mailové adresy kindle@programujte.com.

E-mailová adresa (např. novak@kindle.com):

TIP: Pokud chceš dostávat naše články každé ráno do svého Kindle, koukni do sekce Články do Kindle.

3 názory  —  3 nové  
Hlasování bylo ukončeno    
0 hlasů
Google
(fotka) Jiří ChytilAutor programuje ve VB, zajímá se o elektrotechniku, studuje na SOŠ Elektrotechnické - obor číslicová technika.
Web    

Nové články

Obrázek ke článku NEWTON Media prohledá 200  milionů mediálních zpráv během sekund díky Cisco UCS

NEWTON Media prohledá 200 milionů mediálních zpráv během sekund díky Cisco UCS

Česká společnost NEWTON Media provozuje největší archiv mediálních zpráv ve střední a východní Evropě. Mezi její zákazníky patří například ministerstva, evropské instituce nebo komerční firmy z nejrůznějších oborů. NEWTON Media rozesílá svým zákazníkům každý den monitoring médií podle nastavených klíčových slov a nabízí online službu, kde lze vyhledat mediální výstupy v plném znění od roku 1996.

Reklama
Reklama
Obrázek ke článku Delphi 10.1.2 (Berlin Update 2) – na co se můžeme těšit

Delphi 10.1.2 (Berlin Update 2) – na co se můžeme těšit

Touto roční dobou, kdy je zem pokrytá barevným listím a prsty křehnou v mrazivých ránech, se obvykle těšíme na zbrusu novou verzi RAD Studia. Letos si však ale budeme muset počkat na Godzillu a Linux až do jara. Vezměme tedy za vděk alespoň updatem 2 a jelikož dle vyjádření pánů z Embarcadero se budou nové věci objevovat průběžně, pojďme se na to tedy podívat.

Obrázek ke článku Konference: Moderní datová centra pro byznys dneška se koná už 24. 11.

Konference: Moderní datová centra pro byznys dneška se koná už 24. 11.

Stále rostoucí zájem o cloudové služby i maximální důraz na pružnost, spolehlivost a bezpečnost IT vedou k výrazným inovacím v datových centrech. V infrastruktuře datových center hraje stále významnější roli software a stále častěji se lze setkat s hybridními přístupy k jejich budování i provozu.

Obrázek ke článku Konference: Mobilní technologie mají velký potenciál pro byznys

Konference: Mobilní technologie mají velký potenciál pro byznys

Firmy by se podle analytiků společnosti Gartner měly  rychle přizpůsobit skutečnosti, že mobilní technologie už zdaleka nejsou horkou novinkou, ale standardní součástí byznysu. I přesto - nebo možná právě proto - tu nabízejí velký potenciál. Kde tedy jsou ty největší příležitosti? I tomu se bude věnovat již čtvrtý ročník úspěšné konference Mobilní řešení pro business.

loadingtransparent (function() { var po = document.createElement('script'); po.type = 'text/javascript'; po.async = true; po.src = 'https://apis.google.com/js/plusone.js'; var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(po, s); })();
Hostujeme u Českého hostingu       ISSN 1801-1586       ⇡ Nahoru Webtea.cz logo © 20032016 Programujte.com
Zasadilo a pěstuje Webtea.cz, šéfredaktor Lukáš Churý