C语言分析假频现象（频谱折叠）

yeahamen

已于 2023-02-13 23:00:11 修改

阅读量677

点赞数

文章标签： c语言 java 开发语言

于 2023-01-29 12:44:14 首次发布

本文链接：https://blog.csdn.net/amen_yanwang/article/details/128787616

版权

本文通过C语言实现了一个简单的振幅谱分析程序，展示了如何从正弦波信号中提取频谱信息，并讨论了采样定理的重要性及其对频谱分析的影响。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

这里是用振幅谱代替频谱，前者振幅/后者=1/2*Π。

#include "stdio.h"
#include "math.h"
#include "malloc.h"
#include "stdlib.h"
#define PI 3.1415926
#define f_1 10//T_1 = 1/10
#define f_2 35//T_2 = 1/35
#define f_3 70//T_3 = 1/70
#define fk 1    
#define M 1025   
#define dt 0.01
#define N 100
int main (void)
{
    int i,n;
    FILE *fp1,*fp2;
    float an,bn,cn,t;
    float function[1025];
    fp1 = fopen("function_Cn(A).xls","w");
    fp2 = fopen("function.xls","w");
    //先写原信号至fp2
    for(i=-M/2;i<=M/2;i++)
    {
        t=i*dt;
        function[i+M/2] = sin(2*PI*f_1*t)+sin(2*PI*f_2*t)+sin(2*PI*f_3*t);
        fprintf(fp2,"%f\t%f\n",t,function[i+M/2]);
    }
    //再写振幅谱至fp1
    for(n=-2*N;n<=2*N;n++)
    {
        an=0;//必不可少，否则结果cn是一恒定常数
        bn=0;
        for(i=-M/2;i<=M/2;i++)
        {
            t=i*dt;
            function[i+M/2] = sin(2*PI*f_1*t)+sin(2*PI*f_2*t)+sin(2*PI*f_3*t);
            an = an + function[i+M/2]*cos(2.0*PI*n*fk*t);
            bn = bn + function[i+M/2]*sin(2.0*PI*n*fk*t);
        }
        cn = sqrt(an*an+bn*bn)*dt;
        fprintf(fp1,"%d\t%10.4f\n",n*fk,cn);//算一个写一个
    }
    fclose(fp1);
    fclose(fp2);
    return 0;
}