matlab完整程序，急！大佬帮帮孩子吧

Question

 

Answer 1

% 传递函数建模
clear; clc;
R1 = 5; R2 = 5;
C1 = 20E-6; C2 = 20E-6;
s = tf([1, 0], 1);
C1s = 1 / (s*C1);
C2s = 1 / (s*C2);
G1 = 1 / (1/C1s+1/R1);
G2 = R2 + C2s;
clear R1 R2 C1 C2 C1s C2s;
Gs = G2 / (G1+G2)
clear G1 G2;
 
% 能控、能观
[A, B, C, D] = tf2ss(Gs.Numerator{1, 1}, Gs.Denominator{1, 1});
ctr_matrix = ctrb(A, B);
ctr_order = rank(ctr_matrix);   % 能控矩阵不满秩
obs_matrix = obsv(A, C);
obs_order = rank(obs_matrix);   % 能观矩阵不满秩
 
% 闭环响应
Kp = 3;
Ti = 0.02;
Td = 0.02;
if (Ti==Inf)
    Gc = Kp * (1 + Td*s);
else
    Gc = Kp * (1 + 1/(Ti*s) + Td*s);
end
Gs_close_loop = feedback(Gc*Gs, 1);
 
step(Gs_close_loop, 0.5);        % 阶跃响应
impulse(Gs_close_loop, 0.5);     % 脉冲响应
step(Gc*Gs_close_loop/s, 0.5);   % 斜坡响应
title("Slope Response");

看你提交了两次。。。