地震预测系统项目实现

?整个项目思路即在一组观测数据中，地震专家（即用户）输入观测窗口的最小数量和最大数量，进行预测峰值点

?

?

数据文件如图所示：

#define _CRT_SECURE_NO_WARNINGS
#include<fstream>
#include<string>
#include<iostream>
#include<cmath>

using namespace std;

const double THRESHOLD = 1.5;
//长时间窗口的平方和除个数，与短时间窗口的平方和除个数



//计算短/长时间窗口能量数据的采样值
//统计窗口对应能量值的采样
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;

}

int main() {
	string filename;
	ifstream fin;//从文件读取数据的输入文件流
	int num=0,short_window=0,long_window=0;//定义读取能量数据的长度是11
	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 piont for short_window" << endl;
			cin >> short_window;
			cout << "Enter number of piont for long_window" << endl;
			cin >> long_window;

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

				ratio = short_power / long_power;
				if (ratio > THRESHOLD) {
					cout << "Passible event at "<<time_incr*i<<"seconds\n";
				}
			}

		}
		

		delete[] sensor;
		fin.close();
	}
	
	return 0;


}

运行结果：

?在0.05秒和0.06秒到达地震，

?