13.04.2014 Views

Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...

Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...

Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

zegowej warunek brzegowy spełniany jest poprzez kollokowanie go<br />

w wybranych punktach na brzegu obszaru (rysunek 3.1). Metoda kollokacji<br />

brzegowej jest metodą analityczno-numeryczną, co pozwala na dalszą analizę<br />

otrzymanego rozwiązania, m.in. na jego całkowanie, wyznaczanie wartości<br />

numerycznej rozwiązania dla <strong>do</strong>wolnie wybranego punktu rozważanego obszaru<br />

etc.<br />

Przyjmując definicję problemu jak we wzorze (3.1) rozwiązanie można<br />

przedstawić w następującej postaci:<br />

M<br />

u(x) = u n (x) + c i φ i (x), (3.3)<br />

i=1<br />

gdzie u n (x) oznacza rozwiązanie szczególne równania niejednorodnego (3.1), zaś<br />

φ i (x) oznacza funkcje próbne spełniające w sposób ścisły równanie jednorodne<br />

uzyskiwane z równania (3.1).<br />

Jako że metodę kollokacji brzegowej można traktować jako szczególny<br />

przypadek metody ważonych reszt, wobec tego wymaga się spełnienia<br />

następującego warunku:<br />

v̅j(x)[Bu(x) − g(x)]d∂Ω = 0,<br />

dla j = 1, … M, (3.4)<br />

∂Ω<br />

gdzie v̅j są pewnymi funkcjami wagowymi.<br />

Zakładając, że funkcje wagowe będą w postaci funkcji Diraca:<br />

v̅j(x) = δx j − x, (3.5)<br />

gdzie x 1 , … , x N oznaczają punkty kollokacji na brzegu ∂Ω, wówczas warunek<br />

brzegowy (3.2) jest spełniony w N punktach na brzegu rozważanego obszaru.<br />

Postępując w ten sposób otrzymujemy układ N równań o M niewia<strong>do</strong>mych c i .<br />

W przypadku gdy różniczkowy operator przedstawiony w równaniu (3.1) jest<br />

liniowy powstaje układ równań liniowych o następującej postaci:<br />

gdzie A={aij}, b={bj} oznaczają kolejno:<br />

Ac = b (3.6)<br />

29

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!