2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Обратное преобразование Лапласа.
Сообщение17.12.2009, 22:02 
Нужна Ваша помощь, уважаемые господа учёные. Проблема в следующем - нужно перейти от изображения F(p,$\xi$) к оригиналу f(t,$\xi$). В данном случае t - это время, а $\xi$ - координата. Лучший вариант это получить аналитическую зависимость, но если она не находится ничем, то хотя бы программу, реализующую численное решение, написанную для любой из сред, будь то матлаб, мапл, маткад или математика (хотя последнюю я ещё даже не устанавливал). Трудность в том, что фунция F - очень громоздская.


Итак, для начала введём некие вспомогательные величины:
$$\beta_{1,2} = \sqrt{\frac{p}{2} \left( 1+\epsilon+p\pm \sqrt{(1+\epsilon+p)^2 - 4p} \right) }$$

$$ C_1 = -\frac{\epsilon p_0(1-2\nu)\beta_2sh\beta_2\xi_0}{B(p)E}$$
$$ C_3 = -\frac{\epsilon p_0(1-2\nu)\beta_1sh\beta_1\xi_0}{B(p)E}$$
Где $B(p) = [\beta_1^2 -p(1+\epsilon)]\beta_2sh\beta_2\xi_0ch\beta_1\xi_0 - [\beta_2^2 -p(1+\epsilon)]\beta_1ch\beta_2\xi_0sh\beta_1\xi_0$

Тогда первое изображение, от которого надо
Второе изображение, от которого надо "избавиться":
$$F_2(p,\xi) = \frac{1}{\epsilon p} \left( C_1 [\beta_1^2 -p(1+\epsilon)]ch\beta_1\xi + C_3 [\beta_2^2 -p(1+\epsilon)]ch\beta_2\xi \right)$$

Величины $p_0,\epsilon,\xi_0,\nu,E$ считаются известными.

 
 
 
 Re: Обратное преобразование Лапласа.
Сообщение17.12.2009, 22:49 
Maple хорошо находит оригинальі. Попробуйте сначала убедиться что Маple безсилен.

 
 
 
 Re: Обратное преобразование Лапласа.
Сообщение18.12.2009, 00:49 
Да, пробовал мэплом. Не берёт.

 
 
 
 Re: Обратное преобразование Лапласа.
Сообщение18.12.2009, 02:15 
С помощью вьічетов пробовали?

 
 
 
 Re: Обратное преобразование Лапласа.
Сообщение31.12.2009, 17:47 
Что-то у меня не получается выделить особые точки...

 
 
 
 Re: Обратное преобразование Лапласа.
Сообщение13.02.2010, 12:47 
Разобрался. Если кому нужно, то вот.
Тут надо численно находить корни функции B(p), т.к. только корни функции B(p) создают особенность для F1 и F2. А потом нужно использовать третью теорему разложения из операционного исчисления.

 
 
 [ Сообщений: 6 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group