ПЗ / 1 / Вычисления

%расчёт для задания 1)
clear;
L1=0.5*(10^-9);
L2=5*(10^-9);
dL=0.04*(10^-9);
L=L1:dL:L2;
me=9.10*(10^-31);
h=1.05*(10^-34);
n=2;
m=1;
for i=1:113
E(i)=(1/(2*me))*(((pi*h)/(L(i)))^2)*(n^2-m^2);
end
plot(L,E,'black','Linewidth',2);
xlabel('L [м]','FontSize',14,'FontName','Times New Roman');
ylabel('∆_{E_m}E_n(L) [Дж]','FontSize',14,'FontName','Times New Roman');
%расчёт для задания 2)
clear;
%задание параметров и дискретизация величин
x1=-0.25*(10^-9);
x2=0.75*(10^-9);
dx=0.004*(10^-9);
x=x1:dx:x2;
n=1:1:6;
L=0.5*(10^-9);
Psi=zeros(length(x),length(n));
PSI=zeros(length(x),length(n));
%вычисление значений пси-функций и вероятностей
for j=1:6
for i=1:251
if x(i)<0
Psi(i,j)=0;
PSI(i,j)=0;
elseif x(i)>L
Psi(i,j)=0;
PSI(i,j)=0;
else
Psi(i,j)=sqrt(2/L)*sin((n(j)*pi*x(i))/L);
PSI(i,j)=abs(Psi(i,j))^2;
end
end
end
%графики Ψ-функций
f1=figure;
for j=1:6
subplot(2,3,j)
plot(x,Psi(:,j),'black','Linewidth',2);
title('n=','FontSize',14,'FontName','Times New Roman');
xlabel('x [м]','FontSize',14,'FontName','Times New Roman');
ylabel('Ψ','FontSize',14,'FontName','Times New Roman');
end
%графики вероятностей
f2=figure;
for j=1:6
subplot(2,3,j)
plot(x,PSI(:,j),'black','Linewidth',2);
title('n=','FontSize',14,'FontName','Times New Roman');
xlabel('x [м]','FontSize',14,'FontName','Times New Roman');
ylabel('Ψ^2','FontSize',14,'FontName','Times New Roman');
end