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

UptoLike

ВУЗ: 

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

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

Рубрика: 

13
решить данную СЛАУ и вывести результаты решения, которые представлены
на рисунке 7.
close all
clear all
% Решение уравнения теплопроводности
k=120; % удельная теплоемкость Si (Вт/м/К)
f=1e6; % плотность мощности источника тепла Вт/м^3
Lx=20e-2;
Ly=30e-2;
Tenv=300; % температура (К)
% Создание геометрии
I=15; % число точек по x;
J=17; % число точек по y;
dx=Lx/(I-1);
dy=Ly/(J-1);
X=0:dx:Lx;
Y=0:dy:Ly;
% Решение СЛАУ
N=I*J;
a=zeros(N,N);
b=zeros(1,N);
for i=2:I-1
for j=2:J-1
n=i+(j-1)*I;
a(n,n) =-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;
b(n) =-f;
end
end
% граничные
условия
for j=1:J
i=1;
n=i+(j-1)*I;
a(n,n) = 1;
b(n) = Tenv;
i=I;
n=i+(j-1)*I;
a(n,n) = 1;
b(n) = Tenv;
end
for i=2:I-1
j=1;
n=i+(j-1)*I;
a(n,n) = 1;
b(n) = Tenv;
j=J;
n=i+(j-1)*I;
a(n,n) = 1;

Страницы