Unable to plot a wavefunction

Question

Sergio Manzetti 2018-1-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/375391-unable-to-plot-a-wavefunction

评论： Sergio Manzetti 2018-1-4

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 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

David Goodmanson 2018-1-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/375391-unable-to-plot-a-wavefunction#answer_298641

在 MATLAB Online 中打开

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.