ru.algorithms

 
 - RU.ALGORITHMS ----------------------------------------------------------------
 From : Alexander Shevchenko                 2:5020/1817.26 17 Feb 2003  11:38:15
 To : Evgenij Masherov
 Subject : интеграл
 -------------------------------------------------------------------------------- 
 

 
 14 февраля 2003 10:51, Evgenij Masherov => Alexander Shevchenko:
 
  [skip]
 
  AS>> Тут чего я зашел в тупик. Интеграл считаю методом симпсона.
  AS>> Интеграл "обычный", Гаусовский, аргумент функции гауса 1,4,
  AS>> второй 0.... вообщем это, пожалуй, неважно, а и эха по
  AS>> алгоритмам... вообщем действовать по принципу чем больше, тем
  AS>> лучше не подходит, выдает ответ 0. Похоже он проскакивает
  AS>> единственный выступ (там график с пиком в начале, а после 5 он
  AS>> уже едва виден....). Hу вообщем уменьшаю правую границу...
  AS>> уменьшаю до 21, и точность 10^-9 и получаю более менее правдивый
  AS>> ответ с точностью 8 знаков после запятой... првда считает
  AS>> сравнительно долго (хотя с этим я готов смирится)... воооот...
  AS>> вообщем кто что подскажет? Какой алгоритм более оптимальный в
  AS>> данном случаи?
  EM> А разбить на участки?
 
 Разбил. Имхо наиболее оптимальное разбивание оказалось от 0 до 21 и от 21 до... 
 вообщем считаю как сумму двух интегралов от 0 до 21 и от 21 до 10000. Первый
 интервал считается нармально (в допустимых пределах, правда долго), второй
 считается моментально, но... получаю ответ 2,61201772070907E-9 а маткад говорит 
 что там 1.74066157223855E-12. Разница огромадная... вроде кое-как подгоняю... но
 хочется хоть что бы хотя бы порядок соответствовал, а то кто знает, может с
 другими параметрами вообще будет неверный ответ. Все считаю методом симпсона.
 Значения функции там действительно мизерные (после 21), в районе 40 функция
 равна 2*10^-17.... может другой алгоритм подсчета взять?
 
  EM> Вблизи пика шаг малый, вдали большой?
 
 Hевыходит...
 
  EM>>>>> Требуется дополнительная оценка погрешности. 2. Заменяют
  EM>>>>> переменную, скажем, на y=1/x или y=exp(-x).
  AS>>>> А что это даст?
  EM>>> Бесконечный предел станет конечным.
 
 Значения функции станет бесконечным :)
 
  AS>> Ага, я "работаю" в этом направлении. Правда в ru.math, пока тут
  AS>> не заплюсовали :)
  EM> Тут в принципе может быть смена одной неприятности другой. Вмето
  EM> бесконечного предела - конечный, но с быстро изменяющейся функцией...
 
 Ага, типа того... :)
 
                                                        Alexander Shevchenko
 
 ... Default tagline #2
 ---
  * Origin: !C:\FIDO\GOLDED\JAM\BAD         (2:5020/1817.26)
 
 

Вернуться к списку тем, сортированных по: возрастание даты  уменьшение даты  тема  автор