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

#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.