Апроксимация 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 точек. Не удивительно, что он тормозит.

Новые блоги и статьи Все статьи Все блоги /
Отчёт о затраченных материалах за определенный период с макетом печатной формы Maks 21.04.2026 Отчёт из решения ниже размещён в конфигурации КА2. Задача: разработка отчёта по затраченным материалам за определённый период, с возможностью вывода печатной формы отчёта с шапкой и подвалом. В. . .	Отчёт о спецтехнике находящейся в ремонте Maks 20.04.2026 Отчёт из решения ниже размещен в конфигурации КА2. Задача: отобразить спецтехнику, которая на данный момент находится в ремонте. Есть нетиповой документ "Заявка на ремонт спецтехники" который. . .	Памятка для бота и "визитка" для читателей "Semantic Universe Layer (Слой семантической вселенной)" Hrethgir 19.04.2026 Сгенерировано для краткого описания по случаю сборки и компиляции скелета серверного приложения. И пусть после этого скажут, что статьи сгенерированные AI - туфта и не интересно. И это не реклама -. . .	Запрет удаления строк ТЧ документа при определённом условии Maks 19.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "Аккумуляторы", разработанного в конфигурации КА2. У данного документа есть ТЧ, в которой в зависимости от прав доступа. . .
Модель заражения группы наркоманов alhaos 17.04.2026 Условия задачи сформулированы тут Суть: - Группа наркоманов из 10 человек. - Только один инфицирован ВИЧ. - Колются одной иглой. - Колются раз в день. - Колются последовательно через. . .	Мысли в слух. Про "навсегда". kumehtar 16.04.2026 Подумалось тут, что наверное очень глупо использовать во всяких своих установках понятие "навсегда". Это очень сильное понятие, и я только начинаю понимать край его смысла, не смотря на то что давно. . .	My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .	Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.

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

Опции темы