Спектры элементарных возбуждений в двупериодических одномерных системах

#include <iostream>

#include <math.h>

using namespace std;

#define N 30

#define a 1.0

#define b 1.1

#define d 0.5

// Перегрузка +

Complex Complex::operator+(Complex &fp1)

{

fp1.real = real + fp1.real;

fp1.image = image + fp1.image;

return fp1;

}

// Перегрузка -

Complex Complex::operator-(Complex &fp1)

{

fp1.real = real - fp1.real;

fp1.image = image - fp1.image;

return fp1;

}

// Перегрузка *

Complex Complex::operator*(Complex &fp1)

{

double i, j;

i = real * fp1.real - image * fp1.image;

j = real * fp1.image + fp1.real * image;

fp1.real = i;

fp1.image = j;

return fp1;

}

// Перегрузка /

Complex Complex::operator/(Complex &fp1)

{

double k, i, j;

k = fp1.real * fp1.real + fp1.image * fp1.image;

i = (real * fp1.real + image * fp1.image) / k;

j = (fp1.real * image - real * fp1.image) / k;

fp1.real = i;

fp1.image = j;

return fp1;

}

int main()

{

Complex Gkp;

double m;

int i,j;

for(i=0;i<N;i++)

for(j=0;j<N;j++)

{

Gkp.real=0;

Gkp.image=0;

Gkp.real=Gkp.real+1/(double)N*exp(-1/a*sqrt(pow(i*a-j*b,2)+d*d))*cos(6.28*i-6.28*j);

Gkp.image=Gkp.image-1/(double)N*exp(-1/a*sqrt(pow(i*a-j*b,2)+d*d))*sin(6.28*i-6.28*j);

}

Gkp.real=pow(Gkp.absolute(),2);

cout<<"Gkp"<<" "<<Gkp.real<<"\n";

getchar();

}

Перейти на страницу: 1 2 3 4 

Другие публикации

Усовершенствование материнской платы
Тема дипломной работы – "Усовершенствование материнской платы", являющаяся предметом исследования. Цель работы – выяснить неблагоприятные факторы работы м ...

Устройство импульсного управления исполнительным двигателем постоянного тока
Машины постоянного тока до сих пор активно применяются в качестве двигателей (ДПТ) и генераторов (ГПТ). ДПТ имеют хорошие регулировочные свойства, значительную пере ...

Меню

Copyright @2020, TECHsectors.ru.