Adam Rudziński

Adam Rudziński inżynier elektronik,
fizyk

Temat: interpolacja "na bieżąco"

Krzysztof Łatuszyński:
ad wielomiany, moze to bedzie pomocne:

http://portal.acm.org/citation.cfm?id=1237678

jestem w domu i nie mam dostepu do artykulu, ale googlowanie slow
sequential spline interpolation
cos tam wyrzuca...

Dzięki wielkie! Zwrot-klucz to właśnie "sequential interpolation"
albo "forward interpolation".
Adam Rudziński

Adam Rudziński inżynier elektronik,
fizyk

Temat: interpolacja "na bieżąco"

A jak liczysz błąd?

Wariancja różnicy pomiędzy sygnałem interpolowanym
a dokładnym.
Znaczy całka po całej dziedzinie?

Całka czy suma, po całym czasie trwania jednego
symbolu, czyli przedziale, w którym amplitudy
i fazy (nazwijmy je początkowymi) sinusoid nie
ulegają zmianie. Dobre jest każde coś, co jest
większe, jak przebiegi bardziej się różnią.
Z lenistwa stosowałem funkcję var z Octave.
A mógłbyś podać konkretny przykład? Np. obrazek. I napisz jeszcze jak Ty ten wielomian interpolacyjny wyznaczasz.

Późno jest, nie chce mi się robić obrazka wybuchającej
funkcji sklejanej. :)
A tu wrzuciłem PDFa, w którym m.in. opisałem taki algorytm:
http://home.elka.pw.edu.pl/~arudzins/kwinOFDM.pdf
Adam Rudziński

Adam Rudziński inżynier elektronik,
fizyk

Temat: interpolacja "na bieżąco"

Adam Rudziński:
Adam Rudziński:
Krzysztof Łatuszyński:
(...)

f(t) = (x1 + x2)/2 + ((x2-x1)/2)*sin(-pi/2 + pi*(t-t1)/(t2-t1))

W wolnej chwili zobaczę, co z tego wychodzi.

Sprawdziłem i okazuje się, że ta metoda daje błąd 2 razy
większy, niż klejenie kawałków parabol. Taki sam błąd
jak parabole daje kosinus ze skladowa stala i dowolna faza:
A cos (omega t + phi) + B,
gdzie omega jest wzięte z najszybciej zmiennego przebiegu
w oryginalnym sygnale, a A, B, phi dobierałem na podstawie
wartości trzech sąsiednich próbek. Tyle, że tu trzeba
umieć policzyć arctg i cos. Ale w ogóle pomysł z sinusoidami
wydaje mi się ciekawy - wielkie dzięki za sugestię!

Oczywiste usprawnienie: zamiast
A cos (omega t + phi) + B
zapisać
A cos (omega t) + B + C sin (omega t)
i dobierać A, B, C - to się już da obsłużyć samym mnożeniem
macierzowym. A pewnie jeszcze lepiej wyjdzie jak się
dołoży sinusoidy z "przedostatnią" częstotliwością.

konto usunięte

Temat: interpolacja "na bieżąco"

Adam Rudziński:
Dzięki wielkie! Zwrot-klucz to właśnie "sequential interpolation"
albo "forward interpolation".

I jak rezultaty? Tego szukałeś?
Adam Rudziński

Adam Rudziński inżynier elektronik,
fizyk

Temat: interpolacja "na bieżąco"

Wit Jakuczun:
Adam Rudziński:
Dzięki wielkie! Zwrot-klucz to właśnie "sequential interpolation"
albo "forward interpolation".

I jak rezultaty? Tego szukałeś?

Tegom chciał.

Pierwsza ważna rzecz była już na początku abstraktu - takie algorytmy
mają wrodzoną niestabilność. Czyli to, że mi wybuchało, to normalne.

W artykule twierdzą, że to można jednak wyleczyć, ale jeszcze nie
miałem dłuższej chwili, żeby się w niego lepiej wgryźć i dowiedzieć
jak.

konto usunięte

Temat: interpolacja "na bieżąco"

Adam Rudziński:
Pierwsza ważna rzecz była już na początku abstraktu - takie algorytmy
mają wrodzoną niestabilność. Czyli to, że mi wybuchało, to normalne.
Czytałem. Ciekawe to jest.
W artykule twierdzą, że to można jednak wyleczyć, ale jeszcze
A czy byłbyś w stanie udostępnić ten artykuł? Na priv oczywiście. Nie mam do niego dostępu a temat jest dla mnie bardzo interesujący...



Wyślij zaproszenie do