How do I plot a time-dependent probability density function on MATLAB?

Click For Summary
SUMMARY

The discussion focuses on plotting a time-dependent probability density function (PDF) for a wave function in MATLAB, specifically for a free particle represented by the equation P(x,t)=\frac{\sigma_0}{\mid \alpha \mid^2\sqrt{\pi}}exp(-(\frac{\sigma_0}{\mid \alpha \mid})^4\frac{(x-x_0-p_0t/m)^2}{\sigma_0^2}). The user seeks guidance on how to visualize this PDF over the interval -L/2 to L/2, considering its dependence on time. Suggestions include creating a vector for the spatial domain and selecting specific time points to illustrate key characteristics, or alternatively, developing an animation to depict the time evolution of the PDF.

PREREQUISITES
  • Understanding of wave functions and probability density functions in quantum mechanics
  • Familiarity with MATLAB for plotting and data visualization
  • Knowledge of basic mathematical concepts such as exponential functions and Gaussian distributions
  • Ability to manipulate vectors and matrices in MATLAB
NEXT STEPS
  • Learn how to implement MATLAB's plotting functions for 2D graphs
  • Explore MATLAB's capabilities for creating animations with the 'plot' and 'pause' functions
  • Study the mathematical properties of Gaussian functions and their applications in quantum mechanics
  • Investigate the use of MATLAB's 'meshgrid' function for visualizing multi-dimensional data
USEFUL FOR

Students and researchers in quantum mechanics, MATLAB users interested in data visualization, and anyone looking to understand time-dependent probability density functions in wave mechanics.

Elekko
Messages
14
Reaction score
0

Homework Statement



Given that, in free space the probability density for a wave function (free particle) is \mid \Psi(x,t)\mid^2=P(x,t)=\frac{\sigma_0}{\mid \alpha \mid^2\sqrt{\pi}}exp(-(\frac{\sigma_0}{\mid \alpha \mid})^4\frac{(x-x_0-p_0t/m)^2}{\sigma_0^2})

What is need to be done is to study a wave packet of size \sigma_0=\frac{L}{10} at x_0=0 on the interval -L/2\leq x \leq L/2

I need to plot this probability density function on MATLAB.
Since this PDF is time dependent I actually don't know how I should plot it? I know the interval should go between -L/2 and L/2 but how should this depend on time?

The Attempt at a Solution


Should I just create a vector from -L/2 to L/2 and then just have an arbitrary time?
Basically I cannot "visualize" in mind how it should look like before plotting, as I don't know how it exactly has to be done?
 
Physics news on Phys.org
That's right - you pick a time that shows the important characteristics of the density function.
It's either that or create an animation.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 21 ·
Replies
21
Views
2K
Replies
10
Views
3K
  • · Replies 9 ·
Replies
9
Views
3K
Replies
12
Views
3K
  • · Replies 13 ·
Replies
13
Views
4K
Replies
1
Views
2K