Columns 1 through 7

-0.9346 -0.8120 -0.7002 -0.5987 -0.5071 -0.4249 -0.3519

Columns 8 through 14

-0.2875 -0.2312 -0.1825 -0.1411 -0.1064 -0.0778 -0.0548

Columns 15 through 16

-0.0368 -0.0233

Для реалізації даного алгоритму створимо m-файл PТr.m, що містить опис функції, що повертає координати точок фазової тpaeкторії для кожного свого значення енергії.

% Лістнг 'файлу PTr.m

function z = PTr (EL, Xleft, Xright, gamma, Np)

% Функція, що повертає координати точок фазової траєкторії

% Для кожного свого значення енергії

N=length(EL);

k=I:Np;

for i=I:N

dx=(Xright(i)-Хlеft(i))/(Np-l);

x(k)=Xleft(i)+dx*(k-l);

for j=I:Np

Z(j,I)=x(j);

Z(j,2)=real(p(x(j),EL(i),gamma));

End;

for j=Np+l:2*Np

Z(j,I)=x(2*Np-j+ 1);

Z(j,2)=-Z(2*Nр-j+ 1,2);

End;

Z(2*Np+l,I)=x(I);

Z(2*Np+l,2)=Z(1,2);

if i==l

Z=z;

Else

z=саt(З,z,Z);

End;

End;

function z=p(x,E,gamma)

% Функція, що повертає значення імпульсу

z=gamma*(E­ V(x)). ЛО.5;

functiOll z=V(x)

% Функція, що повертає значення потенціалу

% Лейнарда-Джонсона

z=4*(1./x. л 12-1./x. Л6);

Далі необхідно виконати наступну послідовність команд:

» Хmiп=2Л(1/6); % координата мінімуму потенціалу

» gamma=80; % коефіцієнт, що входить в без вимірне

% рівняння Шредінгера (16.22)

»Emin=-l; % мінімальне значення енергії

» Еmах=-lОЛ-4; % максимальне значення енергії

» z = Action(Emin,Emax,Xmin,gamma); % обчислення залежності S(E)

» [Xleft,Xright,EL] = Spectr(z,Xmin) % обчислення енергетичного

% Спектру коливань

% Двоатомних молекули

» Np=2000; % число точок однієї фазової траєкторії

» Ph=PTr(EL,Xleft,Xright,gamma,Np); % обчислення координат точок

% Фазових траєкторій для кожного

% Значення енергетичного спектру

% Візуалізація фазових траєкторій

» plot(Ph(:,1,1),Ph(:,2,1),'k')

» hold оп

» for k=2:length(EL)

plot(Ph(:,1,k),Ph(:,2,k), 'k');

End;

» hold off

СПИСОК ЛІТЕРАТУРИ

Основна:

1 Поршнев С. В. Комп’ютерне моделювання фізичних процесів в пакеті MATLAB. – М.: Гаряча лінія – Телеком, 2003. – 592 с.

2 Храмов В. Ю. Розрахунок елементів лазерних систем для інформаційних та технологічних комплексів. – СПб: СПбДУ ИТМО, 2008, 79 с.

Додаткова

3.Рейф Ф.Статична фізика М Наука 1972

4. Metropolis N Rosenbluth A W Teller A H Equation of state

calculations for fast computing machines, j, Chem, Phys, 6, 1087, 1953,

5. Bhanot G" Creutz М., Neuberger Н., Microcanonical simulation of Ising system, Nucl, Phys, 235,417,1984,

6. Creutz М" Microcanonical Monte Simulation, Phys, Rev, Letts, 50,1411,1983,

Навчальне видання

Методичні вказівки до виконання лабораторних робіт з дисципліни «Квантова електроніка» для студентів всіх форм навчання напряму 6.050802“Електронні пристрої та системи” та 6.050801 «Мікро та наноелектроніка»,

укладач:

Литвиненко В.А.

Підписано до друку ___________ 2012р.

Формат________Обсяг ______др. арк.

Тираж ______екз. Заказ

51918, м. Дніпродзержинськ,

вул. Дніпробудівська, 2.

Дата добавления: 2015-07-11