function Z=kfb(v)
N=1024;
f3db=40;
f=-f3db:2*f3db/N:f3db;
S=exp(-(1.665*f/f3db).^2);
s=ifft(S(1:N));
s1=0:0.1:1;
y=[0 0.0034 0.0134 0.0302 0.0538 0.0844...
0.1224 0.1682 0.2278 0.3536 1.0000];
p1=polyfit(y,s1,4);
r=polyval(p1,s);
pp=r.^10;
R=fft(r);
RF=abs(R);
a0=RF(1,1);
RF=RF./a0;
HW1=sqrt(RF);
PP=fft(pp);
PPF=abs(PP);
a1=PPF(1,1);
PPF=PPF./a1;
HW2=sqrt(PPF);
v1=randn(1,1024);
v2=randn(1,1024);
v3=randn(1,1024);
v4=randn(1,1024);
v5=randn(1,1024);
v6=randn(1,1024);
v7=randn(1,1024);
v8=randn(1,1024);
v9=randn(1,1024);
v10=randn(1,1024);
v11=randn(1,1024);
v12=randn(1,1024);
V1=fft(v1);
V2=fft(v2);
V3=fft(v3);
V4=fft(v4);
V5=fft(v5);
V6=fft(v6);
V7=fft(v7);
V8=fft(v8);
V9=fft(v9);
V10=fft(v10);
V11=fft(v11);
V12=fft(v12);
W1=HW1.*V1;W2=HW1.*V2;W3=HW1.*V3;W4=HW1.*V4;W5=HW1.*V5;
W6=HW1.*V6;W7=HW1.*V7;W8=HW1.*V8;W9=HW1.*V9;W10=HW1.*V10;
W11=HW2.*V11;W12=HW2.*V12;
w1=ifft(W1);w2=ifft(W2);w3=ifft(W3);w4=ifft(W4);w5=ifft(W5);
w6=ifft(W6);w7=ifft(W7);w8=ifft(W8);w9=ifft(W9);w10=ifft(W10);
w11=ifft(W11);w12=ifft(W12);
w1=w1.^2;w2=w2.^2;w3=w3.^2;w4=w4.^2;w5=w5.^2;
w6=w6.^2;w7=w7.^2;w8=w8.^2;w9=w9.^2;w10=w10.^2;
w11=w11.^2;w12=w12.^2;
Y=w1+w2+w3+w4+w5+w6+w7+w8+w9+w10;
X=w11+w12;
z=X+Y;
z=sqrt(z);
y=xcorr(z);
Z=fft(y(1:N+1));
Z=10*log10(abs(Z(1:N/2+1))/sum(abs(Z)));
plot(Z);
Tme.rar_behindocl_杂波_通讯/手机编程_雷达 杂波
版权申诉
143 浏览量
2022-07-14
06:00:16
上传
评论
收藏 2KB RAR 举报
钱亚锋
- 粉丝: 90
- 资源: 1万+
评论0