![Page 1: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/1.jpg)
AOE 3054. Statistical Analysis of Experimental Data
• Relevant to all experiments, IDLab and uncertainty analysis
• See lab manual, “Basic Concepts in Experiments”
1. Mean, Variance, Histograms2. Probability density functions and the normal
distribution3. Correlation and Regression
![Page 2: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/2.jpg)
Analysis?
Separation
Transition to turbulence
Slowing of flow around stagnation point
Stall zone
Wake closes
![Page 3: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/3.jpg)
Example
USS Jefferson City (SSN 759) navigates past Orote Point as the submarine enters Apra Harbor, Guam.
![Page 4: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/4.jpg)
Wind Tunnel Model
![Page 5: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/5.jpg)
Typical Velocity Signal
![Page 6: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/6.jpg)
Mean, Variance, HistogramsDefinitions
Consider a series of measurements x1, x2, x3 … xN,
Mean
Variance
Standard deviation
For a continuous function x(t), where t could be anything ...
Mean
Variance
Standard deviation
![Page 7: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/7.jpg)
Mean, Variance, HistogramsDefinitions
Histogram: A chart showing the distribution of values in a series of samples
• Divide up range of the measurements into a series of equal intervals (called bins)• Add up the number of measurements falling in each bin• The mean usually falls near the center of the histogram• The standard deviation is typically 1/4 to 1/6th of the spread• Histograms may be used to estimate probability
05
101520253035404550
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
Wave height (ft)N
umbe
r of s
ampl
es
![Page 8: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/8.jpg)
Typical Velocity Signal
![Page 9: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/9.jpg)
Mean Velocity on Center-plane
‘Hull’
‘Con
ning
Tow
er L
eadi
ng E
dge’
![Page 10: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/10.jpg)
Velocity Variance
‘Hull’
‘Con
ning
Tow
er L
eadi
ng E
dge’
![Page 11: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/11.jpg)
Velocity
Num
ber o
f sam
ples
Velocity0-.5U∞
Velocity Histograms From in the Vortex
0-.5U∞
![Page 12: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/12.jpg)
USS Seawolf (SSN 21), completing its initial sea trial
![Page 13: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/13.jpg)
Probability Density FunctionsDefinition
Definition: The area under p(x) between two values x0 and x1 is the probability P that a given sample of x will fall between these values.
Mathematically:
(Add labels)
![Page 14: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/14.jpg)
Normal DistributionDefinition
The probability P that a given sample of x will fall between x0 and x1 is:
Which can be re-written as
Where
0
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0 5 10 15 20 25
p(x)
x
πσ=
σ−
−
2)(
2
2
2)( xx
exp
x
σ
= Matlab Excel
![Page 15: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/15.jpg)
Normal DistributionExample of Use
A sensor is used to detect the flow rate of fuel to a jet engine. The flow rate fluctuates and the following are 21 successive readings readings (in arbitrary units),
(a) Determine the mean and standard deviation.
(b) Estimate the probability that a flow rate will have a value between .5 and .6.
(c) What percentage of the time is the flow rate likely to lie above a value of .8?
(d) What percentage of a large number of readings is likely to lie further than two standard deviations from the mean?
Reading Flow rate1 0.5122 0.4773 0.7944 0.6725 0.7136 0.5887 0.6218 0.7349 0.77110 0.48611 0.55912 0.61413 0.68714 0.72215 0.62716 0.70117 0.57318 0.72119 0.80220 0.55321 0.605
(a) Determine the mean and standard deviation644.0...)794.477.512(.21
1 =+++=x
098.0...))644.477(.)644.512((. 22201 =+−+−=xσ
![Page 16: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/16.jpg)
Normal DistributionExample of Use
(b) Estimate the probability that the flow rate will have a value between .5 and .7.
x = 0.644σ = 0.098
Reading Flow rate1 0.5122 0.4773 0.7944 0.6725 0.7136 0.5887 0.6218 0.7349 0.77110 0.48611 0.55912 0.61413 0.68714 0.72215 0.62716 0.70117 0.57318 0.72119 0.80220 0.55321 0.605
⎟⎟⎠
⎞⎜⎜⎝
⎛ −−⎟⎟
⎠
⎞⎜⎜⎝
⎛ −
=<<
σσxxIxxI
xxxP
01
10 )(
xo x1
![Page 17: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/17.jpg)
Normal DistributionExample of Use
(c) What percentage of the time is the flow rate likely to lie above a value of .8?
x = 0.644σ = 0.098
Reading Flow rate1 0.5122 0.4773 0.7944 0.6725 0.7136 0.5887 0.6218 0.7349 0.77110 0.48611 0.55912 0.61413 0.68714 0.72215 0.62716 0.70117 0.57318 0.72119 0.80220 0.55321 0.605
⎟⎟⎠
⎞⎜⎜⎝
⎛ −−⎟⎟
⎠
⎞⎜⎜⎝
⎛ −
=<<
σσxxIxxI
xxxP
01
10 )(
xo x1
![Page 18: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/18.jpg)
Normal DistributionExample of Use
(d) What percentage of a large number of readings will lie further than two standard deviations from the mean?
x = 0.644σ = 0.098
Reading Flow rate1 0.5122 0.4773 0.7944 0.6725 0.7136 0.5887 0.6218 0.7349 0.77110 0.48611 0.55912 0.61413 0.68714 0.72215 0.62716 0.70117 0.57318 0.72119 0.80220 0.55321 0.605
⎟⎟⎠
⎞⎜⎜⎝
⎛ −−⎟⎟
⎠
⎞⎜⎜⎝
⎛ −
=<<
σσxxIxxI
xxxP
01
10 )(
xo x1
So, on average, there is close to a 95% chance that the value will lie within 2 standard deviations of the mean. This result is important when we come to estimating experimental error.
![Page 19: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/19.jpg)
Correlation (co-variance)Definition
Consider two series of measurements of related quantities x1, x2, x3 … xN, y1, y2, y3 … yN,
Correlation
For continuous functions x(t) and y(t), where t could be anything ...
Correlation
Correlation coefficientVaries from -1 to 1
We also define the Correlation function
where τ is a shift in time, distance of whatever t represents, and the
![Page 20: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/20.jpg)
Correlation (co-variance)Application Example
• GE 90 Aircraft Engine• Eddies from fan tip entering stator row impinge in stators generating noise
![Page 21: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/21.jpg)
∆y
∆z
τ
Zero time-delay correlation function, USpace Time Correlations∆y
/ca
Fan blade wakeFan blade wake
Casing
Fan blade tip vortices
Correlation with this point
Correlation coefficient
![Page 22: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/22.jpg)
-1
-0.5
0
0.5
1
-τUe
/ c a
-0.4-0.3
-0.2-0.1
00.1
0.20.3
0.4
∆z / ca
0
0.1
0.2
0.3
0.4
∆y/ c
a
LSE of coherent structure based onU fluctuation.
Vorticity magnitude isosurface, andstreamtraces.
-τUe/ca
∆x3/ca
∆x2/c
a
Streamtraces
Vorticity isosurfaces
Note dominant elongated vortical structure with axis about 20o from the
time axis (i.e. the flow direction)
Reconstruction of 3D field from correlations
Linear stochastic estimate of flowfield around location 2
![Page 23: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/23.jpg)
Linear RegressionConsider a series of measured points (xi,yi) describing a relationship between two quantities x and y. We wish to find a straight line of the form y = A+Bxthat lies as closely as possible to the data. That is, we wish to choose A and Bto minimize the error
This is done by setting
Which gives us two simultaneous equations for A and B, which may be solved to yield
![Page 24: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/24.jpg)
Linear Regression• Using Excel:
– Under Tools>Data Analysis… (may have to add this feature by going to Tools>Add-ins… and checking Analsys Toolpak
– Alternatively …
• Using Matlab (preferred)
![Page 25: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/25.jpg)
Matlab Example
⎟⎟⎠
⎞⎜⎜⎝
⎛
⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜
⎝
⎛
=
⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜
⎝
⎛
BA
xxx
yyy
MMM3
2
1
3
2
1
111
1 2 3 4 5 6-5
0
5
10
15
20
25
30
35
PointsFit
![Page 26: AOE 3054. Statistical Analysis of Experimental Data](https://reader031.vdocument.in/reader031/viewer/2022012411/616b0e74a46a9a278329b0e3/html5/thumbnails/26.jpg)
Matlab Example – Generalized Regression
⎟⎟⎟
⎠
⎞
⎜⎜⎜
⎝
⎛
⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜
⎝
⎛
=
⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜
⎝
⎛
CBA
xxxxxx
yyy
MMMM
233
222
211
3
2
1
111
1 2 3 4 5 60
5
10
15
20
25
30
35
40PointsFit