prefourier.m

来自「滤波器的设计以及零极点对应频率响应分析」· M 代码 · 共 10 行

function [t,omg,FT,IFT]=prefourier(Trg,N,OMGrg,K)
%用矩阵左乘来实现傅里叶变换前的准备工作
%输入：TRG，时域起止范围，N，时域抽样点数，OMGrg，频域起止范围和频域抽样点数
%输出： 时域和频域抽样点数，傅里叶变换和逆变换
F=Trg(2)-Trg(1); %时域范围
t=linspace(Trg(1),Trg(2)-T/N,N); %生成抽样时间点
OMG=OMGrg(2)-OMGrg(1);          %频域范围
omg=linsapce(OMGrg(1),OMGrg(2)-OMG/K,K); %生成抽样频率点
FT=F/N*exp(-j*kron(omg,t.'));            %构造带有系数的傅里叶变换矩阵
IFT=OMG/2/pi/K*exp(j*kron(t,omg.'));     %构造带有系数的的傅里叶逆变换矩阵