Anonymní profil mori – Programujte.com

Zdrawím,
chtěl bych poradit s následujícím příkladem:

% aproximace korenu polynomu pomoci metody Durand-Kernerovou 

% vstupní parametr je vektor polynomu (pol)

%výstupní parametr (roots) - koreny polynomu



function roots = dkm(pol)

%pocatecni hodnota

sv = 0.4+0.6i;



l = length(pol) - 1;

% vektor korenu

s = 1:1:l;

t = 0+0i; %promena slouzi k zasteveni iteraci, kdyz uz jsou hodnoty 

pol = pol / pol(1); % normalizace polynomu (celý polynom se vydelí prvním koeficientem, aby koeficient byl roven 0)

% nastaveni pocatecnich korenu pro aproximaci 

for i = 2:1:l; 

  s(i) = s(i - 1) * sv;

end

% iteranci promena

i = 0;



while not (t == s(l) && i < 100)% t==s(l) = zastaví cyklus, kdyz hodnata predchozi iterace se opakuji

  t = s(l);

  i = i + 1;

  % algoritnusm výpočtu kořenů

  for n = 1:1:l

    d = 1+0i;

    for m = 1:1:l

      if not (m == n)

        d = d * (s(n) - s(m));

      end

    end

    r = 0+0i;

    for h = 1:1:l

      r = (r + pol(h)) * s(n);

    end

    s(n) = s(n) - (r + pol(l + 1)) / d;

  end

  

end



roots = s;

To rudolf_franek : No dobře. :)

TextBox.Text = string.Format("P = [{0:0}, {0:0}]", x, y);

Ahoj, zabýval se někdo z vás výše zmíněným tématem? Potřeboval bych doporučit pár odkazů k prostudování. Nacházím jen samé hardcore věděcké práce plné záhadných vzorců bez jediného příkladu jejich praktického využití.

Možná by u otázek na používaný jazyk a OS mohly být checkboxy místo radiobuttonů. :)

To arithmetic overflow : Nelze "to dotahnout do maxima, konce, nenechat skulinku, vyzdimat to, pokorit to az na kost, rozparat to, ...." neboť jak praví jeden ze zákonů programování "V každém programu je chyba a v každém programu je instrukce navíc". :) Pointa je dotáhnou to do uspokojivého stavu a mě třetina vteřiny připadá uspokojivá. :)
Ale jestli budete provokovat, vrátím se k původnímu plánu a nechám to počítat grafickou kartu a takovej fofr ještě neviděli. :)

To Qwerty : Ty máš pro ten algoritmus nějaké praktické využití nebo je to jen hraní?

To Qwerty : C#. Jednoduchá jednovláknová aplikace. Nejsem si úplně jist, že u úlohy která běží třetinu vteřiny (ostatní je manipulace se souborem), se vyplatí obsluhovat vlákna. :)

Taky mi vyšlo 96312 a běželo 1.676s (včetně parsování souboru). Spotřeba RAM 64084 KB. Athlon 6000+.

Původně jsem to chtěl řešit na grafické kartě (nvidia cuda), ale vykašlal jsem se na to - aby se dosáhlo značného zvýšení výkonu oproti CPU, tak by se zřejmě muselo jednat náročnější výpočet.