dsp_perio.m

来自「经过调试的平稳随机信号的AR模型」· M 代码 · 共 27 行

clear all
%signal genergate
sdb=20;
w1=2*pi*0.7;
w2=2*pi*0.9;
n=1:1024;
sn1=exp(j*n*w1-j*pi)+exp(j*n*w2-j*0.9*pi);
Vsn=var(sn1);
noise_power=Vsn/10^(sdb/10);
%noise genergate
p=noise_power;
N=1024;
a=sqrt(p);
v=rand(1,N);
v=v/std(v);
v=v-mean(v);
v=a*v
un=sn1+v

N=255;
Unk=fft(un,N)
w=-pi:1/127*pi:pi    
SPER=abs(Unk).^2/N;
plot(w,SPER)
title('SPER method')
ylabel('Sper(w)')
xlabel('w')