dsp_periogram.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

r=xcorr(un,1024,'coeff');
r=r(1025:2048);
SBT=fft(r,255);
w=-pi:1/127*pi:pi;
plot(w,SBT)