地震检测系统

最新推荐文章于 2024-10-07 19:01:55 发布

weixin_54124514

最新推荐文章于 2024-10-07 19:01:55 发布

阅读量481

点赞数

文章标签： c++

本文链接：https://blog.csdn.net/weixin_54124514/article/details/126874983

版权

该程序读取输入文件中的传感器数据，通过计算短/长时间窗口的能量值比来检测地震活动。用户输入短窗口和长窗口的点数，当比值超过阈值1.5时，标记可能的地震事件。程序使用了`power_w`函数来统计能量值。

摘要由CSDN通过智能技术生成

#include<fstream>
#include<iostream>
#include<Windows.h>
#include<string>
#include<cmath>
using namespace std;

const double THRESHOLD = 1.5;

//计算短/长时间窗口能量数据的采样值。
double power_w(double arr[], int length, int n);

int main()
{
   string filename;
   ifstream fin;
   int num =0,short_windows=0,long_windows=0;
   double time_incr =0,*sensor=NULL,short_power=0,long_power=0;
   double ratio;
   cout << "Enter name of input file" << endl;
   cin >> filename;
   fin.open(filename.c_str());
   if (fin.fail())
   {
       cerr << "error opening input file" << endl;
       exit(-1);
   }
   else
   {
       fin >> num >> time_incr;
       cout << "num:" << num << " time_incr:" << time_incr << endl;

       if (num >= 0)
       {
           sensor = new double[num];

           for (int i = 0; i < num; i++)
           {
               fin >> sensor[i];
           }
           cout << "Enter number of points for short-windows" << endl;
           cin >> short_windows;

cout << "Enter number of points for long-windows" << endl;
cin >> long_windows;

           //分析能量数据找出地震时间
           for (int i = long_windows - 1; i < num; i++)
           {
               short_power = power_w(sensor, i, short_windows);
               long_power = power_w(sensor, i, long_windows);
               ratio = short_power / long_power;

               if (ratio> THRESHOLD)
               {
                   cout << "Possible event at"<<time_incr*i << "seconds"<<endl;
               }

}
delete[] sensor;

}

fin.close();
}

   system("pause");
   return 0;
}
//统计短/长事件窗口对应的能量值
double power_w(double arr[], int length, int n)
{
   double xsquare = 0;
   for (int i = 0; i < n; i++)
   {
       xsquare += pow(arr[length - i], 2);
   }
   return xsquare / n;
}