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

UptoLike

ВУЗ: 

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

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

Рубрика: 

100
-1/dy;
else
error('Parameter v4 have incorrect value');
end
end
end
% Определение коэффициентов и свободных членов СЛАУ,
% соответствующих внутренним точкам области
for l=2:s
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)=...
+R(x(i),y(j))*C(x(i),y(j))/dt+...
(K(x(i),y(j))+K(x(i-1),y(j)))/dx^2+...
(K(x(i),y(j))+K(x(i),y(j-1)))/dy^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*i+j)=...
-K(x(i),y(j))/dx^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*(i-2)+j)=...
-K(x(i-1),y(j))/dx^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*(i-1)+j+1)=...
-K(x(i),y(j))/dy^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-1)+m*(i-1)+j-1)=...
-K(x(i),y(j-1))/dy^2;
a(n*m*(l-1)+m*(i-1)+j,n*m*(l-2)+m*(i-1)+j)=...
-R(x(i),y(j))*C(x(i),y(j))/dt;
b(n*m*(l-1)+m*(i-1)+j)=F(x(i),y(j));
end
end
end
% Решение СЛАУ
u=b/a';
% Преобразование вектора-строки значений искомой функции
% в узлах координатной сетки в матрицу размерности n x m,
% удобную для представления результатов
% в графическом виде
for l=1:s
for i=1:n
for j=1:m
T(i,j,l)=u(n*m*(l-1)+m*(i-1)+j);

Страницы