Gaussian Wavepacket Tunnel Through double barrier

[mehdikhann]

Hello,
I was trying to design a movie in Matlab. A Gaussian wave packet moves toward double barrier then some wave reflect and some pass out from the barrier . i design the double barrier in Matlab. but don't know how to make movie in Matlab with Gaussian wavepacket. can you help. Do u have the code in Matlab

 

[Pythagorean]

Hello,
I was trying to design a movie in Matlab. A Gaussian wave packet moves toward double barrier then some wave reflect and some pass out from the barrier . i design the double barrier in Matlab. but don't know how to make movie in Matlab with Gaussian wavepacket. can you help. Do u have the code in Matlab

There is a programming subforum down below.

Short answer: the plot command in a for loop. You plot your function (the gaussian wave packet) and then compute the new shape of the function every time through the for loop, and then plot it.

You can use the pause command (type help pause) to put time between your frame changes.

 

[Pythagorean]

also check out 'help movie'. You can collect each of your plots as a movie frame and build a movie out of it, then once you have the movie array, you can adjust the frame-rate and the number of plays as you want.

 

[mehdikhann]

Hej Pythoagorean,
Thankyou very much for reply. i use the movie frame and it work for some equation. but i was trying move to gaussian wave packet but don't know how implement on MATLAB mathematically . do u have code for Gaussian wave packet

 

[Pythagorean]

Here's the general analytic structure:

http://musr.physics.ubc.ca/~jess/hr/skept/GWP/

Shooting from the hip, I'd use a runge-kutta algorithm to handle the integral, but you may want to see if there's some more experienced replies. I've never coded a gaussian wave packet before.

 

[f95toli]

A split-operator FFT metthod is usually more efficient than solving the integral directly for these kinds of problems.
You should be able to find a description of this method in just about any book on computational physics.

 

[conway]

If I'm understanding these suggestions correctly, then what is being advised is to model the Gaussian wave packet as the superposition of its component frequencies. While this is certainly correct, it's probably overkill for the animation in question. I'm pretty sure you'll get almost the identical picture simply by calculating the reflection/transmission of a monochromatic beam and then just putting a Gaussian envelope over it.