Ознакомьтесь с нашей политикой обработки персональных данных
  • ↓
  • ↑
  • ⇑
 
Записи с темой: численные методы (список заголовков)
11:07 

Добрый день! Нужно минимизировать функционал такого вида (численно)
`min(x) (int_0^tau f(t)-g(t, x)dt)`
Подскажите, пожалуйста, можно ли как-нибудь упростить задачу и какой здесь метод можно применить.

@темы: Численные методы, Методы оптимизации

15:15 

Приведение к каноническому виду

Добрый вечер! Помогите, пожалуйста, привести систему дифференциальный уравнений:
dx/dt=2*x(t)-y(t)+2*z(t);
dy/dt=x(t)+2*z(t);
dz/dt=2*x(t)+y(t)-z(t)
к каноническому виду u'=f(t,u), где u', u, f - векторы, t - скаляр.

@темы: Численные методы

20:24 

Квадратура Гаусса

Построить квадратуру Гаусса с двумя узлами `I(f) = int_0^pi (sin(x)*f(x))dx = c_1*f(x_1)+c_2*f(x_2)`
читать дальше

@темы: Численные методы

19:07 

соль_по_вкусу
Если в многошаговых методах численного решения ОДУ не выполняется условие нормировки коэффициентов, то это влияет на устойчивость? Или вообще может такое быть, чтоб в этом классе методов не выполнялось условие нормировки коэффициентов?

@темы: Численные методы

15:05 

Численные методы

sorata
Чем дороже нам кто-то,тем хуже мы видим,что причиняем боль этому человеку...
Здравствуйте! Извините, что прошу такой глупой помощи, ведь существует гугл и прочие поисковые системы, но в моих попытках поиска искомых материалов в нужной степени подробности я отыскать не сумел(. Не могли бы вы посоветовать литературу, где просто, но доходчиво разъяснен метод обратных итераций для нахождения собственных чисел и собственных векторов матрицы, желательно с разобранным примером?

Заранее благодарю за ответы, даже за отрицательные.

@темы: Матрицы, Численные методы

12:35 

Разностные схемы + конечный ряд Фурье

Ethera
I'll put a gun to your head and pull the fuckin' trigger. (c)
Всем доброго времени суток.
Я решаю одномерное уравнение теплопроводности в дискретном виде, и мне нужно (по заданию) использовать разложение в конечный ряд Фурье. (Не путать с методом Фурье, где используется разделение переменных). Выжав из гугла всё, что можно, обращаюсь к вам за кое-какими разъяснениями :)
По теме (в более-менее понятном виде) я нашла одну-единственную книжку. ССЫЛКА
Вместе покумекав с преподавателями, мы решили, что алгоритм следующий:
1) Взяв начальное распределение температуры, разложить его в конечный ряд Фурье и найти коэффициенты с[k] по формуле на стр.247(255) учебника по ссылке.
При этом начальное распределение у меня просто в виде набора значений в дискретных точках.
2) Зная коэффициенты, найти на каждом временном слое значение температуры по формуле 14 со стр.247.

Как и следовало ожидать, на выходе у меня какие-то бессвязные астрономические цифры.
Если изначально у меня температура всего стержня 20 градусов в каждой точке, а температура на его концах - 0 градусов в каждый момент времени, при этом коэф-т теплопроводности единица, то температура всего стержня должна уменьшаться до нуля. Но этого не происходит. Отсюда у меня несколько вопросов:

1) Не должна ли в алгоритме использоваться собственно разностная схема, представленная на стр. 246 (254) учебника по ссылке? Мы ищем u[m][p], но подставляем их на новые временные слои, не говоря о u[m][p+1], где u - значение температуры, p - индекс по времени, а m - индекс по координате.
2) Второй вопрос, наверное, вытекает из первого. Мне непонятно, как мы собираемся найти значения температуры, забив на температуру на концах стержня? Т.е. я на каждом шаге говорю что она равна нулю, но в формулах для соседних ячеек это никак не фигурирует.

Всем спасибо за внимание :mosk:

@темы: Численные методы, Уравнения мат. физики

19:19 

Уравнение теплопроводности в конечных разностях в цилиндрической СК на оси

Само уравнение теплопроводности имеет вид:
$ \frac{\partial T}{\partial t} = a\left \{\frac{1}{r}\frac{\partial T}{\partial r} + \frac{\partial^{2} T}{\partial r^{2}} + \frac{\partial^{2} T}{\partial z^{2}}\right \}$

Решаю его методом конечных разностей численно в цилиндрической СК. С точками внутри цилиндра и на его границах все вроде понятно. Внутри, например, получается:
$T_{i+1,j,k} &=& a(T_{i,j,k}) \Delta t \left \{\frac{1}{jh} + \frac{T_{i,j+1,k} -T_{i,j,k}}{h} + \frac{T_{i,j-1,k} - 2T_{i,j,k} + T_{i,j+1,k}}{h^{2}} + \frac{T_{i,j,k+1} - 2T_{i,j,k} +T_{i,j,k-1}}{h^{2}}\right \}$

Но как быть с точками на оси, когда r=0, возникает неопределенность. Не знаю, как избавиться от неё (Т.е. как вывести формулу расчета, а она существует....). Текущий радиус 1/r ведь никуда не деть. Как быть?
Заранее спасибо.

@темы: Численные методы, Уравнения мат. физики

Не решается алгебра/высшая математика?.. ПОМОЖЕМ!

главная