matlab在量子力学中的应用,一个关于量子力学中的matlab的问题

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

syms E U x U1 a a1 a2 k k1 k2 k3 b c d e f g h l S

%其中，b c d e f g h l是为计算简明而引入的字符，S为透射系数D与反射系数R之和。

E=input('请输入入射电子能量E=');

U=input('请输入势垒函数U=');

a1=input('请输入势垒起始坐标a1=');

a2=input('请输入势垒终止坐标a2=');

u=9.10908e-031;%u为电子质量

h=1.0545e-34; %h为约化普朗克常数

k1=sqrt(2*u*E)/h;

a=a2-a1; %a为势垒宽度

k=sqrt(2*u*(E-U));

k2=(int(k,a1,a2))/(h*(a2a1));

k3=-1i*k2;

U1=int(U,a1,a2)/(a2-a1); %势垒近似平均值

if numeric(E)>numeric(U1)

b=4*k1^2*k2^2;

c=(k1^2-k2^2)^2*(sin(a*k2))^2;

d=b+c;

D=b/d;

R=c/d;

% D = ( 4 * k 1 ^ 2 * k 2 ^ 2 ) / ( ( k 1 ^ 2 -k2^2)^2*(sin(a*k2))^2+4*k1^2*k2^2);

% R=(((k1^ 2-k 2 ^ 2)^ 2)*(s i n(a*k 2))^ 2))/((k1^ 2-k2^2)^2*(sin(a*k2))^2+4*k1^2*k2^2);

else

if numeric(E)==numeric(U1)

e=a^2*k1^2;

f=4+a^2*k1^2;

D=4/f;

R=e/f;

%D=4/(4+a^2*k1^2);

%R=(k1^2*a^2)/(4+a^2*k1^2);

else

g=4*(k1^2)*(k3^2);

h=((k1^2+k3^2)^2)*(sinh(a*k3))^2;

l=g+h;

D=g/l;

R=h/l; %D=(4*k1^2*k3^2)/((k1^2+k3^2)^2*((sinh(a*k3))^2)+4*k1^2*k3^2; %R=((k1^2+k3^2)^2*(sinh(a*k3))^2)/((k1^2+k3^2)^2*((sinh(a*k3))^2)+4*k1^2*k3^2);

end

end

D=numeric(D);

R=numeric(R);

S=D+R;