fftdrawingdlg.cpp

实现快速傅立叶算法

// FftDrawingDlg.cpp : implementation file
//
#include "windows.h"
#include "stdafx.h"
#include "FftDrawing.h"
#include "FftDrawingDlg.h"

////自己添加.
#include"math.h"
//#include "fstream.h"
#include "stdlib.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()

/////////////////////////////////////////////////////////////////////////////
// CFftDrawingDlg dialog

CFftDrawingDlg::CFftDrawingDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CFftDrawingDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CFftDrawingDlg)
	m_njishu = 0;
	m_caidianpinlv = 0.0;
	m_f1 = 0.0;
	m_f2 = 0.0;
	m_A = 0.0;
	m_B = 0.0;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	

	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);

}




void CFftDrawingDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CFftDrawingDlg)
	DDX_Text(pDX, IDC_JISHU, m_njishu);
	DDV_MinMaxInt(pDX, m_njishu, 1, 9);
	DDX_Text(pDX, IDC_CAIDIANPINLV, m_caidianpinlv);
	DDV_MinMaxDouble(pDX, m_caidianpinlv, 1., 10000.);
	DDX_Text(pDX, IDC_F1, m_f1);
	DDV_MinMaxDouble(pDX, m_f1, 0., 10000.);
	DDX_Text(pDX, IDC_F2, m_f2);
	DDV_MinMaxDouble(pDX, m_f2, 0., 10000.);
	DDX_Text(pDX, IDC_FUZHIA, m_A);
	DDV_MinMaxDouble(pDX, m_A, 1.e-003, 1000.);
	DDX_Text(pDX, IDC_FUZHIB, m_B);
	DDV_MinMaxDouble(pDX, m_B, 1.e-003, 1000.);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CFftDrawingDlg, CDialog)
	//{{AFX_MSG_MAP(CFftDrawingDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_SHIBU, OnShibu)
	ON_BN_CLICKED(IDC_SHUJUDAORU, OnShujudaoru)
	ON_BN_CLICKED(IDC_XUBU, OnXubu)
	ON_BN_CLICKED(IDC_XIANGJIAO, OnXiangjiao)
	ON_BN_CLICKED(IDC_MO, OnMo)
	ON_BN_CLICKED(IDC_GONGLVPU, OnGonglvpu)
	//}}AFX_MSG_MAP
	ON_BN_CLICKED(IDC_IsCaoed, &CFftDrawingDlg::OnBnClickedIscaoed)
	ON_BN_CLICKED(IDC_Cao, &CFftDrawingDlg::OnBnClickedCao)
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CFftDrawingDlg message handlers

BOOL CFftDrawingDlg::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
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CFftDrawingDlg::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 CFftDrawingDlg::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();
	}
    */
   
	CPaintDC dc(this);
	CWnd *pWnd = GetDlgItem(IDC_STATIC);
	pWnd->UpdateWindow();
	CDC *pDC = pWnd->GetDC();
	CPen gridPen(0,0,RGB(250,0,250));
	CPen  *oldPen=pDC->SelectObject(&gridPen);

	CRect rcClient;
	pWnd->GetClientRect(rcClient);
	pDC->Rectangle(rcClient);
  
    int dx=rcClient.Width();
	int dy=rcClient.Height();
    
    //CWnd *pWnd = GetDlgItem(IDC_STATIC);
    //CDC *pDC = pWnd->GetDC();
    //CRect rcClient;
	//pWnd->GetClientRect(rcClient);
	//pDC->Rectangle(rcClient);
	//int dx=rcClient.Width();
	//int dy=rcClient.Height();

    int zhongxian=dy/2;
	int zhongxian2=dx/2;

	pDC->MoveTo(20,dy-zhongxian);
	pDC->LineTo(546,dy-zhongxian);
	pDC->MoveTo(20,dy);
	pDC->LineTo(20,0);
	
    pDC->SelectObject(oldPen);
    
	//////设置画笔属性
	CPen gridPen2(0,0,RGB(0,0,0));
	CPen  *oldPen2=pDC->SelectObject(&gridPen2);
    /////////////
    switch(shuomingwenzi)
	{
	case 1:
		pDC->TextOut(dx-60,dy-30,"实部图");
		break;
	case 2:
		pDC->TextOut(dx-60,dy-30,"虚部图");
		break;
	case 3:
		pDC->TextOut(dx-60,dy-30,"模部图");
		break;
    case 4:
		pDC->TextOut(dx-60,dy-30,"相角图");
		break;
	case 5:
		pDC->TextOut(dx-60,dy-30,"功率谱图");
		break;
	default:
		;
	}

    
    /////以下画极值点的横坐标:
    pDC->SetPixel(20+m_f1*512/m_caidianpinlv+1,zhongxian,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv-1,zhongxian,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv,zhongxian+1,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv,zhongxian-1,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv+1,zhongxian+1,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv+1,zhongxian-1,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv-1,zhongxian+1,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv-1,zhongxian-1,2500);
	pDC->SetPixel(20+m_f1*512/m_caidianpinlv,zhongxian,2500);
   

	      ///////////////////写数据值
	char zifu[100];
	itoa(m_f1,zifu,10);
	pDC->TextOut(20+m_f1*512/m_caidianpinlv+3,zhongxian+3,zifu);
	
    



	pDC->SetPixel(20+m_f2*512/m_caidianpinlv+1,zhongxian,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv-1,zhongxian,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv,zhongxian+1,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv,zhongxian-1,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv+1,zhongxian+1,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv+1,zhongxian-1,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv-1,zhongxian+1,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv-1,zhongxian-1,2500);
	pDC->SetPixel(20+m_f2*512/m_caidianpinlv,zhongxian,2500);
    
        /////////写数据值
	itoa(m_f2,zifu,10);
	pDC->TextOut(20+m_f2*512/m_caidianpinlv+3,zhongxian+3,zifu);


	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv+1,zhongxian,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv-1,zhongxian,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv,zhongxian+1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv,zhongxian-1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv+1,zhongxian+1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv+1,zhongxian-1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv-1,zhongxian+1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv-1,zhongxian-1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv,zhongxian,2500);


	    /////////写数据值
	itoa((m_caidianpinlv-m_f1),zifu,10);
	pDC->TextOut(20+(m_caidianpinlv-m_f1)*512/m_caidianpinlv+3,zhongxian+3,zifu);




	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv+1,zhongxian,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv-1,zhongxian,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv,zhongxian+1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv,zhongxian-1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv+1,zhongxian+1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv+1,zhongxian-1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv-1,zhongxian+1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv-1,zhongxian-1,2500);
	pDC->SetPixel(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv,zhongxian,2500);

	////////写数据值
	itoa((m_caidianpinlv-m_f2),zifu,10);
	pDC->TextOut(20+(m_caidianpinlv-m_f2)*512/m_caidianpinlv+3,zhongxian+3,zifu);



    
	pDC->SetPixel(20+512+1,zhongxian,2500);
	pDC->SetPixel(20+512-1,zhongxian,2500);
	pDC->SetPixel(20+512,zhongxian+1,2500);
	pDC->SetPixel(20+512,zhongxian-1,2500);
	pDC->SetPixel(20+512+1,zhongxian+1,2500);
	pDC->SetPixel(20+512+1,zhongxian-1,2500);
	pDC->SetPixel(20+512-1,zhongxian+1,2500);
	pDC->SetPixel(20+512-1,zhongxian-1,2500);
	pDC->SetPixel(20+512,zhongxian,2500);

    ////////写数据值
	itoa(m_caidianpinlv,zifu,10);
	pDC->TextOut(20+512+3,zhongxian+3,zifu);
	pDC->TextOut(10,zhongxian+3,48);



    ////////////////////////////////////以上画极值点的横坐标


    //////////////////////////////以下画极值点的纵坐标
         ///////画f1处的极值纵坐标
	for(int js_jizhi=0;js_jizhi<N;js_jizhi++)
	{
	    if(js_jizhi==int(m_f1*N/m_caidianpinlv+0.5))
		{
			
			pDC->SetPixel(20,dy-zhongxian-shuju[js_jizhi],2500);
			pDC->SetPixel(20+1,dy-zhongxian-shuju[js_jizhi],2500);
			pDC->SetPixel(20,dy-zhongxian-shuju[js_jizhi]+1,2500);
			pDC->SetPixel(20-1,dy-zhongxian-shuju[js_jizhi],2500);
			pDC->SetPixel(20,dy-zhongxian-shuju[js_jizhi]-1,2500);
			pDC->SetPixel(20+1,dy-zhongxian-shuju[js_jizhi]+1,2500);
			pDC->SetPixel(20+1,dy-zhongxian-shuju[js_jizhi]-1,2500);
			pDC->SetPixel(20-1,dy-zhongxian-shuju[js_jizhi]+1,2500);
			pDC->SetPixel(20-1,dy-zhongxian-shuju[js_jizhi]-1,2500);
		
			
			////////写数据值
			double xieshuzi=shuju[js_jizhi];
			if(fuhao==-1)
			{
				for(int n_xieshuzi=0;n_xieshuzi<cishu;n_xieshuzi++)
					xieshuzi=xieshuzi*2;
				itoa(xieshuzi,zifu,10);
			   pDC->TextOut(20+3,dy-zhongxian-shuju[js_jizhi]+3,zifu);
			}	
			
			else if(fuhao==1)
			{
                
				for(int n_xieshuzi=0;n_xieshuzi<cishu;n_xieshuzi++)
					xieshuzi=xieshuzi/2;
				
				if(xieshuzi<100)
				{
					xieshuzi=xieshuzi*100;
					itoa(xieshuzi,zifu,10);
					char xiaoyu_10[10]={'*','e','-','1','0','0','\0'};
					strcat(zifu,xiaoyu_10);
				}
				else
					itoa(xieshuzi,zifu,10);
				pDC->TextOut(20+3,dy-zhongxian-shuju[js_jizhi]+3,zifu);
			    
			}
			else if(fuhao==0)
			{
				itoa(xieshuzi,zifu,10);	               
			    pDC->TextOut(20+3,dy-zhongxian-shuju[js_jizhi]+3,zifu);
			}
		
			/////////////////以上写数值



		}
	}
         //////画f2处的极值的纵坐标
	for(js_jizhi=0;js_jizhi<N;js_jizhi++)
	{
	    if(js_jizhi==int(m_f2*N/m_caidianpinlv+0.5))
		{
			pDC->SetPixel(20,dy-zhongxian-shuju[js_jizhi],2500);
			pDC->SetPixel(20+1,dy-zhongxian-shuju[js_jizhi],2500);
			pDC->SetPixel(20,dy-zhongxian-shuju[js_jizhi]+1,2500);
			pDC->SetPixel(20-1,dy-zhongxian-shuju[js_jizhi],2500);
			pDC->SetPixel(20,dy-zhongxian-shuju[js_jizhi]-1,2500);
			pDC->SetPixel(20+1,dy-zhongxian-shuju[js_jizhi]+1,2500);
			pDC->SetPixel(20+1,dy-zhongxian-shuju[js_jizhi]-1,2500);
			pDC->SetPixel(20-1,dy-zhongxian-shuju[js_jizhi]+1,2500);
			pDC->SetPixel(20-1,dy-zhongxian-shuju[js_jizhi]-1,2500);
			

			////////写数据值
			double xieshuzi=shuju[js_jizhi];
			if(fuhao==-1)
			{
				for(int n_xieshuzi=0;n_xieshuzi<cishu;n_xieshuzi++)
					xieshuzi=xieshuzi*2;
				itoa(xieshuzi,zifu,10);
			   pDC->TextOut(20+3,dy-zhongxian-shuju[js_jizhi]+3,zifu);
			}