下面是代码,估计脉冲响应与3等效的方法:
clc,清晰
inpSig = (1 1 1);
outSig = (1、4、8、10、8、4、1];
%垫与零输入和输出信号,一个常见的长度(1024个样本)
comLength = 1024;
inpSigPad = [inpSig, 0 (1, comLength -长度(inpSig)));
outSigPad = [outSig, 0 (1, comLength -长度(outSig)));
%方法1
%估计脉冲响应直接使用Matlab
%这个函数使用韦尔奇的平均周期图估计信号的psd
%我使用矩形窗重叠0%(相当于普通fft)
(H1, f) = tfestimate (inpSigPad、outSigPad rectwin (comLength), 0, [], 1,“中心”);
%的方法2
%可以估计十字架和汽车使用Matlab psd
%假设统一的采样频率
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));
%重排样品之间的频率响应(f / 2, fs / 2)
X =[(2:结束),X (1)];
Y =[(2:结束),Y (1)];
H3 = Y (:)。/ X (:);
情节(f(1:30:结束),mag2db (abs (H3(1:30:结束))),“o”);网格小
包含(频率-赫兹的),ylabel (数据库级的)
传奇({“使用tfestimate”,“使用运行cpsd”,“使用FFT”})
