Exponenta.ru

Ребус - это неотформатированное месиво из команд maple в ваших листингах :), а это ещё куда ни шло. В общем, p(T) - это искомая функция в виде десятки с показателем степени b[1] + b[2]/T + b[3] * log10(T). Даны 10 точек экспериментальных данных (насколько они реальны я не спрашивал). Эти самые точки...

Мне показали вот такие данные:

Да в ряд-то разложить каждый может, ещё со времён дедушки Ньютона, он был спец. Да, с заменой это тема, но хотелось как у О.Бендера - сразу. А на счёт численного Якобиана я всё равно попробую. Считать оно может будет и дольше, зато это можно реализовать без символьных вычислений. Над точностью можно...

И ещё немного странный вопрос: А зачем мы Якобиан вычисляем символьно? Не проще ли численно его аппроксимировать?

Постановка задачи могла бы выглядеть так: http://s48.radikal.ru/i120/1210/f2/41d2aea74edd.png Как видно из предыдущего моего поста, для алгебраических кривых (суммы степеней x) эта система очень просто решается методом даже на SMath. Для K = 4 параметров и M = 10 точек - t = 2 секунды работы интерпр...

Ну, не знаю, Алексей Борисович, это конечно красиво смотрится и всё такое, но чем это лучше, скажем, вот этого: http://en.smath.info/forum/yaf_postsm7901_Non-Linear-Solvers--root-finding-methods.aspx#post7901 В плане демонстрации. Я вот от алгебраической регрессии почему-то к нелинейной перейти не м...

Алексей Борисович, я типа не понял, а где моя типа медаль или типа печенька за труды праведные? Если товарищ из Израиля специалист, то он и без меня разберётся, куда уж мне, самокатному... И вообще, инженеры сами по себе не склонны усложнять себе жизнь, составляя системы нелинейных уравнений. Что эт...

Короче говоря, у него есть как-бы нелинейная функция с базисами (1, 1/x, log10(x)). Да, похоже, что суперпозиция базиса с коэффициентами находится в показатели степени десятки. Точки с данными у него - это давления при какой-то температуре. Т.е. иксом является как-бы T, а игриком - p, т.е. f(x) = p(...

Народ начал предлагать свои задачки, но немного не с того конца: http://en.smath.info/forum/yaf_postsm8042_Draghilev-s-method.aspx#post8042 Один участник, увидев нахождение точек экстремума у Химмельблау, решил, что вполне можно сделать тоже самое и для критерия минимума суммы квадратов невязок. В п...

Я немного доработал документ и теперь отчёты у меня влезают на одну страничку:

Draghilev's method. Finding extremes. Himmelblau function version 0.2.pdf

Теперь можно ставить дело на поток. Я автоматизировал работу для размерностей n >=2, как когда-то для Mathcad.

Вот, кажется правильно выписал:

Вот нашёл на просторах функции для тестирования: http://deepblue.lib.umich.edu/bitstream/2027.42/30811/1/0000469.pdf Я вот не знаю, а дробные показатели относят к алгебраичности? Наверное да, но всё равно реальный многомерный пример из химии куда интереснее искусственных.

Да уж как обычно оформляю, ещё со времён техники работы с системой UMP'а. Этим методиками отчёта по решениям уже несколько лет. Там в теме я прикладываю pdf файлы с визуализацией расчётного документа для каждого примера. Этот ихний пример очень примитивен. По сути - две плоскости и одна фигурка вида...

Я ошибся немного, это картинка чуть для другого уравнения. Я в первом уравнении вместо того, чтобы отнять 0,5 - умножил на 1/2. Потом нарисую правильную картинку.

Если взять точку вдоль линии пересечения зелёной и синей поверхности, то будет два корня: