基本波信号に 3 次高調波を追加し、フィルターを使用して基本波を取り戻すこのコードがあります。3 次、5 次、7 次の高調波を追加してフィルターで除去するには、このコードを変更する必要がありますが、フィルターでそれを行う方法がわかりません。
t = [0:199];
A1 = 100;
s1 = A1*sin(2*pi*(1/40)*t); % fundamental
A3 = A1/3;
s3 = A3*sin(2*pi*(1/40)*3*t); %3rd harmonic
s13 = s1 + s3;
% A5 = A1/5;
% s5 = A5*sin(2*pi*(1/40)*5*t); 5th
% A7 = A1/7;
% s7 = A7*sin(2*pi*(1/40)*7*t); 7th
% filter
[b13, a13] = ellip(6,0.5,20,[5.7/40 6.3/40],'stop') %elliptic filter
h13 = impz(b13,a13,length(s13)); %impulse
y13 = filter(b13,a13,s13);