Регуляризация Тихонова в MathCAD Опции

[Ufauser]

Добрый вечер, Коллеги.
В моей работе встал вопрос об использовании регуляризации Тихонова. Интернет не блещет примерами, есть лишь один пример в MathCAD12. Однако, набрав листинг я не получил результата, похожего на разобранный пример.

Функция Y0(x) почему-то пошла вдоль оси абсцисс. Если кто-то сталкивался с подобной проблемой, прошу помощи.

Эскизы прикрепленных изображений

 

[Vallen2006]

Добрый вечер, Коллеги.
В моей работе встал вопрос об использовании регуляризации Тихонова. Интернет не блещет примерами, есть лишь один пример в MathCAD12. Однако, набрав листинг я не получил результата, похожего на разобранный пример.

Функция Y0(x) почему-то пошла вдоль оси абсцисс. Если кто-то сталкивался с подобной проблемой, прошу помощи.

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

[Ufauser]

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

в том-то все и дело, что невязка по учебнику(примеру) идет следующим этапом. Выходит, что приведен или неполный листинг программы, или закралась ошибка прямо в учебнике. Лямбду я проставляю ниже, приравниваю к 1, как в примере. Но после крайней формулы уже программа прекращает решать и выводить графики.

1. Для начала протестируйте код на rhs без шума (в этом случае лямбда надо установить в ноль при поиске решения).
Когда ПР равен нулю метод Тихонова полностью эквивалентен методу максимального правдоподобия (МНК может быть реализацией поиска решения). Если на данных без шума вы сможете получить решение, то можно двигаться дальше.

2. можно запостить хелп либо линку на функцию lsolve (не пользуюсь маткадом)?

2. "lsolve(M, v) Returns the solution x for the linear system of equations M·x = v, using LU decomposition. The BLAS/LAPACK libraries (http://www.intel.com/software/products/mkl/features/lin_alg.htm) from Intel are used.

Arguments:
M is a real or complex matrix. If the matrix is square, it must be non-singular.
v is a real or complex vector or matrix having the same number of rows as M. "

1. к сожалению, я не понимаю что такое RHS без шума. Как раз на примере из учебника хотел понять и решить свое уравнение (Максвелла для E,B,H,D в обратном радиоканале)
"Когда ПР равен нулю метод Тихонова полностью эквивалентен методу максимального правдоподобия (МНК может быть реализацией поиска решения). Если на данных без шума вы сможете получить решение, то можно двигаться дальше." спасибо, попробую. Опишусь, если получится.