fematrixk.m

来自「采用有限元法计算半导体的带隙结构和波函数。」· M 代码 · 共 24 行

%***************************************************************************************
% FEMQWEBS @Finite element method for quantum well electronic band structure calculation
% | ------------ |
% | Description: |
% | ------------ |
% The tool box provides the procedure to calculate all band edge energies 
% and corresponding wavefunctions in sinlge quantum square well using
% Finite Element Method.
% | ------------ |
% | Author:      |
% | ------------ |
% Full name: Le Quang Khai
% Email: [e-mail] or [e-mail]
%***************************************************************************************
function [k,p,q]=feMatrixK(xL,xR)
% Purpose 
% element matrix for [P] + pi^2*[Q] using liner element
eleng=xR-xL;
xav=(xR+xL)/2;
a=1/(mass(xav)*eleng);
b=pi^2*V(xav)*eleng/6;
k=[a+2*b -a+b; -a+b a+2*b];
p=[a -a;-a a];
q=[2*b b; b 2*b];