求助：模糊系统、模糊神经网络及应用程序设计例子什么不对
请高手帮我看看这个程序为什么不对?

#define q1         0.00012
#define q2         0.008
#define M          81
#define SAMPLES    80
#define OK         0

/*定义第二个图象的位置*/
#define graph2_x1  8
#define graph2_x2  632
#define graph2_y1  80
#define graph2_y2  450

/*显示活动条的位置*/
#define act_bar_x1  40
#define act_bar_y1  110
#define act_bar_x2  600
#define act_bar_y2  140

/*显示训练次数和均方误差的位置*/
#define train_x  156
#define train_y  400
#define train_x  448
#define train_y  400

/*显示询问框的位置*/
#define ask_x1  100
#define ask_y1  200
#define ask_x2  540
#define ask_x2  260

/* 定义测试样本的次数 */
#define T      10

#define LEFT   0x02
#define CTRL   0x04
#define ALT    0x08

#include "bios.h"
#include "ctype.h"
#include "math.h"
#include "stdio.h"
#include "conio.h"
#include "stdlib.h"
#include "graphics.h"
#include "iostream.h"
#include "alloc.h"
#include "string.h"

static double        p[M+1][5];
static long double   w[M+1];
static double        a[M+1][5];
static double        b[M+1][5];

static long double   y[M+1];
static long double   xi[5],yd;

static double        t_Y[T];
static double        t_X1[T];
static double        t_X2[T];
static double        t_X3[T];
static double        t_X4[T];

long double (*f[4])(long double,long double,long double);

void far *ptr[4];

int maxx=639,maxy=479;

int train;
double error=1,mean;

/* 样本数据x1x2x3x4以及期望值输出yd */
static long double YD[]=
{
/* 样本数据省略*/
};

static long double X1[]=
{
/* 样本数据省略*/
};

static long double X2[]=
{
/* 样本数据省略*/
};

static long double X3[]=
{
/* 样本数据省略*/
};

static long double X4[]=
{
};

/* 初始化隶属度函数的参数a,b*/
static long double al_before[4]{0,15,0,7.5};
static long double bl_before[4]{3,3,3,1.5};

/* 隶属度函数*/
inline long double uPL(long double x, long double a, long double b)
{
if(x>15) return(1);
else return(exp(-(x-a)*(x-a)/b));
}

inline long double uNL(long double x, long double a, long double b)
{
if(x>=0) return(1);
else return(1);
}

inline long double uZO(long double x, long double a, long double b)
{
 return(exp(-(x-a)*(x-a)/b));
}

/* 计算y[1]到y[81]的值*/

 int yL()
{int i;
for(i=1;i<=M;i++)
y[i]=p[i][0]+p[i][1]*xi[1]+p[i][2]*xi[2]+p[i][3]*xi[3]+p[i][4]*xi[4];
return OK;}

/* 计算w[1]到w[81]的值*/

 int wL()
{
int i,j,k,l,m;
for(i=1;i<=3;i++)
 for(j=1;j<=3;j++)
  for(k=1;k<=3;k++)
   for(l=1;l<=3;l++)
    {m=27*(i-1)+9*(j-1)+3*(k-1)+1;
     w[m]=((*f[i])(xi[1],a[m][1],b[m][1]))*((*f[i])(xi[2],a[m][2],b[m][2]))*
     ((*f[i])(xi[3],a[m][3],b[m][3]))*((*f[i])(xi[4],a[m][4],b[m][4]));
}
   return OK;}

/* 计算BP神经网络的输出值*/
long double Y_effect()
{int i;
 long double WY=0,W=0;
 yL();
 wL();

 for(i=1;i<=M;i++)
{WY=WY+w[i]*y[i];
W=W+w[i];}

return(WY/W);}

int border(int x1,int y1,int x2,int y2,int color)
{sectcolor(color);
rectangle(x1,y1,x2,y2);
return OK;}

/* 在X轴上画刻度*/
int ke_du_x(int x1,int x2,int y,int range,int color,int flag)
{
float i,j;
int n=3;
float rule;
char string[10];
rule=(x2-x1)/range;
setcolor(color);
for(i=x1;i<=x2;i=i+rule)

{if(flag) ltoa(((i-x1)/rule+1),string,10);
else ltoa(((i-x1)/rule+1),string,10);

if(((on++)%3==0)&&(flag))
 outtextxy(i-3,y+3,string);
else if(!flag)
 outtextxy(i-3,y+3,string);

 line(i,y,i,y-6);

if((i>x1)&&(!flag))
 for(j=i-rule;j<=i;j=j+rule/5)
 line(i,y,i,y-3);
return OK;}

/* 在y轴上画刻度*/
int ke_du_y(int x,int y1,int y2,int range,int color)
{
float i,j;
float rule;
char string[10];
rule=(y2-y1)/range;
settextjustify(0,2);
setcolor(color);

for(i=(y2-rule);i>=y1;i=i-rule)

{ltoa((y2-i)/rule),string,10);
 outtextxy(x-7,i-3,string);
 line(x,i,x+6,i);
 for(j=(i+rule);j>=i;j=j-rule/5)
 line(x,j,x+3,j);
}
return OK;}

/* 图形模型初始化 */
initgra()
{
int gdriver=DETECT,gmode,errorcode;
if(registerbgidriver(EGAVGA_driver)<0) goto err;
if(registerbgifont(triplex_font)<0) goto err;
if(registerbgifont(small_font)<0) goto err;

initgraph(&gdriver,&gmode,"c:\\tc");

err:
errorcode=graphresult();

if(errorcode != grOK)
{
 printf("graphics error:%s\n",grapherrormsg(errorcode));
 printf("press any key to halt:");
  getch();
  exit(1);
}
return OK;}

/* 画第一幅图象*/
int draw_graph1()
{
long double i,PL;
cleardriver();
setbkcolor(WHITE);
border(0,0,639,479,LIGHTBLUE);
settextstyle(0,0,0);
for(i=0;i<=15;i=i+0.005)
{putpixel(20+40*i,(230-200*uPL(i,al_before[1],bl_before[1])),DARKGRAY);
putpixel(20+40*i,(230-200*uNL(i,al_before[2],bl_before[2])),LIGHTBLUE);
putpixel(20+40*i,(230-200*uZO(i,al_before[3],bl_before[3])),LIGHTRED);

}

outtextxy(25+40*1,80,"uPL");
outtextxy(20+40*7.2,80,"uZO");
outtextxy(15+40*13.5,80,"uNL");

border(20,25,620,230,DARKGRAY);
ke_du_x(20,620,230,15,LIGHTBLUE,0);
ke_du_y(20,25,230,1,LIGHTBLUE);

setfillstyle(l,LIGHTBLUE);
bar(0,0,639,20);
settextjustify(CENTER_TEXT,CENTER_TEXT);
setcolor(WHITE);
setfillstyle(DEFULT_FONT,0,1);
outtextxy(320,10,"GRAPH1");

for(i=0;i<=15;i=i+0.005)
{putpixel(20+40*i,(455-200*uPL(i,a[1][4],b[1][4]+1.2))),DARKGRAY);
putpixel(20+40*i,(455-200*uNL(i,a[2][4],b[2][4]+1.5))),LIGHTBLUE);
putpixel(20+40*i,(455-200*uZO(i,a[3][4],b[3][4]+1.2))),LIGHTRED);
}

setcolor(LIGHTBLUE);
settextjustify(0,2);
outtextxy(30+40*1,320,"uPL");
outtextxy(20+40*7.2,320,"uZO");
outtextxy(20+40*13.2,320,"uNL");
}
{border(20,250,620,455,DARKGRAY);
ke_du_x(20,620,455,15,LIGHTBLUE,0);
ke_du_y(20,250,455,1,LIGHTBLUE);}

return OK;}

 

/* 画第二幅图象*/
int draw_graph2()
{
long double i,PL;
float rule=(graph2_x2-graph2_x1)/(SAMPLES-1);
cleardriver();
setbkcolor(WHITE);
border(0,0,639,479,LIGHTBLUE);

border(graph2_x1,graph2_y1-20,graph2_x2,graph2_y2,DARKGRAY);
ke_du_x(graph2_x1,graph2_x2,graph2_y2,(SAMPLES-1),LIGHTBLUE,1);
ke_du_y(graph2_x1,graph2_y1,graph2_y2,1,LIGHTBLUE);

setfillstyle(1,LIGHTBLUE);
bar(1,1,638,29);
border(1,1,638,29,DARKGRAY);
settextjustify(CENTER_TEXT,CENTER_TEXT);
setcolor(WHITE);
settextstyle(DEFULT_FONT,0,1);
outtextxy(320,15,"GRAPH2");
setfillstyle(1,YELLOW);
bar(245,80,395,115);
border(245-1,80-1,395+1,115+1,DARKGRAY);
settextstyle(DEFULT_FONT,0,0);
settextjustify(LEFT_TEXT,CENTER_TEXT);
outtextxy(295,90,"PREDICTED");
setlinestyle(DOTTED_LINE,1,1);
line(265,90,285,90);
setfillstyle(SOLIDD_FILL,LIGHTRED);
fillellipse(275,90,3,3);
outtextxy(295,105,"RAEL");
line(265,105,285,105);
setcolor(BLUE);
circle(275,105,4);

moveto(graph2_x1,graph2_y2-(YD[0]*400));
for(i=0;i<=SAMPLES-1+T;i++)
{
if(i<=SAMPLES-1)
{
setcolor(LIGHTRED);
lineto(graph2_x1+rule*i,graph2_y2-(YD[i]*400));

setcolor(BLUE);
circle(graph2_x1+rule*i,graph2_y2-(YD[i]*400),4);
}

else
{
x[1]=t_X1[i-SAMPLES];
x[2]=t_X2[i-SAMPLES];
x[3]=t_X3[i-SAMPLES];
x[4]=t_X4[i-SAMPLES];
t_Y[i-SAMPLES]=(0.472*((1+sin(M_PI*(xi[1]-0.5)))/2)
 +0.607*((1+sin(M_PI*(xi[2]-0.5)))/2)
 +0.779*((1+sin(M_PI*(xi[3]-0.5)))/2)
 +((1+sin(M_PI*(xi[4]-0.5)))/2)/2.858;

setcolor(LIGHTRED);
lineto(graph2_x1+rule*i,graph2_y2-(t_Y[i-SAMPLES]*400));

setcolor(BLUE);
circle(graph2_x1+rule*i,graph2_y2-(t_Y[i-SAMPLES]*400),4);
}
}

moveto(graph2_x1,graph2_y2-(Y_effect())*400;
for(i=0;i<=SAMPLES-1+T;i++)
{
if(i<=SAMPLES-1)
{
xi[1]=X1[i];
xi[2]=X2[i];
xi[3]=X3[i];
xi[4]=X4[i];
setcolor(BLUE);
lineto(graph2_x1+rule*i,graph2_y2-(Y_effect())*400);

setfillstyle(SOLID_FILL,LIGHTRED);
fillellipse(graph2_x1+rule*i,graph2_y2-(Y_effect())*400),3,3);
}

else
{
xi[1]=t_X1[i-SAMPLES];
xi[2]=t_X2[i-SAMPLES];
xi[3]=t_X3[i-SAMPLES];
xi[4]=t_X4[i-SAMPLES];
setcolor(BLUE);
lineto(graph2_x1+rule*i,graph2_y2-(Y_effect())*400);

setfillstyle(SOLID_FILL,LIGHTRED);
fillellipse(graph2_x1+rule*i,graph2_y2-(Y_effect())*400),3,3);
}
}

return OK;}

/* 保存屏幕图象函数*/
void save_screen(void far *buf[4])
{
unsigned size;
int ystart=0,yend,yincr,block;

yincr=(maxy+1)/4;
yend=yincr;
size=imagesize(0,ystart,maxx,yend);

for(block=0;block<=3;block++)
{
if((buf[block]=farmalloc(size))==NULL)
{
closegraph();
printf("Error:not enough heap space in save_screen().\n"
exit(1);
}

getimage(0,ystart,maxx,yend,buf[block]);
ystart=yend+1;
yend+=yincr+1;
}
}

/* 恢复屏幕图象函数*/
void restore_screen(void far *buf[4])
{
int ystart=0,yend,yincr,block;

yincr=(maxy+1)/4;
yend=yincr;

for(block=0;block<=3;block++)
{
putimage(0,ystart,buf[block],COPY_PUT);
farfree(buf[block]);
ystart=yend+1;
yend+=yincr+1;
}
}

/*定期义用户界面*/
int interfce()
{
setbkcolor(LIGHTBLUE);
boeder(0,0,639,479,WHITE);
settextstyle(0,0,0);
settextjustify(1,1);
outtextxy(Train_x,Train_y,"Train time:");
outtextxy(error_x,error_y,"Mean error:");

outtextxy(Train_x-30,Train_y-80,"X1=:");
outtextxy(error_x,Train_y-80,"X2=:");
outtextxy(Train_x-30,Train_y-60,"X3=:");
outtextxy(error_x,error_y-60,"X4=:");

outtextxy(act_bar_x1-5,act_bar_y1-5,act_bar_x2+5,act_bar_y2+5);

setfillstyle(1,WHITE);
bar(act_bar_x1,act_bar_y1,act_bar_x2,act_bar_y2);
settextjustify(1,1);
settextjustify(TRIPLEX_FONT,0,4);
{
setcolor(WHITE);
outtextxy(319,24,"FUZZY LOGIC NETWORK MODELLING");
setcolor(DARKGRAY);
outtextxy(321,26,"FUZZY LOGIC NETWORK MODELLING");
setcolor(LIGHTGRAY);
outtextxy(320,25,"FUZZY LOGIC NETWORK MODELLING");
settextstyle(0,0,0);
outtextxy(320,470,"Press shife+lt+ctrl to see the graph......");
}

setcolor(WHITE);
line(0,50,639,50);
line(0,460,639,460);
line(0,Train_y-8,639,Train_y-8);
line(0,Train_y+8,639,Train_y+8);
line(320,Train_y-8,320,Train_y+8);
line(0,Train_y-90,639,Train_y-90);
line(0,Train_y-50,639,Train_y-50);
line(0,Train_y-70,639,Train_y-70);
line(320,Train_y-50,320,Train_y-90);

return OK;
}

/*交互窗口函数*/
int ask()
{
char key;
char *str1="Press ESC to exit system";
char *str2="Press ENTER to contiue system";
char *str3="to Go on system";
border(ask_x1,ask_y1,ask_x2,ask_y2,LIGHTBLUE);
setfillstyle(SOLIDD_FILL,LIGHTRED);
bar3d(ask_x1+1,ask_y1+1,ask_x2-1,ask_y2,2,1);
setcolor(BLUE);
settextstyle(2,0,6);
outtextxy(ask_x1+20,ask_y1+15,str1);
outtextxy(ask_x1+20,ask_y1+30,str2);
outtextxy(ask_x1+20,ask_y1+45,str3);

re_input:
key=getch();
if(key==13) return(1);
else if(key==27) return(0);
}

/* 保存p[][],a[][],b[][]数组的值及误差和训练次数*/
int save_pab()
{
FILE *fpp,*fpa,*fpb,*fpt,*fpe;
int i,j;

fpp=fopen("p.dat","wt");
fpa=fopen("a.dat","wt");
fpb=fopen("b.dat","wt");
fpt=fopen("Train.dat","wt");
fpe=fopen("error.dat","wt");
fprintf(fpp,"%s","    P0    P1    P2    P3    P4\n";
for(i=1;i<=M;i++)
    {for(j=1;j<=4;j++)
        {if(j==0)
            fprintf(fpp,"%d",i);
        if(p[i][j]<0)
        {
        if(j>=10)
            fprintf(fpp,"%-3.8f",p[i][j]);
        else
            fprintf(fpp,"%-3.8f",p[i][j]);
            }
         else
            {if(i>=10)
            fprintf(fpp,"%-3.8f",p[i][j]);
               else
            fprintf(fpp,"%-3.8f",p[i][j]);
              }       }
            fprintf(fpp,"\n");
              }
       for(i=1;i<=M;i++)
         {for(j=1;j<=4;j++)
            {fprintf(fpa,"%-3.6f",a[i][j]);
            fprintf(fpb,"%-3.6f",b[i][j]);
              }
            fprintf(fpa,"\n");
            fprintf(fpb,"\n");
              }
            fprintf(fpt,"%s%d,","THE Train time is:",Train);
            fprintf(fpe,"%s%f,","THE Mean error is:",error);

            fcloss(fpp);
            fcloss(fpa);
            fcloss(fpb);
            fcloss(fpt);
            fcloss(fpe);
return OK;
}

/* 生成测试样本的随机数*/
int create_rand()
{
int i;
randomize();
for(i=0;i<=T-1;i++)
{
t_xi[1]=(double) random(1500)/100;
t_xi[2]=(double) random(1500)/100;
t_xi[3]=(double) random(1500)/100;
t_xi[4]=(double) random(1500)/100;
}
return OK;
}

/* 用测试样本检验建模的效果并生成测试文件*/
test_sample()
{
FILE *fpt;
int i;
double y,x1,x2,x3,x4;
double Y;
double error,mean=0;

fpt=fopen("test.dat","wt");
create_rand();
i=0;
fprintf(fpt,"%s","    x1    x2    x3    x4        Predict        Real        Error";
While(i<T)
{
xi[1]=t_xi[1];
xi[2]=t_xi[2];
xi[3]=t_xi[3];
xi[4]=t_xi[4];

Y=Y_effect();

t_Y[i]=(0.472*((1+sin(M_PI*(xi[1]-0.5)))/2)
 +0.607*((1+sin(M_PI*(xi[2]-0.5)))/2)
 +0.779*((1+sin(M_PI*(xi[3]-0.5)))/2)
 +((1+sin(M_PI*(xi[4]-0.5)))/2)/2.858;

error=fabs(Y-t_Y[i]);
if(i<(T-1))
{fprintf(fpt,"\n%-d  %-5.2f    %-5.2f  %-5.2f  %-5.2f  %-5.6f  %-5.6f  %-5.6f","\n",
(i+1),t_X1[i],t_X2[i],t_X3[i],t_X4[i],Y,t_Y[i],error);}

else fprintf(fpt,"\n%-d  %-5.2f    %-5.2f  %-5.2f  %-5.2f  %-5.6f  %-5.6f  %-5.6f","\n",
(i+1),t_X1[i],t_X2[i],t_X3[i],t_X4[i],Y,t_Y[i],error);}

mean=mean+error;
i++;
}
fprintf(fpt,"\n Mean error is%-3.6f",mean/10);
fclose(fpt);
return OK;
}

/* 在浮点数字符串的位置插入小数点*/
insert_dot(char *str,int dec)
{
char c;
int i,j,len;
len=strlen(str);
if(dec>0)
{for(i=0;i<=(len-dec);i++)
 *(str+len-i+1)= *(str+len-i)
 *(str+dec)=".";
}

else if(i<0)
{j=abs(dec);
for(i=1;i<=len;i++)
 *(str+len-i+1+j)= *(str+len-i-1);
 *str="0",
 *(str+1)=".";

for(i=1;i<=j;i++)
 *(str+i+1)="0";
}

else
{for(i=0;i<=(len-dec);i++)
 *(str+len-i+2)= *(str+len-i)
 *(str+1)=".";
 *str="0";
}
return 0;

}

/*初始化用户界面的所有参数*/
init_graph_all_date(int count,int ndig,char *str)
{
int dign,dec;
setcolor(LIGHTRED);
settextjustify(0,2);
itoa(train,str,10);
outtextxy(train_x+50,train_y-4,str);
mean=1,0;
str=fcvt(mean,ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+50,error_y-4,str);
itoa(count,str,10);
outtextxy(act_barx2,act_bar_y2+10,str);
setcolor(YELLOW);

str=fcvt(X1[0],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-84,str);

str=fcvt(X2[0],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-84,str);

str=fcvt(X3[0],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-64,str);

str=fcvt(X4[0],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-64,str);
return OK;
}

/*show1_x1_x4和show2_x1_x4在用户界面中动态显示*/
/*样本数据X1,X2,X3,X4的值*/
int show1_x1_x4(int count,int ndig,char *str)
{
int dec,sign;
if(count>=2)
{
settextstyle(0,0,0);
srtcolor(getbkcolor());
srttextjustify(0,2);

str=fcvt(X1[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-84,str);

str=fcvt(X2[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-84,str);

str=fcvt(X3[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-64,str);

str=fcvt(X4[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-64,str);
}

return OK;
}

int show2_x1_x4(int count,int ndig,char *str)
{
int dec,sign;
settextstyle(0,0,0);
srttextjustify(0,2);
if((train>=1)&&(count==2))
{setcolor(getbkcolor());

str=fcvt(X1[79],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-84,str);

str=fcvt(X2[79],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-84,str);

str=fcvt(X3[79],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-64,str);

str=fcvt(X4[79],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-64,str);
}

setcolor(YELLOW);
str=fcvt(X1[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-84,str);

str=fcvt(X2[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-84,str);

str=fcvt(X3[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(train_x-10,train_y-64,str);

str=fcvt(X4[count-2],ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+20,error_y-64,str);

return OK;
}

/*显示活动条的当前运行状态*/
int show_act_bar(int count,float setp)
{
settextjustify(0,2);
setfillstyle(1,LIGHTGREEN);
bar(act_bar_x1,act_bar_y1,act_bar_x1+((count-1)*setp),act_bar_y2);
return OK;
}

/*show_mean1和show_mean2显示当前误差状态*/
int show_mean1(int ndig,char *str)
{
int dec,sign;
setcolor(getbkcolor());

str=fcvt(mean,ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+50,error_y-4,str);
return OK;
}

int show_mean2(int ndig,char *str)
{
int sign,dec;
settextstyle(0,0,0);

str=fcvt(mean,ndig,&dec,&sign);
insert_dot(str,dec);
outtextxy(error_x+50,error_y-4,str);
return OK;
}

/*在用户界面中动态显示样本组数*/
int show_sample_no(int count,char *str)
{settextstyle(0,0,0);
setcolor(getbkcolor());

itoa((count-2),str,10);
outtextxy(act_bar_x2,act_bar_y2+10,str);
if((train>=1)&&(count==2))
outtextxy(act_bar_x2,act_bar_y2+10,"80");
setcolor(LIGHTRED);

itoa((count-1),str,10);
outtextxy(act_bar_x2,act_bar_y2+10,str);
return OK;
}

/*初始化用户界面的活动条*/
int init_bar()
{
setfillstyle(1,WHITE);
bar(act_bar_x1,act_bar_y1,act_bar_x2,act_bar_y2);
return OK;
}

/*用户界面动态显示训练次数*/
int show_train(char *str)
{
setfillstyle(0,0,0);
setcolor(getbkcolor());

itoa(train,str,10);
outtextxy(train_x+50,train_y-4,str);
setcolor(LIGHTRED);

itoa((train+1),str,10);
outtextxy(train_x+50,train_y-4,str);
setcolor(LIGHTRED);
return OK;
}

/* 主函数*/
main()
{
unsigned i,j,k,l,m;
int count=0;
int modifiers;
int contin;
char *str;
int dec,sign,ndig=6;
float setp=(act_bar_x2-act_bar_x1)/SAMPLES;

long double ww;
long double wwy;

long double Y;
long double J;

long double temp_ai[5];

f[1]=uPL;
f[2]=uNL;
F[3]=uZo;

str=(char *) malloc(20);
}

/*初始化a[]b[]*/
for(i=0;i<=3;i++)
for(j=0;j<=3;j++)
for(k=0;k<=3;k++)
for(l=0;l<=3;l++)
{m=27*(i-1)+9*(j-1)+3*(k-1)+l;
a[m][1]=al_before[i];
a[m][2]=al_before[j];
a[m][3]=al_before[k];
a[m][4]=al_before[l];
b[m][1]=bl_before[i];
b[m][2]=bl_before[j];
b[m][3]=bl_before[k];
b[m][4]=bl_before[l];
}

initgra();
interface();
init_graph_all_data(count,ndig,str);
create_rand();

/*调整平均误差*/
while(mean>0.005)
{
xi[1]=X1[count=0];
xi[2]=X2[count];
xi[3]=X3[count];
xi[4]=X4[count];
yd[1]=YD[count++];

Y=Y_effect();

/*对80组样本调整误差*/
while(count<SAMPLES+1)
{
while(((fabs(Y-yd))*(fabs(Y-yd)))>0.000015)
{
ww=0;
wwy=0;
for(j=1;j<=M;j++)
{ww=ww+w[j];
wwy=wwy+w[j]*y[j];
}

for(i=1;i<=M;i++)
p[i][0]=p[i][0]-(q1*(Y-yd)*w[i])/ww;
p[i][1]=p[i][1]-(q1*(Y-yd)*w[i])/ww*xi[1];
p[i][2]=p[i][2]-(q1*(Y-yd)*w[i])/ww*xi[2];
p[i][3]=p[i][3]-(q1*(Y-yd)*w[i])/ww*xi[3];
p[i][4]=p[i][4]-(q1*(Y-yd)*w[i])/ww*xi[4];

for(j=1;j<=4;j++)
a[i][j]=a[i][j]-(q2*(Y-yd)*(y[i]*ww-wwy)*2*(xi[j]-(temp_ai[j]=a[i][j]))*w[i])/(b[i][j]*pow

(ww,2));
for(j=1;j<=4;j++)
b[i][j]=b[i][j]-(q2*(Y-yd)*(y[i]*ww-wwy)*pow((xi[j]-temp_ai[j]),2)*w[i])/(pow(b[i][j]*ww,2));

if((modifiers&LEFT)&&(modifiers&CTRL)&&(modifiers&ALT))
{
save_screen(ptr);
fflush(stdin);
 goto FINISH;}
}

BB:
Y=Y_effect();
}
show1_x1_x4(count,ndig,str)
xi[1]=X1[count];
xi[2]=X2[count];
xi[3]=X3[count];
xi[4]=X4[count];
yd=Yd[count++];
show2_x1_x4(count,ndig,str);
show_cat_bar(count,step);
show_sample_no(count,str);
}
show_mean1(ndig,str);
for(i=0;i<=SAMPLES;I++)
{
xi[1]=X1[i];
xi[2]=X2[i];
xi[3]=X3[i];
xi[4]=X4[i];
yd=Yd[i];
J=(yd-Y_effect());
mean=mean+fabs(J);
}

show_train(str);
mean=mean/SAMPLES;
show_mean2(ndig,str);
init_bar();
train++;
}

FINISH:
draw_graph1();
getch();
fflush(stdin);

draw_graph2();
getch();
fflush(stdin);

setlinestyle(SOLID_LINE,0,BLUE);
contin=ask();
if(contin==l)
{cleardevice();
setbkcolor(LIGHTBLUE);
restore_screen(ptr);
 goto BB;
}

closegraph();
save_pad();
test_sample();
cout<<"THe program is finised";

return OK;
}