Ищу велосипеды

Feb 09, 2016 02:28

Хочу несколько вещей реализовать, и по работе, и для своего ScanCombine. Найти их реализованными пока не смог, но может быть плохо искал, не по тем ключевым словам ( Read more... )

imagegraph2txt, Программки

Leave a comment

Comments 7

axsmyth February 9 2016, 00:41:03 UTC
но зачастую начинаются "высокочастотные пульсации", почти всегда вблизи нуля, иногда мне настолько "везет", что излучение некоторых спектров дает мне отрицательную освещенность!

Очень вероятно, что я многое упускаю, но не в теореме ли Котельникова тут дело?

Reply

nabbla1 February 9 2016, 00:51:28 UTC
Скорее нечто такое: https://ru.wikipedia.org/wiki/%D0%A4%D0%B5%D0%BD%D0%BE%D0%BC%D0%B5%D0%BD_%D0%A0%D1%83%D0%BD%D0%B3%D0%B5

И можно на видеоролике выше в записи увидеть, как иногда вместо гладкой кривой получаем выбросы вверх-вниз, и чуть-чуть меняя положение точки можно с этой ситуацией справиться, и человеку это как-то на удивление легко дается, а как это запрограммировать - не знаю пока.

Reply

axsmyth February 9 2016, 01:06:45 UTC
Да, вы скорей всего правы.
Если обращать внимание на последовательность ваших ручных действий, то бросается в глаза, что сильнее всего влияет на появление этих высокочастотных искажений положение не краёв регулируемого отрезка, а закраевых точек. То есть, если есть кусок интерпляции А-B-C-D, то на точность интерполяции на отрезке B-C (а в данном случае речь не совсем об интерполяции, скорее уж об алгоритмическом подгоне) влияют (может в порядке очерёдности, но не факт):
1. абсцисса точек В и С
2. ордината точек A и D

Это то, что удалось заметить. Как это алгоритмизировать мне пока не ясно.

ЗЫ. Хотя играет ещё роль угол наклона большей части линии фактического изображения. Если он будет меньше 45 градусов, то, возможно, абсцисса и ордината поменяются местами.

Reply


thexder666 February 9 2016, 10:49:00 UTC
Так надо СНАЧАЛА построить интерполяцию, а потом брать интеграл произведения!

Reply

nabbla1 February 9 2016, 11:54:33 UTC
Определенный интеграл (особенно по всей области определения) посчитать вообще не проблема - на каждом отрезке у нас даны величины X0[k], Y0[k], b[k], c[k], d[k], так что на интервале X[k]..X[k+1] функция задается как Y(x) = Y0[k]+b[k]*(x-X0[k])+c[k]*(x-X0[k])^2+d[k]*(x-X0[k])^3 - проинтегрировать ее по отрезку раз плюнуть, выходит Y0[k]*dX+b[k]*dX^2/2+c[k]*dX^3/3+d[k]*dX^4/4, где dX=X[k+1]-X[k ( ... )

Reply

thexder666 February 9 2016, 12:04:38 UTC
Насколько ЧМы помню, надо брать корни полинома Чебышева

Reply

nabbla1 February 11 2016, 11:20:41 UTC
Это для интерполяции полиномом, когда у нас N точек данных и мы хотим через них провести полином степени N-1.

А тут у нас сплайны, и хотя коэффициенты каждого сплайна определяются в конечном счете всеми точками, они все-таки куда более локальны, поэтому если мы середину графика возьмем, окажется, на нее не влияет, как мы крайние точки расставили.

Reply


Leave a comment

Up