Matlab - Vypocet polynomu – Matematika – Fórum – Programujte.com
 x   TIP: Přetáhni ikonu na hlavní panel pro připnutí webu
Reklama
Reklama

Matlab - Vypocet polynomu – Matematika – Fórum – Programujte.comMatlab - Vypocet polynomu – Matematika – Fórum – Programujte.com

 

Hledá se programátor! Plat 1 800 € + bonusy (firma Boxmol.com)
mori
~ Anonymní uživatel
2 příspěvky
28. 8. 2009   #1
-
0
-

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;

musím to předělat pomocí funkcí varargin a varargout. A také jak přidělat vstup a výstup do souboru.

Nahlásit jako SPAM
IP: 213.191.119.–
Reklama
Reklama
Zjistit počet nových příspěvků

Přidej příspěvek

Toto téma je starší jak čtvrt roku – přidej svůj příspěvek jen tehdy, máš-li k tématu opravdu co říct!

Ano, opravdu chci reagovat → zobrazí formulář pro přidání příspěvku

×Vložení zdrojáku

×Vložení obrázku

Vložit URL obrázku Vybrat obrázek na disku
Vlož URL adresu obrázku:
Klikni a vyber obrázek z počítače:

×Vložení videa

Aktuálně jsou podporována videa ze serverů YouTube, Vimeo a Dailymotion.
×
 
Podporujeme Gravatara.
Zadej URL adresu Avatara (40 x 40 px) nebo emailovou adresu pro použití Gravatara.
Email nikam neukládáme, po získání Gravatara je zahozen.
-
Pravidla pro psaní příspěvků, používej diakritiku. ENTER pro nový odstavec, SHIFT + ENTER pro nový řádek.
Sledovat nové příspěvky (pouze pro přihlášené)
Sleduj vlákno a v případě přidání nového příspěvku o tom budeš vědět mezi prvními.
Reaguješ na příspěvek:

Uživatelé prohlížející si toto vlákno

Uživatelé on-line: 0 registrovaných, 6 hostů

Podobná vlákna

Java- Řešení polynomů — založil SIMICIBI

Program na ratanie polynomu — založil Rendo

 

Hostujeme u Českého hostingu       ISSN 1801-1586       ⇡ Nahoru Webtea.cz logo © 20032016 Programujte.com
Zasadilo a pěstuje Webtea.cz, šéfredaktor Lukáš Churý