开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 行业应用 > 能源行业(电力石油煤炭) > 查看源码
ewaveDlg.cpp
资源名称：ewave.rar　[点击查看]
上传用户：mayisishi
上传日期：2007-08-08
资源大小：80k
文件大小：13k
源码类别：

能源行业(电力石油煤炭)

开发平台：

Visual C++

ewaveDlg.cpp：源码内容
  1. // ewaveDlg.cpp : implementation file
  2. //
  4. #include "stdafx.h"
  5. #include "ewave.h"
  6. #include "ewaveDlg.h"
  7. #include "math.h"
  8. #include "fstream.h"
  10. double adAg[2000],adt[2000],amax,tmax,dOpx=65,dOpy=220,adAgtemp[2000];
  13. #ifdef _DEBUG
  14. #define new DEBUG_NEW
  15. #undef THIS_FILE
  16. static char THIS_FILE[] = __FILE__;
  17. #endif
  19. /////////////////////////////////////////////////////////////////////////////
  20. // CAboutDlg dialog used for App About
  22. class CAboutDlg : public CDialog
  23. {
  24. public:
  25. CAboutDlg();
  27. // Dialog Data
  28. //{{AFX_DATA(CAboutDlg)
  29. enum { IDD = IDD_ABOUTBOX };
  30. //}}AFX_DATA
  32. // ClassWizard generated virtual function overrides
  33. //{{AFX_VIRTUAL(CAboutDlg)
  34. protected:
  35. virtual void DoDataExchange(CDataExchange* pDX);    // DDX/DDV support
  36. //}}AFX_VIRTUAL
  38. // Implementation
  39. protected:
  40. //{{AFX_MSG(CAboutDlg)
  41. //}}AFX_MSG
  42. DECLARE_MESSAGE_MAP()
  43. };
  45. CAboutDlg::CAboutDlg() : CDialog(CAboutDlg::IDD)
  46. {
  47. //{{AFX_DATA_INIT(CAboutDlg)
  48. //}}AFX_DATA_INIT
  49. }
  51. void CAboutDlg::DoDataExchange(CDataExchange* pDX)
  52. {
  53. CDialog::DoDataExchange(pDX);
  54. //{{AFX_DATA_MAP(CAboutDlg)
  55. //}}AFX_DATA_MAP
  56. }
  58. BEGIN_MESSAGE_MAP(CAboutDlg, CDialog)
  59. //{{AFX_MSG_MAP(CAboutDlg)
  60. // No message handlers
  61. //}}AFX_MSG_MAP
  62. END_MESSAGE_MAP()
  64. /////////////////////////////////////////////////////////////////////////////
  65. // CEwaveDlg dialog
  67. CEwaveDlg::CEwaveDlg(CWnd* pParent /*=NULL*/)
  68. : CDialog(CEwaveDlg::IDD, pParent)
  69. {
  70. //{{AFX_DATA_INIT(CEwaveDlg)
  71. m_dAmax = 0.08;
  72. m_dTg = 0.65;
  73. m_pd = 35;
  74. //}}AFX_DATA_INIT
  75. // Note that LoadIcon does not require a subsequent DestroyIcon in Win32
  76. m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
  77. }
  79. void CEwaveDlg::DoDataExchange(CDataExchange* pDX)
  80. {
  81. CDialog::DoDataExchange(pDX);
  82. //{{AFX_DATA_MAP(CEwaveDlg)
  83. DDX_Text(pDX, IDC_AMAX, m_dAmax);
  84. DDX_Text(pDX, IDC_TG, m_dTg);
  85. DDX_Text(pDX, IDC_PD, m_pd);
  86. //}}AFX_DATA_MAP
  87. }
  89. BEGIN_MESSAGE_MAP(CEwaveDlg, CDialog)
  90. //{{AFX_MSG_MAP(CEwaveDlg)
  91. ON_WM_SYSCOMMAND()
  92. ON_WM_PAINT()
  93. ON_WM_QUERYDRAGICON()
  94. ON_BN_CLICKED(IDOK, OnSTART)
  95. ON_BN_CLICKED(IDC_OUTPUT, OnOutput)
  96. ON_BN_CLICKED(IDC_INPUT, OnInput)
  97. //}}AFX_MSG_MAP
  98. END_MESSAGE_MAP()
  100. /////////////////////////////////////////////////////////////////////////////
  101. // CEwaveDlg message handlers
  103. BOOL CEwaveDlg::OnInitDialog()
  104. {
  105. CDialog::OnInitDialog();
  107. // Add "About..." menu item to system menu.
  109. // IDM_ABOUTBOX must be in the system command range.
  110. ASSERT((IDM_ABOUTBOX & 0xFFF0) == IDM_ABOUTBOX);
  111. ASSERT(IDM_ABOUTBOX < 0xF000);
  113. CMenu* pSysMenu = GetSystemMenu(FALSE);
  114. if (pSysMenu != NULL)
  115. {
  116. CString strAboutMenu;
  117. strAboutMenu.LoadString(IDS_ABOUTBOX);
  118. if (!strAboutMenu.IsEmpty())
  119. {
  120. pSysMenu->AppendMenu(MF_SEPARATOR);
  121. pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);
  122. }
  123. }
  125. // Set the icon for this dialog.  The framework does this automatically
  126. //  when the application's main window is not a dialog
  127. SetIcon(m_hIcon, TRUE); // Set big icon
  128. SetIcon(m_hIcon, FALSE); // Set small icon
  130. // TODO: Add extra initialization here
  132. return TRUE;  // return TRUE  unless you set the focus to a control
  133. }
  135. void CEwaveDlg::OnSysCommand(UINT nID, LPARAM lParam)
  136. {
  137. if ((nID & 0xFFF0) == IDM_ABOUTBOX)
  138. {
  139. CAboutDlg dlgAbout;
  140. dlgAbout.DoModal();
  141. }
  142. else
  143. {
  144. CDialog::OnSysCommand(nID, lParam);
  145. }
  146. }
  148. // If you add a minimize button to your dialog, you will need the code below
  149. //  to draw the icon.  For MFC applications using the document/view model,
  150. //  this is automatically done for you by the framework.
  152. void CEwaveDlg::OnPaint() 
  153. {
  154. if (IsIconic())
  155. {
  156. CPaintDC dc(this); // device context for painting
  158. SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);
  160. // Center icon in client rectangle
  161. int cxIcon = GetSystemMetrics(SM_CXICON);
  162. int cyIcon = GetSystemMetrics(SM_CYICON);
  163. CRect rect;
  164. GetClientRect(&rect);
  165. int x = (rect.Width() - cxIcon + 1) / 2;
  166. int y = (rect.Height() - cyIcon + 1) / 2;
  168. // Draw the icon
  169. dc.DrawIcon(x, y, m_hIcon);
  170. }
  171. else
  172. {
  173. if(adt[0]==0)
  174. {
  175.     CPaintDC dcb(this);
  176.      CBrush br(RGB(0,0,0));
  177.      CBrush * pOldBr;
  178.      pOldBr=dcb.SelectObject(&br);
  179.      dcb.Rectangle(dOpx-50,dOpy-70,dOpx+420,dOpy+70);
  180.      dcb.SelectObject(pOldBr);
  182.      CClientDC dc(this);
  183.         CPen pen(PS_SOLID,1,RGB(0,255,0));
  184.          CPen *pOldPen;
  185.         pOldPen=dc.SelectObject(&pen);
  186.          dc.MoveTo(dOpx,dOpy-40);
  187.         dc.LineTo(dOpx,dOpy+40);
  188.          dc.MoveTo(dOpx+100,dOpy-40);
  189.         dc.LineTo(dOpx+100,dOpy+40);
  190.          dc.MoveTo(dOpx+200,dOpy-40);
  191.         dc.LineTo(dOpx+200,dOpy+40);
  192.          dc.MoveTo(dOpx+300,dOpy-40);
  193.         dc.LineTo(dOpx+300,dOpy+40);
  194.      dc.MoveTo(dOpx+400,dOpy-40);
  195.         dc.LineTo(dOpx+400,dOpy+40);
  196.      dc.MoveTo(dOpx,dOpy-40);
  197.         dc.LineTo(dOpx+400,dOpy-40);
  198.      dc.MoveTo(dOpx,dOpy);
  199.         dc.LineTo(dOpx+400,dOpy);
  200.      dc.MoveTo(dOpx,dOpy+40);
  201.         dc.LineTo(dOpx+400,dOpy+40);
  202.         dc.SelectObject(pOldPen);
  204. CClientDC dcf(this);
  205. CFont font,font1;
  206. font.CreateFont(15,0,0,0,560,0,0,0,ANSI_CHARSET,OUT_DEFAULT_PRECIS,
  207. OUT_DEFAULT_PRECIS,PROOF_QUALITY,FF_ROMAN,"a");
  208. font1.CreateFont(17,0,0,0,0,0,0,0,ANSI_CHARSET,OUT_DEFAULT_PRECIS,
  209.         OUT_DEFAULT_PRECIS,PROOF_QUALITY,FF_ROMAN,"a1");
  211. CFont * pOldFont;
  212.     pOldFont=dcf.SelectObject(&font);
  213. dcf.SetBkMode(TRANSPARENT);
  214. dcf.SetTextColor(RGB(255,255,0));
  215. dcf.TextOut(dOpx+5,dOpy+40,"0");
  216. dcf.TextOut(dOpx+95,dOpy+40,"10");
  217. dcf.TextOut(dOpx+195,dOpy+40,"20");
  218. dcf.TextOut(dOpx+295,dOpy+40,"30");
  219. dcf.TextOut(dOpx+395,dOpy+40,"40");
  220.             dcf.TextOut(dOpx-10,dOpy-6,"0");
  221. dcf.TextOut(dOpx+193,dOpy+53,"t(s)");
  222. dcf.TextOut(dOpx-45,dOpy-6,"a(gal)");
  223. pOldFont=dcf.SelectObject(&font1);
  224.         dcf.SetBkMode(TRANSPARENT);
  225.         dcf.SetTextColor(RGB(255,255,0));
  226.         dcf.TextOut(dOpx+105,dOpy-63,"人工地震波加速度时程曲线");
  227. dcf.SelectObject(pOldFont);
  228. }
  229. else
  230. {
  231. Wavedrawing();
  232. }
  235. CDialog::OnPaint();
  236. }
  237. }
  239. // The system calls this to obtain the cursor to display while the user drags
  240. //  the minimized window.
  241. HCURSOR CEwaveDlg::OnQueryDragIcon()
  242. {
  243. return (HCURSOR) m_hIcon;
  244. }
  246. void CEwaveDlg::Wavegener()
  247. {
  248. UpdateData(true);
  249. double dAmax=m_dAmax*9.8;
  250.     double adGpu[10000],adFpu[10000],ddOmiga=0.02,adOmiga[10000];
  251. double Omiga1=2*pi/(5*m_dTg),Omiga2=2*pi/m_dTg,Omiga3=20*pi;
  253. /////////////////////////////////////////////////////////////////////
  254. ///////得到反应谱
  255.     for(int loop=0;loop<10000;loop++)
  256. {
  257. adOmiga[loop]=(loop+1-0.5)*ddOmiga;
  258. if(adOmiga[loop]<=Omiga1)
  259. adFpu[loop]=pow(0.2,0.9)*dAmax-0.02*dAmax*(2*pi/adOmiga[loop]-5*m_dTg);
  260. else
  261. if(adOmiga[loop]>=Omiga1 && adOmiga[loop]<Omiga2)
  262. adFpu[loop]=pow(m_dTg*adOmiga[loop]/(2*pi),0.9)*dAmax;
  263. else
  264. if(adOmiga[loop]>=Omiga2 && adOmiga[loop]<Omiga3)
  265. adFpu[loop]=dAmax;
  266. else
  267. adFpu[loop]=11*dAmax*pi/adOmiga[loop]+0.45*dAmax;
  268. }
  270.     /////////////////////////////////////////////////////////////////////
  271. ///////得到功率谱
  272. for(loop=0;loop<10000;loop++)
  273. {
  274. if(adFpu[loop]<0)
  275. adFpu[loop]=0.0;
  276. adGpu[loop]=2*0.05*pow(adFpu[loop],2)/(-2*(pi*adOmiga[loop])*log(-pi*log(0.85)/(adOmiga[loop]*40)));
  277. }
  279. /////////////////////////////////////////////////////////////////////
  280. ///////生成地震波
  281. double adFi[10000];
  282. int i;
  283.     srand((unsigned)time( NULL ));
  284. for(loop=0;loop<10000;loop++)
  285. {
  286.         int temp=rand();
  287. temp=temp%6282+1;
  288.     adFi[loop]=double(temp)/1000;         //随机相位角
  289. }
  290. for(loop=0;loop<2000;loop++)
  291. adt[loop]=0.02*(loop+1);
  292. for(loop=0;loop<2000;loop++)
  293. {
  294. adAg[loop]=0;
  295. i=0;
  296. while(i<10000)
  297. {
  298. adAg[loop]+=pow(0.08*adGpu[i],0.5)*cos(adOmiga[i]*adt[loop]+adFi[i]);   //人工波地面加速度
  299. i++;
  300. }
  301. }
  302. for(loop=0;loop<200;loop++)
  303. adAg[loop]*=pow(adt[loop]/4,2);
  304. for(loop=750;loop<2000;loop++)
  305. adAg[loop]*=pow(2.71828,-0.15*(adt[loop]-15));  //乘包络线
  306. for(loop=0;loop<2000;loop++)
  307. adAg[loop]*=100;                            //转换单位
  309.     for(loop=0;loop<2000;loop++)
  310. {
  311. if(adAg[loop]<0)
  312. adAgtemp[loop]=-adAg[loop];
  313. else
  314. adAgtemp[loop]=adAg[loop];
  315. }
  316. amax=0;
  317. for(loop=0;loop<2000;loop++)
  318. {
  319. if(adAgtemp[loop]>amax)
  320. {amax=adAgtemp[loop];
  321. tmax=adt[loop];
  322. }
  323. }                                              //求峰值
  324. for(loop=0;loop<2000;loop++)
  325. adAg[loop]*=m_pd/amax;                     //修正峰值
  327. amax=0;
  328.  for(loop=0;loop<2000;loop++)
  329. {
  330. if(adAg[loop]<0)
  331. adAgtemp[loop]=-adAg[loop];
  332. else
  333. adAgtemp[loop]=adAg[loop];
  334. }
  336. for(loop=0;loop<2000;loop++)
  337. {
  338. if(adAgtemp[loop]>amax)
  339. {amax=adAgtemp[loop];
  340. tmax=adt[loop];
  341. }
  342. }
  344.     UpdateData(false);
  345. }
  346.     /////////////////////////////////////////////////////////////////////
  347. ///////绘制地震波加速度时程曲线
  348. void CEwaveDlg::Wavedrawing()
  349. {
  350. Invalidate();
  351. CPaintDC dcb(this);
  352. CBrush br(RGB(0,0,0));
  353. CBrush * pOldBr;
  354. pOldBr=dcb.SelectObject(&br);
  355. dcb.Rectangle(dOpx-50,dOpy-70,dOpx+420,dOpy+70);
  356. dcb.SelectObject(pOldBr);
  359. CClientDC dc(this);
  360. CPen pen(PS_SOLID,1,RGB(255,0,0));
  361. CPen *pOldPen;
  362. CPen pen1(PS_SOLID,1,RGB(0,255,0));
  363.     
  364. pOldPen=dc.SelectObject(&pen);
  365. dc.MoveTo(dOpx,dOpy);
  366. for(int loop=0;loop<2000;loop++)
  367. dc.LineTo(dOpx+adt[loop]*10,dOpy-adAg[loop]*40/amax);
  369. pOldPen=dc.SelectObject(&pen1);
  370.     dc.MoveTo(dOpx,dOpy-40);
  371.     dc.LineTo(dOpx,dOpy+40);    
  372. dc.MoveTo(dOpx+100,dOpy-40);
  373. dc.LineTo(dOpx+100,dOpy+40);
  374.     dc.MoveTo(dOpx+200,dOpy-40);
  375. dc.LineTo(dOpx+200,dOpy+40);
  376.     dc.MoveTo(dOpx+300,dOpy-40);
  377.     dc.LineTo(dOpx+300,dOpy+40);
  378. dc.MoveTo(dOpx+400,dOpy-40);
  379.     dc.LineTo(dOpx+400,dOpy+40);
  380. dc.MoveTo(dOpx,dOpy-40);
  381.     dc.LineTo(dOpx+400,dOpy-40);
  382. dc.MoveTo(dOpx,dOpy);
  383.     dc.LineTo(dOpx+400,dOpy);
  384. dc.MoveTo(dOpx,dOpy+40);
  385.     dc.LineTo(dOpx+400,dOpy+40);
  386. dc.SelectObject(pOldPen);
  388. CString str,str1;
  389. str.Format("%3.0f",amax);
  390.     str1.Format("%3.0f",-amax);
  391. CClientDC dcf(this);
  392. CFont font,font1;
  393. font.CreateFont(15,0,0,0,560,0,0,0,ANSI_CHARSET,OUT_DEFAULT_PRECIS,
  394. OUT_DEFAULT_PRECIS,PROOF_QUALITY,FF_ROMAN,"a");
  395. font1.CreateFont(17,0,0,0,0,0,0,0,ANSI_CHARSET,OUT_DEFAULT_PRECIS,
  396. OUT_DEFAULT_PRECIS,PROOF_QUALITY,FF_ROMAN,"a1");
  398. CFont * pOldFont;
  399.     pOldFont=dcf.SelectObject(&font);
  400. dcf.SetBkMode(TRANSPARENT);
  401. dcf.SetTextColor(RGB(255,255,0));
  402. dcf.TextOut(dOpx+5,dOpy+40,"0");
  403. dcf.TextOut(dOpx+95,dOpy+40,"10");
  404. dcf.TextOut(dOpx+195,dOpy+40,"20");
  405. dcf.TextOut(dOpx+295,dOpy+40,"30");
  406. dcf.TextOut(dOpx+395,dOpy+40,"40");
  407.     dcf.TextOut(dOpx-10,dOpy-6,"0");
  408. dcf.TextOut(dOpx+193,dOpy+53,"t(s)");
  409. dcf.TextOut(dOpx-45,dOpy-6,"a(gal)");
  410. dcf.TextOut(dOpx-25,dOpy-45,str);
  411.     dcf.TextOut(dOpx-25,dOpy+35,str1);
  412. pOldFont=dcf.SelectObject(&font1);
  413. dcf.SetBkMode(TRANSPARENT);
  414. dcf.SetTextColor(RGB(255,255,0));
  415. dcf.TextOut(dOpx+105,dOpy-63,"人工地震波加速度时程曲线");
  417.     dcf.SelectObject(pOldFont);
  418.     UpdateWindow();
  420. }
  422. void CEwaveDlg::OnSTART() 
  423. {
  424. // TODO: Add your control notification handler code here
  425. Wavegener();
  426. Wavedrawing();
  427. }
  429.     /////////////////////////////////////////////////////////////////////
  430. ///////数据输出
  431. void CEwaveDlg::OnOutput() 
  432. {
  433. // TODO: Add your control notification handler code here
  434. CString sFileName;
  435. ofstream fout;
  436. CFileDialog dlg(false);
  437. char sFilter[]="All file(*.*)*.*";
  438. dlg.m_ofn.lpstrFilter=sFilter;
  440. if(dlg.DoModal()==IDOK)
  441. {
  442. sFileName=dlg.GetPathName();
  443. if(sFileName.GetLength()==0)return;
  444. fout.open(sFileName);
  445. }
  446. fout<<"t"<<" "<<tmax<<"    "<<"a"<<" "<<amax<<endl;
  447. for(int loop=0;loop<2000;loop++)
  448. fout<<adt[loop]<<" "<<adAg[loop]<<endl;
  449. fout.close();
  450. }
  453.     /////////////////////////////////////////////////////////////////////
  454. ///////数据输入
  455. void CEwaveDlg::OnInput() 
  456. {
  457. // TODO: Add your control notification handler code here
  458. CString sFileName;
  459. ifstream fin;
  460. char tempt,tempa;
  461. double temptmax,tempamax;
  462. CFileDialog dlg(true);
  463. char sFilter[]="All file(*.*)*.*";
  464. dlg.m_ofn.lpstrFilter=sFilter;
  466. if(dlg.DoModal()==IDOK)
  467. {
  468. sFileName=dlg.GetPathName();
  469. if(sFileName.GetLength()==0)return;
  470. fin.open(sFileName);
  471. }
  472.     fin>>tempt>>temptmax>>tempa>>tempamax;
  473. for(int loop=0;loop<2000;loop++)
  474. fin>>adt[loop]>>adAg[loop];
  475. fin.close();
  477.  for(loop=0;loop<2000;loop++)
  478. {
  479. if(adAg[loop]<0)
  480. adAgtemp[loop]=-adAg[loop];
  481. else
  482. adAgtemp[loop]=adAg[loop];
  483. }
  484.  
  485.  amax=0;
  486. for(loop=0;loop<2000;loop++)
  487. {
  488. if(adAgtemp[loop]>amax)
  489. {amax=adAgtemp[loop];
  490. tmax=adt[loop];
  491. }
  492. }
  493. Wavedrawing();
  494. UpdateData(false);
  495. }