How to plot the Free particle wave function

5 次查看（过去 30 天）
Pouyan Msgn
Pouyan Msgn 2021-4-20
评论： Pouyan Msgn 2021-4-20
Hello!
I have the function:
And a plot from my book is :
As you see in the picture, as a function of x/a is ploted. I tried to program this in Matlab with optional parameters but I dont get the correct picture
My code is:
clc
clear all
a=10;
k=-100:1:100;
c1=pi*sqrt(2*a);
h=(6.626/(2*pi))*10^-34;
m=9.11*10^-31;
N=length(k);
xb=60;
h=2*xb/(N-1);
x=-xb:h:xb;
t=m*a*a/h;
V=zeros(N,1);
S=0;
f=@(k,x)(1/c1)*((sin(k*a)./k)*exp(i*(k.*x-(h/(2*m))*t.*(k.^2))));
for i=1:length(k)-1
c=f(k(i),x(i))+f(k(i+1),x(i));
S=S+c*0.5;
V(i)=S;
end
Ps=abs(V).^2;
plot(x/a,a*Ps)
grid on
The plot is :
How can I get the correct plot?!

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

采纳的回答

David Hill
David Hill 2021-4-20
a=10;
c1=pi*sqrt(2*a);
h=(6.626/(2*pi))*10^-34;
m=9.11*10^-31;
X=-60:.1:60;
t=m*a^2/h;
V=zeros(1,201);
for j=1:length(X)
f=@(k)((sin(k*a)./k).*exp(1i*(k.*X(j)-(h/(2*m))*t.*(k.^2))));
V(j)=integral(f,-100,100);
end
plot(X/a,a*abs(V).^2);
  2 个评论
Pouyan Msgn
Pouyan Msgn 2021-4-20
This is the only message I get :
The integral may not exist, or it may be difficult to approximate numerically to the requested accuracy.
> In integralCalc/iterateScalarValued (line 372)
In integralCalc/vadapt (line 132)
In integralCalc (line 75)
In integral (line 88)
In QM3 (line 10)
Pouyan Msgn
Pouyan Msgn 2021-4-20
It worked for me in intervall -10 to 10! thank you

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by