修改现有的代码对于FFT函数,您可以使用下面的代码:
函数Xk = myfft (x, n)
N =长度(x);% DFT的基本周期
如果N <
x = [x, 0 (1, n n)];
N = N;
结束
n = 0:1: n - 1;% n行向量
k = 0:1: n - 1;% k行向量
WN = exp (1 i2pi / N);% Wn因素
nk = n ' * k;%创建一个N×N矩阵nk的价值观
WNnk = WN。^ nk;% DFT矩阵
Xk = x * WNnk;%的行向量DFT系数
修改包括' n '一个额外的参数传递给函数指定所需的FFT输出的大小,和检查如果输入“x”的长度小于“n”,在这种情况下进行补零长度等于“n”。
同样,为传输线函数修改现有的代码,您可以使用下面的代码:
函数x_hat = myifft (ftx, n)
N =长度(ftx);% DFT的基本周期
如果N <
ftx = [ftx, 0 (1, n n)];
N = N;
结束
n = 0:1: n - 1;% n行向量
k = 0:1: n - 1;% k行向量
WN = exp (1 i2pi / N);% Wn因素
nk = n ' * k;%创建一个N×N矩阵nk的价值观
WNnk = WN。^ (nk);%,台湾矩阵
x_hat = (ftx * WNnk) / N;%,台湾值的行向量
x_hat =实际(x_hat);
