foremddlg.cpp

vc代用matcom编程

// ForEmdDlg.cpp : implementation file
//

#include "stdafx.h"
#include "ForEmd.h"
#include "ForEmdDlg.h"
#include "matlib.h"
#include "math.h"
#include "emd.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()

/////////////////////////////////////////////////////////////////////////////
// CForEmdDlg dialog

CForEmdDlg::CForEmdDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CForEmdDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CForEmdDlg)
	m_edit = 0;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
	mydata=Alloc2DArr(4,2048);
}

void CForEmdDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CForEmdDlg)
	DDX_Text(pDX, IDC_EDIT1, m_edit);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CForEmdDlg, CDialog)
	//{{AFX_MSG_MAP(CForEmdDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(ID_FOR_EMD, OnForEmd)
	ON_WM_DESTROY()
	ON_BN_CLICKED(IDC_BUTTON1, OnButton1)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CForEmdDlg message handlers

BOOL CForEmdDlg::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
	
	// TODO: Add extra initialization here
	
	CRect Rect;
	GetClientRect(&Rect);
	Rect.right = Rect.left+495;
	Rect.bottom = Rect.top +325;
	Rect.left = Rect.left+15  ;
	Rect.top = Rect.top + 15;
	//Allocate space for series .
	if ( !m_Chart.AllocSerie(4096*2) ) {
		AfxMessageBox("Error allocating chart serie") ;
        return 0;
	}
    m_Chart.SetChartTitle("拟合后幅值谱图");
	m_Chart.SetGridXYNumber(5,2);
	m_Chart.SetAxisStyle(0);
	m_Chart.SetRange(0,100,0,50);
	m_Chart.SetAxisStyle(0);
	m_Chart.mpSerie[0].m_plotColor = RGB(255,255,255);
	m_Chart.SetGridXYNumber(10,10);
	//m_Chart2d.bLogScale = TRUE ;
    m_Chart.Create(WS_CHILD|WS_VISIBLE,Rect,this,IDC_CHART1);
	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CForEmdDlg::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 CForEmdDlg::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 CForEmdDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CForEmdDlg::OnForEmd() 
{
	// TODO: Add your control notification handler code here
    

	initM(MATCOM_VERSION);//初始化matcom库
	Mm x=zeros(1,2048);
	Mm a,b;
	dMm(y);
	for(int i=1;i<=2048;i++)
	{
		x.r(1,i)=10*sin(10*i/1024+1);//给x阵赋值
	}
    emd(x,i_o,y,a,b);
	for(i=1;i<=4;i++)
	for(int j=1;j<=2048;j++)
        mydata[i-1][j-1]=y.r(i,j);
	exitM();

}
double **CForEmdDlg::Alloc2DArr(int row, int col)
{
	//分配二维数组空间
	int i;
	double *x,**y;
	
	x=(double*)calloc(row*col,sizeof(double));	
	for(i=0;i<row*col;i++)	x[i]=0;	//赋初值	
	y=(double**)calloc(row,sizeof(double*));
	for( i=0;i<row;i++)	y[i]=&x[col*i];	
	return y;
}
void CForEmdDlg::Free2DArr(double **arr)
{
	//收回二维数组空间
	if(arr==NULL) return;
	free(arr[0]);
	free(arr);
}
void CForEmdDlg::OnDestroy() 
{
	CDialog::OnDestroy();
	
	// TODO: Add your message handler code here
    Free2DArr(mydata);
	m_Chart.ClearChart();
    

}

void CForEmdDlg::OnButton1() 
{
	// TODO: Add your control notification handler code here
	UpdateData();
	int i;
	m_Chart.SetRange(0,2048,-20,20);
	if(m_edit<4)
	{for(i=0;i<2048;i++)
		m_Chart.SetXYValueA(i,mydata[m_edit][i],i,0);
    	m_Chart.Invalidate();
	}
	
}