sin sindlg.cpp

使用VC2005编写的自己生成不同频率的数据

// sin sinDlg.cpp : implementation file
//

#include "stdafx.h"
#include "sin sin.h"
#include "sin sinDlg.h"
#include "math.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#endif


// CAboutDlg dialog used for App About

class CAboutDlg : public CDialog
{
public:
	CAboutDlg();

	// Dialog Data
	enum { IDD = IDD_ABOUTBOX };

protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support

	// Implementation
protected:
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD) 

{
}

void CAboutDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
}

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
END_MESSAGE_MAP()


// CsinsinDlg dialog




CsinsinDlg::CsinsinDlg(CWnd* pParent /*=NULL*/)
: CDialog(CsinsinDlg::IDD, pParent)
, m_fd(20)
, m_pl(1)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CsinsinDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT2, m_fd);
	DDX_Text(pDX, IDC_EDIT1, m_pl);
}

BEGIN_MESSAGE_MAP(CsinsinDlg, CDialog)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	//}}AFX_MSG_MAP
	ON_BN_CLICKED(IDC_BUTTON1, &CsinsinDlg::DrawSin)
	ON_BN_CLICKED(IDC_BUTTON2, &CsinsinDlg::cancel)

END_MESSAGE_MAP()


// CsinsinDlg message handlers

BOOL CsinsinDlg::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

	ShowWindow(SW_MINIMIZE);

	// TODO: Add extra initialization here

	return TRUE;  // return TRUE  unless you set the focus to a control
}

void CsinsinDlg::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 CsinsinDlg::OnPaint()
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, reinterpret_cast<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 function to obtain the cursor to display while the user drags
//  the minimized window.
HCURSOR CsinsinDlg::OnQueryDragIcon()
{
	return static_cast<HCURSOR>(m_hIcon);
}


void CsinsinDlg::DrawSin()
{	
	Complex *A,U,W,T;
	const float PI=3.1415926;


	int M=10;
	int N=1024;
	int LE,LE1,I,J,IP; 

	A = new Complex[2000];
	//void ReverseOrder(int ,Complex *); 

	CDC *pDC;
	pDC=GetDC();

	CRect rect;
	GetClientRect(&rect);

	int i;

	UpdateData(TRUE);
	pDC->MoveTo(100,100);

	for(i=0;i<360;i++)
	{
		pDC->LineTo( 100+i,100 - (m_fd)*sin((double)((m_pl)*i*3.1415926/180.0)) );
		//pDC->LineTo(100+i,100+(m_fd)*sin( (double)((m_pl)*i*180/3.1415926)) + 
		//					(m_fd)*sin( (double)((m_pl*20)*i*180/3.1415926)) );
	}
	
	for(i=0;i<1024;i++)
	{ 
		A[i].Re=(m_fd)*sin( (double)((m_pl)*i*3.1415926/180.0));
		A[i].Im=0;
	} 
	ReverseOrder(N,A); 
	int L=1; 
	while(L<=M) 
	{ 
		LE=(int)pow(2.0,L); 
		LE1=LE/2; 
		U.Re=1.0f; 
		U.Im=0.0f; 
		//计算W算子的值 
		W.Re=(float)cos(PI/(1.0*LE1));  
		W.Im=(float)-1.0*sin(PI/(1.0*LE1)); 
		if(abs(W.Re)<1.0e-12) 
			W.Re=0.0f; 
		if(abs(W.Im)< 1.0e-12) 
			W.Im=0.0f; 
		J=1; 
		while(J<=LE1) 
		{ 
			I=J; 
			while(I<=N) 
			{ 
				IP=I+LE1; 
				//复数运算A×U 
				T.Re=(float)A[IP-1].Re*U.Re-A[IP-1].Im*U.Im;  

				T.Im=(float)A[IP-1].Re*U.Im+A[IP-1].Im*U.Re;  

				//复数运算A－T 
				A[IP-1].Re=(float)A[I-1].Re-T.Re;  

				A[IP-1].Im=(float)A[I-1].Im-T.Im;  

				//复数运算A＋T 
				A[I-1].Re+=T.Re; 
				A[I-1].Im+=T.Im; 
				I+=LE; 
			} 
			float temp=U.Re; 
			//复数运算U×W 
			U.Re=(float)U.Re*W.Re-U.Im*W.Im;  

			U.Im=(float)temp*W.Im+U.Im*W.Re; 

			J++; 
		} 
		L++; 
	} 
	
	pDC->MoveTo(100,300-100);
	pDC->LineTo(100,300+100);

	pDC->MoveTo(100,300);
	pDC->LineTo(100+1000,300);
	CPen Red;

	Red.CreatePen(PS_SOLID,1,RGB(255,0,0));;
	pDC->SelectObject(&Red);
	pDC->MoveTo(100,300);
	for(i=0;i<1000;i++)
	{
		pDC->LineTo(100+i,300-20*log(sqrt(A[i].Re * A[i].Re + A[i].Im * A[i].Im)));
	}	

	// TODO: Add your control notification handler code here

}
void CsinsinDlg::cancel()
{	
	OnOK();
	// TODO: Add your control notification handler code here
}



void CsinsinDlg::ReverseOrder(int N, Complex *A)
{

	int NV2=N/2; 
	int NM1=N-1; 
	int I,J,K=0; 
	//用于中介的复数变量T 
	Complex T; 
	I=J=1; 
	while(I<=NM1) 
	{ 
		if(I< J) 
		{ 
			//借助于中间变量T，将A[J－1]的内容和A[I－1]的内容互换 
			T=A[J-1]; 
			A[J-1]=A[I-1]; 
			A[I-1]=T; 
		} 
		K=NV2; 
		while(K< J) 
		{ 
			J-=K; 
			K/=2; 
		} 
		J+=K; 
		I++; 
	} 

}