oct3dsgn.m

来自「车辆工程平顺性方面的1/3倍中心频率的有关程序」· M 代码 · 共 13 行

function [B,A] = oct3dsgn(Fc,Fs,N)
  
pi = 3.14159265358979;
f1 = Fc/(2^(1/6)); 
f2 = Fc*(2^(1/6)); 
Qr = Fc/(f2-f1); 
Qd = (pi/2/N)/(sin(pi/2/N))*Qr;
alpha = (1 + sqrt(1+4*Qd^2))/2/Qd; 
W1 = Fc/(Fs/2)/alpha; 
W2 = Fc/(Fs/2)*alpha;
[B,A] = butter(N,[W1,W2]);