crdlg.cpp

这是计算海上双船间碰撞危险度的程序

// CRDlg.cpp : implementation file
//

#include "stdafx.h"
#include "CR.h"
#include "CRDlg.h"
#include "Math.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif




	double DLA;//
    double k1,k2;
	double F1;
	double DE;
	double F2;
	double F3;
	double d1,d2;
	double Arena;	
	double t1,t2;
    double crdcpa,crtcpa;
	double aT;

/////////////////////////////////////////////////////////////////////////////
// CAboutDlg dialog used for App About

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

// Dialog Data
	//{{AFX_DATA(CAboutDlg)
	enum { IDD = IDD_ABOUTBOX };
	//}}AFX_DATA

	// ClassWizard generated virtual function overrides
	//{{AFX_VIRTUAL(CAboutDlg)
	protected:
	virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
	//}}AFX_VIRTUAL

// Implementation
protected:
	//{{AFX_MSG(CAboutDlg)
	//}}AFX_MSG
	DECLARE_MESSAGE_MAP()
};

CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
{
	//{{AFX_DATA_INIT(CAboutDlg)
	//}}AFX_DATA_INIT
}

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

BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
	//{{AFX_MSG_MAP(CAboutDlg)
		// No message handlers
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CCRDlg dialog

CCRDlg::CCRDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CCRDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CCRDlg)
	m_dcr = 0.0;
	m_dlong = 0.0;
	m_dvelocity1 = 0.0;
	m_dvelocity2 = 0.0;
	m_dcourse1 = 0.0;
	m_dcourse2 = 0.0;
	m_dXT = 0.0;
	m_dRT = 0.0;
	m_Chuiyu = _T("");
	m_Cnengjiandu = _T("");
	m_Cshuiyu = _T("");
	m_nT = 0;
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CCRDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CCRDlg)
	DDX_Text(pDX, IDC_CRZ, m_dcr);
	DDX_Text(pDX, IDC_SHIPLONG, m_dlong);
	DDV_MinMaxDouble(pDX, m_dlong, 0., 10000.);
	DDX_Text(pDX, IDC_VELOCITY1, m_dvelocity1);
	DDV_MinMaxDouble(pDX, m_dvelocity1, 0., 60.);
	DDX_Text(pDX, IDC_VELOCITY2, m_dvelocity2);
	DDV_MinMaxDouble(pDX, m_dvelocity2, 0., 60.);
	DDX_Text(pDX, IDC_COURSE1, m_dcourse1);
	DDV_MinMaxDouble(pDX, m_dcourse1, 0., 360.);
	DDX_Text(pDX, IDC_COURSE2, m_dcourse2);
	DDV_MinMaxDouble(pDX, m_dcourse2, 0., 360.);
	DDX_Text(pDX, IDC_XORIENTATION, m_dXT);
	DDV_MinMaxDouble(pDX, m_dXT, -180., 180.);
	DDX_Text(pDX, IDC_DISTANCE, m_dRT);
	DDX_LBString(pDX, IDC_HUIYU, m_Chuiyu);
	DDX_CBString(pDX, IDC_NENGJIANDU, m_Cnengjiandu);
	DDX_CBString(pDX, IDC_SHUIYU, m_Cshuiyu);
	DDX_Text(pDX, IDC_SHIJIAN, m_nT);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CCRDlg, CDialog)
	//{{AFX_MSG_MAP(CCRDlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_TDCALCULATE, OnTdcalculate)
	ON_BN_CLICKED(IDC_CRCALCULATE, OnCrcalculate)
	ON_LBN_SELCHANGE(IDC_HUIYU, OnSelchangeHuiyu)
	ON_CBN_SELCHANGE(IDC_SHUIYU, OnSelchangeShuiyu)
	ON_CBN_SELCHANGE(IDC_NENGJIANDU, OnSelchangeNengjiandu)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CCRDlg message handlers

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

	// Add "About..." menu item to system menu.
	CListBox* pLB=(CListBox*)GetDlgItem(IDC_HUIYU);
	pLB->InsertString(0,"对遇");
	pLB->InsertString(1,"追越");
	pLB->InsertString(2,"右舷交叉相遇");
    pLB->InsertString(3,"被追越");
    pLB->InsertString(4,"左舷交叉相遇");

    CComboBox* pLC = (CComboBox*)GetDlgItem(IDC_SHUIYU);
	pLC->InsertString(0,"宽深水域");
	pLC->InsertString(1,"浅窄水域");

    CComboBox* pLD = (CComboBox*)GetDlgItem(IDC_NENGJIANDU);
	pLD->InsertString(0,"良好");
	pLD->InsertString(1,"不良");





	// 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
	
	// TODO: Add extra initialization here
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

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

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

void CCRDlg::OnTdcalculate() 
{
	// TODO: Add your control notification handler code here
	UpdateData(TRUE);
	double vx0;
	double vy0;
	double vxt;
	double vyt;
	double vxr,vyr,vxr2,vyr2,vr2;
	int a;
	double Cvr;
	double change=2*3.1415926/360;
    char szBuf[100];


	vx0=m_dvelocity1*sin(m_dcourse1*change);
	vy0=m_dvelocity1*cos(m_dcourse1*change);//本船在x,y上的速度分量

	vxt=m_dvelocity2*sin(m_dcourse2*change);
	vyt=m_dvelocity2*cos(m_dcourse2*change);//他船在x,y上的速度分量

	vxr=vxt-vx0;
	vyr=vyt-vy0;//相对速度分量
     
	vxr2=pow(vxr,2);
	vyr2=pow(vyr,2);
	vr2=vxr2+vyr2;
	vr=pow(vr2,0.5);//vr=(vxr*vxr+vyr*vyr)1/2,计算相对速度大小

	if(vxr>=0)
	{
		if( vyr>=0)
		    a=0;
		else 
			a=180;
	}
	else 
	{
		if(vyr>=0)
			a=360;
		else
			a=180;
	}

   
	Cvr=atan(vxr/vyr)/change+a;//计算相对航向


	if(m_dXT>=0)
	{
		if(m_dXT+m_dcourse1>=360)
			aT=m_dXT+m_dcourse1-360;
		else
			aT=m_dXT+m_dcourse1;
	}
	else
	{
		if(m_dXT+m_dcourse1<=0)
			aT=m_dXT+m_dcourse1+360;
		else
            aT=m_dXT+m_dcourse1;
	}


	//aT=m_dXT+m_dcourse1;//计算真方位

	m_ddcpa=m_dRT*sin((Cvr-aT-180)*change);
	m_dtcpa=60*m_dRT*cos((Cvr-aT-180)*change)/vr;//

//----------------------------将数值变为字符型显示------------------------------------------
	sprintf(szBuf, "%03.1f", m_ddcpa);
	SetDlgItemText(IDC_DCPA, szBuf);

	sprintf(szBuf, "%03.1f", m_dtcpa);
	SetDlgItemText(IDC_TCPA, szBuf);
//-------------------------------------------------------------------------------------
	

	UpdateData(FALSE);

}



void CCRDlg::OnSelchangeShuiyu() 
{
	// TODO: Add your control notification handler code here
    CComboBox* pLB = (CComboBox*)GetDlgItem(IDC_SHUIYU);
	int Num = pLB->GetCurSel();

    
	if(Num==0||Num==-1)
		F3=1.0;
	else
		F3=0.75;//有待考虑
	
}

void CCRDlg::OnSelchangeNengjiandu() 
{
	// TODO: Add your control notification handler code here
CComboBox* pLB = (CComboBox*)GetDlgItem(IDC_NENGJIANDU);
	int Num = pLB->GetCurSel();

    if(Num==0||Num==-1)
		F2=1.0;
	else
		F2=1.25;
	
}

void CCRDlg::OnSelchangeHuiyu() 
{
	// TODO: Add your control notification handler code here
	CListBox* pLB=(CListBox*)GetDlgItem(IDC_HUIYU);
    int Num = pLB->GetCurSel();

	if(Num==0|| Num==1||Num==2)
		k1=1.0;
	else
		k1=0.8;
	
}
void CCRDlg::OnCrcalculate() 
{
	// TODO: Add your control notification handler code here
    UpdateData(TRUE);


	DLA=12*m_dlong/1852;//

	if(m_ddcpa<=0.5)
		DE=aT;//没有完整
	

	if(DE<=180)
		F1=(360-DE)/360;
	else
		F1=DE/360;

	

	d1=(k1*DLA+m_nT/3600*vr*F1)*F2*F3;

    k2=2;//可以调整
	d2=k2*d1;
//---------------------------------------------------------------------------------------//
	if(m_ddcpa<=d1)
		crdcpa=1;
	else
	{
		if(m_ddcpa<=d2 && m_ddcpa>d1)
			crdcpa=0.5-0.5*sin(3.14/(d2-d1)*(m_ddcpa-(d1+d2)/2));
		else
			crdcpa=0;
	}
//--------------------------------------------------------------------------------------//
    Arena=d2+vr*m_nT/3600;

	if(m_ddcpa<=DLA)
		t1=pow((pow(DLA,2)-pow(m_ddcpa,2)),0.5)/vr*60;
	else
		t1=(DLA-m_ddcpa)/vr*60;


	//t2=pow((pow(Arena,2)-pow(m_ddcpa,2)),0.5)/vr*60;
    t2=pow((pow(12,2)-pow(m_ddcpa,2)),0.5)/vr*60;

	if(m_dtcpa<=t1)
		crtcpa=1;
	else
	{
		if(m_dtcpa<=t2 && m_dtcpa>t1)
			crtcpa=(t2-m_dtcpa)/(t2-t1);
		else
			crtcpa=0;
	}
//---------------------------------------------------------------------------------------//
	
	if(crdcpa==0)
		m_dcr=0;
	else
	{
		if(crtcpa==0)
			m_dcr=0;
		else
            m_dcr=0.5*crdcpa+0.5*crtcpa;//有待考虑
	}


    UpdateData(FALSE);

	
}