(MATLAB)fvtool(FilterObject)和fvtool(FilterObject.SOSMatrix)之间的差异
我正在设计一个数字滤波器来放大低频,我使用IIR的任意大小。代码如下所示:(MATLAB)fvtool(FilterObject)和fvtool(FilterObject.SOSMatrix)之间的差异,matlab,filtering,Matlab,Filtering,我正在设计一个数字滤波器来放大低频,我使用IIR的任意大小。代码如下所示: % Acquisition Frequency FreqAq = 1024; % Low Frequency Gain LFGain = 10^(14.4/20); % Frequency Intervals FreqIntv = [2 6 600]/FreqAq; % Numerator Order NumOrder = 40; % Denominator Order DenOrder = 45; % Fre
% Acquisition Frequency
FreqAq = 1024;
% Low Frequency Gain
LFGain = 10^(14.4/20);
% Frequency Intervals
FreqIntv = [2 6 600]/FreqAq;
% Numerator Order
NumOrder = 40;
% Denominator Order
DenOrder = 45;
% Frequency Spectrum
Freq = linspace(0,FreqAq,FreqAq + 1)/FreqAq;
% Amplitude Design
bLF = log(1/LFGain)/(FreqIntv(1) - FreqIntv(2));
aLF = LFGain*exp(bLF*FreqIntv(1));
A = [LFGain*ones(1,find(Freq == FreqIntv(1))),...
aLF*exp(-bLF*Freq(find(Freq == FreqIntv(1)) + 1:find(Freq == FreqIntv(2)))),...
ones(1,find(Freq == FreqIntv(3)) - find(Freq == FreqIntv(2))),...
ones(1,length(Freq) - find(Freq == FreqIntv(3)))];
% Arbitrary response magnitude filter specification object
d = fdesign.arbmag('Nb,Na,F,A',NumOrder,DenOrder,Freq,A);
% Filter Weights
W = [ones(1,find(Freq == FreqIntv(1)))*1 ...
ones(1,length(Freq(find(Freq == FreqIntv(1)) + 1:find(Freq == FreqIntv(2)))))*.1 ...
ones(1,find(Freq == 10/FreqAq,1,'first') - find(Freq == FreqIntv(2)))*1 ones(1,find(Freq == FreqIntv(3)) - find(Freq == 10/FreqAq,1,'first'))...
ones(1,length(Freq) - find(Freq == FreqIntv(3)))*1];
% Filter Design
Hd = design(d, 'iirlpnorm', 'Weights', W, 'Norm', 2, 'DensityFactor',30,...
'SystemObject',true);
过滤数据时,我没有得到我期望的结果
DataFiltered = sosfilt(Hd.SOSMatrix,Data)
为了验证错误,我绘制了:
fvtool(Hd, 'MagnitudeDisplay','Magnitude (dB)', ...
'NormalizedFrequency','On','Color','White');
fvtool(Hd.SOSMatrix, 'MagnitudeDisplay','Magnitude (dB)', ...
'NormalizedFrequency','On','Color','White');
注意到第二个偏移了+35dB,这是用于滤波器的偏移
如果我过滤数据的方式不正确,或者大小的计算方式有任何差异,我似乎不明白我的代码出了什么问题。任何信息都将不胜感激
先谢谢你