Lecture10-FromAnalogToDigitalControllers.pdf

(418 KB) Pobierz
Lecture 10:
A Digital Control System
From Analog to Digital Controllers
( )
y t
(
)
u t
·
The sampling theorem
t
t
·
Aliasing and presampling lters
(
)
( )
u t
y t
Process
Unsampling signals
·
Hold
Sampler
y k
u k
·
Discrete time approximation of continuous time controller
u k
D-A
Computer
A-D
y k
Frequency domain
State space
t
t
The Sampling Theorem (Shannon 1949)
Proof idea
A continuous-time signal with Fourier transform
F D
w
E 0
for
F (
)
F s (
)
w
w
h
0 , is uniquely dened by its sampled values, provided
the sampling frequency
w
h
w
w
2
p
h satises
w
2
w
0 .
s
s
-
2
w
0
2
w
-
2
w
- w
0
w
2
w
- w
w
N
N
N
N
N
N
N
N
The signal can be reconstructed by the interpolation formula
Let
F D
w
E be the Fourier transform of
and expand the peri-
f
X
[
f D kh E sin
w
D t kh E 2
odic function
F s D
w
E as a Fourier series
s
f D t E
w
D t kh E 2
s
X
[
X
[
k [
E 1
h
C k e ikh w
F s D
F D
k
E
w
w
w
s
k [
k [
It is straightforward to verify that
C k f D kh E
k 0
1
2
,
,
, . . .
The frequency
w
2
is denoted
w
N and called the Nyquist
s
frequency .
Hence M f D kh EN0 F s D
w
E 0 F D
w
E 0 f D t E
&
Dept. of Automatic Control, Lund
1
775864840.551.png 775864840.639.png 775864840.650.png 775864840.661.png 775864840.001.png 775864840.012.png 775864840.023.png 775864840.034.png 775864840.045.png 775864840.056.png 775864840.067.png 775864840.078.png 775864840.089.png 775864840.100.png 775864840.111.png 775864840.122.png 775864840.133.png 775864840.144.png 775864840.155.png 775864840.166.png 775864840.177.png 775864840.188.png 775864840.199.png 775864840.210.png 775864840.221.png 775864840.232.png 775864840.243.png 775864840.254.png 775864840.265.png 775864840.276.png 775864840.287.png 775864840.298.png 775864840.309.png 775864840.320.png 775864840.331.png 775864840.342.png 775864840.353.png 775864840.364.png 775864840.375.png 775864840.386.png 775864840.397.png 775864840.408.png 775864840.419.png 775864840.430.png 775864840.441.png 775864840.452.png 775864840.463.png 775864840.474.png 775864840.485.png 775864840.496.png 775864840.507.png 775864840.518.png 775864840.529.png 775864840.540.png 775864840.552.png 775864840.563.png 775864840.574.png 775864840.585.png 775864840.596.png 775864840.607.png 775864840.618.png 775864840.629.png 775864840.637.png 775864840.638.png 775864840.640.png 775864840.641.png 775864840.642.png 775864840.643.png 775864840.644.png 775864840.645.png 775864840.646.png 775864840.647.png 775864840.648.png 775864840.649.png 775864840.651.png 775864840.652.png 775864840.653.png 775864840.654.png 775864840.655.png 775864840.656.png
Aliasing and frequency folding
1
0
Aliasing and presampling lters
0
5
10
Time
F (w )
F s (w )
-2w N
- w N
0
w N
2w N
-2w N
- w N
0
w N
2w N
w
w
2
Nyquist frequency
N
s
New frequencies since the system is NOT time-invariant
w
w
n
w
n 0
,
1
,
2
, . . .
sampl ed
s
Mini-problem
Example — Feedwater heating in a ship boiler
A wheel, rotating 90 turns per second, is lmed using a
camera with sampling time 10ms. What will the wheel rotation
look like on the lm?
38 min
Pressure
Steam
2 min
Valve
Feed
water
To boiler
Pump
Condensed
water
Temperature
2.11 min
Time
&
Dept. of Automatic Control, Lund
2
775864840.657.png 775864840.658.png 775864840.659.png 775864840.660.png 775864840.662.png 775864840.663.png 775864840.664.png 775864840.665.png 775864840.666.png 775864840.667.png 775864840.668.png 775864840.669.png 775864840.670.png 775864840.671.png 775864840.002.png 775864840.003.png 775864840.004.png 775864840.005.png 775864840.006.png 775864840.007.png 775864840.008.png 775864840.009.png 775864840.010.png 775864840.011.png 775864840.013.png 775864840.014.png 775864840.015.png 775864840.016.png 775864840.017.png 775864840.018.png 775864840.019.png 775864840.020.png 775864840.021.png 775864840.022.png 775864840.024.png 775864840.025.png 775864840.026.png 775864840.027.png 775864840.028.png 775864840.029.png 775864840.030.png 775864840.031.png 775864840.032.png 775864840.033.png 775864840.035.png 775864840.036.png 775864840.037.png 775864840.038.png 775864840.039.png 775864840.040.png 775864840.041.png 775864840.042.png 775864840.043.png 775864840.044.png 775864840.046.png 775864840.047.png 775864840.048.png 775864840.049.png 775864840.050.png 775864840.051.png 775864840.052.png 775864840.053.png 775864840.054.png 775864840.055.png 775864840.057.png 775864840.058.png 775864840.059.png 775864840.060.png 775864840.061.png 775864840.062.png 775864840.063.png 775864840.064.png 775864840.065.png 775864840.066.png 775864840.068.png 775864840.069.png 775864840.070.png 775864840.071.png 775864840.072.png 775864840.073.png 775864840.074.png 775864840.075.png 775864840.076.png 775864840.077.png 775864840.079.png 775864840.080.png 775864840.081.png 775864840.082.png 775864840.083.png 775864840.084.png 775864840.085.png 775864840.086.png 775864840.087.png 775864840.088.png 775864840.090.png 775864840.091.png 775864840.092.png 775864840.093.png 775864840.094.png 775864840.095.png 775864840.096.png 775864840.097.png 775864840.098.png 775864840.099.png 775864840.101.png 775864840.102.png 775864840.103.png 775864840.104.png 775864840.105.png 775864840.106.png 775864840.107.png 775864840.108.png 775864840.109.png 775864840.110.png 775864840.112.png 775864840.113.png 775864840.114.png 775864840.115.png 775864840.116.png 775864840.117.png 775864840.118.png 775864840.119.png 775864840.120.png 775864840.121.png 775864840.123.png 775864840.124.png 775864840.125.png 775864840.126.png 775864840.127.png 775864840.128.png 775864840.129.png 775864840.130.png 775864840.131.png 775864840.132.png 775864840.134.png 775864840.135.png 775864840.136.png 775864840.137.png 775864840.138.png 775864840.139.png 775864840.140.png 775864840.141.png 775864840.142.png 775864840.143.png 775864840.145.png 775864840.146.png 775864840.147.png 775864840.148.png 775864840.149.png 775864840.150.png 775864840.151.png 775864840.152.png 775864840.153.png 775864840.154.png 775864840.156.png 775864840.157.png 775864840.158.png 775864840.159.png 775864840.160.png 775864840.161.png 775864840.162.png 775864840.163.png 775864840.164.png 775864840.165.png 775864840.167.png 775864840.168.png 775864840.169.png 775864840.170.png 775864840.171.png 775864840.172.png
Pre- and postsampling lters
Example – Preltering
(a)
(b)
Typical control problem:
1
1
"Decrease the inuence of a low-frequency process distur-
bance despite high-frequency measurement noise."
0
0
·
Frequency separation
0
10
20
30
0
10
20
30
·
Prelter
w
N
(c)
(d)
1
1
Bessel, Butterworth, ITAE lters
Postsampling lters
·
0
0
Avoid exciting mechanical resonances
Higher order hold
0
10
20
30
0
10
20
30
Time
Time
w
0
.
9 ,
w
0
.
5 ,
w
0
.
1
(6th order Bessel lter)
d
N
alias
Consequence of using prelter
Unsampling signals
Pre- and postlter should be included in the process model
Exception: Fast sampling
·
Shannon
X
[
f D kh E sin D
w
D t kh E 2 E
s
f D t E
A Bessellter can be approximated with a delay:
w
D t kh E2
s
k [
1
·
Zero order hold
0.01
·
First order hold
0.1
1
10
·
Predictive rst order hold
0
0.1
1
10
Frequency, rad/s
6th order Bessel (solid line), time delay (dashed line)
&
Dept. of Automatic Control, Lund
3
775864840.173.png 775864840.174.png 775864840.175.png 775864840.176.png 775864840.178.png 775864840.179.png 775864840.180.png 775864840.181.png 775864840.182.png 775864840.183.png 775864840.184.png 775864840.185.png 775864840.186.png 775864840.187.png 775864840.189.png 775864840.190.png 775864840.191.png 775864840.192.png 775864840.193.png 775864840.194.png 775864840.195.png 775864840.196.png 775864840.197.png 775864840.198.png 775864840.200.png 775864840.201.png 775864840.202.png 775864840.203.png 775864840.204.png 775864840.205.png 775864840.206.png 775864840.207.png 775864840.208.png 775864840.209.png 775864840.211.png 775864840.212.png 775864840.213.png 775864840.214.png 775864840.215.png 775864840.216.png 775864840.217.png 775864840.218.png 775864840.219.png 775864840.220.png 775864840.222.png 775864840.223.png 775864840.224.png 775864840.225.png 775864840.226.png 775864840.227.png 775864840.228.png 775864840.229.png 775864840.230.png 775864840.231.png 775864840.233.png 775864840.234.png 775864840.235.png 775864840.236.png 775864840.237.png 775864840.238.png 775864840.239.png 775864840.240.png 775864840.241.png 775864840.242.png 775864840.244.png 775864840.245.png 775864840.246.png 775864840.247.png 775864840.248.png 775864840.249.png 775864840.250.png 775864840.251.png 775864840.252.png 775864840.253.png 775864840.255.png 775864840.256.png 775864840.257.png 775864840.258.png 775864840.259.png 775864840.260.png 775864840.261.png 775864840.262.png 775864840.263.png 775864840.264.png 775864840.266.png 775864840.267.png 775864840.268.png 775864840.269.png 775864840.270.png 775864840.271.png 775864840.272.png 775864840.273.png 775864840.274.png 775864840.275.png 775864840.277.png 775864840.278.png 775864840.279.png 775864840.280.png 775864840.281.png 775864840.282.png 775864840.283.png 775864840.284.png 775864840.285.png 775864840.286.png 775864840.288.png 775864840.289.png 775864840.290.png 775864840.291.png 775864840.292.png 775864840.293.png 775864840.294.png 775864840.295.png 775864840.296.png 775864840.297.png 775864840.299.png 775864840.300.png 775864840.301.png 775864840.302.png 775864840.303.png 775864840.304.png 775864840.305.png 775864840.306.png 775864840.307.png 775864840.308.png 775864840.310.png 775864840.311.png 775864840.312.png 775864840.313.png 775864840.314.png 775864840.315.png 775864840.316.png 775864840.317.png 775864840.318.png 775864840.319.png 775864840.321.png 775864840.322.png 775864840.323.png 775864840.324.png 775864840.325.png 775864840.326.png 775864840.327.png 775864840.328.png 775864840.329.png 775864840.330.png 775864840.332.png 775864840.333.png 775864840.334.png 775864840.335.png 775864840.336.png 775864840.337.png 775864840.338.png 775864840.339.png 775864840.340.png 775864840.341.png 775864840.343.png 775864840.344.png 775864840.345.png 775864840.346.png 775864840.347.png 775864840.348.png 775864840.349.png 775864840.350.png 775864840.351.png 775864840.352.png 775864840.354.png 775864840.355.png 775864840.356.png 775864840.357.png 775864840.358.png 775864840.359.png 775864840.360.png 775864840.361.png 775864840.362.png 775864840.363.png 775864840.365.png 775864840.366.png 775864840.367.png 775864840.368.png 775864840.369.png 775864840.370.png 775864840.371.png 775864840.372.png 775864840.373.png 775864840.374.png 775864840.376.png 775864840.377.png 775864840.378.png 775864840.379.png 775864840.380.png
Shannon reconstruction
Zero and rst order hold
X
[
f D kh E sin D
w
D t kh E 2 E
s
f D t E
w
D t kh E2
s
k [
The exact formula is not causal !
Approximation with information from
d
time steps ahead:
t 0
t 2
t 6
t 1
t 3
t 4
t 5
Time
E n d
X
w
t
2
E sin
f D nh
s
t
f D kh E h D nh
t
kh E
h D
t
w
t
2
s
k [
1
p
sin
t
p
t
0
t 0
t 2
t 6
Time
t 1
t 3
t 4
t 5
0
10
Time
Predictive rst order hold
Sinusoidal signal with h 1 and h 0
.
5
h=1
h=0.5
·
Use forward difference instead of backward difference
1
1
t h
u D kh
t
E u D kh E
D u D kh h E u D kh EE
0
0
0
5
10
0
5
10
·
Use model of the controller.
1
1
0
0
0
5
10
0
5
10
1
1
0
0
t 0
t 1
t 2
t 3
t 4
t 5
t 6
Time
0
5
10
0
5
10
Time
Time
&
Dept. of Automatic Control, Lund
4
775864840.381.png 775864840.382.png 775864840.383.png 775864840.384.png 775864840.385.png 775864840.387.png 775864840.388.png 775864840.389.png 775864840.390.png 775864840.391.png 775864840.392.png 775864840.393.png 775864840.394.png 775864840.395.png 775864840.396.png 775864840.398.png 775864840.399.png 775864840.400.png 775864840.401.png 775864840.402.png 775864840.403.png 775864840.404.png 775864840.405.png 775864840.406.png 775864840.407.png 775864840.409.png 775864840.410.png 775864840.411.png 775864840.412.png 775864840.413.png 775864840.414.png 775864840.415.png 775864840.416.png 775864840.417.png 775864840.418.png 775864840.420.png 775864840.421.png 775864840.422.png 775864840.423.png 775864840.424.png 775864840.425.png 775864840.426.png 775864840.427.png 775864840.428.png 775864840.429.png 775864840.431.png 775864840.432.png 775864840.433.png 775864840.434.png 775864840.435.png 775864840.436.png 775864840.437.png 775864840.438.png 775864840.439.png 775864840.440.png 775864840.442.png 775864840.443.png 775864840.444.png 775864840.445.png 775864840.446.png 775864840.447.png 775864840.448.png 775864840.449.png 775864840.450.png 775864840.451.png 775864840.453.png 775864840.454.png 775864840.455.png 775864840.456.png 775864840.457.png 775864840.458.png 775864840.459.png 775864840.460.png 775864840.461.png 775864840.462.png 775864840.464.png 775864840.465.png 775864840.466.png 775864840.467.png 775864840.468.png 775864840.469.png 775864840.470.png 775864840.471.png 775864840.472.png 775864840.473.png 775864840.475.png 775864840.476.png 775864840.477.png 775864840.478.png 775864840.479.png 775864840.480.png 775864840.481.png 775864840.482.png 775864840.483.png 775864840.484.png 775864840.486.png 775864840.487.png 775864840.488.png 775864840.489.png 775864840.490.png 775864840.491.png 775864840.492.png 775864840.493.png 775864840.494.png 775864840.495.png 775864840.497.png 775864840.498.png 775864840.499.png 775864840.500.png 775864840.501.png 775864840.502.png 775864840.503.png 775864840.504.png 775864840.505.png 775864840.506.png 775864840.508.png 775864840.509.png 775864840.510.png 775864840.511.png 775864840.512.png 775864840.513.png 775864840.514.png 775864840.515.png 775864840.516.png 775864840.517.png 775864840.519.png 775864840.520.png 775864840.521.png 775864840.522.png 775864840.523.png 775864840.524.png 775864840.525.png 775864840.526.png 775864840.527.png 775864840.528.png 775864840.530.png 775864840.531.png 775864840.532.png 775864840.533.png 775864840.534.png 775864840.535.png 775864840.536.png 775864840.537.png 775864840.538.png 775864840.539.png 775864840.541.png 775864840.542.png 775864840.543.png 775864840.544.png 775864840.545.png 775864840.546.png 775864840.547.png 775864840.548.png 775864840.549.png 775864840.550.png 775864840.553.png 775864840.554.png 775864840.555.png 775864840.556.png 775864840.557.png 775864840.558.png 775864840.559.png 775864840.560.png 775864840.561.png 775864840.562.png 775864840.564.png 775864840.565.png 775864840.566.png 775864840.567.png 775864840.568.png 775864840.569.png 775864840.570.png 775864840.571.png 775864840.572.png 775864840.573.png 775864840.575.png 775864840.576.png 775864840.577.png 775864840.578.png 775864840.579.png 775864840.580.png 775864840.581.png 775864840.582.png 775864840.583.png 775864840.584.png 775864840.586.png 775864840.587.png 775864840.588.png 775864840.589.png 775864840.590.png 775864840.591.png 775864840.592.png 775864840.593.png 775864840.594.png 775864840.595.png 775864840.597.png 775864840.598.png 775864840.599.png 775864840.600.png 775864840.601.png 775864840.602.png 775864840.603.png 775864840.604.png 775864840.605.png
Implementing a controller using a computer
H ( z)
»
G ( s)
u(t)
u kh
(
)
y kh
(
)
y (t)
{
}
{
}
A-D
Algorithm
D-A
Discrete time approximation of
continuous time controller
Clock
— Frequency domain
Want to get
A/D + Algorithm + D/A G D s E
Methods:
·
Approximate s i.e
G D s E 0 H D z E
·
Pole-zero matching
·
Use sample/hold to get
H D z E from
G D s E
Matlab
Approximation methods
Forward difference (Euler’s method)
SYSD = C2D(SYSC,TS,METHOD) converts the continuous
system SYSC to a discrete-time system SYSD with
sample time TS. The string METHOD selects the
discretization method among the following:
'zoh'
dx D t E
dt
x D t h E x D t E
h
q 1
h
x D t E
Backward difference
dx D t E
dt
Zero-order hold on the inputs.
'foh'
Linear interpolation of inputs
(triangle appx.)
1 q 1
h
x D t E x D t h E
h
x D t E
'tustin'
Bilinear (Tustin) approximation.
Trapezoidal method (Tustin, bilinear)
'prewarp'
Tustin approximation with frequency
prewarping.
The critical frequency Wc is specified
last as in C2D(SysC,Ts,'prewarp',Wc)
x D t h E x D t E
2
x D t h E x D t E
h
'matched'
Matched pole-zero method
(for SISO systems only).
&
Dept. of Automatic Control, Lund
5
775864840.606.png 775864840.608.png 775864840.609.png 775864840.610.png 775864840.611.png 775864840.612.png 775864840.613.png 775864840.614.png 775864840.615.png 775864840.616.png 775864840.617.png 775864840.619.png 775864840.620.png 775864840.621.png 775864840.622.png 775864840.623.png 775864840.624.png 775864840.625.png 775864840.626.png 775864840.627.png 775864840.628.png 775864840.630.png 775864840.631.png 775864840.632.png 775864840.633.png 775864840.634.png 775864840.635.png 775864840.636.png
Zgłoś jeśli naruszono regulamin