Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
ПЗ / 2 / 2.docx
Скачиваний:
6
Добавлен:
09.06.2023
Размер:
123.99 Кб
Скачать

2 Решение системы (13)

Дана система

Для начала её решения найдём производные , :

т.к. , и , и , и выражение является сложной функцией от x, то её производная равна:

Производная находится аналогично , поэтому:

Производная находится аналогично предыдущим, поэтому:

Теперь с учётом найденных производных и (9) преобразуем систему (13), получим:

На время решения системы (13.4) обозначим выражение как , получим:

Разделим все уравнения в (13.5) на C21, получим:

Сделаем ещё одни замены:

Произведём замену в (13.5) на сделанные замены (13.6—13.10), получим:

Совершая последовательные преобразования над системой (13.11), мы должны выразить величину t… и получить:

Т.к. t определяется выражением (13.8), а коэффициент прохождения T (11), то:

где .

Приложение б

Код программы необходимый для вычислений:

clear;

%данные

me=9.10*(10^-31);

h=1.05*(10^-34);

l=0.5*(10^-9);

w1=0;

w2=-1.282*10^(-18);

w3=-6.409*10^(-19);

w0=0*10^(-19):0.01*10^(-19):300*10^(-19);

k1=zeros();

k2=zeros();

k3=zeros();

T1=zeros();

T2=zeros();

%задание значений k1

for i=1:size(w0,2)

k1(i)=sqrt(2*me*(w0(i)-w1))/h;

end

%задание значений k2

for i=1:size(w0,2)

k2(i)=sqrt(2*me*(w0(i)-w2))/h;

end

%задание значений k3

for i=1:size(w0,2)

k3(i)=sqrt(2*me*(w0(i)-w3))/h;

end

%задание значений T для несим. ямы

for i=1:size(w0,2)

T1(i)=abs((4*k1(i)*k2(i))/(((k1(i)-k2(i))*(k2(i)-k3(i))*exp(-1i*l*(k3(i)-k2(i))))-(k1(i)+k2(i))*(k2(i)+k3(i))*exp(-1i*l*(k3(i)+k2(i)))))^2;

end

%переопределение значений k3 для условия k3=k1

for i=1:size(w0,2)

k3(i)=k1(i);

end

%задание значений T для сим. ямы

for i=1:size(w0,2)

T2(i)=abs((4*k1(i)*k2(i))/(((k1(i)-k2(i))*(k2(i)-k3(i))*exp(-1i*l*(k3(i)-k2(i))))-(k1(i)+k2(i))*(k2(i)+k3(i))*exp(-1i*l*(k3(i)+k2(i)))))^2;

end

%график T для (не)сим. ямы

plot(w0,T1,w0,T2,'Linewidth',2);

legend('Несимметричная яма','Симметричная яма','FontSize',14,'FontName','Times New Roman')

xlabel('E [J]','FontSize',14,'FontName','Times New Roman');

ylabel('T1(E), T2(E)','FontSize',14,'FontName','Times New Roman');

Таганрог 2022

Соседние файлы в папке 2