基于偏微分方程离散化计算的地下换热器建模与温度检测matlab仿真

1.课题概述

2.系统仿真结果

3.核心程序与模型

4.系统原理简介

4.1地下换热器的建模

4.2温度检测技术

5.完整工程文件

1.课题概述

基于偏微分方程离散化计算的地下换热器建模与温度检测，模拟这个不锈钢圆桶中土壤的温度场和湿度场。

2.系统仿真结果

3.核心程序与模型

版本：MATLAB2022a

................................................................
    Signal_Check(Index0) = 0;
    MM_updata            = Model_Matrix;
    Model_Matrix         = Signal_Check;
    ij2                  = mod(ij-1,time2)+1;
    Signal_Checks(Index0)= 0;
    MM_updatas           = Model_Matrixs;
    Model_Matrixs        = Signal_Checks;
 
    %矩阵赋值，乘以一个增益因子25.
    if ij2<=time4 
       if check1
          Model_Matrix(index) = abs(waves1(:,ij2)');
       else
          waves1(:,ij2)=25*Model_Matrix(index)';
       end
    end
 
    if mod(ij,time2)==(time4+1)
       check1=~check1; 
       waves1 =-fliplr(waves1);
    end


    %矩阵赋值，乘以一个增益因子25.
    if ij2<=time4 
       if check1
          Model_Matrixs(index) = abs(waves1s(:,ij2)');
       else
          waves1s(:,ij2)=25*Model_Matrixs(index)';
       end
    end
 
    if mod(ij,time2)==(time4+1)
       check1=~check1; 
       waves1s =-fliplr(waves1s);
    end   
 

    if mod(ij,25) == 0
       set(figmatrix,'CData',Model_Matrix);%减去初始温度，使得显示效果明显
       axis equal;
       drawnow;
       axis([0,dx*(LenX-1),0,dx*(LenY-1)]);
       pause(1);
    end
 
    %定义不同深度下的温度
    SX      = LenX/2;
    SY      = LenY/2;
    tmp1s = [tmp1s,mean(mean(Model_Matrix(SX-20:SX+20,SY-20:SY+20)))+initial_tmp];
    %检测位置1
    tmp2s = [tmp2s,mean(mean(Model_Matrix(SX-100:SX-60,SY-100:SY-60)))+initial_tmp];
    %检测位置2
    tmp3s = [tmp3s,mean(mean(Model_Matrix(SX+30:SX+80,SY+30:SY+80)))+initial_tmp];
    %检测位置3
    tmp4s = [tmp4s,mean(mean(Model_Matrix(SX-150:SX-40,SY+50:SY+150)))+initial_tmp];


    %定义不同分层下的湿度
    [XS,YS] = find(abs(Model_Matrixs)>0.1);
    Model_Matrixs1 = Model_Matrixs;
    for jj = 1:length(XS)
        Model_Matrixs1(XS(jj),YS(jj))= 0;
   end
    tmp1ss = [tmp1ss,mean(mean(Model_Matrixs1))+initial_sd];

    [XS,YS] = find(abs(Model_Matrixs)>0.2);
    Model_Matrixs2 = Model_Matrixs;
    for jj = 1:length(XS)
        Model_Matrixs2(XS(jj),YS(jj))= 0;
   end
    tmp2ss = [tmp2ss,mean(mean(Model_Matrixs2))+initial_sd];

    [XS,YS] = find(abs(Model_Matrixs)>0.3);
    Model_Matrixs3 = Model_Matrixs;
    for jj = 1:length(XS)
        Model_Matrixs3(XS(jj),YS(jj))= 0;
   end
    tmp3ss = [tmp3ss,mean(mean(Model_Matrixs3))+initial_sd];

    [XS,YS] = find(abs(Model_Matrixs)>0.4);
    Model_Matrixs4 = Model_Matrixs;
    for jj = 1:length(XS)
        Model_Matrixs4(XS(jj),YS(jj))= 0;
   end
    tmp4ss = [tmp4ss,mean(mean(Model_Matrixs4))+initial_sd];
end

figure;
subplot(121);
plot(tmp1s,'r','linewidth',2);
hold on
plot(tmp2s,'b','linewidth',2);
hold on
plot(tmp3s,'m','linewidth',2);
hold on
plot(tmp4s,'g','linewidth',2);
hold on
xlabel('time(hour)');
ylabel('土壤温度');
grid on
legend('检测点1','检测点2','检测点3','检测点4');
axis([0,150,14,35]);


subplot(122);
plot(tmp1ss,'r','linewidth',2);
hold on
plot(tmp2ss,'b','linewidth',2);
hold on
plot(tmp3ss,'m','linewidth',2);
hold on
plot(tmp4ss,'g','linewidth',2);
hold on
xlabel('time(hour)');
ylabel('土壤湿度');
grid on
legend('检测点1','检测点2','检测点3','检测点4');

axis([0,150,0.188,0.20]);
02_028m