地震道振幅统计代码c语言

在C语言中编写地震道振幅统计的代码通常涉及信号处理和数据分析的部分。这里是一个简单的示例，用于计算一组地震数据中振幅的最大值、最小值和平均值：

#include <stdio.h>
#include <stdlib.h>

#define MAX_LENGTH 1000 // 你可以根据需要调整数组长度

// 结构体来存储地震道的数据
typedef struct {
    double amplitude[MAX_LENGTH];
    int length;
} SeismicTrace;

// 功能：读取地震数据并计算振幅统计
void calculate_stats(SeismicTrace *trace) {
    FILE *fp = fopen("earthquake_data.txt", "r"); // 假设数据文件名是earthquake_data.txt
    if (fp == NULL) {
        printf("Error opening file.\n");
        return;
    }

    // 从文件中读取数据，填充结构体
    for (int i = 0; i < trace->length; i++) {
        fscanf(fp, "%lf", &trace->amplitude[i]);
    }
    fclose(fp);

    // 计算并打印振幅统计
    double max_amp = *max_element(trace->amplitude, trace->amplitude + trace->length);
    double min_amp = *min_element(trace->amplitude, trace->amplitude + trace->length);
    double avg_amp = (double)sum(trace->amplitude, trace->length) / trace->length; // 假设sum函数计算总和

    printf("Max amplitude: %.2f\n", max_amp);
    printf("Min amplitude: %.2f\n", min_amp);
    printf("Average amplitude: %.2f\n", avg_amp);
}

int main() {
    SeismicTrace trace;
    trace.length = MAX_LENGTH;

    calculate_stats(&trace);

    return 0;
}