-
Path: news-archive.icm.edu.pl!newsfeed.gazeta.pl!newsfeed.pionier.net.pl!news.glorb.c
om!postnews.google.com!t33g2000yqe.googlegroups.com!not-for-mail
From: bartekltg <b...@g...com>
Newsgroups: pl.comp.programming
Subject: Re: Pocedura całkowania
Date: Tue, 7 Jul 2009 06:54:08 -0700 (PDT)
Organization: http://groups.google.com
Lines: 146
Message-ID: <8...@t...googlegroups.com>
References: <h2t00t$4h$1@atlantis.news.neostrada.pl>
<6...@q...googlegroups.com>
<h2t7he$gdo$1@atlantis.news.neostrada.pl>
<3...@a...googlegroups.com>
<h2umtd$qs5$1@atlantis.news.neostrada.pl>
<9...@c...googlegroups.com>
<h2vibs$13o$1@atlantis.news.neostrada.pl>
NNTP-Posting-Host: 82.210.189.188
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-2
Content-Transfer-Encoding: quoted-printable
X-Trace: posting.google.com 1246974848 30682 127.0.0.1 (7 Jul 2009 13:54:08 GMT)
X-Complaints-To: g...@g...com
NNTP-Posting-Date: Tue, 7 Jul 2009 13:54:08 +0000 (UTC)
Complaints-To: g...@g...com
Injection-Info: t33g2000yqe.googlegroups.com; posting-host=82.210.189.188;
posting-account=CvUQzQoAAABvVQmR58QmR6N4Cev1qhAS
User-Agent: G2/1.0
X-HTTP-UserAgent: Mozilla/4.0 (compatible; MSIE 8.0; Windows NT 5.1; Trident/4.0;
.NET CLR 3.0.04506.648),gzip(gfe),gzip(gfe)
Xref: news-archive.icm.edu.pl pl.comp.programming:182614
[ ukryj nagłówki ]On 7 Lip, 15:23, "slawek" <s...@h...pl> wrote:
> Użytkownik "bartekltg" <b...@g...com> napisał w wiadomości grup
> dyskusyjnych:98a4aeba-51da-4e56-8f7f-edbf2e378...@c1
g2000yqi.googlegroups.com...
Teraz odpowiem wstepnie na szybko
> Teoretycznie tak. W praktyce - weź sobie metodę która łyka na raz 23 węzły
> (dlaczego tyle? nieważne, to tylko przykład). Podziel sobie 10000 punktów w
> następujące sposoby:
> 23+23+..., 1+23+23+..., 2+23+23+...,...,22+23+23+... (oczywiście, dla
> punktów na początku i na końcu jakoś trzeba dać sobie radę - ale tego będzie
> mniej więcej 50 punktów na 10000 czyli 0.5%). Uśrednij te wszystkie wyniki.
> Co ci wyjdzie? Surprise, surprise - wzór trapezów (z małymi zwichrowaniami
> na końcach, ale to 0.5% przypadków, do olania).
Pokazales, jak zerpsuc metode. Brawo. Nie wierzysz, ze kwadratura
rzedu
23 bedzie szybceij zbiegac niz trapezy? To ciezko mi bedzie dalej
rozmawiac:(
> Trochę to śmieszne. Ale robienie tych wszystkich tabelek i współczynników
> nie ma sensu.
Czemu smieszna? jakich tabelek, w programie nic nie robie.
>
> Natomiast dla spline'ów wrzucasz całą krzywą na raz do procedury (wysłałem
> ją do wątku, więc możesz ją sobie przetestować, to uproszczony wariant dla
To sprawy implementacji. WIadomo, ze duzo rzeczy mozna zrobic zle:)
> W zasadzie dla wielomianu n-tego stopnia przechodzącego przez n punktów -
> też tak jest. Cała krzywa idzie od razu. Ale wielomian stopnia 10 tysięcy?
> Brrr... to nie zadziała.
Robisz tak jak w splajnach - _kwadratury zlozone_ pierwsze 23 punkty,
pozniej
kolejne (22 nowe) na kazdym przedziale dopasowujesz wielomian.
Splajny robia dokladnie to samo, tylko dla n punktow bierzesz wyzszy
stopien wielomiany i rzadasz odpowiedniej ciaglosci pochodnych
na granicach.
> > Po pierwsze jednego nie dodasz, dla parabol musisz dodac 2, dla
> > simpspona 3.
>
> Tak, po prostu dodajesz po jednym punkcie - i patrzysz jak wzór ślizga się
> po sobie.
Co to znaczy slizga?
Mowilem juz, nie dodasz tylko jednego punktu. A jak dodales bez
zastanowienia,
to moze i dadala Ci sie liczba postaci 0xFFFFF8 :)
> Spline jest dokładniejszy. Sprawdzone - połowienie kroku całkowania daje
> wartości szybciej zbieżne.
Jaki spline? kubiczny?
Szybiej od czego? trapezow, parabol - oczywiscie. Simpsona- wielomiany
sa
i tu i tu stopnia 3, wiec ok. Wyniki porownania moga byc ceikawe. Ale
nie
beda szybsze od romberga czy gaussa oaprtego na 23 wezlach:)
> > z 'wypuklosci'. Usredniajac to spowrotem dostajesz trapezy - tracisz
> > aproksymowana informacje f'.
>
> Im więcej uśredniania, tym mniej arbitralny jest wynik.
To jest _nieprawda_.
Jakby tak bylo, splajny kubiczne czy simpson wcale ine byly by lepsze
od trapezow, a oczywiscie sa.
> A to już jest ciekawe. Vide procedurka - pętla jest po n, więc byłyby to
> unrolling loopsa po n, gdzie n jest nieznane.
NIe o to chodzilo.
> > Ale po co? Jedyne zadanie, jakie przed Toba stoi to zwiac taka
> > kwadrature,
> > aby jej wezly pokrywalty sie z x[k]. No, chyba, ze x[k] nie sa rozno
> > oddalone,
> > ale o to sie pytalem w pierwszym poscie;p
>
> Są. Na razie. Ale jest ich 10000. A to oznacza, że kwadratura
> "równouprawniająca" te punkty byłaby stopnia 10000. Co jest do zrobienia
To jest nieprawda. Kwadratura zlozona prostokotow tez rownouprawnia
punkty. Zreszta, skad pomysl na "rownouprawnienie" punktow.
Splajny nie maja tej wlasnosci! Zadajac tej wlasnosci tracisz
informacje
o pochodnych i psujesz zbieznosc kwadrtatur.
Chyba rzeczywiscie nie wierzysz, ze prostokaty zbiegaja wolniej niz
trapezy,
trapezy niz parabole, a te wolniej niz simpson..;]
> (Integrate na InterpolatingPolynomial)... ale kompletnie niepraktyczne - ze
> wzrostem stopnia wielomianu rosą szanse na to że owszem będzie interpolował,
Tak, interpolowanie wielomianem tak wysokiego stopnia na
rownooddalonych
wezlach bylo by kiepskim pomyslem.
> ale będzie też złą aproksymacją w sensie L2. Jak chcesz to ci mogę taką
> kwadraturę wysłać :) Powinno być mniej niż 1000 linijek.
I w normie max, tec... takie przyklady rzuca sie na podstawowym kursie
numerkow
pomiedzy interpolacja wielomianami o dowolnych wezlach a wspomnieniem
o
tym, zeby wybrac ladne wezly, np miejsca zerowe weilomianow
czebyszewa.
[o, nawet ladny wykresik dali http://en.wikipedia.org/wiki/Runge_function
http://www.mste.uiuc.edu/exner/ncsa/runge/ ]
Pamietac trzeba , ze w pewien sposob dotyczy to tez splajnow wykokiego
rzedu:)
pozdr
bartekltg
Następne wpisy z tego wątku
- 07.07.09 13:58 slawek
- 07.07.09 14:12 slawek
- 07.07.09 14:22 slawek
- 07.07.09 14:24 bartekltg
- 07.07.09 16:15 slawek
- 07.07.09 16:31 slawek
- 07.07.09 16:36 Wit Jakuczun
- 07.07.09 16:53 slawek
- 07.07.09 17:11 Michoo
- 07.07.09 20:44 Mariusz Marszałkowski
- 08.07.09 04:21 Wit Jakuczun
- 08.07.09 09:45 slawek
- 08.07.09 17:17 Mariusz Marszałkowski
- 08.07.09 20:05 Wit Jakuczun
- 08.07.09 20:34 Mariusz Marszałkowski
Najnowsze wątki z tej grupy
- TCL - problem z escape ostatniego \ w nawiasach {}
- Nauka i Praca Programisty C++ w III Rzeczy (pospolitej)
- testy-wyd-sort - Podsumowanie
- Tworzenie Programów Nieuprzywilejowanych Opartych Na Wtyczkach
- Do czego nadaje się QDockWidget z bibl. Qt?
- Bibl. Qt jest sztucznie ograniczona - jest nieprzydatna do celów komercyjnych
- Co sciaga kretynow
- AEiC 2024 - Ada-Europe conference - Deadlines Approaching
- Jakie są dobre zasady programowania programów opartych na wtyczkach?
- sprawdzanie słów kluczowych dot. zła
- Re: W czym sie teraz pisze programy??
- Re: (PDF) Surgical Pathology of Non-neoplastic Gastrointestinal Diseases by Lizhi Zhang
- CfC 28th Ada-Europe Int. Conf. Reliable Software Technologies
- Młodzi programiści i tajna policja
- Ada 2022 Language Reference Manual to be Published by Springer
Najnowsze wątki
- 2024-10-04 Warszawa => QA Engineer <=
- 2024-10-04 Gdańsk => Specjalista ds. Sprzedaży <=
- 2024-10-04 Warszawa => Senior PHP Laravel Developer (e-commerce) <=
- 2024-10-04 Warszawa => Data Scientist / Data Engineer (predictive modelling) <=
- 2024-10-03 Nieparzyste dmuchanie
- 2024-10-03 Prognozowanie zużycia energii przez PGE?
- 2024-10-03 Re: Drugi ekran na Androidzie
- 2024-10-03 sprawiedliwosc nierychliwa
- 2024-10-03 zloto
- 2024-10-03 Odkurzacz mnie bije :(
- 2024-10-03 Gdańsk => Technical Lead ( (Java Background)) <=
- 2024-10-03 Warszawa => Mid IT Recruiter <=
- 2024-10-03 Olsztyn => Sales Specialist <=
- 2024-10-03 Leszczyna nie zna prawa?
- 2024-10-03 Warszawa => OpenText ECM Specialist <=