编写C语言主函数测试程序，计算出uN(n)的功率谱，并用图示显示出功率谱；（为完成图示显示，可将功率谱数据结果以文件形式保存，然后在MATLAB中导入并进行图示）

抱歉，我无法直接为您编写完整的程序。但是，我可以提供一些思路和参考代码：

1. 首先，需要编写一个函数来计算uN(n)的功率谱。可以使用傅里叶变换实现：

#include <stdio.h>
#include <math.h>

#define PI 3.14159265358979323846

void power_spectrum(double *u, int N, double *P) {
    // 计算uN(n)的功率谱
    double re, im;
    for (int k = 0; k < N; k++) {
        re = im = 0.0;
        for (int n = 0; n < N; n++) {
            re += u[n] * cos(2 * PI * n * k / N);
            im -= u[n] * sin(2 * PI * n * k / N);
        }
        P[k] = re * re + im * im;
    }
}

其中，u是输入的信号，N是信号的长度，P是输出的功率谱。

2. 接下来，在主函数中读取信号，调用功率谱函数计算功率谱，并将结果保存到文件中：

int main() {
    // 读取信号数据
    int N = 100;  // 假设信号长度为100
    double u[N];  // 存储信号数据
    FILE *fp = fopen("signal.txt", "r");
    for (int i = 0; i < N; i++) {
        fscanf(fp, "%lf", &u[i]);
    }
    fclose(fp);

    // 计算功率谱
    double P[N/2+1];  // 存储功率谱
    power_spectrum(u, N, P);

    // 将结果保存到文件
    fp = fopen("power_spectrum.txt", "w");
    for (int k = 0; k <= N/2; k++) {
        fprintf(fp, "%d\t%lf\n", k, P[k]);
    }
    fclose(fp);

    return 0;
}

3. 最后，在MATLAB中导入功率谱数据，并用plot函数绘制功率谱图：

data = load('power_spectrum.txt');
plot(data(:,1), data(:,2));
xlabel('频率');
ylabel('功率');

希望这些代码可以为您提供一些帮助。请注意，这只是一个简单的参考示例，您可能需要根据实际情况进行修改和优化。