conio.cpp

DSP中输入信号的生成过程。 conio.cpp实现X(n)信号

/*
以下是输入信号X(n)的生成过程。X(n)中有两个正弦信号，
频率分别是140Hz和70Hz,同时还含有高斯白噪声。
该信号共产生2000个点，这些点都存储在 "x.txt"文本中。
*/
#include "stdio.h"
#include "sinwn.c"
main()
{int i,m,n;
long seed;
double fs,snr,x[2000];
static double a[2]={1,1};
static double f[2]={70,140};
static double ph[3]={45,88};
FILE *fp;
m=2;
n=2000;
seed=13579l;
fs=2800;
snr=5.0;
sinwn(a,f,ph,m,fs,snr,seed,x,n);
printf("\nTwo sinusoidal signals plus gauss white noise\n");
for(i=0;i<2000;i++)
{printf("%10.7lf",x[i]);
 if(i%4==3)
 printf("\n");
}
fp=fopen("x.txt","w");
for(i=0;i<n;i++)
{fprintf(fp,"%12.7lf\n",x[i]);}
fclose(fp);
}