firdlg.cpp

数字信号处理算法离散卷积FIR

// FIRDlg.cpp : 实现文件
//

#include "stdafx.h"
#include "FIR.h"
#include "FIRDlg.h"
#include ".\firdlg.h"
#include <cmath>
#include "ListOper.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#endif

float pi = (float)3.1415926;
// CFIRDlg 对话框



CFIRDlg::CFIRDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CFIRDlg::IDD, pParent)
	, m_N(100)
	, m_Wc(0.2)
{
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CFIRDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDTN, m_N);
	DDV_MinMaxLong(pDX, m_N, 10, 100000);
	DDX_Text(pDX, IDC_EDTF, m_Wc);
	DDV_MinMaxDouble(pDX, m_Wc, 0, 1);
}

BEGIN_MESSAGE_MAP(CFIRDlg, CDialog)
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	//}}AFX_MSG_MAP
	ON_BN_CLICKED(IDC_BTNCAL, OnBnClickedBtncal)
	ON_BN_CLICKED(IDC_BTNCLOSE, OnBnClickedBtnclose)
END_MESSAGE_MAP()


// CFIRDlg 消息处理程序

BOOL CFIRDlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	// 设置此对话框的图标。当应用程序主窗口不是对话框时，框架将自动
	//  执行此操作
	SetIcon(m_hIcon, TRUE);			// 设置大图标
	SetIcon(m_hIcon, FALSE);		// 设置小图标

	// TODO: 在此添加额外的初始化代码
	this->m_N = 100;
	this->m_Wc = 0.025;
	this->UpdateData(FALSE);
	return TRUE;  // 除非设置了控件的焦点，否则返回 TRUE
}

// 如果向对话框添加最小化按钮，则需要下面的代码
//  来绘制该图标。对于使用文档/视图模型的 MFC 应用程序，
//  这将由框架自动完成。

void CFIRDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // 用于绘制的设备上下文

		SendMessage(WM_ICONERASEBKGND, reinterpret_cast<WPARAM>(dc.GetSafeHdc()), 0);

		// 使图标在工作矩形中居中
		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;

		// 绘制图标
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

//当用户拖动最小化窗口时系统调用此函数取得光标显示。
HCURSOR CFIRDlg::OnQueryDragIcon()
{
	return static_cast<HCURSOR>(m_hIcon);
}

/****************************************************
*   Name    : FuncHd
*	Function: Hd()--Required frequency response function
*
*
*****************************************************/
COMP FuncHd(double LowLimit,double UpperLimit,COMP x)
{
	if(x.real()>UpperLimit||x.real()<LowLimit)
		return 0;
	else
		return 1;
	
}
void Integral(double LowLimit,double UpperLimit,int N,COMP Hn[])
{
	int M = 2*N;
	for(int i=0;i<N;i++)
	{
		Hn[i] = COMP(0,0);
		for(int k=0;k<M;k++)
		{
			COMP C = COMP(cos(2*pi*i*k/(double)M),sin(2*pi*i*k/(double)M));
			Hn[i] += C*FuncHd(LowLimit,UpperLimit,COMP(cos(2*pi*k/(double)M),sin(2*pi*k/(double)M)));

		}
		Hn[i] = Hn[i]*COMP(1/(double)M,0);
	}

}

void CFIRDlg::OnBnClickedBtncal()
{
	// TODO: 在此添加控件通知处理程序代码
	int r = 0;
	r=(0?1:0);
	this->UpdateData(TRUE);
	int nN = this->m_N;
	double Limit = this->m_Wc*pi; 
	COMP* Hn = new COMP[nN];
	Integral(-Limit,Limit,nN,Hn);

	CListCtrl* pList=(CListCtrl*) this->GetDlgItem(IDC_LISTRET);
	CListOper oper;
	oper.FillList(*pList,nN,Hn);

	delete[] Hn;
}

void CFIRDlg::OnBnClickedBtnclose()
{
	// TODO: 在此添加控件通知处理程序代码
	this->DestroyWindow();
}

void CFIRDlg::OnOK()
{
	// TODO: 在此添加专用代码和/或调用基类

	//CDialog::OnOK();
}