请给我编写一个三容液位系统matlab仿真曲线

以下是一个简单的三容液位系统的Matlab仿真曲线，你可以根据需要进行修改和调整：

% 设置系统参数
A = [0.8 0.1 0.1; 0.2 0.7 0.1; 0.3 0.3 0.4];
B = [0.2 0.5; 0.1 0.1; 0.7 0.4];
C = [1 0 0; 0 1 0; 0 0 1];
D = [0 0; 0 0; 0 0];
x0 = [0.5; 0.5; 0];

% 设置控制器参数
Q = diag([1 1 1]);
R = diag([1 1]);
[K,S,e] = lqr(A,B,Q,R);

% 设置仿真时长和时间步长
t = 0:0.1:50;
n = length(t);

% 仿真系统响应
x = zeros(3,n);
y = zeros(3,n);
u = zeros(2,n);
x(:,1) = x0;
for i = 1:n-1
    u(:,i) = -K*x(:,i);
    x(:,i+1) = A*x(:,i) + B*u(:,i);
    y(:,i) = C*x(:,i);
end
y(:,n) = C*x(:,n);

% 绘制仿真曲线
subplot(2,1,1)
plot(t,y(1,:),'r',t,y(2,:),'g',t,y(3,:),'b')
xlabel('Time (s)')
ylabel('Liquid level (m)')
title('Three-tank liquid level system')
legend('Tank 1','Tank 2','Tank 3')

subplot(2,1,2)
plot(t,u(1,:),'r',t,u(2,:),'g')
xlabel('Time (s)')
ylabel('Input flow rate (m^3/s)')
title('Control input')
legend('Pump 1','Pump 2')

这段代码中首先设置了系统的参数，包括状态空间模型矩阵A、B、C、D以及初始状态x0。然后通过LQR方法设计了控制器，并设置了仿真时长和时间步长。最后进行仿真，得到了液位和控制输入的仿真曲线，并通过subplot函数将其分别绘制在两个子图中。