patternsearchview.cpp

非线性规划中求解无约束问题中的一个有效方法--模矢法(步长加速法),并用图表示出来步骤.

// PatternSearchView.cpp : CPatternSearchView 类的实现
//

#include "stdafx.h"
#include "PatternSearch.h"
#include "math.h"

#include "PatternSearchDoc.h"
#include "PatternSearchView.h"
#include ".\patternsearchview.h"

#define LX	15
#define LY	15
#define RX	335
#define RY	335
#define SPAN	32
#define MAX	2
#define JINGDU	0.1

#ifdef _DEBUG
#define new DEBUG_NEW
#endif

float JiDian[MAX];	//存储基点数据
float OldJiDian[MAX];
float ShiDian[MAX];
float OldShiDian[MAX];
float tempMatrix[MAX];
float tempJiDian[MAX];
float tempShiDian[MAX];
float BuChang;		//步长
BOOL JustOnce;
BOOL IsShiDian;
BOOL GetIt;


// CPatternSearchView

IMPLEMENT_DYNCREATE(CPatternSearchView, CFormView)

BEGIN_MESSAGE_MAP(CPatternSearchView, CFormView)
	ON_BN_CLICKED(IDC_BUTTON_CALCULATE, OnBnClickedButtonCalculate)
END_MESSAGE_MAP()

// CPatternSearchView 构造/析构

CPatternSearchView::CPatternSearchView()
	: CFormView(CPatternSearchView::IDD)
	, m_cResult(_T(""))
	, m_bStep(FALSE)
	, m_cJiDian(_T(""))
	, m_fBuChang(0)
{
	// TODO: 在此处添加构造代码
	for(int i=0;i<MAX;i++)
	{
		JiDian[i] = 0;
		OldJiDian[i] = 0;
	}
	
	m_fBuChang = 0.5;
	m_cJiDian = "-3,1";
	m_cResult.Format("\r\n步长加速法\r\n"
		"--------------------------\r\n"
		"模矢法"
		"(Pattern Search)\r\n"
		"Hooke,Jeeves1961提出");
	m_bStep = FALSE;
	JustOnce = TRUE;
	IsShiDian = FALSE;
	GetIt = TRUE;
}

CPatternSearchView::~CPatternSearchView()
{
}

void CPatternSearchView::DoDataExchange(CDataExchange* pDX)
{
	CFormView::DoDataExchange(pDX);
	DDX_Text(pDX, IDC_EDIT_RESULT, m_cResult);
	DDX_Check(pDX, IDC_CHECK_ZHUBU, m_bStep);
	DDX_Text(pDX, IDC_EDIT_JIDIAN, m_cJiDian);
	DDX_Text(pDX, IDC_EDIT_BUCHANG, m_fBuChang);
}

BOOL CPatternSearchView::PreCreateWindow(CREATESTRUCT& cs)
{
	// TODO: 在此处通过修改 CREATESTRUCT cs 来修改窗口类或
	// 样式
	
	return CFormView::PreCreateWindow(cs);
}

void CPatternSearchView::OnInitialUpdate()
{
	CFormView::OnInitialUpdate();
	GetParentFrame()->RecalcLayout();
	ResizeParentToFit();

}

void CPatternSearchView::OnDraw(CDC *pDC)
{


	CPen Pen;
	Pen.CreatePen(PS_SOLID,1,RGB(0,0,0));

	CPen dPen;
	dPen.CreatePen(PS_DOT,1,RGB(255,255,255));

	CBrush Brush;
	Brush.CreateSolidBrush(RGB(255,255,255));

	CBrush *pOldBrush;
	pOldBrush = pDC->SelectObject( &Brush );
	pDC->SetBkColor(RGB(0,0,0));
	pDC->Rectangle(LX,LY,RX,RY);
	pDC->SelectObject( pOldBrush ); 
	Brush.DeleteObject();

	
	CPen *pOldPen;
	//绘辅助线
	pOldPen = pDC->SelectObject( &dPen );

	for(int i=1;i<(RY-LY)/SPAN;i++)
	{
		pDC->MoveTo( LX,LY+i*SPAN);
		pDC->LineTo( RX,LY+i*SPAN);
		pDC->MoveTo( LX+i*SPAN,LY);
		pDC->LineTo( LX+i*SPAN,RY);
	}

	pDC->SelectObject( pOldPen );
	dPen.DeleteObject();

	//绘坐标轴
	pOldPen = pDC->SelectObject( &Pen);
	pDC->MoveTo( LX,LY+(RY-LY)/2 );
	pDC->LineTo( RX,LY+(RY-LY)/2 );
	pDC->MoveTo( LX+(RX-LX)/2,LY );
	pDC->LineTo( LX+(RX-LX)/2,RY );
	pDC->SelectObject( pOldPen );
	Pen.DeleteObject();
}

// CPatternSearchView 诊断

#ifdef _DEBUG
void CPatternSearchView::AssertValid() const
{
	CFormView::AssertValid();
}

void CPatternSearchView::Dump(CDumpContext& dc) const
{
	CFormView::Dump(dc);
}

CPatternSearchDoc* CPatternSearchView::GetDocument() const // 非调试版本是内联的
{
	ASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CPatternSearchDoc)));
	return (CPatternSearchDoc*)m_pDocument;
}
#endif //_DEBUG


// CPatternSearchView 消息处理程序

void CPatternSearchView::OnBnClickedButtonCalculate()
{
	// TODO: 在此添加控件通知处理程序代码
	UpdateData();
	if(JustOnce)	//第一次执行才运行，初始化数据
	{
		InitializeParameter();
		CopyMatrix(ShiDian,JiDian);
		m_cResult = "";
	}
	JustOnce = FALSE;
	if(m_bStep)
	{
		GetDlgItem(IDC_BUTTON_CALCULATE)->SetWindowText("继 续");
	}
	else
	{
		GetDlgItem(IDC_BUTTON_CALCULATE)->SetWindowText("计 算");
	}

	//---------------------------------------------------------

	//step one
	CopyMatrix(OldJiDian,JiDian);
	if(!IsShiDian)
	{
		CopyMatrix(ShiDian,JiDian);
	}
	IsShiDian = FALSE;

	//step two
	for(int i=0;i<MAX;i++)
		SheDong(ShiDian,ShiDian,BuChang,i);
	CopyMatrix(OldShiDian,ShiDian);

	//step three
	for(int i=0;i<MAX;i++)
		SheDong(ShiDian,ShiDian,BuChang,i);
	CopyMatrix(JiDian,ShiDian);

	//step four
	CopyMatrix(ShiDian,MatrixOperationTwo( OldJiDian,JiDian ));

	if(TargetFunValue(ShiDian)<TargetFunValue(JiDian))
	{
		CopyMatrix(tempShiDian,ShiDian);
		for(int i=0;i<MAX;i++)
			SheDong(tempShiDian,tempShiDian,BuChang,i);
		if(!GetIt)
		{
			CopyMatrix(JiDian,ShiDian);
		}
		else
		{
			IsShiDian = TRUE;
		}
		Recursion();
		if(BuChang<JINGDU)
			return;
		if(!m_bStep)
			OnBnClickedButtonCalculate();	//结束one
	}
	else
	{
		CopyMatrix(tempShiDian,ShiDian);
		CopyMatrix(tempJiDian,JiDian);
		GetIt = FALSE;
		for(int i=0;i<MAX;i++)
			SheDong(tempShiDian,tempJiDian,BuChang,i);
		if(!GetIt)
		{
			CopyMatrix(tempJiDian,JiDian);
			for(int i=0;i<MAX;i++)
				SheDong(tempJiDian,tempJiDian,BuChang,i);
			if(GetIt)
			{
				Recursion();	
				if(BuChang<JINGDU)
					return;
				if(!m_bStep)
					OnBnClickedButtonCalculate();	//结束two
			}
			else
			{
				BuChang = BuChang/2;
				Recursion();
				if(BuChang<JINGDU)
					return;
				if(!m_bStep)
					OnBnClickedButtonCalculate();//结束three
			}
		}
		if(BuChang<JINGDU)
			return;
	}

	//---------------------------------------------------------
	
}

// 目标函数求值
float CPatternSearchView::TargetFunValue(float fVariable[])
{
	float result=0;
	result = fVariable[0]*fVariable[0]+2*fVariable[1]*fVariable[1]-4*fVariable[0]-2*fVariable[0]*fVariable[1];
	return result;
}

// 初始化手工输入的参数
void CPatternSearchView::InitializeParameter()
{

	LPTSTR str = new TCHAR[m_cJiDian.GetLength()+1];
	_tcscpy(str, m_cJiDian);
	int i=0;
	char * pch;
	pch = strtok (str,",");
	while (pch != NULL)
	{
		JiDian[i]=atof(pch);
		i++;
		pch = strtok (NULL, ",");
	}
	BuChang = m_fBuChang;
}

// 对独立变量进行摄动
void CPatternSearchView::SheDong(float MatrixA[],float MatrixB[],float BuChang,int i)
{
	if( TargetFunValue(MatrixOperation( MatrixA,BuChang,i,1))<TargetFunValue(MatrixB))
	{
		CopyMatrix(MatrixA,MatrixOperation( MatrixA,BuChang,i,1));
		GetIt = TRUE;
	}
	else if( TargetFunValue(MatrixOperation( MatrixA,BuChang,i,-1))<TargetFunValue(MatrixB) &&  TargetFunValue(MatrixOperation( MatrixA,BuChang,i,1))>=TargetFunValue(MatrixB))
	{	
		CopyMatrix(MatrixA,MatrixOperation( MatrixA,BuChang,i,-1));
		GetIt = TRUE;
	}
	else
	{
		GetIt = FALSE;
	}

}

// 矩阵运算（可用运算符重载实现）
float * CPatternSearchView::MatrixOperation(float ShiDian[],float BuChang,int i,int Symbol)
{
	CopyMatrix(tempMatrix,ShiDian);
	tempMatrix[i] = tempMatrix[i]+Symbol*BuChang;
	return tempMatrix;
}

// 6.86
float * CPatternSearchView::MatrixOperationTwo(float MatrixA[],float MatrixB[])
{
	for(int i=0;i<MAX;i++)
		tempMatrix[i] = 2*MatrixB[i]-MatrixA[i];
	return tempMatrix;
}

float CPatternSearchView::CopyMatrix(float MatrixA[],float MatrixB[])
{
	for(int i=0;i<MAX;i++)
		MatrixA[i] = MatrixB[i];
	return 0;
}

void  CPatternSearchView::DrawPic()
{
	DrawPoint(OldJiDian);
	DrawPoint(JiDian);
	DrawPoint(OldShiDian);
	DrawPoint(ShiDian);
	DrawLine(OldJiDian,JiDian);
}

void CPatternSearchView::DrawPoint(float Matrix[])
{
	float Point[2] = { LX+(RX-LX)/2+Matrix[0]*SPAN,LY+(RY-LY)/2-Matrix[1]*SPAN };
	float PointR = 2;
	
	CClientDC dc(this);
	CRect rClient;
	GetClientRect(rClient);
	
	CBrush Brush;
	Brush.CreateSolidBrush(RGB(0,0,0));

	CBrush *pOldBrush;
	pOldBrush = dc.SelectObject( &Brush );
	dc.SetBkColor(RGB(0,0,0));
	dc.Ellipse( Point[0]-PointR,Point[1]-PointR,Point[0]+PointR,Point[1]+PointR );

	dc.SelectObject( pOldBrush ); 
	Brush.DeleteObject();
}

void CPatternSearchView::DrawLine(float MatrixA[],float MatrixB[])
{
	float PointOne[2] = { LX+(RX-LX)/2+MatrixA[0]*SPAN,LY+(RY-LY)/2-MatrixA[1]*SPAN };
	float PointTwo[2] = { LX+(RX-LX)/2+MatrixB[0]*SPAN,LY+(RY-LY)/2-MatrixB[1]*SPAN };

	CClientDC dc(this);
	CRect rClient;
	GetClientRect(rClient);

	CPen Pen;
	Pen.CreatePen(PS_SOLID,1,RGB(255,0,0));
	CPen *pOldPen;
	pOldPen = dc.SelectObject( &Pen );
	dc.MoveTo( PointOne[0],PointOne[1] );
	dc.LineTo( PointTwo[0],PointTwo[1] );

	dc.SelectObject( pOldPen );
	Pen.DeleteObject();	
}


void CPatternSearchView::DisplayZB(void)
{
	CString tempResult;
	tempResult.Format("(%f,%f)\r\nf(X)=%f\r\n",JiDian[0],JiDian[1],TargetFunValue(JiDian));
	m_cResult = tempResult + "\r\n" + m_cResult;

	//状态栏坐标显示
	CPatternSearchDoc * pDoc = GetDocument();
	ASSERT_VALID(pDoc);

	CString strPointX,strPointY;
	CStatusBar * pStatus = (CStatusBar *)AfxGetApp()->m_pMainWnd->GetDescendantWindow(ID_VIEW_STATUS_BAR);
	strPointX.Format("x:%f",JiDian[0]);
	strPointY.Format("y:%f",JiDian[1]);
	pStatus->SetPaneText(1,strPointX);
	pStatus->SetPaneText(2,strPointY);
}

void CPatternSearchView::Recursion(void)
{
	DisplayZB();
	DrawPic();
	UpdateData(FALSE);
}