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

D=numeric(D);

R=numeric(R);

S=D+R;

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

pandas使用isna函数和any函数判断dataframe中的每一个数据列中是否包含缺失值
pandas使用isna函数和any函数判断dataframe中的每一个数据列中是否包含缺失值(check if column contains any missing values in dataf ...
在Java中如何从一个多层嵌套循环中直接跳出？
一.问题描述最近发现一个有趣的问题,在Java中如何从一个多层嵌套循环中退出呢? 例如下面,有两层循环,break只能退出第一个for循环,不能直接跳出第二个for循环外面. public stat ...
C++ 在一个cpp文件中使用另一个cpp文件中定义的函数
C++ 在一个cpp文件中使用另一个cpp文件中定义的函数建立一个console项目头文件 cpp文件 main.cpp 建立一个console项目以dev为例,在[文件][新建][项目],新 ...
python 在一个py文件中调用另一个py文件中的变量以及一个有趣的问题
python 在一个py文件中调用另一个py文件中的变量以及一个有趣的问题 1.在一个py文件中调用另一个py文件中的变量新建一个py文件,记为t1.py # t1.py m = 2 n = 1 a ...
php中goto另一个页面,php中goto操作符的详细介绍
php goto操作符可以用来跳转到程序的另一个位置.该位置可以用目标名称加冒号来标记,而跳转指令是在goto之后接上目标的位置标记对已goto操作符使用的一些限制目标位置只能位于同一个文件和作 ...
layui中从上一个js模块中取参数_layui的引用js踩坑
前言: 今天因为项目需要,需要使用layui,因为本身不多的前端经验,以为layui的用法和其他的前端框架,例如jquery,bootstrap等等,只需要直接引入layui.js,和layui.cs ...
matlab与origin关联,Matlab与Origin软件在“电路及磁路”教学中的应用
根据"电路及磁路"教学过程中出现的正弦交流电波形图难以绘制的问题,提出用Matlab与Origin软件相结合的方法来描述它.通过教程中实例的分析,逐步阐述了Matlab与Origi ...
标签的for循环和if_SO面试题08：如何从一个多层嵌套循环中直接跳出？
# 问题 Java中如何从一个多层嵌套循环中退出? 例如,下面有两个循环,break只能退出一个for循环,不能直接跳过第二个for循环 for (Type type : types) { for ...
如何在Excel中将工作表复制或移动到另一个工作簿中
There may be times when you need to copy or move a worksheet into another workbook in Excel or make ...

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

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

最新文章

热门文章