hfilter.cpp

波形显示程序,高通滤波和低通滤波..........波形显示程序,高通滤波和低通滤波

// HFilter.cpp : implementation file
//

#include "stdafx.h"
#include "EEGProcessing.h"
#include "HFilter.h"
#include "math.h"
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// CHFilter dialog


CHFilter::CHFilter(CWnd* pParent /*=NULL*/)
	: CDialog(CHFilter::IDD, pParent)
{
	//{{AFX_DATA_INIT(CHFilter)
		// NOTE: the ClassWizard will add member initialization here
	//}}AFX_DATA_INIT
}


void CHFilter::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CHFilter)
		// NOTE: the ClassWizard will add DDX and DDV calls here
	//}}AFX_DATA_MAP
}


BEGIN_MESSAGE_MAP(CHFilter, CDialog)
	//{{AFX_MSG_MAP(CHFilter)
	ON_WM_PAINT()
	ON_WM_TIMER()
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CHFilter message handlers

BOOL CHFilter::OnInitDialog() 
{
	CDialog::OnInitDialog();
	
	// TODO: Add extra initialization here
	m_Low = 0;
	m_High = 1024;
	m_now =0;
	SetTimer(1,100,NULL);
	return TRUE;  // return TRUE unless you set the focus to a control
	              // EXCEPTION: OCX Property Pages should return FALSE
}

void CHFilter::OnPaint() 
{
	CPaintDC dc(this); // device context for painting
	
	// TODO: Add your message handler code here
	
	// Do not call CDialog::OnPaint() for painting messages
}
float CHFilter::hfilter(float *data ,long len)
{
	int i;
	float *data_t;
// 输入，输出临时缓冲区；
	float xv[5], yv[5];
//滤波器增益
	float GAIN=1.0162;
//输出采样值；
	double out;
//初始化输出缓冲区；
	for(i=0;i<5;i++)
		yv[i]=0;
		data_t=data;
//读取当前采样值
		xv[4] = *data_t / GAIN;
		xv[0] = xv[4]; xv[1] = xv[4];
		xv[2] = xv[4]; xv[3] = xv[4];
	for(i=0;i<len;i++)
		{
			data_t++;
			xv[0] = xv[1]; xv[1] = xv[2];
			xv[2] = xv[3]; xv[3] = xv[4];
			xv[4] = float(*data_t) / GAIN;
//更新状态；
			yv[0] = yv[1]; yv[1] = yv[2];
			yv[2] = yv[3]; yv[3] = yv[4];
//计算输出采样值；
			out =(xv[0] + xv[4]) - 4 * (xv[1] + xv[3]) + 6 * xv[2]+(-0.9684 * yv[0])+ ( 3.9048 * yv[1]) +( -5.9043 * yv[2]) + ( 3.9679 * yv[3]);
			yv[4] =(float)out;
			*data_t=yv[4];
		}
return *data_t;
}
void CHFilter::DrawWave1(CDC *pDC)
{
	CRect rect;
	
	CString str;
	int i;
	int m_left,m_top,m_right,m_bottom;
	
    int m_Interval = (m_High - m_Low)/10;
	if (m_Interval < 1)  m_Interval = 1;
	
	// 获取绘制坐标的文本框
	CWnd* pWnd = GetDlgItem(IDC_COORD1);
	
	pWnd->GetClientRect(&rect);
	pDC->Rectangle(&rect);
	
	
    m_left = rect.left+10;
	m_top = rect.top+10;
	m_right = rect.right-10;
	m_bottom = rect.bottom-20;
	
	int m_IntervalPan = (m_right - m_left)/11;
	if (m_IntervalPan < 1 ) m_IntervalPan =1;
	
	// 创建画笔对象
	CPen* pPenRed = new CPen;
	
	// 红色画笔
	//pPenRed->CreatePen(PS_SOLID,1,RGB(255,0,0));
	
	// 创建画笔对象
	CPen* pPenBlue = new CPen;
	
	// 蓝色画笔
	pPenBlue->CreatePen(PS_SOLID,1,RGB(0,0, 255));
	
	// 创建画笔对象
	CPen* pPenGreen = new CPen;
	
	// 绿色画笔
	pPenGreen->CreatePen(PS_DOT,1,RGB(0,255,0));
	
	// 选中当前红色画笔，并保存以前的画笔
	CGdiObject* pOldPen = pDC->SelectObject(pPenRed);
	
	// 绘制坐标轴
	pDC->MoveTo(m_left,m_top);
	
	// 垂直轴
	pDC->LineTo(m_left,m_bottom);
	
	// 水平轴
	pDC->LineTo(m_right,m_bottom);
	
	
	// 写X轴刻度值
	for(i=0;i<10;i++)
	{
		//str.Format(_T("%d"),m_Min+i*m_Interval);
		str.Format(_T("%d"),m_Low+i*m_Interval);
		pDC->TextOut(m_left+i*m_IntervalPan,m_bottom+3,str);		
	}
	//str.Format(_T("%d"),m_Max);
	str.Format(_T("%d"),m_High);
	pDC->TextOut(m_left+10*m_IntervalPan,m_bottom+3,str);		
	
	
	// 绘制X轴刻度
	for (i = m_left; i < m_right-20; i += 5)
	{
		if ((i & 1) == 0)
		{
			// 10的倍数
			pDC->MoveTo(i + 10, m_bottom);
			pDC->LineTo(i + 10, m_bottom+4);
		}
		else
		{
			// 10的倍数
			pDC->MoveTo(i + 10, m_bottom);
			pDC->LineTo(i + 10, m_bottom+2);
		}
	}
	
	// 绘制Y轴箭头
	pDC->MoveTo(m_right-5,m_bottom-5);
	pDC->LineTo(m_right,m_bottom);
	pDC->LineTo(m_right-5,m_bottom+5);
	
	// 绘制X轴箭头	
	pDC->MoveTo(m_left-5,m_top+5);
	pDC->LineTo(m_left,m_top);
	pDC->LineTo(m_left+5,m_top+5);
	

	
	// 绘制Y轴网格 选择绿色画笔
    pDC->SelectObject(pPenGreen);	
    
	int iTemp = (m_bottom  - m_top)/5;
    for (i = 1 ;i <= 5 ;i++)
	{
		pDC->MoveTo(m_left,m_bottom - i*iTemp);
		pDC->LineTo(m_right,m_bottom - i*iTemp);
	}

    //数组赋值
	if(m_now<1024)
	{
	    m_now++;
	}
	else
	{
		m_now = 0;
	}
	for(i = m_Low;i<m_High;i++)
	{ 
		//m_lCount[i]=(rand()%10+rand()%100+rand()%1000)%1024;
        m_lCount[i] = /*(int)*/(sin(i-m_now)*13); 
	//	/*m_lCount[i] =*/hfilter(m_lCount,1024);
	}
	int xTemp;
	int yTemp;
	// 更改成蓝色画笔
	pDC->SelectObject(pPenBlue);
	for (i = m_Low; i <= m_High; i=i+3)
	{		
		xTemp = m_left+(i-m_Low)*m_IntervalPan/m_Interval;
		yTemp = m_bottom/2 - /*(int) */(hfilter(m_lCount,1024) * m_bottom / 1024);
		if (yTemp < m_top) yTemp = m_top;
		if((xTemp >=m_left)&&(xTemp <=m_right))
		{
			pDC->MoveTo(xTemp, m_bottom/2);
			pDC->LineTo(xTemp, yTemp);
		}
	}
	// 恢复以前的画笔
	pDC->SelectObject(pOldPen);	
	delete pPenRed;
	delete pPenBlue;
	delete pPenGreen;
	return;
}
void CHFilter::OnTimer(UINT nIDEvent) 
{
	// TODO: Add your message handler code here and/or call default
	CRect rect;
	
	// 获取绘制坐标的文本框
	CWnd* pWnd = GetDlgItem(IDC_COORD1);                
	
	//获得对话框上的picture的窗口句柄
	pWnd->GetClientRect(&rect);
	// 指针
	pDC = pWnd->GetDC();
		pWnd->Invalidate();
	pWnd->UpdateWindow();
	
	//内存绘图
	CBitmap memBitmap;
	CBitmap* pOldBmp = NULL;
	
	//创建内存绘图设备
	memDC.CreateCompatibleDC(pDC);               
	memBitmap.CreateCompatibleBitmap(pDC,rect.right,rect.bottom);
	pOldBmp = memDC.SelectObject(&memBitmap);
	memDC.BitBlt(rect.left,rect.top,rect.right,rect.bottom,pDC,0,0,SRCCOPY);
	
	//自定义绘图函数，
	DrawWave1(&memDC);     
	
	//把内存绘图拷贝到屏幕
	pDC->BitBlt(rect.left,rect.top,rect.right,rect.bottom,&memDC,0,0,SRCCOPY); 
	
	memDC.SelectObject(pOldBmp);
	memDC.DeleteDC();
	memBitmap.DeleteObject();
	CDialog::OnTimer(nIDEvent);
}