|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
产生对称的频率之后,再取离散正傅里叶变换,能回实域?
有本书上,把频域相应,直接往F轴的正方向对称。再FFT得出了其实域相应。困惑。
function[impw]=jakes_filter(fd)
n=512;nn=2*n;
fs=0:fd/64:fd;
H=zeros(1,n);
for k=1:(n/8+1)
jpsd(k)=1/((1-((fs(k))/fd)^2+eps)^0.5);
if(jpsd(k))>1000
jpsd(k)=1000;
end
H(k)=jpsd(k)^0.5;
end
for k=1:n
H(n+k)=H(n+1-k);
end
[inv,time]=linear_fft(H,nn,fd/64);
imp=real(inv(450:577));
impw=imp.*hanning(128)';
energy=sum(impw.^2);
impw=impw/(energy^0.5);
function[fftx,freq]=linear_fft(x,n,ts)
y=zeros(1,n);
for k=1:n
freq(k)=(k-1-(n/2))/(n*ts);
y(k)=x(k)*((-1.0)^(k+1));
end
fftx=ifft(y)/n; |
|