Poisson distribution (radioactive decay)

Click For Summary
SUMMARY

The discussion focuses on fitting experimental data of beta particle counts to a Poisson distribution using MATLAB. The user encountered issues with the probability equation, specifically the behavior of the function P(k) = e^-λ * λ^K/K!. It was clarified that the exponential decay component (e^-λ) is crucial for the distribution to curve downwards for large K values, and the user was advised to check for coding errors that may have led to incorrect graphing results. Ultimately, the problem was identified as a coding mistake, emphasizing the importance of careful implementation in statistical modeling.

PREREQUISITES
  • Understanding of Poisson distribution and its application in radioactive decay.
  • Proficiency in MATLAB for data analysis and visualization.
  • Familiarity with Excel for initial data processing and chart creation.
  • Basic knowledge of probability theory and statistical functions.
NEXT STEPS
  • Learn how to implement Poisson distribution fitting in MATLAB.
  • Explore the use of MATLAB's built-in functions for statistical analysis.
  • Study the significance of the exponential decay factor in probability distributions.
  • Investigate common coding errors in MATLAB that affect data visualization.
USEFUL FOR

Physics students, data analysts, and anyone involved in statistical modeling of experimental data, particularly in the context of radioactive decay and Poisson distributions.

steamyoshi
Messages
5
Reaction score
0

Homework Statement


I am a freshman in physics, just done a lab about radioactive decay.
I've measured the # of beta particles per second 400 times and got the frequency of each number K using Excel.
I'm supposed to take the data and fit it to a puason distribution in MATlab.
The data points themselves seem to be on a nice curve.
The problem I'm having is that the equation for the probability doesn't seem to curve at all.

Homework Equations


Puason distribution in radioactive decay: the chance that K beta particles will be detected in 1 second is
P(k) = e^-λ * λ^K/K!


The Attempt at a Solution


I've made an excel chart which shows λ^K/K! for many different K values, and then tweaked the lambda value.
For all the values of λ I've tried, λ^K/K! always increases for increasing values of K, meaning the graph never curves back down. What am I doing wrong?
 
Last edited:
Physics news on Phys.org
Puason distribution in radioactive decay: the chance that K beta particles will be detected in 1 second is
P(k) = e^-λ * λ^K/K!
... that would be poisson distribution, like the fish.
Careful to keep the same variables through your expressions.

$$P(k)=\frac{\lambda^k e^{-k}}{k!}$$

I've made an excel chart which shows λ^K/K! for many different K values, and then tweaked the lambda value.
For all the values of λ I've tried, λ^K/K! always increases for increasing values of K, meaning the graph never curves back down. What am I doing wrong?
##\lambda^k/k!## does increase for positive k, if ##\lambda > 1##but the poisson distribution function has a negative exponential in it which makes it converge for large k.
 
Last edited:
I'm sorry, I don't understand why e is ^-k, all the formulae I've found have e^-λ, which is a constant
thanks for the reply
 
That's because I'm an idiot... I should know better than to answer questions at 2am.
concentrating:

When I plot P vs k, I get a decreasing exponential for small values of lambda, and an approximate gaussian for large values of lambda.

You probably have a mistake in your code.
Check - sounds like a misplaced minus sign.

Time for bed.
 
Thanks, it WAS a code problem.
 
No worries.
And I got to demonstrate not to be afraid of making dumb mistakes too :)
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
7
Views
1K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 6 ·
Replies
6
Views
6K
  • · Replies 7 ·
Replies
7
Views
3K