Forum Forum 1 Grupy Ćwiczeniowej Strona Główna

Forum 1 Grupy Ćwiczeniowej
Forum studentów informatyki Politechniki Białostockiej
 

Metody Numeryczne - Domówki

 
Napisz nowy temat   Odpowiedz do tematu    Forum Forum 1 Grupy Ćwiczeniowej Strona Główna -> 3 semestr
Zobacz poprzedni temat :: Zobacz następny temat  
Autor Wiadomość
boro



Dołączył: 27 Lut 2007
Posty: 286
Przeczytał: 0 tematów

Ostrzeżeń: 0/5

PostWysłany: Wto 17:11, 20 Lis 2007    Temat postu: Metody Numeryczne - Domówki

Ma ktoś może zrobione zadanie z MN na jutro i jest chętny do podzielenia sie nim /ps1/ ?

mam nadzieje ze to nie bedzie zakwalifikowane jako spam :]


Ostatnio zmieniony przez boro dnia Nie 21:58, 17 Lut 2008, w całości zmieniany 2 razy
Powrót do góry
Zobacz profil autora
czart



Dołączył: 02 Mar 2007
Posty: 168
Przeczytał: 0 tematów

Ostrzeżeń: 0/5
Skąd: Z lasu

PostWysłany: Wto 18:11, 20 Lis 2007    Temat postu:

Ja bym prosił o zadania z 2-3 ostatnich zajęć bo troche odpuściłem sobie i juz teraz nie wiem nawet co trzeba przyniesc ;d
Powrót do góry
Zobacz profil autora
Cheester



Dołączył: 02 Mar 2007
Posty: 44
Przeczytał: 0 tematów

Ostrzeżeń: 0/5

PostWysłany: Wto 18:25, 20 Lis 2007    Temat postu:

wynik Newtona jest bardzo niedokładny, może ktoś ulepszy:

program newka;
uses crt;
var
licznik,licznik2,licznik3:integer;
a,b,x,fa,fb:real;

procedure bisekcja(var x:real; a,b:real;z:integer);
var
fa,fb,fx,eps:real;
begin
fa:=sin(a)+5*cos(5*a);
fb:=sin(b)+5*cos(5*b);
repeat
if (z=1) then licznik:=licznik+1;
if (z=2) then licznik2:=licznik2+1;
eps:=b-a;
if (z=1) then x:=a+(b-a)/2; if (z=2) then x:=a-fa*(b-a)/(fb-fa);
fx:=sin(x)+5*cos(5*x);
if (fa*fx<0) then
begin
fb:=fx;
b:=x;
end
else
begin
fa:=fx;
a:=x;
end;
until eps<0.001;
end;



procedure Newton( xb,xc:real;var w:real;var i:integer);
var xa:real;
eps:real;
begin
eps:=0.001;
xa:=xc;
if (cos(xb)-25*sin(5*xb))=0 then exit;
xb:=xa-(sin(xa)+5*cos(5*xa))/(cos(xa)-25*sin(5*xa));
i:=1;
while abs(xa-xb)>Eps do
begin
if (cos(xb)-25*sin(5*xb))=0 then exit;
xc:=xb-(sin(xb)+5*cos(5*xb))/(cos(xb)-25*sin(5*xb));
xa:=xb;
xb:=xc;
inc(i);
w:=xa;
end;
end;


begin
a:=0;
b:=1;
clrscr;
while (b<3.14*4) do
begin
licznik:=0;
licznik2:=0;
fa:=sin(a)+5*cos(5*a);
fb:=sin(b)+5*cos(5*b);
if fa*fb<0 then
begin
newton(a,b,x,licznik3);
write(x,' '); write(licznik3,' ');
bisekcja(x,a,b,1); write(x,' '); write(licznik,' ');
bisekcja(x,a,b,2); write(x,' '); writeln(licznik2,' ');
end;
a:=a+1;
b:=b+1;
end;
readln;
end.
Powrót do góry
Zobacz profil autora
boro



Dołączył: 27 Lut 2007
Posty: 286
Przeczytał: 0 tematów

Ostrzeżeń: 0/5

PostWysłany: Wto 18:47, 20 Lis 2007    Temat postu:

Dzięki Cheester. Wysłałem ci Czarcie na maila kilka programów. Może coś ci się przyda.
Powrót do góry
Zobacz profil autora
wicher



Dołączył: 07 Mar 2007
Posty: 70
Przeczytał: 0 tematów

Ostrzeżeń: 0/5
Skąd: z 13 posterunku

PostWysłany: Wto 19:31, 20 Lis 2007    Temat postu:

Ja bym też prosił o programy z dwóch (a nawet trzech) ostatnich zajęć bo za bardzo nie wiem co trzeba było zrobić.
Powrót do góry
Zobacz profil autora
fala (aka tomek)



Dołączył: 03 Lis 2007
Posty: 67
Przeczytał: 0 tematów

Ostrzeżeń: 0/5
Skąd: Łapy

PostWysłany: Wto 22:04, 20 Lis 2007    Temat postu:

tylko jak ktoś by chciał się podzielić tymi programami to najlepiej może niech zamieści na forum kody Wink tak żeby większa rzesza mogła z nich skorzystać, bo np ja też nie pamiętam już co trzeba napisać Razz
Powrót do góry
Zobacz profil autora
boro



Dołączył: 27 Lut 2007
Posty: 286
Przeczytał: 0 tematów

Ostrzeżeń: 0/5

PostWysłany: Wto 22:19, 20 Lis 2007    Temat postu:

[link widoczny dla zalogowanych] w katalogu '2 rok' jest kilka programow
Powrót do góry
Zobacz profil autora
czart



Dołączył: 02 Mar 2007
Posty: 168
Przeczytał: 0 tematów

Ostrzeżeń: 0/5
Skąd: Z lasu

PostWysłany: Sob 10:51, 24 Lis 2007    Temat postu:

Co trzeba zrobić z tą ostatnią funkcją (x^2-2) pamięta ktoś ??
Powrót do góry
Zobacz profil autora
eewkaa



Dołączył: 01 Mar 2007
Posty: 50
Przeczytał: 0 tematów

Ostrzeżeń: 0/5
Skąd: z mieszkania

PostWysłany: Nie 16:50, 25 Lis 2007    Temat postu:

Rozwiązać metodą iteracji.
Wypisać rozwiązania (kolejne przybliżenia?) do pliku i w Excelu narysować wykresy.

Albo coś w tym stylu...
Jeśli ktoś pamięta coś więcej niech mnie poprawi.
Powrót do góry
Zobacz profil autora
boro



Dołączył: 27 Lut 2007
Posty: 286
Przeczytał: 0 tematów

Ostrzeżeń: 0/5

PostWysłany: Czw 23:42, 29 Lis 2007    Temat postu:

PS1

wie ktoś co trzeba zrobić z metod numerycznych na środę? przynam się ze nie słuchałem faceta na zajęciach a bardzo chciałbym oddać program w nast.tygodniu :)
Powrót do góry
Zobacz profil autora
vbazyl



Dołączył: 06 Mar 2007
Posty: 94
Przeczytał: 0 tematów

Ostrzeżeń: 0/5
Skąd: Podlasie

PostWysłany: Pią 1:06, 30 Lis 2007    Temat postu:

Z tego co pamiętam to obliczyć całkę metodą prostokątów i jakiś układ równań metodą gaussa. Całka jakby ktoś chciał:

Kod:

float f(float x)
{
  return x*sin(x*x)+5;
};

float calka(int N)
{
  float xp=0;
  float xk=2*3.14;
  float sz=(xk-xp)/N;
  xk=(xp+sz);
  float suma=0;
  float sr=(xp+xk)/2;
  for(int i=0;i<N;i++)
  {
    suma+=sz*f(sr);
    xp=xk;
    xk=xp+sz;
    sr=(xp+xk)/2;
  };
  return suma;
};

int main(int argc, char *argv[])
{
    calka(4);
    setprecision(3);
    cout<<"    n     |    I      |    I(n)-I(n+1)"<<endl;
    cout<<"    2     |  "<<calka(2)<<"  |    "<<calka(2)-calka(4)<<endl;
    cout<<"    4     |  "<<calka(4)<<"  |    "<<calka(4)-calka(8)<<endl;
    cout<<"    8     |  "<<calka(8)<<"   |    "<<calka(8)-calka(16)<<endl;
    cout<<"    16    |  "<<calka(16)<<"  |    "<<calka(16)-calka(32)<<endl;
    cout<<"    32    |  "<<calka(32)<<"  |    "<<calka(32)-calka(64)<<endl;
    system("PAUSE");
    return EXIT_SUCCESS;
}
Powrót do góry
Zobacz profil autora
Wyświetl posty z ostatnich:   
Napisz nowy temat   Odpowiedz do tematu    Forum Forum 1 Grupy Ćwiczeniowej Strona Główna -> 3 semestr Wszystkie czasy w strefie EET (Europa)
Strona 1 z 1

 
Skocz do:  
Możesz pisać nowe tematy
Możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach


fora.pl - załóż własne forum dyskusyjne za darmo
Powered by phpBB © 2001, 2005 phpBB Group
deoxBlue v1.0 // Theme created by Sopel stylerbb.net & programosy.pl

Regulamin