随机信号做谱估计(1)-周期图法——C语言
随机信号模型
这里我们用X(n)=s(n)+w(n)来表示所分析的随机信号,其中s(n)是具有两个频率分量的正弦波,可表示为是s(n)=A1sin(2pif1+pi/3)+A2sin(2pif2+pi/4)。w(n)是随机高斯白噪声,它的谱线是一个直线,用它代表信号在传播中所夹杂的噪声不仅可以体现随机过程而且对于功率谱的求解也更为简单。
估计质量的评价
可以用均分误差以及偏差来评价估计质量的好坏,通过画出均分误差和均值随SNR变化的图来观察估计质量。
原理
实验代码
使用VS进行C语言编程,得到一系列数据。然后将这些数据保存到TXT文件中,再将文件类型改为.dat文件。用gnuplot调用这些数据,画出均分误差和均值随SNR变化的图。
代码(C语言):
#include <stdlib.h>
#include <time.h>
#include <math.h>
#include<cstdio>
#define pi 3.1415916
#define N 64 //样本点数
#define M 100 //信噪比取的个数
#define C 50 //相同信噪比进行次数
#define f1 0.1
#define f2 0.3
void xcorr(double* number, double* sn);
void awgn(double* num);
void dft(double* xn, double* xk);
int findPeaks(double* xn, int* address);
void rank(double* mse1, double* mse2, double* expect1, double* expect2, double* snr);
void prin(double* snr, double* mse1, double* mse2, double* approximate1, double* approximate2);
struct C_complex_double
{
double real;
double image;
}; //复数结构体
int main(int argc, char* argv[])
{
double xn[N] = {
.0 };
double Pn = .0;
double Pnoise[M] = {
.0 };
double noise[M] = {
.0 };
double Pn1 = 0.0;
double rk[N] = {
0.0 };
int address[N / 2] = {
0 };
double rk_flag[N / 2] = {
0.0 };
double x2n[N] = {
.0 };
double mse_f1[M] = {
.0 };
double mse_f2[M] = {
.0 };
double num[N] = {
.0 };
double approximate1[M] = {
.0 };
double approximate2[M] = {
.0 };
srand