Здравствуйте, гость ( Вход | Регистрация )
| nqks |
5.10.2010, 22:51
Сообщение
#1
|
|
Школьник ![]() Группа: Продвинутые Сообщений: 10 Регистрация: 18.4.2010 Город: Москва |
см пост номер 5(в этом сообщении приведена исходная задача,которая была упрощена)
предмет-защита предприятий и гражданская оборона,но знания требуются в области тервера.(не обращайте внимание на формулировки,важны лишь формульные значения). собственно,в результате решения задачи о получении модели поля концентрации для мгновенного линейного источника,расположенного перпендикулярно ветру,при равномерном распространении вещества. исходное-dc(x,y,z,t)=dQf(x,t)u(y-y')v(z) при решении данного уравнения получился результат: с(x,y,z,t)=(1/4yo)*[erf((y+yo)/(sqrt(2)*sigmay))-erf((y-yo)/(sqrt(2)sigmay))]; (1) yo-пределы интегрирования в ходе решения. sigmay-среднее квадратическое отклонение. собственно,вопрос:как проинтегрировать выражение (1) от минус бесконечности до плюс бесконечности? *разумеется,в ответе должна получиться единица* по поводу своих идей:была пока только единственная идея-перейти в сферическую систему координат,но она особым успехом не увенчалась. |
![]() ![]() |
| malkolm |
6.10.2010, 15:22
Сообщение
#2
|
|
Старший преподаватель ![]() ![]() ![]() ![]() ![]() Группа: Преподаватели Сообщений: 2 167 Регистрация: 14.6.2008 Город: Н-ск Вы: преподаватель |
с(x,y,z,t)=(1/4yo)*[erf((y+yo)/(sqrt(2)*sigmay))-erf((y-yo)/(sqrt(2)sigmay))]; (1) собственно,вопрос:как проинтегрировать выражение (1) от минус бесконечности до плюс бесконечности? *разумеется,в ответе должна получиться единица* по поводу своих идей:была пока только единственная идея-перейти в сферическую систему координат,но она особым успехом не увенчалась. Запишите каждый erf через интеграл по t и измените порядок интегрирования: вместо интегрирования внитри по t, снаружи по y, сделайте интегрирование внутри по y, а снаружи по t. Единица должна вылезти, из свойств erf на бесконечности. Вообще, если F(x) есть произвольная функция распределения, то функция (F(x-a)-F(x-b )) / (b-a) является плотностью распределения. А именно, распределения суммы двух независимых случайных величин - с распределением F и с равномерным на [a, b ] распределением. Так что если константы у Вас там правильные, то функция обязана быть плотностью, т.е. быть нормированной. Разумеется, чтобы к этому виду привести, нужно erf заменить на функцию распределения любым преобразованием сдвига и масштаба, например erf(x) = 2Ф(x*sqrt{2})-1, где Ф(x) - функция стандартного нормального распределения. |
| nqks |
6.10.2010, 15:26
Сообщение
#3
|
|
Школьник ![]() Группа: Продвинутые Сообщений: 10 Регистрация: 18.4.2010 Город: Москва |
Запишите каждый erf через интеграл по t и измените порядок интегрирования: вместо интегрирования внитри по t, снаружи по y, сделайте интегрирование внутри по y, а снаружи по t. Единица должна вылезти, из свойств erf на бесконечности. Вообще, если F(x) есть произвольная функция распределения, то функция (F(x-a)-F(x-b )) / (b-a) является плотностью распределения. А именно, распределения суммы двух независимых случайных величин - с распределением F и с равномерным на [a, b ] распределением. Так что если константы у Вас там правильные, то функция обязана быть плотностью, т.е. быть нормированной. Разумеется, чтобы к этому виду привести, нужно erf заменить на функцию распределения любым преобразованием сдвига и масштаба, например erf(x) = 2Ф(x*sqrt{2})-1, где Ф(x) - функция стандартного нормального распределения. дада,это понятно,но нужно обойтись без функции стандартного нормального распределения. в общем,я решил для 2-мерного пространства,осталось для 3-мерного.чуть выше я привёл интеграл,в котором происходит загвоздка. |
| malkolm |
6.10.2010, 20:28
Сообщение
#4
|
|
Старший преподаватель ![]() ![]() ![]() ![]() ![]() Группа: Преподаватели Сообщений: 2 167 Регистрация: 14.6.2008 Город: Н-ск Вы: преподаватель |
дада,это понятно,но нужно обойтись без функции стандартного нормального распределения. Не понял: Вам ехать или шашечки? То, что данная функция имеет единичную площадь подграфика, получается элементарной сменой порядка интегрирования. Впрочем, хозяин - барин (IMG:style_emoticons/default/smile.gif) |
nqks интеграл от erf 5.10.2010, 22:51
tig81 Простите, запамятовала, наверное, а что такое erf? 6.10.2010, 12:45
nqks функция ошибок.я чуть позже напишу сегодняшние нар... 6.10.2010, 14:36
tig81
функция ошибок.
Спасибо. 6.10.2010, 14:37
nqks собственно говоря,задача упростилась.теперь всё св... 6.10.2010, 15:16
nqks я только-только начал тервер учить,нету у меня еще... 6.10.2010, 21:04
malkolm Первоначальная задача вроде бы иначе выглядела: ... 7.10.2010, 16:33![]() ![]() |
|
Текстовая версия | Сейчас: 19.4.2026, 19:56 |
Зеркало сайта Решебник.Ру - reshebnik.org.ru