Exponenta.ru

Пример 7 . Уже всем известный поиск экстремумов функции Химмельблау : f(x,y)=(x^2+y-11)^2+(x+y^2-7)^2 4*x^3+4*x*y-42*x+2*y^2-14=0 2*x^2+4*x*y+4*y^3-26*y-22=0 Примечания : 1. Здесь я немного упростил программный блок решателя. 2. Чёрная точка - начало движения по кривой. https://dl.dropbox.com/u/683...

Пример 6 . tg( x * y + 0.1 ) - x^2 = 0 x^2 + 2 * y^2 - 1 = 0 https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example6/Draghilev%27s%20method.%20Finding%20roots.%20System%205.%20Mathcad15%202D.jpg https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example6/Draghilev%27s%20method.%20Find...

Пример 5 . x^2 + y^2 - 8 = 0 sin( x * y ) * sin( exp( x * y ) ) = 0 https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example5/Draghilev%27s%20method.%20Finding%20roots.%20System%204.%20Mathcad15%202D.jpg https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example5/Draghil...

Пример 4 . ( x^2 + y^2 )^2 - 2 * 1^2 * ( x^2 - y^2 ) = 0 sin( x + y + π / 12 ) - 1.1 * x - ( √6 - √2 ) / 4 = 0 https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example4/Draghilev%27s%20method.%20Finding%20roots.%20System%203.%20Mathcad15%202D.jpg ...

Пример 3 . sin( x + y ) - 1.1 * x - 0.1 = 0 x^2 + y^2 - 1 = 0 https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example3/Draghilev%27s%20method.%20Finding%20roots.%20System%201.1.%20Mathcad15%202D.jpg Документы: 1. Draghilev's method. Finding roots. System 1.1. MC15_M010_20110622.pdf 2. D...

Пример 2 . x^2 + y^2 + z^2 - 2 = 0 2 * x + y - z = 0 x - 2 * y - z = 0 https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example2/Draghilev%27s%20method.%20Finding%20roots.%20System%202.%20Mathcad15%203D.jpg Документы: 1. Draghilev's method. Finding roots. System 2. MC15_M010_20110622.pdf 2. Drag...

Освежим, так сказать, тему. Старые примеры с обновлённым решателем для MC15. Пример 1 . x^2 - y * x + y^2 = 1 sin( 5 * x^2 ) + sin( 4 * y^2 ) = 0 https://dl.dropbox.com/u/68387680/Mathcad/Draghilev/Example1/Draghilev%27s%20method.%20Finding%20roots.%20System%201.%20Mathcad15%202D.jpg...

Доверительный интервал - это статистическая мера, насколько я помню, а я описал задачу аналитически. Вот причём тут дуализм. И мы движемся по кривой, если используем метод Драгилева. В данном примере из 4-х мерного пространства. У нас 4 параметра: b1, b2, b3, v. Если повторить эксперимент ещё раз, т...

Странный пример, вообще говоря. Я понял только одно, что здесь есть один минимум и вполне конкретный без всяких там доверительных интервалов. http://i023.radikal.ru/1210/4c/f7e59c2574bb.jpg Точка b[1] - это минимум. Там вектор производных обращается в ноль. В точке b[2] вектор производных совсем в н...

Нет, Алексей Борисович, увы, это не подходит как точное решение исходной постановки задачи как минимума суммы квадратов. Желательно бы получить более точное решение: b1 = 28.4502559095434933; b2 = -2720.12633981294994; b3 = -7.01208324395984839. С ними сумма квадратов разностей равна: ~= 20,9. Как я...

Логарифмирование действительно упрощает задачу. Во всяком случае в "классическом" Mathcad оно позволяет применить функцию linfit. Задача упрощается, но точность не повышается http://s017.radikal.ru/i424/1210/6d/21f45d8e004d.jpg with(Statistics): T := Vector( [ 236.3, 253.4, 26...

Алексей Борисович, этот знак появился у вас, ко мне какие вопросы? Я лишь скопировал ваш пример и показал, что вы не правильно переписали.

Жду указания на ошибки…

Вот я и указал.

Вот так надо:


Или прологарифмировать заданные r[i]. Видимо тогда всё гораздо упростится. Как это мне раньше в голову не пришло.

T - те же самые, там 273 справа прибавлено к вектору. http://s52.radikal.ru/i138/1210/d7/52fb100c9640.png with(Statistics): T := Vector( [ 236.3, 253.4, 261.5, 270.4, 280.6, 288.4, 299.1, 315.2, 333.6, 353.1 ], datatype = float ): p := Vector( [ 1, 5, 10, 20, 40, 60, 100, 200, 40...

Решения вот такие: