wn = [0.105254918858264, 0.153593077243090];
fs = 1/0.003;
bw = 0.0002;
我们= wn /(2 *π)/ (fs / 2);
[num1, den1] = iirnotch(我们(1),bw);
[num2, den2] = iirnotch(禾(2),bw);
f1 =特遣部队(num1, den1, 0.003);
f2 =特遣部队(num2, den2, 0.003);
fcomb = f1 * f2;
fnorm = fcomb / dcgain (fcomb);
disp ([DC模式1:num2str (dcgain (f1))))
disp ([DC模式2:num2str (dcgain (f2))))
disp (['DC for combined: 'num2str (dcgain (fcomb))))
disp (['DC for normalized: 'num2str (dcgain (fnorm))))