?? 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);
}
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -