cronig.m

来自「Resolution of the Cronig-Penney s predat」· M 代码 · 共 43 行

clear all;
clc;

global P;
global Rg;
global a;
global abg;

P = 13.1;
Rg = 5.76;  % meV
a = 10;     % nm
abg = 9.97; % nm

N = 2000;

x = [-pi : 2*pi/(N-1) : pi];
y = [0.001 : (5*pi - 0.001)/(N-1) : 5*pi];

E = zeros(4, N);

for j = 1 : length(x)
    
    Fw = Ffun(x(j), y);
    [nn, zi] = findzeros(Fw);

    zx = x(zi);
    zy = y(zi);
    
    for k = 1 : length(zy)
        E(k, j) = zy(k)^2 * (abg/a)^2 * Rg;
    end
    
    %plot(x, Fw, 'b', zx, zy, 'ro');
    %hold on;
end
    
figure;

plot(x, E([1 : 4], :), '.');
xlim([-pi pi]);