Методы математической физики. Куликова И.В - 23 стр.

UptoLike

ВУЗ: 

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

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

Рубрика: 

23
решить СЛАУ (19) и вывести результаты решения, которые приведены на
рисунке 12 и 13.
close all
clear all
% Решение уравнения теплопроводности
k=120; % коэффициент теплопроводности Si (Вт/м/К)
ro=2330; % плотность (кг/м^3)
c=800; % удельная теплоемкость Дж/(кг*К)
f=-1e6; % плотность источника тепла Вт/м^3
Lx=20e-2;
Ly=30e-2;
Lt=300;
Tenv=300; % температура (К)
% Создание геометрии
I=11; % число точек по x;
J=12; % число точек по y;
S=17; % число точек по t;
dx=Lx/(I-1);
dy=Ly/(J-1);
dt=Lt/(S-1);
X=0:dx:Lx;
Y=0:dy:Ly;
time=0:dt:Lt;
% Заполнение матриц коэффициентов
N=I*J*S;
a=zeros(N,N);
b=zeros(1,N);
% Начальные условия
s=1;
for i=1:I
for j=1:J
n=i+(j-1)*I+(s-1)*I*J;
a(n,n)=1;
b(n)=Tenv;
end
end
for s=2:S
for i=2:I-1
for j=2:J-1
n=i+(j-1)*I+(s-1)*I*J;
a(n,n) =ro*c/dt+k*2/dx^2+k*2/dy^2;
a(n,n+1)=-k/dx^2;
a(n,n-1)=-k/dx^2;
a(n,n+I)=-k/dy^2;
a(n,n-I)=-k/dy^2;
a(n,n-I*J)=-ro*c/dt;
b(n) =-f;
end
end
% граничные условия
for j=1:J
i=1;

Страницы