hc471.m

来自「这是数字信号处理——matlab释义与实现一书的源代码」· M 代码 · 共 16 行

%《数字信号处理教程——MATLAB释义与实现》第四章例4.7.1程序hc471
% 求矩形窗频谱能量的程序hc471
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
T=0.01;N=16384; a=2;    % 输入原始数据
D=2*pi/(N*T);          % 频率分辨率
M=floor(a/T);           % 求边界点下标
x=[ones(1,2*M+1)];      % 给出时间序列
X=fft(x,N);            % 求傅立叶变换FFT
k1=input('k1=（书上取零） '),       % 求能量的起点频率下标
k2=input('k2=（书上取floor(0.5*pi/D)）'),  % 求能量的终点频率下标
k=k1:k2;                % 求能量频段的下标序列
E=T/N*sum(abs(X(k+1)).^2)   % 求正频率区间能量
E1=2*T/N*sum(abs(X(k+1)).^2) % 求全频率区间能量
E1/4.01*100                 % 占总能量百分比