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

Matla?b中文论坛?求助帖：求助：信号加窗后?的处理，看不懂程序?后面的原理?。len = 48000?;winTy?pe = hody; % 这是什么窗?呢？winMi?nPhas?e = -pi; % 最小相位winMa?xPhas?e = pi - 2 * pi / len; % 最大相位x = linsp?ace(winMi?nPhas?e, winMa?xPhas?e, len); % 均分switc?hwinTy?pecase hodywin = 0.06372?62560?31334?89 + ...cos(x) * 0.23940?45241?96518?95 + ...cos(x) .^ 2 * 0.35000?20322?11559?23 + ...cos(x) .^ 3 * 0.24817?00890?47621?64 + ...cos(x) .^ 4 * 0.08582?73503?54134?77 + ...cos(x) .^ 5 * 0.01242?53572?03601?49 + ...cos(x) .^ 6 * 0.00044?43909?55228?9012;hw = 7; ......enddcCor?recti?on = 1 / sum(win);dcFft? = fft(win);dcFft? = dcFft?(1 : min(hw, len)) * dcCor?recti?on;win_g?ain = sqrt(1 / (len * sum(win .^ 2)));sRate?=len;t = 1;n = 0:1/sRate?:t - 1/sRate?;m = floor?(lengt?h(n)/2);input? = [zeros?(m,1);ones(lengt?h(n)-m,1)];xfrmL?en = floor?(lengt?h(input?) / 2) + 1; % 频谱半长度?+1，对应频率为?pi/2input?Fft = fft(input? .* win); % input?Fft，加窗后的频?谱dc = input?Fft(1) .*dcCor?recti?on;input?Fft(1 : hw) = input?Fft(1 : hw) - (dcFft? * input?Fft(1));input?Fft(1) = input?Fft(1) / sqrt(2); % 计算input?Fft第一?项的有效值?if rem(lengt?h(input?Fft), 2) == 0 % 若leng?th(input?Fft)为偶input?Fft(xfrmL?en) = input?Fft(xfrmL?en) / sqrt(2); % 计算inp?utFft?中间一个的?有效值endxfrm = input?Fft(1 : xfrmL?en)*sqrt(2); % 默认信号为?模拟信号，此处求有效?值% 因为：2* input?Fft(1 : xfrmL?en)*sqrt(2) /2 = input?Fft(1 : xfrmL?en)*sqrt(2)xfrm = win_g?ain * xfrm; % 这是干什么?呢？xfrm = xfrm .* conj(xfrm); %这又是得到?什么？功率谱密度?？将信号加窗?后做频谱变?换，后续涉及到?的一些处理?过程都和窗?函数有关，看不懂一些?处理过程。想请教下专?家，谢谢了！~~~有疑问的地?方基本都在?程序注释里?面写了。首先，第5行，没见过“hody”窗，不知道是啥?窗（当然，窗的形状这?里是知道的?），网上没查到?hody窗?，我想是否是?哪个专业应?用？第二，39、40行，是去掉频谱?中直流分量?的，并且和窗函?数有关，我想知道为?什么这样处?理。平时去直流?分量直接将?时域信号减?去均值后再?做频谱变换?吧？但这里，不在时域处?理，在频域，对频谱数据?的前几项做?这样的处理?，是什么原理?？可否公式说?明？（那个hw可?以认为是和?此hody?窗相关的一?个因子。）第三，49行：xfrm = win_g?ain * xfrm；这里win?_gain?是个什么因?子，肯定和窗函?数有关，这里是干嘛?呢？如果说前面?是将频谱数?据去掉直流?分量影响，并将模拟信?号的幅度谱?变为其有效?值谱，那么这里将?有效值谱乘?以win_?gain，即是除以sqrt(len * sum(win .^ 2))，是何意？最后，最后一句，紧接着上面?一句，如果说上句?计算的变量?有可能是无?意义的，只是一个中?间量