《matlab论坛求助：信号加窗后的处理,看不懂程序后面的原理.docxVIP

Matlab中文论坛求助帖：求助：信号加窗后的处理，看不懂程序后面的原理。len = 48000;winType = hody; % 这是什么窗呢？winMinPhase = -pi; % 最小相位winMaxPhase = pi - 2 * pi / len; % 最大相位x = linspace(winMinPhase, winMaxPhase, len); % 均分switch winType case hody win = 0.06372625603133489 + ... cos(x) * 0.23940452419651895 + ... cos(x) .^ 2 * 0.35000203221155923 + ... cos(x) .^ 3 * 0.24817008904762164 + ... cos(x) .^ 4 * 0.08582735035413477 + ... cos(x) .^ 5 * 0.01242535720360149 + ... cos(x) .^ 6 * 0.0004443909552289012; hw = 7; ......enddcCorrection = 1 / sum(win);dcFft = fft(win);dcFft = dcFft(1 : min(hw, len)) * dcCorrection;win_gain = sqrt(1 / (len * sum(win .^ 2)));sRate=len;t = 1;n = 0:1/sRate:t - 1/sRate;m = floor(length(n)/2);input = [zeros(m,1);ones(length(n)-m,1)];xfrmLen = floor(length(input) / 2) + 1; % 频谱半长度+1，对应频率为pi/2inputFft = fft(input .* win); % inputFft，加窗后的频谱dc = inputFft(1) .*dcCorrection;inputFft(1 : hw) = inputFft(1 : hw) - (dcFft * inputFft(1));inputFft(1) = inputFft(1) / sqrt(2); % 计算 inputFft第一项的有效值if rem(length(inputFft), 2) == 0 % 若length(inputFft)为偶 inputFft(xfrmLen) = inputFft(xfrmLen) / sqrt(2); % 计算inputFft中间一个的有效值endxfrm = inputFft(1 : xfrmLen)*sqrt(2); % 默认信号为模拟信号，此处求有效值% 因为：2* inputFft(1 : xfrmLen)*sqrt(2) /2 = inputFft(1 : xfrmLen)*sqrt(2)xfrm = win_gain * xfrm; % 这是干什么呢？xfrm = xfrm .* conj(xfrm); %这又是得到什么？功率谱密度？将信号加窗后做频谱变换，后续涉及到的一些处理过程都和窗函数有关，看不懂一些处理过程。想请教下专家，谢谢了！~~~有疑问的地方基本都在程序注释里面写了。首先，第5行，没见过“hody”窗，不知道是啥窗（当然，窗的形状这里是知道的），网上没查到hody窗，我想是否是哪个专业应用？第二，39、40行，是去掉频谱中直流分量的，并且和窗函数有关，我想知道为什么这样处理。平时去直流分量直接将时域信号减去均值后再做频谱变换吧？但这里，不在时域处理，在频域，对频谱数据的前几项做这样的处理，是什么原理？可否公式说明？ （那个hw可以认为是和此hody窗相关的一个因子。）第三，49行：xfrm = win_gain * xfrm；这里win_gain 是个什么因子，肯定和窗函数有关，这里是干嘛呢？如果说前面是将频谱数据去掉直流分量影响，并将模拟信号的幅度谱变为其有效值谱，那么这里将有效值谱乘以win_gain ，即是除以 sqrt(len * sum(win .^ 2))，是何意？最后，最后一句，紧接着上面一句，如果说上句计算的变量有可能是无意义的，只是一个中间量，那么这句是程序最后一句，得到的量应该是有意义的了。为了保证变量都是有意义的，我们从原去掉了直流成分影响的频谱数据inputFft ，推导到此处xfrm，具体是：xfrm = (win_gain * sqrt(2)) * inputFft(1 :xfrmLen).*conj((wi_gain *