Решение задач математической физики в системе MatLab. Рындин Е.А - 24 стр.

UptoLike

Составители: 

Рубрика: 

% области х = х(1) (1 - ГУ Дирихле, 2 - ГУ Неймана);
% g1 - функция в правой части граничного условия на первой границе,
% задаваемая строкой символов, заключенных
% в одинарные кавычки;
% v2 - параметр, значение которого определяет
% тип граничного условия на второй границе
% области х = х(n) (1 - ГУ Дирихле, 2 - ГУ Неймана);
% g2 - функция в правой части граничного условия на второй границе,
% задаваемая строкой символов, заключенных
% в одинарные кавычки;
% v3 - параметр, значение которого определяет
% тип граничного условия на третьей границе
% области y = y(1) (1 - ГУ Дирихле, 2 - ГУ Неймана);
% g3 - функция в правой части граничного условия на третьей границе,
% задаваемая строкой символов, заключенных
% в одинарные кавычки;
% v4 - параметр, значение которого определяет
% тип граничного условия на четвертой границе
% области y = y(m) (1 - ГУ Дирихле, 2 - ГУ Неймана);
% g4 - функция в правой части граничного условия на четвертой границе,
% задаваемая строкой символов, заключенных одинарные кавычки.
% Выходные параметры:
% х – вектор-строка координатной сетки по оси х размерности 1 х n;
% y - вектор-строка координатной сетки по оси y размерности 1 х m;
% t - вектор-строка сетки по оси времени размерностью 1 х s;
% U - матрица значений результирующей функции
% в узлах координатной сетки размерностью n х m x s.
% Функции и переменные по умолчанию
if exist('t0')==0
t0=0;
end
if exist('ts')==0
ts=0.2;
end
if exist('s')==0
s=6;
end
if exist('x0')==0
x0=-1;
end
if exist('xn')==0
xn=1;
end
if exist('n')==0
n=18;
end
if exist('y0')==0
y0=-1;
end
if exist('ym')==0
ym=1;
end
24