Pid Control For Boiler
Pid Control For Boiler
Pid Control For Boiler
01
Backward differencing method:
𝑈(𝑠) 𝑠 + 𝑎
𝐶(𝑠) = =
𝐸(𝑠) 𝑠 + 𝑏
Hence it is proved that, the backward difference method results in digital implementation of
compensator.
Question no. 02
Continuous time representation of system using
Part (a): Forward differencing implementation
The forward differencing implementation of equation
1 1 1 + 𝑎𝑇
𝑢(𝑘 + 1) = 𝑢(𝑘) − 𝑒(𝑘) + 𝑒(𝑘 + 1)
(1 + 𝑏𝑇) (1 + 𝑏𝑇) (1 + 𝑏𝑇)
is done in MATLAB Simulink and the is shown in figure 1:
The step response of figure 1 with sample time of 0.001 𝑠𝑒𝑐, 0.01 sec and 157 𝑚𝑠𝑒𝑐 is shown in figure
2 to figure 4 respectively.
Figure 2: Step response of figure 1 with Ts=0.001 s:
Question 03
Problem A-6-19
Part (1): Open loop transfer function
10(𝑠 + 𝑎)
𝐺(𝑠) =
𝑠(𝑠 + 1)(𝑠 + 8)
Part (2):
Closed loop characteristics equation:
𝑠 3 + 9𝑠 2 + 19𝑠 + 10𝑎 = 0
Part (3):
Simulink diagram:
CL1=feedback(G1*K1,1);
figure;
rlocus(CL1)
grid on
𝑎=2
figure;
plot(t,st,'linewidth',1.5)
hold on
plot(t,sr6_19,'linewidth',1.5)
grid on
title('Step response')
legend('Reference','Output')
Part (2):
Closed loop transfer function:
10
𝐺𝑐 (𝑠) =
𝑠(𝑠 + 1 + 10𝑘)
10
1+
𝑠(𝑠 + 1 + 10𝑘)
After simplification:
1
𝐺𝑐(𝑠) =
𝑠2 + 𝑠(1 + 10𝑘) + 10
Part (3):
Closed loop characteristics equation:
𝑠 2 + 𝑠 + 10𝑘𝑠 + 10 = 0
Now, rearranging above equation:
10𝑘𝑠
1+ =0
𝑠2 + 𝑠 + 10
Taking 𝐾 = 10𝑘, we get:
𝐾𝑠
1+ =0
𝑠2 + 𝑠 + 10
Part (4)-(6):
Simulink diagram:
Part (5)-(6):
%% problem 6-20
k=1;
sys=tf(10,[1 1]);
CL20a=feedback(sys,k);
CL20=feedback(CL20a*tf(1,[1 0]),1);
figure;
rlocus(CL20)
grid on
Figure 13: Root locus of figure 12
Step Response:
figure;
plot(t,st20,'linewidth',1.5)
hold on
plot(t,sr6_20,'linewidth',1.5)
grid on
title('Step response')
legend('Reference','Output')
ylim([0 1.2])
Figure 14: Step response of figure 12