Апроксимация spline1dfitpenalized из ALGLIB

@ToLexx · Регистрация: 21.10.2013

Студворк — интернет-сервис помощи студентам

Здравствуйте!
Необходимо было почистить сигнал аппроксимацией, выбрал spline1dfitpenalized.
Всё прогоняется нормуль, вот только ну очень долго.
при проверке оказалось, что дольше всего проводится сама аппроксимация

C#
1
alglib.spline1dfitpenalized(Time, mas, m, rho, out info, out s, out rep);

Внутрь заношу массив из всего-лишь 50'000 значений, при параметрах rho=1 и m =100, расчёт длиться чуть ли не по 20-30 сек.
Сам метод в который и передаю массивы вот такой:

C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
public List<double> ALGaproc(double[] Time, List<double> Vin, int m, double rho, int FramesADC)
        {
            alglib.spline1dinterpolant s;
            alglib.spline1dfitreport rep;
            int info;
            int len = Vin.Count;
            double[] mas = new double[len];
            for (int i = 0; i < len; i++)
                mas[i] = Vin[i];
            alglib.spline1dfitpenalized(Time, mas, m, rho, out info, out s, out rep);
            var Vout = new List<double>();
            for (int i = 0; i < FramesADC; i++)
                Vout.Add(alglib.spline1dcalc(s, Time[i]));
            return Vout;
        }

на прогоне 10 строки всё встаёт на 20-30 сек(зависит от размер m, при m=750 до 1мин и наоборот при m=20 ->3сек).

Если честно уже неделю сижу, периодически возвращаясь к проблеме. Слишком маленькое значение m=20 взять не могу, аппроксимация проходит дерьмово, если не сказать хуже.

Добавлено через 3 минуты
Забыл сказать) Друг мой также использует эту функцию и при прогоне массива в почти 25000 точек,rho=0.1, m=250, рассчитывается всё практически моментально.

Добавлено через 2 часа 48 минут
Ну если никто не ответил буду пока считать, что слишком долгий расчёт это даже нормально.(

@Storm23 · 20.11.2015, 16:39

ToLexx,
Во-первых, не мешало объяснить что такое spline1dfitpenalized, потому что тут не форум alglib, и скорее всего никто не знает, что это за функция и что она делает.
Во-вторых, что это за загадочный параметр m? И почему он должен быть аж 750 ? Если это число точек для расчета сплайна, то для кубического сплайна достаточно 4 точки. Абсолютно не понятно, зачем вам аж 750 точек, а при меньшем значении "аппроксимация проходит дерьмово". Как вы определили "дерьмовость" аппрокимации третьего порядка? На глаз чтоли? Я вообще не понимаю, как аппроксимация может выгладеть дерьмово. Может у вас и не аппроксимация а интерполяция?
В-третьих не мешало бы видеть исходные данные.

Ну и в любом случае, если m это число точек для каждого сплайна, то алгоритму нужно решить систему с 750 уравнениями для каждой(!) из 50'000 точек. Не удивительно, что он тормозит.

Новые блоги и статьи Все статьи Все блоги /
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .
Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка библиотек: SDL3, Box2D, FreeType, SDL3_ttf, SDL3_mixer и SDL3_image из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .

@Storm23 10427 / 5157 / 1825 Регистрация: 11.01.2015 Сообщений: 6,226 Записей в блоге: 34
	20.11.2015, 16:39
	ToLexx, Во-первых, не мешало объяснить что такое spline1dfitpenalized, потому что тут не форум alglib, и скорее всего никто не знает, что это за функция и что она делает. Во-вторых, что это за загадочный параметр m? И почему он должен быть аж 750 ? Если это число точек для расчета сплайна, то для кубического сплайна достаточно 4 точки. Абсолютно не понятно, зачем вам аж 750 точек, а при меньшем значении "аппроксимация проходит дерьмово". Как вы определили "дерьмовость" аппрокимации третьего порядка? На глаз чтоли? Я вообще не понимаю, как аппроксимация может выгладеть дерьмово. Может у вас и не аппроксимация а интерполяция? В-третьих не мешало бы видеть исходные данные. Ну и в любом случае, если m это число точек для каждого сплайна, то алгоритму нужно решить систему с 750 уравнениями для каждой(!) из 50'000 точек. Не удивительно, что он тормозит. 1