該樓層疑似違規已被系統折疊?隱藏此樓查看此樓

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;

總結

以上是生活随笔為你收集整理的量子力学问题matlab求解,一个关于量子力学中的matlab的问题的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。