Построение функции принадлежности нечеткого множества

#include <windows.h>
#include <stdio.h>
#define _USE_MATH_DEFINES
#include <math.h>
#include "resource.h"
 
HINSTANCE hInst;      // handle программы
const int DataLength=1000; // Длина массива
double x[DataLength],y[DataLength]; // сами массивы
 
void Graph(HWND hwnd) // функция построения графика
{
    HWND hwndHist=GetDlgItem(hwnd,IDC_GRAF);
    HDC hdc=GetDC(hwndHist);      // Контекст окна, на котором рисуюем график
    RECT rc;
    const int padd=35; // отступ от края окна
    GetClientRect(hwndHist,&rc); // размер окна
    IntersectClipRect(hdc,0,0,rc.right,rc.bottom); // ограничим область рисования этими размерами
    // Рисуем оси координат
    Rectangle(hdc,0,0,rc.right,rc.bottom);
    MoveToEx(hdc,padd,rc.top+padd,NULL);
    LineTo(hdc,padd,rc.bottom-padd);
    LineTo(hdc,rc.right-padd,rc.bottom-padd);
 
    SetTextAlign(hdc,TA_LEFT | TA_BASELINE);
    TextOut(hdc,padd/2, padd/2,"Y(X)",4);
// для определения масштаба найдем экстремумы
    int i;
    double dMaxx=x[0];
    double dMinx=x[0];
    double dMaxy=y[0];
    double dMiny=y[0];
    for(i=1;i<DataLength;i++)
    {
        dMaxx=max(dMaxx,x[i]);
        dMaxy=max(dMaxy,y[i]);
        dMinx=min(dMinx,x[i]);
        dMiny=min(dMinx,y[i]);
    }
// вывод меток по осям  
    const int COLLS=5;
    SetTextAlign(hdc,TA_LEFT | TA_BASELINE);
    int labWidth=(rc.right-2*padd)/COLLS;
    int labHeight=(rc.bottom-2*padd)/COLLS;
    double Deltay=(dMaxy-dMiny)/COLLS;
    for(i=0;i<COLLS;i++)
    {
        double t=dMiny+Deltay*(i+1);
        char buf[16];
        sprintf(buf,"%4lg\0",t);
        int yp=rc.bottom-padd-(i+1)*labHeight;
        TextOut(hdc,0,yp,buf,strlen(buf));
        MoveToEx(hdc,padd-3,yp,NULL);
        LineTo(hdc,padd+3,yp);
    }
 
    SetTextAlign(hdc,TA_RIGHT | TA_BASELINE);
    double Deltax=(dMaxx-dMinx)/COLLS;
    for(i=0;i<COLLS;i++)
    {
        double t=dMinx+Deltax*(i+1);
        char buf[16];
        sprintf(buf,"%4lg\0",t);
        int xp=padd+labWidth*(i+1);
        TextOut(hdc,xp,rc.bottom-padd/2,buf,strlen(buf));
        MoveToEx(hdc,xp, rc.bottom-padd,NULL);
        LineTo(hdc,xp, rc.bottom-padd+5);
    }
    SetTextAlign(hdc,TA_LEFT | TA_BASELINE);
    TextOut(hdc,rc.right-padd+1,rc.bottom-padd+5,"X",1);
 
    HPEN p1=CreatePen(PS_SOLID,1,RGB(0,128,0)); // рисуем пером зеленого цвета
    HPEN pOld=(HPEN)SelectObject(hdc,p1);
    int xp,yp;
    if(dMaxx==dMinx)
        xp=padd;
    else
        xp=padd+int((x[0]-dMinx)/(dMaxx-dMinx)*(rc.right-2*padd)); // пересчитали начальное значение x
    if(dMaxy==dMiny)
        yp=rc.bottom-padd;
    else
        yp=rc.bottom-padd-int((y[0]-dMiny)/(dMaxy-dMiny)*(rc.bottom-2*padd));// пересчитали начальное значение y
    MoveToEx(hdc,xp,yp,NULL);
    for(i=0;i<DataLength;i++)
    {
        if(dMaxx==dMinx)
            xp=padd;
        else
            xp=padd+int((x[i]-dMinx)/(dMaxx-dMinx)*(rc.right-2*padd));// пересчитали очередное значение x
        if(dMaxy==dMiny)
            yp=rc.bottom-padd;
        else
            yp=rc.bottom-padd-int((y[i]-dMiny)/(dMaxy-dMiny)*(rc.bottom-2*padd));// пересчитали очередное значение y
 
        LineTo(hdc,xp,yp);           //Провели линию от i до i+1й точек
    }
    SelectObject(hdc,pOld);
    DeleteObject(p1);
    ReleaseDC(hwndHist,hdc);
}
 
// заполнение массивов x,y
void ReadDataFile(HWND hwnd)
{
 
    char buf[32];
    GetDlgItemText(hwnd,IDC_MAXX,buf,32);
    double max_x;sscanf(buf,"%lf",&max_x);  // получили значение max_x из поля IDC_MAXX диалога
    GetDlgItemText(hwnd,IDC_A,buf,32);
    double a;sscanf(buf,"%lf",&a);                 // получили значение a из поля IDC_A диалога
    GetDlgItemText(hwnd,IDC_B,buf,32);
    double b;sscanf(buf,"%lf",&b);
    for(int i=0;i<DataLength;i++)
    {
// вычисление x[i] и y[i]
        x[i]=double(i)/DataLength*max_x; 
        if(x[i]<a)
        {
            y[i]=1.0;
        }else if(x[i]<b)
        {
            y[i]=0.5*(1.0+cos((x[i]-a)/(b-a)*M_PI));
        }else
        {
            y[i]=0;
        }
    }
}
 
//  функция обработчик сообщений окна
BOOL WINAPI fnObr(HWND hwnd,WORD msg,WPARAM wp,LPARAM lp)
{
    switch(msg)
    {
    case WM_COMMAND: // нажатие кнопок и пунктов меню
    {
        switch(LOWORD(wp))
        {
        case ID_HELP_ABOUT: // меню about
            MessageBox(hwnd,
            "Программа для вывода функции принадлежности\n(С) Задорожный С.С. 2014",
            "О программе",MB_OK);
            return 1;
        case ID_CREATE: // кнопка Построить
        {
            ReadDataFile(hwnd); // заполнить X,y
            InvalidateRect(hwnd,NULL,true); // перерисовать окно
            return 1;
        }
        case ID_FILE_EXIT: // кнопка Выход
            SendMessage(hwnd,WM_CLOSE,0,0);
            return 1;
        }// end switch WM_COMMAND
        return 0; // остальные команды не обрабатываем
    }
    case WM_INITDIALOG: // открытие окна
        SetClassLong(hwnd,GCL_HICON,(LONG)LoadIcon(hInst,"STAT")); // установить логотип программы
        SetDlgItemText(hwnd,IDC_A," 0.5"); // заполнить параметры стандартными значениями
        SetDlgItemText(hwnd,IDC_B," 0.75");
        SetDlgItemText(hwnd,IDC_MAXX," 1.0");
        ReadDataFile(hwnd); // заполнить x,y
        return 1;
    case WM_SIZE: // изменение размеров окна
        InvalidateRect(hwnd,NULL,true); // перерисовать
        return 1;
    case WM_PAINT: // сообщение о необходимости перерисовки окна
    {
        PAINTSTRUCT ps;
        HDC hdc=BeginPaint(hwnd,&ps);
        Graph(hwnd);   //  нарисовать график
        EndPaint(hwnd,&ps);
        return 1;
 
    }
    case WM_CLOSE: // закрытие окна
        EndDialog(hwnd,1);
        return 1;
    } // end switch
    return 0;
}
 
 
// точка входа в программу
int WINAPI WinMain(HINSTANCE hi,HINSTANCE,LPSTR,int)
{
    hInst=hi;
    return DialogBox(hi,"STAT32",0,(DLGPROC)fnObr); // открываем диалог "STAT32"
}