算法设计: 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;
}