使用数据文件中的一组地震检波器测量值确定可能的地震事件的位置。

算法设计： 1) 读取文件头并分配内存；
2) 从数据文件读取地震数据，从键盘读取计算能量的短时间和长时间窗口测量 值的数目;
3) 计算各个时间点上的短时间窗口和长时间窗口的能量值，打印出可能的地震 事件时间

#include<fstream>
#include<string>
#include<iostream>
#include<Windows.h>
#define THRESHOULD 1.5
using namespace std;
double power_w(double *arr,int length,int n); 
//功能：计算短/长时间窗口能量的采样值
//length从第几个开始 n表数量
int main(){
 string filename;
 ifstream fin;
 double time;
 double *sensor;
 int num,short_count=0,long_count=0;
 double short_pow,long_pow;
 cout << "请输入文件名：";
 cin >> filename;
 fin.open(filename.c_str());
 if(fin.fail()){
  cout<<"打开失败"<<endl;
  exit(-1);
 }else{fin>>num>>time;
  if(num>0){
   sensor=new double[num];
   for(int i=0;i<num;i++){
    fin>>sensor[i];
   }
   cout<<"请输入短时间和长时间测量的数目：";
   cin>>short_count>>long_count;
   for(int i=long_count;i<=num;i++){
    short_pow=power_w(sensor,i,short_count);
    long_pow=power_w(sensor,i,long_count);
    if((short_pow/long_pow)>THRESHOULD){
     cout<<"possible event at"<<time*i<<"seconds\n";
    }
   }
   delete[] sensor;
  }
  fin.close();
 }
 system("pause");
 return 0;
}
double power_w(double *arr,int length,int n){
 double avr=0;
 for(int i=length-1;i>=length-n;i--){
  avr+=arr[i]*arr[i];
 }
 return avr/n;
}