linear_fft.asv

来自「谢谢站长！谢谢给我这次下载的机会！非常感谢！现在上传一些我仿真的多普勒功率普的程」· ASV 代码 · 共 12 行

%file:linear_fft.m
% this function takes n (must be even)time domain samples(real or complex)
% and finds the psd by taking (fft/n)^2
% the two sided spectrum is produced by shifting the psd.
% the array freq provides the appropriate frequency values for plotting purposes
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=fft(y)/n;