wbfb.m

来自「这是关于雷达杂波的仿真的matlab程序」· M 代码 · 共 42 行

function Z=wbfb(p,q)
N=1024;
f3db=40;
f=-f3db:2*f3db/N:f3db;
S=exp(-(1.665*f/f3db).^2);
s=ifft(S(1:N));

r=sqrt(s*(gamma(3)-gamma(2).^2)/gamma(2).^2);
R=fft(r);
RF=abs(R);
a0=RF(1,1);
RF=RF./a0;
HW=sqrt(RF);


v1=rands(1,N);
V1=fft(v1);
W1=HW.*V1;

w1=ifft(W1);

v2=rands(1,N);
V2=fft(v2);
W2=HW.*V2;

w2=ifft(W2);
 
w1=w1*sqrt((q.^p)/2);
w2=w2*sqrt((q.^p)/2);
u=w1.^2+w2.^2;

z=u.^(1/p);
y=xcorr(z);

Z=fft(y(1:N+1));

Z=10*log10(abs(Z(1:N/2+1))/sum(abs(Z)));
plot(Z);