MObl L03


ÿþPolitechnika Zwitokrzyska WydziaB Elektrotechniki, Automatyki i Informatyki Katedra ZastosowaD Informatyki Metody obliczeniowe  laboratorium OpracowaB: dr in|. Grzegorz SBoD Instrukcja laboratoryjna nr 3: Aproksymacja Data: 30.03.2012 r. Aproksymacja  wersja robocza 1. Wprowadzenie Aproksymacja polega na znalezieniu funkcji o z góry zaBo|onej postaci, której przebieg bdzie mo|liwie najbli|szy punktom o zadanych wspóBrzdnych (obrazujcym np. wyniki pomiarów). Powszechnie stosowan metod aproksymacji jest metoda najmniejszych kwadratów, u której podstaw le|y d|enie do minimalizacji sumarycznego bBdu kwadratowego, obliczanego wg nastpujcej zale|no[ci: n 2 Á = ( yi - f (xi )) (1) " i = 1 gdzie: n  liczba punktów próbkowania przebiegu; (x , y )  wspóBrzdne i-tego punktu; i i f(x )  warto[ szukanej funkcji w i-tym punkcie. i Ogólnie rzecz ujmujc, celem aproksymacji jest znalezienie wielomianu o postaci: m Qm (x) = ai Å" xi , (2) " i= 0 a elementami szukanymi s wspóBczynniki a . i Dla ich znalezienia nale|y rozwiza ukBad równaD: a0S0 + a1S1 + ‹à + amSm = T0 ñø ôø a0S1 + a1S2 + ‹à + amSm+ 1 = T1 ôø ôø ôø‹à (3) òø a0Sk + a1Sk + 1 + ‹à + amSk + m = Tk ôø ôø‹à ôø ôø a0Sm + a1Sm+ 1 + ‹à + amS2m = Tm óø gdzie: m  zakBadany rzd wielomianu; n Sk = xi k , k = 0, 1, 2, ..., 2m; (4) " i= 0 n Tk = xi k yi , k = 0, 1, 2, ..., m. (5) " i= 0 - 1 - Katedra ZastosowaD Informatyki, Metody obliczeniowe  laboratorium Instrukcja laboratoryjna nr 3: Aproksymacja UkBad równaD (3) mo|na przedstawi w dogodniejszej do rozwizywania formie: S0 S1 S2 ‹à Sm a0 T0 îø ùø îø ùø îø ùø ïø S1 S2 S3 ‹à Sm+ 1 úø ïø a1 úø ïø T1 úø ïø úø ïø úø ïø úø ïø úø ïø úø ïø úø S2 S3 S4 ‹à Sm+ 2 a2 = T2 (6) ïø úø ïø úø ïø î" î" î" ñ" î" î" î"úø ïø úø ïø úø ïø úø ïø Sm Sm+ 1 Sm+ 2 ‹à S2m úø ïø am úø ïøTm úø ðø ûø ðø ûø ðø ûø S0 S1 S2 ‹à Sm a0 T0 îø ùø îø ùø îø ùø ïø S1 S2 S3 ‹à Sm+ 1 úø ïø a1 úø ïø T1 úø ïø úø ïø úø ïø úø úø ïø úø ïø úø S2 S3 S4 ‹à Sm+ 2 A = a2 T = T2 Oznaczajc: S = ïø , , , uzyskuje si równanie: ïø úø ïø úø ïø î" î" î" ñ" î" î" î"úø ïø úø ïø úø ïø úø ïø Sm Sm+ 1 Sm+ 2 ‹à S2m úø ïø am úø ïøTm úø ðø ûø ðø ûø ðø ûø S Å" A = T , (7) którego rozwizaniem jest: (8) A = S- 1 Å" T 2. PrzykBad Dane pomiarowe (wspóBrzdne punktów): i 0 1 2 3 4 x -1,8 -1 0,2 1,2 2 y 5,68 1,99 1,24 9,81 14,5 Obliczenia dla m = 1 i x y k 0 1 2 k 0 1 i i k k 0 -1,8 5,68 x 1 -1,8 3,24 x y 5,68 -10 0 0 0 k k 1 -1 1,99 x 1 -1 1 x y 1,99 -2 1 1 1 k k 2 0,2 1,24 x 1 0,2 0,04 x y 1,24 0,25 2 2 2 k k 3 1,2 9,81 x 1 1,2 1,44 x y 9,81 11,8 3 3 3 k k 4 2 14,5 x 1 2 4 x y 14,5 28,9 4 4 4 S 5 0,6 9,72 T 33,2 28,7 k k 5 0,6 33,2 îø ùø îø ùø S = T = ïø ïø 0,6 9,72úø 28,7úø ðø ûø ðø ûø 9,72 - 0,6 îø ùø S = Macierz dopeBnieD: ïø úø - 0,6 5 ðø ûø 9,72 - 0,6 îø ùø ST = Macierz dopeBnieD transponowana: ïø úø - 0,6 5 ðø ûø Wyznacznik: |S| = 48,2 0,201 - 0,012 îø ùø Macierz odwrotna: S- 1 = ïø úø - 0,012 0,104 ðø ûø - 2 - Katedra ZastosowaD Informatyki, Metody obliczeniowe  laboratorium Instrukcja laboratoryjna nr 3: Aproksymacja 6,33 îø ùø f (x) = 6,33 + 2,56x A = S- 1 Å" T = , czyli szukana funkcja aproksymujca to: ïø 2,56úø ðø ûø 16 14 12 10 8 6 4 2 0 -2 -1,5 -1 -0,5 0 0,5 1 1,5 2 x Rys. 1. Wynik obliczeD dla m = 1. CigB lini oznaczono funkcj aproksymujc. Czarne punkty oznaczaj poBo|enie punktów pomiarowych. Obliczenia dla m = 2 i x y k 0 1 2 3 4 k 0 1 2 i i k k 0 -1,8 5,68 x 1 -1,8 3,24 -5,8 10,5 x y 5,68 -10 18,4 0 0 0 k k 1 -1 1,99 x 1 -1 1 -1 1 x y 1,99 -2 1,99 1 1 1 k k 2 0,2 1,24 x 1 0,2 0,04 0,01 0 x y 1,24 0,25 0,05 2 2 2 k k 3 1,2 9,81 x 1 1,2 1,44 1,73 2,07 x y 9,81 11,8 14,1 3 3 3 k k 4 2 14,5 x 1 2 4 8 16 x y 14,5 28,9 57,8 4 4 4 S 5 0,6 9,72 2,9 29,6 T 33,2 28,7 92,4 k k 5 0,6 9,72 33,2 îø ùø îø ùø ïø úø ïø S = 0,6 9,72 2,9 T = 28,7úø ïø úø ïø úø ïø úø ïø úø 9,72 2,9 29,6ûø 92,4ûø ðø ðø 279 10,5 - 93 îø ùø ïø10,5 53,4 - 8,7úø S = Macierz dopeBnieD: ïø úø ïø - 93 - 8,7 48,2ûø úø ðø 279 10,5 - 93 îø ùø ïø10,5 53,4 - 8,7úø ST = Macierz dopeBnieD transponowana: ïø úø ïø - 93 - 8,7 48,2ûø úø ðø Wyznacznik: |S| = 500 0,558 0,021 - 0,185 îø ùø Macierz odwrotna: S- 1 = ïø 0,021 0,107 - 0,017úø ïø úø ïø - 0,185 - 0,017 0,096 úø ðø ûø - 3 - y Katedra ZastosowaD Informatyki, Metody obliczeniowe  laboratorium Instrukcja laboratoryjna nr 3: Aproksymacja 1,98 îø ùø ïø 2 f ( x) = 1,98 + 2,16x + 2,26x A = S- 1 Å" T = 2,16úø , czyli szukana funkcja aproksymujca to: ïø úø ïø úø 2,26ûø ðø 18 16 14 12 10 8 6 4 2 0 -2 -1,5 -1 -0,5 0 0,5 1 1,5 2 x Rys. 2. Wynik obliczeD dla m = 2. CigB lini oznaczono funkcj aproksymujc. Czarne punkty oznaczaj poBo|enie punktów pomiarowych. 3. Zadania do samodzielnego rozwizania 3.1. Znalez funkcje aproksymujce nastpujce zbiory punktów: i 0 1 2 3 a) x 1,0 2,0 4,0 7,0 y 1,4 2,1 2,2 3,5 i 0 1 2 3 b) x -3,0 1,0 2,0 5,0 y 15,0 -0,5 1,5 12,0 i 0 1 2 3 c) x -3,0 -1,0 2,0 4,0 y 55,0 1,0 10,0 -8,0 Ka|dy przypadek rozwa|y dla kilku (co najmniej trzech) rzdów wielomianu (2). 3.2. Wyznaczy bBdy aproksymacji dla funkcji znalezionych w punkcie 3.1. Oceni przydatno[ zwikszania stopnia wielomianu aproksymujcego z punktu widzenia poprawy dokBadno[ci aproksymacji. 3.3. Napisa aplikacj (w wybranym jzyku programowania), której zadaniem bdzie wyznaczanie wielomianu aproksymujcego (2) Aplikacja powinna umo|liwia Batwe wprowadzanie wspóBrzdnych punktów odniesienia oraz dobór rzdu wielomianu (w granicach od 1 do 5). - 4 - y

Wyszukiwarka

Podobne podstrony:
MObl L03
MObl L02
r3 l03
MObl L09
MObl L01
RR L03 pthread create attr
l03102013 z09
MObl L10
powt mobl 2
l03$102013 z09
MObl L07 uklady
l03112013 z09
MObl L02 interp
MObl L08 mzero
MObl L05
l03112013 z09
MObl L04 matlab
l03
4 6 m3 L03

więcej podobnych podstron