下面是用3种等效方法估计脉冲响应的代码:
clc,清晰
inpSig = [1 1 1];
outSig = [1,4,8,10,8,4,1];
将带零的输入和输出信号填充到一个公共长度(1024个样本)
comLength = 1024;
inpSig, 0 (1, comLength - length(inpSig))];
outSigPad = [outSig, zero (1, comLength - length(outSig))];
%方法1
%直接用Matlab估计脉冲响应
这个函数使用韦尔奇平均周期图来估计信号的psd
我在这里使用矩形窗口，重叠为0%(相当于常规fft)
[H1, f] = tfestimate(inpSigPad, outSigPad, rectwin(comLength)， 0， []， 1，“中心”）;
%方法2
您可以估计交叉和自动psd使用Matlab
%假设采样频率是单一的
pxx = cpsd(inpSigPad, inpSigPad, rectwin(comLength)， 0， []， 1，“中心”）;
pxy = cpsd(inpSigPad, outSigPad, rectwin(comLength)， 0， []， 1，“中心”）;
传递函数的%定义
H2 = pxy./pxx;
数字
情节(f, mag2db (abs (H1)));
持有在
情节(f, mag2db (abs (H2)),“——”）
%方法3
使用FFT估计输入和输出的PSD
X = fftshift(fft(inpSigPad));
Y = fftshift(fft(outSigPad));
%重新排序样本使频率响应在(-fs/2, fs/2]之间
X = [X(2:end)， X(1)];
Y = [Y(2:结束)，Y(1)];
H3 = y (:)./ x (:);
情节(f(1:30:结束),mag2db (abs (H3(1:30:结束))),“o”）;网格小
包含(“频率-赫兹”), ylabel (“幅度- dB”）
传奇({“使用tfestimate”，“使用运行cpsd”，“使用FFT”}）
