帮我写一个恒温控制系统matlab仿真代码，包括加了智能控制和去掉智能控制的图像对比

好的，以下是一个简单的恒温控制系统的MATLAB仿真代码，包括加了智能控制和去掉智能控制的图像对比：

% 温度控制系统MATLAB仿真代码
% 包括加了智能控制和去掉智能控制的图像对比

% 清空环境变量
clear all;
close all;
clc;

% 设定初始温度和恒温值
T0 = 25;    % 初始温度
T_set = 37; % 恒温值

% 设定系统参数
Kp = 1;     % 比例系数
Ki = 0.5;   % 积分系数
Kd = 0.2;   % 微分系数
K = 0.1;    % 智能控制系数

% 设定仿真时间
tend = 100;
dt = 0.01;
t = 0:dt:tend;

% 初始化变量
T = T0*ones(1,length(t));
e = zeros(1,length(t));
de = zeros(1,length(t));
ie = zeros(1,length(t));
u = zeros(1,length(t));
T_int = T0;

% 模拟PID控制系统
for i = 2:length(t)
    e(i) = T_set - T(i-1);
    de(i) = (e(i) - e(i-1))/dt;
    ie(i) = ie(i-1) + e(i)*dt;
    u(i) = Kp*e(i) + Ki*ie(i) + Kd*de(i);
    T_int = T_int + u(i)*dt;
    T(i) = T_int;
end

% 模拟加入智能控制系统
T_int2 = T0;
for i = 2:length(t)
    e(i) = T_set - T(i-1);
    de(i) = (e(i) - e(i-1))/dt;
    ie(i) = ie(i-1) + e(i)*dt;
    u(i) = Kp*e(i) + Ki*ie(i) + Kd*de(i) + K*(T_set - T(i-1));
    T_int2 = T_int2 + u(i)*dt;
    T(i) = T_int2;
end

% 绘制图像对比
figure;
subplot(2,1,1);
plot(t,T);
hold on;
plot(t,T_set*ones(1,length(t)),'r--');
xlabel('时间 (s)');
ylabel('温度 (℃)');
title('没有智能控制');

subplot(2,1,2);
plot(t,T);
hold on;
plot(t,T_set*ones(1,length(t)),'r--');
xlabel('时间 (s)');
ylabel('温度 (℃)');
title('加入智能控制');

这段代码会模拟一个恒温控制系统，包括PID控制和加入智能控制两种情况。绘制出来的图像可以对比两种情况的温度变化情况。