Unable to plot a wavefunction
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Sergio Manzetti
am 3 Jan. 2018
Kommentiert: Sergio Manzetti
am 4 Jan. 2018
Hi, I wonder if any of you can try to plot this and see if they also get a weird bar?
syms h g x C t m e p c
h = 1.0545718E-34
g = 5.344285879E-28
E = 1
Z = - (exp(-(x*(g*1i + (- 2*g^2 + E)^(1/2)))/h)*(g*1i - (- 2*g^2 + E)^(1/2)))/(2*(E - 2*g^2)^(1/2)) + (exp(-(x*(g*1i - (- 2*g^2 + E)^(1/2)))/h)*(g*1i + (- 2*g^2 + E)^(1/2)))/(2*(E - 2*g^2)^(1/2))
fplot(real(Z), [ -4, 4 ])
figure(1)
Thanks!
If so, why?
0 Kommentare
Akzeptierte Antwort
David Goodmanson
am 4 Jan. 2018
Hi Sergio,
This function appears to have scaling issues. If you plot
figure(1)
fplot(real(Z), [-2*h, 2*h])
figure(21)
fplot(imag(Z), [-2*h, 2*h])
you can see the function. I believe this addresses the question, but there are other things going on too. Possibly you were using E=1 just for test purposes, but compared to that value, 2*g^2 is infinitesimal. Also the function itself is an increasing exponential which is not good for a wave function unless there is some kind of hard boundary.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Line Plots finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!