📄 fcmadlg.cpp
字号:
// FCMADlg.cpp : implementation file
//
#include "stdafx.h"
#include "FCMA.h"
#include "FCMADlg.h"
#include "stdlib.h" //rand() requie
#include "math.h"
#include "windowsx.h"
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif
/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About
class CAboutDlg : public CDialog
{
public:
CAboutDlg();
// Dialog Data
//{{AFX_DATA(CAboutDlg)
enum { IDD = IDD_ABOUTBOX };
//}}AFX_DATA
// ClassWizard generated virtual function overrides
//{{AFX_VIRTUAL(CAboutDlg)
protected:
virtual void DoDataExchange(CDataExchange* pDX); // DDX/DDV support
//}}AFX_VIRTUAL
// Implementation
protected:
//{{AFX_MSG(CAboutDlg)
//}}AFX_MSG
DECLARE_MESSAGE_MAP()
};
CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
//{{AFX_DATA_INIT(CAboutDlg)
//}}AFX_DATA_INIT
}
void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
CDialog::DoDataExchange(pDX);
//{{AFX_DATA_MAP(CAboutDlg)
//}}AFX_DATA_MAP
}
BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
//{{AFX_MSG_MAP(CAboutDlg)
// No message handlers
//}}AFX_MSG_MAP
END_MESSAGE_MAP()
/////////////////////////////////////////////////////////////////////////////
// CFCMADlg dialog
CFCMADlg::CFCMADlg(CWnd* pParent /*=NULL*/)
: CDialog(CFCMADlg::IDD, pParent)
{
//{{AFX_DATA_INIT(CFCMADlg)
m_IFilepath = _T("");
m_QFilepath = _T("");
m_Umatrix = _T("");
m_corematrix = _T("");
m_process = _T("");
m = 2.0f;
threshold = 1.0f;
NodeNum = 500;
repeatNum = 30;
m_result = _T("当前测试类数: 0");
m_MaxNodeNum = 4000;
m_MaxCoreNum = 256;
//}}AFX_DATA_INIT
// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}
void CFCMADlg::DoDataExchange(CDataExchange* pDX)
{
CDialog::DoDataExchange(pDX);
//{{AFX_DATA_MAP(CFCMADlg)
DDX_Text(pDX, IDC_EDIT1, m_IFilepath);
DDX_Text(pDX, IDC_EDIT2, m_QFilepath);
DDX_Text(pDX, IDC_EDIT3, m_Umatrix);
DDX_Text(pDX, IDC_EDIT4, m_corematrix);
DDX_Text(pDX, IDC_STATIC2, m_process);
DDX_Text(pDX, IDC_EDIT5, m);
DDX_Text(pDX, IDC_EDIT6, threshold);
DDX_Text(pDX, IDC_EDIT7, NodeNum);
DDX_Text(pDX, IDC_EDIT9, repeatNum);
DDX_Text(pDX, IDC_STATIC12, m_result);
DDX_Text(pDX, IDC_EDIT8, m_MaxNodeNum);
DDX_Text(pDX, IDC_EDIT10, m_MaxCoreNum);
//}}AFX_DATA_MAP
}
BEGIN_MESSAGE_MAP(CFCMADlg, CDialog)
//{{AFX_MSG_MAP(CFCMADlg)
ON_WM_SYSCOMMAND()
ON_WM_PAINT()
ON_WM_QUERYDRAGICON()
ON_BN_CLICKED(IDC_BUTTON1, OnGetIdata)
ON_BN_CLICKED(IDC_BUTTON2, OnGetQdata)
ON_BN_CLICKED(IDOK, OnStart)
//}}AFX_MSG_MAP
END_MESSAGE_MAP()
/////////////////////////////////////////////////////////////////////////////
// CFCMADlg message handlers
BOOL CFCMADlg::OnInitDialog()
{
CDialog::OnInitDialog();
// Add "About..." menu item to system menu.
// IDM_ABOUTBOX must be in the system command range.
ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
ASSERT(IDM_ABOUTBOX < 0xF000);
CMenu* pSysMenu = GetSystemMenu(FALSE);
if (pSysMenu != NULL)
{
CString strAboutMenu;
strAboutMenu.LoadString(IDS_ABOUTBOX);
if (!strAboutMenu.IsEmpty())
{
pSysMenu->AppendMenu(MF_SEPARATOR);
pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
}
}
// Set the icon for this dialog. The framework does this automatically
// when the application's main window is not a dialog
SetIcon(m_hIcon, TRUE); // Set big icon
SetIcon(m_hIcon, FALSE); // Set small icon
//===============>分配内存块<================//**0909
UpdateData(true);
Idata=(float*)GlobalAllocPtr(GHND,m_MaxNodeNum*sizeof(float));
Qdata=(float*)GlobalAllocPtr(GHND,m_MaxNodeNum*sizeof(float));
nodedata=(Node *)GlobalAllocPtr(GHND,m_MaxNodeNum*sizeof(Node));
coredata=(Core *)GlobalAllocPtr(GHND,m_MaxNodeNum*sizeof(Core));
Umatrix=(float*)GlobalAllocPtr(GHND,m_MaxNodeNum*m_MaxCoreNum*sizeof(Node));
CorematrixSave=(Core *)GlobalAllocPtr(GHND,m_MaxNodeNum*sizeof(Core));
//UmatrixSave1=(float*)GlobalAllocPtr(GHND,NodeNum*m_MaxCoreNum*sizeof(Node));
//--------------------------------------------|
return TRUE; // return TRUE unless you set the focus to a control
}
void CFCMADlg::OnSysCommand(UINT nID, LPARAM lParam)
{
if ((nID & 0xFFF0) == IDM_ABOUTBOX)
{
CAboutDlg dlgAbout;
dlgAbout.DoModal();
}
else
{
CDialog::OnSysCommand(nID, lParam);
}
}
// If you add a minimize button to your dialog, you will need the code below
// to draw the icon. For MFC applications using the document/view model,
// this is automatically done for you by the framework.
void CFCMADlg::OnPaint()
{
if (IsIconic())
{
CPaintDC dc(this); // device context for painting
SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);
// Center icon in client rectangle
int cxIcon = GetSystemMetrics(SM_CXICON);
int cyIcon = GetSystemMetrics(SM_CYICON);
CRect rect;
GetClientRect(&rect);
int x = (rect.Width() - cxIcon + 1) / 2;
int y = (rect.Height() - cyIcon + 1) / 2;
// Draw the icon
dc.DrawIcon(x, y, m_hIcon);
}
else
{
CDialog::OnPaint();
}
}
// The system calls this to obtain the cursor to display while the user drags
// the minimized window.
HCURSOR CFCMADlg::OnQueryDragIcon()
{
return (HCURSOR) m_hIcon;
}
void CFCMADlg::OnGetIdata() //获得I路数据文件地址 : I、Q路数据均为float型!
{
read=false;
CFileDialog Idatadlg(true);
if(Idatadlg.DoModal()==IDOK)
{
m_IFilepath=Idatadlg.GetPathName();
UpdateData(false);
}
else
return;
}
void CFCMADlg::OnGetQdata() //获得Q路数据文件地址
{
read=false;
CFileDialog Qdatadlg(true);
if(Qdatadlg.DoModal()==IDOK)
{
m_QFilepath=Qdatadlg.GetPathName();
UpdateData(false);
}
else
return ;
}
void CFCMADlg::OnStart() //=========>聚类过程开始
{
// TODO: 一些参数、数据的初始化 **
C=2;
K=-1.0;
Km.maxK=0.0;
Km.c=C;
savemax=0.0;
firstcutpeak=false;
float limit=0.0; //阀值
float dij=0.0; //(Xi,Vj)间距离
float dkj=0.0;
float sum=0.0;
float dij2=0.0; //计算J时
int repcount=0;
CString coreinitsave="";
CString wholeinitcore=""; //0909
CString save=""; //0909
UpdateData(true); //得到参数设置
/*===============>分配内存块<================//**0909
Idata1=(float*)GlobalAllocPtr(GHND,NodeNum*sizeof(float));
Qdata1=(float*)GlobalAllocPtr(GHND,NodeNum*sizeof(float));
nodedata1=(Node *)GlobalAllocPtr(GHND,NodeNum*sizeof(Node));
coredata1=(Core *)GlobalAllocPtr(GHND,NodeNum*sizeof(Core));
Umatrix1=(float*)GlobalAllocPtr(GHND,NodeNum*m_MaxCoreNum*sizeof(Node));
UmatrixSave1=(float*)GlobalAllocPtr(GHND,NodeNum*m_MaxCoreNum*sizeof(Node));
CorematrixSave1=(Core *)GlobalAllocPtr(GHND,NodeNum*sizeof(Core));*/
/*for(int a=0;a<MAXCoreNum;a++) //**U初始化
{
for(int b=0;b<NodeNum;b++)
{
Umatrix[a][b]=0.0;
}
}*/
for(int a=0;a<m_MaxCoreNum;a++) //**U初始化0909
{
for(int b=0;b<NodeNum;b++)
{
*(Umatrix+(a*NodeNum+b))=0.0f;
}
}
//===============>参数检查<==================//
if(m_IFilepath==""||m_QFilepath=="")
{
AfxMessageBox("请输入数据!");
return ;
}
if(m<1.5||m>2.5)
{
AfxMessageBox("权重指数m范围为[1.5，2.5]");
return;
}
if(repeatNum==NULL||threshold==NULL)
{
AfxMessageBox("参数未设置完全!请重试!");
return;
}
if(NodeNum>m_MaxNodeNum) //**0909
{
CString a="";
a.Format("最大数据数为：%d,请修改参数!",m_MaxNodeNum);
AfxMessageBox(a);
return;
}
//-------------------------------------------|
//===============>读取文件数据<===============//
if(read==false)
{
//-----I路数据----
int Ifilelength;
CFile idatafile;
idatafile.Open(m_IFilepath,CFile::modeRead,NULL);
Ifilelength=idatafile.GetLength();//得到文件总字节数 //47976
if((int)(Ifilelength/sizeof(float))<NodeNum)
{
AfxMessageBox("文件长度不够大!"); //数据样点数不够
return;
}
int Inum;
idatafile.SeekToBegin();
Inum=idatafile.Read(Idata,NodeNum*sizeof(float));//**0909
idatafile.Close();
//----------------|
//-----Q路数据-----
int Qfilelength;
CFile qdatafile;
qdatafile.Open(m_QFilepath,CFile::modeRead,NULL);
Qfilelength=qdatafile.GetLength();//得到文件总字节数
if((int)(Qfilelength/sizeof(float))<NodeNum)
{
AfxMessageBox("文件长度不够大!"); //数据样点数不够
return;
}
int Qnum;
qdatafile.SeekToBegin();
Qnum=qdatafile.Read(Qdata,NodeNum*sizeof(float));//**0909
qdatafile.Close();
read=true;
}
//---------------------------------------------------|
m_Umatrix="";
m_corematrix="";
m_process.Format("正在进行聚类，请等待...");
UpdateData(false);
//==============>给各数据点赋值<================//**0909
for(int p=0;p<NodeNum;p++)
{
(*(nodedata+p)).x=*(Idata+p);
(*(nodedata+p)).y=*(Qdata+p);
(*(nodedata+p)).id=-1; // -1: 初始时不属于任何类
(*(nodedata+p)).density=0.0; //0.0: 密度指标初始化
}
CString coreinit;
coreinit.Format("**********************\r\n初始随机聚类中心数：%d\r\n中心坐标:\r\n",C);
CString temp="",temp0="",temp3="",temp4="";
while(C<=m_MaxCoreNum) //**0909
{
temp="";
temp0="";
repcount=0;
//e=0.0;
m_result="";
coreinit="**************\r\n初始随机聚类中心坐标\r\n";
save.Format("*********C=%d 时的初始中心坐标:\r\n",C);//0909
wholeinitcore+=save;
m_result.Format("当前测试类数： %d",C);
//sss=m_result;
UpdateData(false);
//*============>减法聚类初始化中心<===============//
if(C==2)
firstcutpeak=true;
else
firstcutpeak=false;
CutPeak(C,firstcutpeak); //提取C个初始中心
for(int i3=0;i3<C;i3++) //将原来保存的密度最大的数据点传给中心
{
(*(coredata+i3)).x=coreXsave[i3];
(*(coredata+i3)).y=coreYsave[i3];
}
for(int i2=0;i2<C;i2++)
{
temp.Format("%f %f\r\n",(*(coredata+i2)).x,(*(coredata+i2)).y);
temp0+=temp;
}
//------------------------------------------------|
coreinit+=temp0;
wholeinitcore+=temp0; //0909+
//--------------------------------------------------|
do{
limit=0.0; //每次迭代都要归0
//==================>计算C时的U阵<====================//
for(int i=0;i<C;i++)
{
for(int j=0;j<NodeNum;j++)
{
dij=Caldistance(0,i,j,0,0,0); //dij:是距离的平方值，便于在计算中无需开方
//=====>计算dkj<======//
for(int k=0;k<C;k++)
{
if(dij==0.0) //如果数据点就是中心点，不用计算dkj
{
for(int l=0;l<C;l++)
{
if(l==i)
(*(Umatrix+(l*NodeNum+j)))=1.0f;
else
(*(Umatrix+(l*NodeNum+j)))=0.0f;
}
break;
}
else //如果不是
{
dkj=Caldistance(0,k,j,0,0,0);
sum+=(float)pow((double)(dij/dkj),(double)(1/(m-1)));
}
}//k
//--------------------|
if(dij!=0.0)
(*(Umatrix+(i*NodeNum+j)))=(float)(1/sum);
sum=0.0f;
}
}
//--------------------------------------------------|ok!
//===============>计算c个新的聚类中心<================//
if(CalCorepos()==false)
return ;
//==================>计算目标函数<====================//
for(int a=0;a<C;a++)
{
for(int b=0;b<NodeNum;b++)
{
dij2=Caldistance(0,a,b,0,0,0);
limit+=(float)(pow((double)(*(Umatrix+(a*NodeNum+b))),(double)m)*dij2);
}
}
preJ=limit;
//==================>计算目标函数差<==================//
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -