📄 dsfb.m
字号:
function Z=dsfb(Rc,Uc)
N=1024;
f3db=30;
f=-f3db:2*f3db/N:f3db;
S=exp(-(1.665*f/f3db).^2);
s=ifft(S(1:N));
r=log(1+s*(exp(Rc)-1))/Rc;
R=fft(r);
RF=abs(R);
a0=RF(1,1);
RF=RF./a0;
RF=sqrt(RF);
HW=RF;
v=rands(1,N);
V=fft(v);
W=HW.*V;
w=ifft(W);
u=sqrt(Rc)*w+log(Uc);
u=abs(u);
z=exp(u);
y=xcorr(z);
Z=fft(y(1:N+1));
Z=10*log10(abs(Z(1:N/2+1))/sum(abs(Z)));
plot(Z);
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -