Приложение В. Листинг функций расчета интегральной фазовой расстройки

(обязательное)

//Plot phase mismatch

void MainWindow::calc_phase()

{

setGraphData();

const int n = 2000;

double *n0r = new double [n];

double *n1r = new double [n];

double *n0l = new double [n];

double *n1l = new double [n];

double *N0r = new double [n];

double *N1r = new double [n];

double *N0l = new double [n];

double *N1l = new double [n];

double *k0r = new double [n];

double *k1r = new double [n];

double *k0l = new double [n];

double *k1l = new double [n];

double *Kr = new double [n];

double *Kl = new double [n];

for (int i=0; i<n; ++i)

{

Q[i]=1.571-Q[i];

n0r[i]=(ne*no)/sqrt((no*no)*(sin(Q[i]+teta)*sin(Q[i]+teta))+(ne*ne)*(cos(Q[i]+teta)*cos(Q[i]+teta)));

n1r[i]=(ne*no)/sqrt((no*no)*(sin(Q[i]-teta)*sin(Q[i]-teta))+(ne*ne)*(cos(Q[i]-teta)*cos(Q[i]-teta)));

n0l[i]=(ne*no)/sqrt((no*no)*(sin(3.14-Q[i]+teta)*sin(3.14-Q[i]+teta))+(ne*ne)*(cos(3.14-Q[i]+teta)*cos(3.14-Q[i]+teta)));

n1l[i]=(ne*no)/sqrt((no*no)*(sin(3.14-Q[i]-teta)*sin(3.14-Q[i]-teta))+(ne*ne)*(cos(3.14-Q[i]-teta)*cos(3.14-Q[i]-teta)));

N0r[i]=asin((n0r[0]/n0r[i])*sin(teta));

N1r[i]=asin((n1r[0]/n1r[i])*sin(teta));

N0l[i]=asin((n0l[0]/n0l[i])*sin(teta));

N1l[i]=asin((n1l[0]/n1l[i])*sin(teta));

k0r[i]=2*(3.14/lambda)*n0r[i]*cos(N0r[i]);

k1r[i]=2*(3.14/lambda)*n1r[i]*cos(N1r[i]);

k0l[i]=2*(3.14/lambda)*n0l[i]*cos(N0l[i]);

k1l[i]=2*(3.14/lambda)*n1l[i]*cos(N1l[i]);

Kr[i]=k0r[i]-k1r[i];

Kl[i]=k0l[i]-k1l[i];

}

for (int i=0; i<n; ++i)

{

Qr[i]=0;

Ql[i]=0;

for (int j=0; j<=i; ++j)

{

Qr[i]=Qr[i]+Kr[j];

Ql[i]=Ql[i]+Kl[j];

}

}

for (int i=0; i<n;++i)

{

Qr[i]=Qr[i]*(d/n);

Ql[i]=Ql[i]*(d/n);

}

delete n0r;

delete n1r;

delete n0l;

delete n1l;

delete N0r;

delete N1r;

delete N0l;

delete N1l;

delete k0r;

delete k1r;

delete k0l;

delete k1l;

delete Kr;

delete Kl;

}

void MainWindow::plot_phase()

{

graph = 2;

const int n = 2000;

double *y1 = new double [n];

ui->qwtPlot->setAxisScale(QwtPlot::xBottom, (-d/2)*1000000, (d/2)*1000000);

ui->qwtPlot->setAxisScale(QwtPlot::yLeft, -4, 4);

for (int i=0; i<n;++i)

y1[i]=(-d/2+i*(d/n))*1000000;

//ui->doubleSpinBox->setValue(fabs(n1r[1000]));

//curve1->attach(ui->qwtPlot);

ui->x_label->setText("d, um");

ui->y_label->setText("Phase mismatch, abs. units");

zoomer = new ScrollZoomer(ui->qwtPlot->canvas());

curve->setSamples(y1,Ql,n);

curve1->setSamples(y1,Qr,n);

curve->attach(ui->qwtPlot);

curve1->attach(ui->qwtPlot);

curve2->detach();

ui->qwtPlot->replot();

deletey1;

}

Приложение Г. Листинг функций расчета дифракционных характеристик

(обязательное)

//Plot diffraction characteristics

void MainWindow::plot_dif()

{

QFile file("data.csv");

graph = 3;

setGraphData();

calc_phase();

const int n = 2000;

double *delta = new double [n];

ui->qwtPlot->setAxisScale(QwtPlot::xBottom, -20, 20);

ui->qwtPlot->setAxisScale(QwtPlot::yLeft, 0, 1);

std::complex<double> *cIdr = new std::complex<double> [n];

std::complex<double> *cIdl = new std::complex<double> [n];

std::complex<double> *cIds = new std::complex<double> [n];

std::complex<double> temp;

double *Idr = new double [n];

double *Idl = new double [n];

double *Ids = new double [n];

double *y1 = new double [n];

for (int i=0; i<n;++i)

y1[i]=(0+i*(d/(n-1)));

for (int i=0; i<n;++i)

{

delta[i]=(-20+i*(2*20)/(n-1));

}

for (int i=0; i<n; ++i)

{

cIdr[i]=std::complex<double> (0,0);

cIdl[i]=std::complex<double> (0,0);

cIds[i]=std::complex<double> (0,0);

for (int j=0; j<n; ++j)

{

temp=std::complex<double> (0,delta[i]*(y1[j]/d)+Qr[j]);

cIdr[i]=cIdr[i]+(1.0)*exp(temp);

temp=std::complex<double> (0,delta[i]*(y1[j]/d)+Ql[j]);

cIdl[i]=cIdl[i]+(1.0)*exp(temp);

cIds[i]=(cIdr[i]+cIdl[i])/2.0;

}

}

for (int i=0; i<n; ++i)

{

Idr[i]=pow(abs(cIdr[i])/n,2);

Idl[i]=pow(abs(cIdl[i])/n,2);

Ids[i]=pow(abs(cIds[i])/n,2);

}

//Recording data to file

//***************

QTextStream out(&file);

if(file.open(QIODevice::WriteOnly|QIODevice::Truncate))

{

out<<"Diff. eff. positive"<<';';

out<<"Diff. eff. negative"<<';';

out<<"Diff. eff. summary"<<';';

out<<"Relative mismatch"<<';'<<"\r\n";

for (int i=0; i<n; ++i)

{

out.setRealNumberNotation(QTextStream::ScientificNotation);

out<<Idr[i]<<';';

out<<Idl[i]<<';';

out<<Ids[i]<<';';

out<<delta[i]<<';'<<"\r\n";

}

file.close();

}

//***************

ui->x_label->setText("delta, abs. units");

ui->y_label->setText("Diffraction efficiency, abs. units");

zoomer = new ScrollZoomer(ui->qwtPlot->canvas());

curve->setSamples(delta,Idl,n);

curve1->setSamples(delta,Idr,n);

curve2->setSamples(delta,Ids,n);

curve->attach(ui->qwtPlot);

curve1->attach(ui->qwtPlot);

curve2->attach(ui->qwtPlot);

ui->qwtPlot->replot();

delete delta;

delete Idr;

delete Idl;

delete y1;

delete cIdr;

delete cIdl;

}

Приложение Д. Листинг программы в MathCAD

(обязательное)