Measurement and control technology and instruments - APS
1 .Singal
1.1 Type
- Deterministic signal
- periodic signal
- aperiodic signal [ continuous frequency spectrum ]
- quasi-periodic signal : 2 and more frequency cosw0t+cos√2w0t
- transient signal : x(t)=x0e−atsin(w0t+φ)
- indeterministic [ random signal ]
- stationary random signal
- nonstationary
continuous s
discrete s
analog s : 独变+幅值为连续
digital s : 离散信号 , 幅值连续
energy s : ∫∞−∞x2(t)dt<∞
- power s : ∫∞−∞x2(t)dt→∞, 1t2−t1∫t2t1x2(t)dt<∞
1.2 Time-domain & frequency domain
- time domain : change with time
- frequency domain : frequency composition, amplitude, phase angle
1.3 Periodic square wave
1.4 Fourier series
use trigonometric function to approximate any signal
1.4.1 One side
x(t)=a0+∑∞n=1(ancosnw0t+bnsinw+0t)=a0+∑∞1sin(nw0t+φn)n=1,2,3,…
a0=1T∫T2−T2x(t)dt\an=2T∫T2−T2x(t)cosnw0tdt\bn=2T∫T2−T2x(t)sinnw0tdt
An=√a2n+b2ntanφn=anbn
1.4.2 Euler equation
e±jwt=coswt±jsinwt\coswt=12(ejwt+e−jwt)\sinwt=12j(ejwt−e−jwt)
1.4.3 Both side
x(t)=∑∞n=−∞cnejnw0t=∑∞−∞|cn|ej(nw0t+φ)n=0,±1,±2,…
- |cn|=An=√a2n+b2nvarphin=arctananbn
- cn=1T∫T2−T2x(t)e−jnw0tdt
- \int^a_af(x)dx=\left{\begin{array}{c,l}2\int^a_0f(x)dx,&偶\0,&奇\end{array}\right.
1.5 Periodic Function
1.5.1 Feature
- discrete
- have a base frequency
- amplitude is less when frequency goes higher
1.5.2 Strength
- average : μx=1T0∫T00x(t)dt
- root mean square value : xrms=√1T0∫T00x2(t)dt
- square value: Paj=x2rms=1T0∫T00x2(t)dt
1.6 Aperiodic
T0→∞ , frequency continus
1.7 Fourier transform
X(f)=∫∞−∞x(t)e−j2πftdt=12π∫∞−∞x(t)ejwtdt
x(t)=∫∞−∞X(f)ej2πftdf=∫∞−∞X(w)ejwtdw
- X(f)=|X(f)|ejφ(t)
- |X(f)| : continuous amplitude spectrum
- φ(f) continuous phase spectrum
1.7.1 Window function
w(t)=\left{\begin{array}{cl}1,&|t|<\dfrac{T}{2}\0,&|t|>\dfrac{T}{2}\end{array}\right.\Leftrightarrow W(f)=T\dfrac{sin\pi fT}{\pi f T}=Tsinc(\pi fT)
frequency spectrum : W(f)=∫∞−∞w(t)e−j2πftdt=1j2πf(ejπfT−ejπfT)Downarrowsin(πfT)=12j(ejπfT−e−πfT)\=TsinπfTπfT=Tsinc(πfT)
- |W(f)|=T|sinc(πfT)|
- only real part no virtual
- amplitude spectrum : |W(f)|=T|sinc(2πfT)|
1.7.2 F-transform princple
Superposition | ax(t)+by(t)↔aX(f)+bY(f) |
---|---|
scale change | x(kt)↔1kX(fk) |
time change | x(t−t0)↔X(f)e−j2πft0 |
frequency change | x(t)e∓j2πf0t↔X(f±f0) |
time domain convolution | x1(t)∗x2(t)↔X1(f)X2(f) |
frequency domain convolution | x1(t)x2(t)↔X1(f)∗X2(f) |
- convolution : ∫∞−∞x1(τ)x2(t−τ)dτ
1.8 δ function
Window function : \varepsilon\rightarrow 0\quad\Rightarrow\quad\delta (t)=\left{\begin{array}{cl}\infty,&t=0\0,&t\ne0\end{array}\right.
∫∞−∞δ(t)dt=limε→0∫∞−∞Sε(t)dt=1
1.8.1 Princple
sample:
∫∞−∞δ(t−t0)f(t)dt=∫∞−∞δ(t−t0)f(t0)dt=f(t0)∫∞−∞δ(t−t0)dt\=f(t0)
convolution
x(t)∗δ(t)=∫∞−∞x(τ)δ(t−τ)dτDownarrowevenfunction\=∫∞−∞x(τ)δ(t−τ)dτ=x(t)
x(t)∗δ(t±t0)=∫∞−∞x(τ)δ(t±t0−τ)dτ=x(t±t0)
1.8.2 Frequency specturm
Δ(f)=∫∞−∞δ(t)e−j2πftdt=e0=1delta(t)=∫∞−∞1ej2πftdf
Time D | Frequency D |
---|---|
δ(t) | 1 |
1 | δ(f) |
δ(t−t0) | e−j2πft0 |
ej2πf0t | δ(f−f0) |
1.9 Sine and cosine
sin2πf0t↔12j[δ(f−f0)−δ(f+f0)]cos2πf0t↔12[δ(f−f0)+δ(f+f0)]
- sin2πf0t=12j(ej2πf0t−e−j2πf0t)cos2πf0t=12(ej2πf0t+e−j2πf0t)
1.10 Comb function
comb(t,Ts)=∑∞−∞δ(t−nts)n=0,±1,±2…Updownarrow\comb(f,fs)=1Ts∑∞−∞δ(f−kfs)=1Ts∑∞−∞δ(f−kTs)
1.11 Random
x(t)=x1(t),x2(t).…xi(t)…
random progress
- stationary : characteristic parameter dont change with time
- nonstationary
characteristic parameter
Constant - average
μx=limT→∞1T∫T0x(t)dt
fluctuation - variance
σ2x=limT→∞1T∫T0[x(t)−μx]2dt
Strength - mean square value
ψ2x=limT→∞1T∫T0x2(t)dt
σ2x=ψ2x−μ2x
Probability density function
p(x)=limΔx→0Pr[x<x(t)≤x+Δx]Δx\=limΔx→0limT→∞TxTΔx
- Tx=n∑1Δti
- 正弦
autocorrelation function
power spectral density function
2. Testing device
2.1 Properities
Statics
device error
contain stable featuure
- Dynamic
- system parameter constant
- linear
- first=0,Laplace Transform H(s),Fourier Transform H(jw),h(t)=L−1[H(s)]
2.2 Static Feature
2.2.1 Linearity
definition : the difference between output and ideal
线性误差=ΔmaxYmax−Ymin×100%
2.2.2 Sensitivity
Definition : unit input change cause the output change →signal-to-noise ratio
灵敏度=ΔYΔx
2.2.3 hysterisis error
when increase and decrease the difference at same time
2.2.4 Resolution
smallest tangible change
分辨力=ΔXmax−Xmin×100%
2.3.5 Zero wander & sensitivity wander
2.3.6 Precision
- FS精度 : Δ误差Max Full Scale
- real precision : Δ误差测量值
2.4 Dynamic performance
- Transfer function H(s) [ complex domain ] →t=0sin function stimulate
- Frequency response function H(jw) [ frequency domain ] →steady state output
- impulse response function h(t) [ time domain ]
2.4.1 Series and parallel
- Series : H(s)=H1(s)H2(s)\A(w)=n∏Ai(w)varphi(s)=n∑1φi(s)
- parallel :H(w)=∑Hi(s)
2.4.2 First-order system
- H(s)=1τs+1
- H(jw)=1jτw+1{A(w)=1√1+(τw)2varphi(w)=−arctan(τw)
- h(t)=1τe−t2
- τ重要
2.4.3 Second-order system
- H(s)=w2ns2+2ζwns+w2n
- H(jw)=11−(wwn)2+2ζjwwn{A(w)=1√[1−(wwn)2]2+4ζ2(wwn)2varphi(w)=−arctan2ζ(wwn)1−(wwn)2
- h(t)=wn√1−ζ2e−ζwntsin√1−ζ2wnt0<ζ<1
- ζ,wn
2.4.4 Distortionless condition
satisfy y(t)=A0x(t−t0)
↓F
H(w)=A(w)ejφ(w)=Y(w)X(w)=A0e−jt0w
Amplitude distortion:A(w)=A0=C
Phase distortion : φ(w)=−t0w
2.4.5 Dynamic measure
the frequency response method : x(t)=X0sin2πft
- 一阶 : A(w)=1√1+(τw)2varphi(w)=−arctan(τw)
Step response methodu(t)
3. Sensor
3.1 Type
- mechanical
- resistance, capacitance , inductance
- Magnetoelectricity , piezoelectricity, thermoelectricity
- laser
3.2 Resistance
measure
strain, stress
force, displacement, pressure, accerleration
3.2.1 Rheostat
R=ρcA
Precision : S=dRdx
3.2.2 Resistance strain
3.2.2.1 Metal strain
deformation cause resistance valve change
material : constantan
dRR=ε(1+2υ+λE)≈(1+2υ)ε
- υ: Possion ratio
Sg=dR/Rdl/l=1+2υ
3.2.2.2 Semi-conductor gauge
piezoresistance : when force the resistance valve change
dRR=λEε§g=dR/Rdl/l=λE
sensitive, but affected by temperature
3.3 Capacitance
C=ε0εAδ
- ε: relative dielectric constant,aps ε=1
- ε0: dielectric constant in vaccum,ε0=8.85×10−12F/m
3.3.1 distance change
dC=−εε0A1δ2dδ§=dCdδ=−εε0A1δ2
unlinear
3.3.2 Area change
A=αr22\C=εε0αr22δ§=dCdα=ε0εr22δ
linear
3.4 Inductance
Variable Reluctance, self- inductance
L=N2Rm=N2μ0A02δ§=N2μ0A02δ2
- N: number of turns
- Rm=2δμ0A0 : sum
- μ0: magnetic permeability,4πtimes10−3H/m
3.5 Magnetoelectricity
Moving-coil
- speed : e=NBlvsinθ
- angle speed : e=kNBAw
3.6 Piezoelectricity
Measure: pressure, stress, acceleration
Reversible : mechanical ↔ electric
3.6.1 Principle
- piezoelectric effect : piezoelectric material generate electric field when pressed
- inverse piezoelectric effect : in electric field size change
3.6.2 Material
- piezoelectric monocrystal
- piezoceramics
3.6.3 Sensitive coefficient
C=εε0Aδ
- ε=4.5F/m,磺
- ε=1200F/m,钛酸铝
3.7 Thermoelectricity
Thermocouple: temperature difference cause electromotive force
thermal resistance : resistance value change with T
3.8 Photoelectricity
- Outside : light on , electronic out
- Inside : light , R change
3.9 Semiconductor
Magneto-dependent sensor :
Hall effect
thermosensitive
3.10 Choose
4. Signal Conditioning
4.1 Bridge
4.1.1 Direct Current Brigde
U0=(R1R1+R2−R4R3+R4)Ue
Signal arm : R1
U0=(R1+ΔRR1+R2+ΔR−R4R3+R4)Uedownarrow等阻\=ΔR2(2R+ΔR)Ue≈ΔR4RUe
Half bridge : R1+R2
U0=(R1+ΔRR1+R2−R4R1+R4)Ue→等阻=ΔR2RUe
Full bridge : ∑Ri
U0=(R1+ΔRR1+R2−R4−ΔRR3+R4)Ue→等阻=ΔRRUe
sensitivityS=U0ΔR/R
- signal arm : Ue4
- half bridge : Ue2
- Full bridge : Ue
4.1.2 AC bridge
4 arm can + L/R/C
Z01Z03=Z02Z04varphi1+φ3=φ2+φ4
4.2 Modulation & demodulation
Modulation : use low frequency signal to control amplitude or frequency of oscillator signal
demodulation : recover the original signal from the modulated signal
Utilize
- denoising
- long distance transmission
4.2.1 Amplitude modulation
调幅=高频载波⋅被测信号\m(t)=x(t)⋅cos2πf0tdownarrowFfrac12X(f+f0)+12X(f−f0)
- f0muss>max(x(t)),or 不重叠
4.2.2 Amplitude demodulation[Detection]
synchronously demodulation
x(t)cos2πf0tcos2πf0t=x(t)2+12x(t)cos4πf0t
Envelope detection
4.2.3 Frequency modulation
low amplitude change with the high frequency signal
LC oscillating circuit
f0=12π√LC0 x(t)=Acos[w0t+k∫x(t)dt+θ0]
Demodulation Frequency discrimination : high pass filter + envelope detection
4.3 Filter
- 低通Low-Pass
- 高通High-Pass
- 带通Band-Pass=H+L
- 陷波/带阻Band-Stop/Notch=H//L
4.3.1 Parameter
Ideal
Real
- Cutoff frequency fc1,fc2 : half power point,A=A0√2,[−3dB=20log(1√2)]
- Bandwidth B=fc2−fc1 : [-3dB bandwidth],B↓,discrimination↑
- range±δ : small the best best
- Quality coefficient : Q=f0B,Q↑,choosing thebz best
- 倍频程选择性 :
- up : |A(fc2)−A(2fc2)|
- down : |A(fc1)−A(fc12)|
- fast the best
- filter coeffecient : λ=B−60dBB−3dB
- ideal=1,normal 1-5
4.3.2 Real filter circuit
low -pass
|H(f)|=1√1+(f/fc)2phi(f)=−arctan(ffc)
- fc=12πRC
high pass
|H(f)|=(f/fc)√1+(f/fc)2phi(f)=π2−arctan(ffc)
4.3.3 Band pass
constant bandwidth ratio
Bifoi=fc2i−fc1if0i=C\fc2i=2nfc1i
- n octave : n=1,octave,n=1/3,1/3 time octave
- center : foi=√fc1ifc2i
constant bandwidth
B=fc2i−fc1i=C
5. Signal processing
5.1 Sampling theorem
sampling frequency is twice time of fh
fs>2fh(3∼4)
no mix and overlap between the signal
5.2 Quantization x(t)s(t)
use one finite level to similarity the real
A/D transfer
5.3 Window functionx(t)s(t)w(t)
- cut off : signal ⋅window function(time)
- give away : W(f) infinite bandwidth sinc function,x(t)带限信号→截断→无限带宽
5.4 Frequency Sampling[x(t)s(t)w(t)]∗d(t)
pulse D(f)⋅signal frequency spectrum→时域窗内信号,窗外周期延拓
5.6 Correlation analysis
5.6.1 correlated coefficient
ρxy=E[(x−μx)(y−μy)]σxσy=E(XY)−E(X)E(Y)\=±12[D(X±Y)−D(X)−D(Y)]
- σ : 标准差
5.6.2 Autocorrelation function
Rx(τ)=limT→∞1T∫T0x(t)x(t+τ)dt
- Autocorrelation coefficient : ρx(τ)=Rx(τ)−μ2xσ2x
- μ2x−σ2x≤Rx(τ)≤μ2x+σ2x
- Rx(τ)max=Rx(0)=ψ2x
- τ→∞,ρx(τ)→0/μ2x,x(t)和x(t+tau)无内部联系
- Rx(τ)=Rx(−τ), 偶函数
- 周期函数的自相关函数仍为同频率周期函数,
[Exp]x(t)=x0sin(wt+φ)→Rx(τ)=x202coswτ
5.6.3 Cross-correlation function
Rxy(τ)=limT→∞1T∫T0x(t)y(t+τ)dt
- τ→∞,ρxy→0,Rxy(τ)→μxμy,x(t),y(t)不相关
- μxμy−σxσy≤Rxy(τ)≤μxμy+σxσy
- 非偶
[Exp]x(t)=x0sin(wt+φ)\y(t)=y0sint(wt+θ)→Rxy(τ)=12x0y0cos(wt+φ−θ)
5.7 Power spectrum
5.7.1 Auto power spectrum density
Sx(f)−=∫∞∞Rx(τ)e−j2πfτdτ\Rx(τ)=∫∞∞Sx(f)ej2πfτdf
5.7.2 Parseval theorem
Energy equation : ∫∞∞x2(t)dt=∫∞∞|X(f)|2df
[Exp]
- Sx(f)=limT→∞1T|X(f)|2df
- Y(f)=H(f)X(f)§y(f)=|H(f)|2Sx(f)§xy(f)=H(f)Sx(f)
5.7.3 Cross power spectrum
Sxy(f)=∫∞∞Rxy(τ)e−j2πfτdτupdownarrow\Rxy(τ)=∫∞∞Sxy(f)ej2πfτdf
Sxy(f)=H(f)Sx(f)
7. Vibration
7.1 Methode
- Light : 振动量to光信号
- 光学读数显微镜测振
- 激光干涉法测振
- Electric : 振动量to电量
- 频率范围宽,动态范围广,测量灵敏
7.2 Excitation
- Sin: 广
- random : 带宽,白噪声
- transient : 宽频带
- fast sin scan
- impact hammer
- step excitation
7.2.1 Exciter
excitation → object → forced vibration
- hammer
- electric exciter
7.2.2 Sensor
Vibration →electric quantity
- 惯性式
- 相对式
- touch
- untouched
- 涡流位移
- 电容加速度
- 磁电加速度
- piezoelectric accelerometer
- impedance head
7.3 Signal analysis
- 振动仪
- frequency analyzer
- 频率特性分析仪
- digital signal processing system