Engineering Probability Class 13 Thurs 2021-03-11
Table of contents::
1 Matlab
Matlab, Mathematica, and Maple all will help you do problems too big to do by hand. Sometime I'll demo one or the other.
-
Matlab
-
Major functions:
cdf(dist,X,A,...) pdf(dist,X,A,...)
-
Common cases of dist (there are many others):
'Binomial' 'Exponential' 'Poisson' 'Normal' 'Geometric' 'Uniform' 'Discrete Uniform'
-
Examples:
pdf('Normal',-2:2,0,1) cdf('Normal',-2:2,0,1) p=0.2 n=10 k=0:10 bp=pdf('Binomial',k,n,p) bar(k,bp) grid on bc=cdf('Binomial',k,n,p) bar(k,bc) grid on x=-3:.2:3 np=pdf('Normal',x,0,1) plot(x,np)
Interactive GUI to explore distributions: disttool
-
Random numbers:
rand(3) rand(1,5) randn(1,10) randn(1,10)*100+500 randi(100,4)
Interactive GUI to explore random numbers: randtool
-
Plotting two things at once:
x=-3:.2:3 n1=pdf('Normal',x,0,1) n2=pdf('Normal',x,0,2) plot(x,n1,n2) plot(x,n1,x,n2) plot(x,n1,'--r',x,n2,'.g')
-
Use Matlab to compute a geometric pdf w/o using the builtin function.
-
Review. Which of the following do you prefer to use?
Matlab
Maple
Mathematica
Paper. It was good enough for Bernoulli and Gauss; it's good enough for me.
Something else (please email about it me after the class).
1.1 My opinion
This is my opinion of Matlab.
-
Advantages
Excellent quality numerical routines.
Free at RPI.
Many toolkits available.
Uses parallel computers and GPUs.
Interactive - you type commands and immediately see results.
No need to compile programs.
-
Disadvantages
Very expensive outside RPI.
Once you start using Matlab, you can't easily move away when their prices rise.
You must force your data structures to look like arrays.
Long programs must still be developed offline.
Hard to write in Matlab's style.
Programs are hard to read.
-
Alternatives
Free clones like Octave are not very good
The excellent math routines in Matlab are also available free in C++ librarues
With C++ libraries using template metaprogramming, your code looks like Matlab.
They compile slowly.
Error messages are inscrutable.
Executables run very quickly.
2 Chapter 4 ctd
Text 4.2 p 148 pdf
Simple continuous r.v. examples: uniform, exponential.
-
The exponential distribution complements the Poisson distribution. The Poisson describes the number of arrivals per unit time. The exponential describes the distribution of the times between consecutive arrivals.
The exponential is the continuous analog to the geometric. If the random variable is the integral number of seconds, use geometric. If the r.v. is the real number time, use exponential.
Ex 4.7 p 150: exponential r.v.
-
Properties
Memoryless.
\(f(x) = \lambda e^{-\lambda x}\) if \(x\ge0\), 0 otherwise.
Example: time for a radioactive atom to decay.
Skip 4.2.1 for now.
The most common continuous distribution is the normal distribution.
4.2.2 p 152. Conditional probabilities work the same with continuous distributions as with discrete distributions.
-
p 154. Gaussian r.v.
\(f(x) = \frac{1}{\sqrt{2\pi} \cdot \sigma} e^{\frac{-(x-\mu)^2}{2\sigma^2}}\)
cdf often called \(\Psi(x)\)
-
cdf complement:
\(Q(x)=1-\Psi(x) = \int_x^\infty \frac{1}{\sqrt{2\pi} \cdot \sigma} e^{\frac{-(t-\mu)^2}{2\sigma^2}} dt\)
-
E.g., if \(\mu=500, \sigma=100\),
P[x>400]=0.66
P[x>500]=0.5
P[x>600]=0.16
P[x>700]=0.02
P[x>800]=0.001
Text 4.3 p 156 Expected value
Skip the other distributions (for now?).
3 Examples
4.11, p153.
4 4.3.2 Variance
p160
6 4.4.3 Normal (Gaussian) dist
p 167.
Show that the pdf integrates to 1.
Lots of different notations:
Generally, F(x) = P(X<=x).
For normal: that is called $\Psi(x)$ .
$Q(x) = 1-\Psi(x)$ .
Example 4.22 page 169.
7 4.4.4 Gamma r.v.
2 parameters
Has several useful special cases, e.g., chi-squared and m-Erlang.
The sum of m exponential r.v. has the m-Erlang dist.
Example 4.24 page 172.