blueghost5dlg.cpp

信道基本特性的仿真程序

// blueghost5Dlg.cpp : implementation file
//

#include "stdafx.h"
#include "blueghost5.h"
#include "blueghost5Dlg.h"
#include "math.h"

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

/////////////////////////////////////////////////////////////////////////////
// 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()

/////////////////////////////////////////////////////////////////////////////
// CBlueghost5Dlg dialog

CBlueghost5Dlg::CBlueghost5Dlg(CWnd* pParent /*=NULL*/)
	: CDialog(CBlueghost5Dlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CBlueghost5Dlg)
		// NOTE: the ClassWizard will add member initialization here
	//}}AFX_DATA_INIT
	// Note that LoadIcon does not require a subsequent DestroyIcon in Win32
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

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

BEGIN_MESSAGE_MAP(CBlueghost5Dlg, CDialog)
	//{{AFX_MSG_MAP(CBlueghost5Dlg)
	ON_WM_SYSCOMMAND()
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(start, Onstart)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CBlueghost5Dlg message handlers

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

	// Add "About..." menu item to system menu.

	// 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 CBlueghost5Dlg::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 CBlueghost5Dlg::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 CBlueghost5Dlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}
int rand_ra();
void CBlueghost5Dlg::Onstart() 
{
	// TODO: Add your control notification handler code here
	int loop1,loop2;
	double FDT, NORM1, NORM2;
	double fading_I, fading_Q;
	double fading1_I, fading1_Q;
	double fading2_I, fading2_Q;
	double ALPHA1[17], PHI1[17], THETA1[17];
	double ALPHA2, PHI2, THETA2;
	double T1_I[17], T1_Q[17];

	double rnd1,rnd2,rnd3,rnd4;
	double PI=3.1415926;
	double K;
	double KK[3]={1,10,1000};
	double fading_abs,fading_log;  
	int N,k;
	srand( (unsigned)time( NULL ) );	
	CClientDC dc(this);

// 竖线						流程图(1)的部分
	for (loop1 = 0; loop1 <=10; loop1=loop1+1)
		{dc.MoveTo (50+50*loop1,50);	dc.LineTo (50+50*loop1,350);}
// 横线
	for (loop1 = -10; loop1 <=50; loop1=loop1+10)
		{dc.MoveTo (50,100+loop1*5);	dc.LineTo (550,100+loop1*5);}
// 初始設定
	for(int i=0;i<3;i++)
	{
	K=KK[i];	

/* 为使散射波与直射波的功率总和为１的归一化设定*/
	/* NORM1:散射波归一化变量 */
	/* NORM2:直接波归一化变量*/
   	NORM1=sqrt(1./16.)*sqrt(1./(K+1));	//详细的请参照后面的实验说明
	NORM2=sqrt(K/(K+1));  		

/* 初始设定衰落的角度?相位 */		
	FDT= 0.0030864;
	rnd1=rand_ra();
	THETA1[1]=(double)(rnd1/30000*2*PI/16);
	for(N=2;N<=16;N++)
	{
		THETA1[N]=2*PI/16 +THETA1[N-1];
	}
	for(N=1;N<=16;N++)
	{
		rnd2=(double)rand_ra()/30000;
		PHI1[N]=(double)(2*PI*rnd2);
		ALPHA1[N]=2*PI*FDT*cos(THETA1[N]);
	}
	THETA2=(double)rand_ra()/30000*2*PI;
	PHI2=(double)rand_ra()/30000*2*PI;
	ALPHA2=2*PI*FDT*cos(THETA2);
for (k=0; k<=1000; k++)  // 以下进入循环．
{
/* 衰落仿真器 */			//流程图(3)的部分

/* 初始設定 */
	fading1_I = 0;		fading1_Q = 0;
	fading2_I = 0;		fading2_Q = 0;
//↓↓↓从这开始追加 
	for(N=1;N<=16;N++)
		{
			T1_I[N]=cos(ALPHA1[N]*k+PHI1[N]);
			T1_Q[N]=sin(ALPHA1[N]*k+PHI1[N]);
			fading1_I=T1_I[N]+fading1_I;
			fading1_Q=T1_Q[N]+fading1_Q;
		}
		fading1_I=fading1_I*NORM1;
		fading1_Q=fading1_Q*NORM1;
		fading2_I=NORM2*cos(ALPHA2*k+PHI2);
		fading2_Q=NORM2*sin(ALPHA2*k+PHI2);
		//fading1_Q=0;fading1_I=0;
		fading_I=fading1_I+fading2_I;
		fading_Q=fading1_Q+fading2_Q;
 
//↑↑↑到此为止

/* 衰落振幅的計算 */			//流程图(4)的部分
	fading_abs = sqrt(pow(fading_I,2)+pow(fading_Q,2));
	fading_log = 20*log10(fading_abs);

/* 做图部分 */
	for(loop2 =0; loop2 <=500000; loop2=loop2+1) {};
	dc.SetPixel (k*0.5+50,100-fading_log*5,RGB(255,0,0));
}
	UpdateData(FALSE);
	}
}
int rand_ra()	
{
	int r,i=0;	//定义发生1～30000的随机計算中如使用随机数，就调用这函数．
	while(i == 0){
		r = rand() - 2767;
		if (r >= 0) break; 
	}
	return r;
}