Квантовый обрыв с питоном - PullRequest
Новая
       29

Квантовый обрыв с питоном

0 голосов
/ 29 ноября 2018

Я пытаюсь построить волновую функцию, которая разбросана по шаговой функции.Я упомянул проблему квантовой механики Гриффита 2-е издание 2.35.здесь энергия E = V / 3 и я устанавливаю V = 3.форма потенциала выглядит как обрыв, левая сторона обрыва выше, чем правая.

Затем я пытаюсь найти коэффициент отражения для этого, и я просто сравнил высоту волновой функции;слева относительно ступенчатой ​​функции 2, а справа 1,7, и я получил разумное значение вероятности отражения;0,15, но точный ответ - 0,11111

Мой вопрос заключается в том, является ли этот код правильным или я просто случайно получил правильный ответ ??? 

    import matplotlib.pyplot as plt
    %matplotlib inline
    from scipy.integrate import odeint
    import numpy as np

    def ps(y,t):
        if t<10:
            u=0
        else:
            u=-3
        psi=y[0]
        dpsi=y[1]
        dp=dpsi
        dp2=(-1+u)*psi
        return [dp,dp2]
    yo=[2,0]
    t=np.linspace(0,100,400)
    psis=odeint(ps,yo,t)    
    #plt.ylim(-3,3)
    plt.xlim(0,20)
    plt.plot(t,psis[:,0],'r')

    plt.plot([0,10,10,40],[0,0,-3,-3],'g')
...