程序如下:
clear;%清除工作区变量
clc;%清屏
close all;%关闭所有图形窗口
z=[];
for h=8:0.01:9 %确定空气交换系数范围%% 材料参数输入
m1=6;m2=60;m3=36;m4=60;% 分别对四种介质分割m=m1+m2+m3+m4;% 介质分割和
n=6000;% 对时间分割
t=6000;% 总时长l1=0.6/1000;l2=6/1000;l3=3.6/1000;l4=6/1000;% 四种材料厚度lam_1=0.082;lam_2=0.37;lam_3=0.045;lam_4=0.028;% 四种材料的热传导率de_1=300;de_2=862;de_3=74.2;de_4=1.18;% 四种材料的密度
c1=1377;c2=2100;c3=1726;c4=1005;% 四种材料的比热容%% 计算热扩散率
a1=lam_1/(c1*de_1);% I层材料的热扩散率
a2=lam_2/(c2*de_2);% II层材料的热扩散率
a3=lam_3/(c3*de_3);% III层材料的热扩散率
a4=lam_4/(c4*de_4);% IV层材料的热扩散率%% 材料长度分割和时间步长分割求解
derta_x1=l1/m1;% I层材料的分割长度
derta_x2=l2/m2;% II层材料的分割长度
derta_x3=l3/m3;% III层材料的分割长度
derta_x4=l4/m4;% IV层材料的分割长度
derta_t=t/n;% 时间步长分割%% 计算各层介质剖分的步长比
r1=derta_t/derta_x1^2*a1;% 第I层介质剖分的步长比r2=derta_t/derta_x2^2*a2;% 第II层介质剖分的步长比r3=derta_t/derta_x3^2*a3;% 第III层介质剖分的步长比r4=derta_t/derta_x4^2*a4;% 第IV层介质剖分的步长比u=zeros(m+1,n+1);% 定义四层耦合介质温度分布矩阵%% 初始条件和边界条件
u(:,1)=37;%初始条件
u(1,:)=75;%边界条件*****************************************************************有问题,改成变化的温度
%% 差分格式的系数矩阵的构造
省略……
%% 构造右端项
省略……
%% 追赶法求解
bb=diag(A)';
aa=[0,diag(A,-1)'];
c=diag(A1)';
N=length(bb);
L=zeros(N);
uu0=0;y0=0;aa(1)=0;
L(1)=bb(1)-aa(1)*uu0;y(1)=(b(1)-y0*aa(1))/L(1);
uu(1)=c(1)/L(1);
for i=2:(N-1)
L(i)=bb(i)-aa(i)*uu(i-1);
y(i)=(b(i)-y(i-1)*aa(i))/L(i);
uu(i)=c(i)/L(i);
end
L(N)=bb(N)-aa(N)*uu(N-1);
y(N)=(b(N)-y(N-1)*aa(N))/L(N);
x(N)=y(N);
for i=(N-1):-1:1x(i)=y(i)-uu(i)*x(i+1);
end
u(2:m+1,k)=x';
end
求问如何把其中的边界条件改成变量,如下图所示