Нейрокомпьютерное моделирование

#include <vcl.h>
#pragma hdrstop
#include <math.h>
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
int a;
const float Pi=3.14;
 
 
 
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
 
float fx1( float am, float dx1,float x)
{
float s;
s=am*sin(x*Pi/(2*dx1));
return s;
 
}
float fx2( float y0, float k,float x)
{
float s;
s=y0+k*x;
return s;
 
}
float f1(float ta, float am,float tn, float ty,float t )
{
float s;
if (t<=ta)
s=0;
if ((t>=ta)&&(t<tn+ta))
s=am*sin((t-ta)*Pi/(2*tn));
if ((t>=ta+tn)&&(t<tn+ta+ty))
s=am-am/(ty)*(t-ta-tn);
if (t>=tn+ta+ty)
s=0;
return s;
 
}
 
 
 
 
void __fastcall TForm1::BitBtn2Click(TObject *Sender)
{
a=random(15);
LabeledEdit1->Text=IntToStr(10*a);
}
//---------------------------------------------------------------------------
void __fastcall TForm1::BitBtn1Click(TObject *Sender)
{
int h1,h2,h3,h4,x0,dx,ta1,ta2,ta3;
//uv1,uv2,uv3,uv4 -амплитуда входного импульса
//ty1, ty2, ty3 - время угасания двухскатного потенциала торможения
float a1,a2,a3,a4,uv1,uv2,uv3,uv4,umax,ty1,ty2,ty3,t,dt;
//tn1, tn2, tn3 - время нарастания двухскатного потенциала торможения
int tn1,tn2,tn3,am1,am2,am3;
am1=20;am2=30,am3=40;
x0=50;
h1=70;
h2=h1+70;
h3=h2+70;
h4=h3+70;
uv1=100;
uv2=110;uv3=120;
uv4=StrToFloat(LabeledEdit1->Text);
umax=200.0;
a1=h1/(umax)*uv1;
a2=(h2-h1)/(umax)*uv2;
a3=(h3-h2)/(umax)*uv3;
a4=(h4-h3)/(umax)*uv4;
Image1->Canvas->Pen->Color=clWhite;
Image1->Canvas->Brush->Color=clWhite;
Image1->Canvas->Rectangle(ClientRect);
Image1->Canvas->Pen->Color=clBlack;
Image1->Canvas->MoveTo(x0,h1);
Image1->Canvas->LineTo(Image1->Width-x0,h1);
Image1->Canvas->MoveTo(x0,h2);
Image1->Canvas->LineTo(Image1->Width-x0,h2);
Image1->Canvas->MoveTo(x0,h3);
Image1->Canvas->LineTo(Image1->Width-x0,h3);
Image1->Canvas->MoveTo(x0,h4);
Image1->Canvas->LineTo(Image1->Width-x0,h4);
dx=10;
if (uv4>=uv1)
{
//1
Image1->Canvas->TextOutA(1,h1-a1-20,"Uпн=100 мВ");
Image1->Canvas->TextOutA(x0+50,h1-20,"Сигнал торможения с выхода первого слоя");
Image1->Canvas->MoveTo(x0,h1);
Image1->Canvas->LineTo(x0,h1-a1);
Image1->Canvas->MoveTo(x0,h1-a1);
Image1->Canvas->LineTo(x0+dx,h1-a1);
Image1->Canvas->MoveTo(x0+dx,h1-a1);
Image1->Canvas->LineTo(x0+dx,h1);
}
if (uv4>=uv2)
{
//2
Image1->Canvas->TextOutA(1,h2-a2-20,"Uпс=110 мВ");
Image1->Canvas->TextOutA(x0+50,h2-20,"Сигнал торможения с выхода второго слоя");
Image1->Canvas->MoveTo(x0,h2);
Image1->Canvas->LineTo(x0,h2-a2);
Image1->Canvas->MoveTo(x0,h2-a2);
Image1->Canvas->LineTo(x0+dx,h2-a2);
Image1->Canvas->MoveTo(x0+dx,h2-a2);
Image1->Canvas->LineTo(x0+dx,h2);
}
if (uv4>=uv3)
{
//3
Image1->Canvas->TextOutA(1,h3-a3-20,"Uпв=120 мВ");
Image1->Canvas->TextOutA(x0+50,h3-20,"Сигнал торможения с выхода третьего слоя");
Image1->Canvas->MoveTo(x0,h3);
Image1->Canvas->LineTo(x0,h3-a3);
Image1->Canvas->MoveTo(x0,h3-a3);
Image1->Canvas->LineTo(x0+dx,h3-a3);
Image1->Canvas->MoveTo(x0+dx,h3-a3);
Image1->Canvas->LineTo(x0+dx,h3);
}
//4
Image1->Canvas->TextOutA(1,h4-a4-20,"Uвх="+FloatToStr(uv4)+"мВ");
Image1->Canvas->TextOutA(x0+50,h4-20,"Суммарный сигнал");
Image1->Canvas->MoveTo(x0,h4);
Image1->Canvas->LineTo(x0,h4-a4);
Image1->Canvas->MoveTo(x0,h4-a4);
Image1->Canvas->LineTo(x0+dx,h4-a4);
Image1->Canvas->MoveTo(x0+dx,h4-a4);
Image1->Canvas->LineTo(x0+dx,h4);
dt=0.01;
t=0;
ta1=StrToInt(LabeledEdit2->Text);
ta2=StrToInt(LabeledEdit3->Text);
ta3=StrToInt(LabeledEdit4->Text);
tn1=StrToInt(LabeledEdit5->Text);
tn2=StrToInt(LabeledEdit6->Text);
tn3=StrToInt(LabeledEdit7->Text);
ty1=StrToFloat(LabeledEdit8->Text);
ty2=StrToFloat(LabeledEdit9->Text);
ty3=StrToFloat(LabeledEdit10->Text);
 
Image1->Canvas->MoveTo(x0+dx,h1);
Image1->Canvas->LineTo(x0+dx+ta1,h1);
if (uv4>=uv1)
{
while (t<tn1)
{
 Image1->Canvas->Pixels[x0+dx+ta1+t][h1+fx1(am1,tn1,t)]=clBlue;
 
 t+=dt;
}
t=0;
while (t<ty1)
{
 Image1->Canvas->Pixels[x0+dx+ta1+tn1+t][h1+fx2(am1,-am1/(ty1),t)]=clBlue;
 
 t+=dt;
}
}
if (uv4>=uv2)
{
t=0;
Image1->Canvas->MoveTo(x0+dx,h2);
Image1->Canvas->LineTo(x0+dx+ta2,h2);
while (t<tn2)
{
 Image1->Canvas->Pixels[x0+dx+ta2+t][h2+fx1(am2,tn2,t)]=clBlue;
 
 t+=dt;
}
t=0;
while (t<ty2)
{
 Image1->Canvas->Pixels[x0+dx+ta2+tn2+t][h2+fx2(am2,-am2/(ty2),t)]=clBlue;
 
 t+=dt;
}
}
if (uv4>=uv3)
{
t=0;
Image1->Canvas->MoveTo(x0+dx,h3);
Image1->Canvas->LineTo(x0+dx+ta3,h3);
while (t<tn3)
{
 Image1->Canvas->Pixels[x0+dx+ta3+t][h3+fx1(am3,tn3,t)]=clBlue;
 
 t+=dt;
}
 
t=0;
while (t<ty3)
{
 Image1->Canvas->Pixels[x0+dx+ta3+tn3+t][h3+fx2(am3,-am3/(ty3),t)]=clBlue;
 
 t+=dt;
}
}
t=0;
if ((uv4>=uv1)&&(uv4<uv2))
{
while (t<Image1->Width-2*x0-dx)
{
 Image1->Canvas->Pixels[x0+dx+t][h4+f1(ta1,am1,tn1,ty1,t)]=clBlue;
 
 t+=dt;
}
}
if ((uv4>=uv2)&&(uv4<uv3))
{
while (t<Image1->Width-2*x0-dx)
{
 Image1->Canvas->Pixels[x0+dx+t][h4+f1(ta1,am1,tn1,ty1,t)+f1(ta2,am2,tn2,ty2,t)]=clBlue;
 
 t+=dt;
}
}
if (uv4>=uv3)
{
while (t<Image1->Width-2*x0-dx)
{
 Image1->Canvas->Pixels[x0+dx+t][h4+f1(ta1,am1,tn1,ty1,t)+f1(ta2,am2,tn2,ty2,t)+f1(ta3,am3,tn3,ty3,t)]=clBlue;
 
 t+=dt;
}
}
}