Методы решения задач математической физики. Рындин Е.А. - 114 стр.

UptoLike

ВУЗ: 

ТТИ ЮФУ | Таганрог

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

Рубрика: 

114
% Определение коэффициентов и свободных членов СЛАУ,
% соответствующих внутренним точкам области
for l=2:s-1
for i=2:n-1
for j=2:m-1
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*(i-1)+j)=...
-2/dt^2+2/dx^2+2/dy^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*i+j)=-1/dx^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*(i-2)+j)=...
-1/dx^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*(i-1)+j+1)=...
-1/dy^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*(i-1)+j-1)=...
-1/dy^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*l+m*(i-1)+j)=1/dt^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-2)+m*(i-1)+j)=...
1/dt^2;
end
end
end
% Решение СЛАУ
u=b/a';
% Преобразование вектора-строки значений искомой функции
% в узлах координатной сетки в матрицу
% размерности n x m x s, удобную для представления
% результатов в графическом виде
for l=1:s
for i=1:n
for j=1:m
U(i,j,l)=u(n*m*(l-1)+m*(i-1)+j);
end
end
end
% Построение графика искомой функции U(x,y,t)
for l=1:s
figure
surf(y,x,U(:,:,l))

Страницы