
% EQ mit Cos-Fenster-Spektrum und Yulewalkeralgorithmus
%=====================================================
clear

fs=40000;  % Abtastfrequenz (eigentlich 40000Hz)
f=0:20000; % fg (eigentlich 20000Hz)


% linearer Spektrum
%===========================
H(1:length(f))=1; % 0Hz-20kHz
H(1:20) = 0;      % 20Hz-20kHz


% Band 1
%===========================
gain=2;
fm=50;
fu=20;
fo=70;
f1=fu:fo;
H1=1+ gain.*(cos((f1-fm)./((fo-fu)/2) .*(pi/2)));
H(fu:fo)=H(fu:fo).* H1;
%===========================
% Band 1 Ende


% Band 2
%===========================
gain2=2;
fm2=100;
fu2=70;
fo2=130;
f2=fu2:fo2;
H2=1+ gain2.*(cos((f2-fm2)./((fo2-fu2)/2) .*(pi/2)));
H(fu2:fo2)=H(fu2:fo2).* H2;
%===========================
% Band 2 Ende


% Band 3
%===========================
gain3=2;
fm3=200;
fu3=130;
fo3=270;
f3=fu3:fo3;
H3=1+ gain3.*(cos((f3-fm3)./((fo3-fu3)/2) .*(pi/2)));
H(fu3:fo3)=H(fu3:fo3).* H3;
%===========================
% Band 3 Ende


% Band 4
%===========================
gain4=2;
fm4=400;
fu4=270;
fo4=530;
f4=fu4:fo4;
H4=1+ gain4.*(cos((f4-fm4)./((fo4-fu4)/2) .*(pi/2)));
H(fu4:fo4)=H(fu4:fo4).* H4;
%===========================
% Band 4 Ende

% Band 5
%===========================
gain5=2;
fm5=800;
fu5=530;
fo5=1070;
f5=fu5:fo5;
H5=1+ gain5.*(cos((f5-fm5)./((fo5-fu5)/2) .*(pi/2)));
H(fu5:fo5)=H(fu5:fo5).* H5;
%===========================
% Band 5 Ende

% Band 6
%===========================
gain6=2;
fm6=1600;
fu6=1070;
fo6=2130;
f6=fu6:fo6;
H6=1+ gain6.*(cos((f6-fm6)./((fo6-fu6)/2) .*(pi/2)));
H(fu6:fo6)=H(fu6:fo6).* H6;
%===========================
% Band 6 Ende


% Band 7
%===========================
gain7=2;
fm7=3200;
fu7=2130;
fo7=4270;
f7=fu7:fo7;
H7=1+ gain7.*(cos((f7-fm7)./((fo7-fu7)/2) .*(pi/2)));
H(fu7:fo7)=H(fu7:fo7).* H7;
%===========================
% Band 7 Ende


% Band 8
%===========================
gain8=-1;
fm8=5700;
fu8=4270;
fo8=7130;
f8=fu8:fo8;
H8=1+ gain8.*(cos((f8-fm8)./((fo8-fu8)/2) .*(pi/2)));
H(fu8:fo8)=H(fu8:fo8).* H8;
%===========================
% Band 8 Ende


% Band 9
%===========================
gain9=2;
fm9=8565;
fu9=7130;
fo9=10000;
f9=fu9:fo9;
H9=1+ gain9.*(cos((f9-fm9)./((fo9-fu9)/2) .*(pi/2)));
H(fu9:fo9)=H(fu9:fo9).* H9;
%===========================
% Band 9 Ende


% Band 10
%===========================
gain10=2;
fm10=12000;
fu10=10000;
fo10=14000;
f10=fu10:fo10;
H10=1+ gain10.*(cos((f10-fm10)./((fo10-fu10)/2) .*(pi/2)));
H(fu10:fo10)=H(fu10:fo10).* H10;
%===========================
% Band 10 Ende

% Band 11
%===========================
gain11=2;
fm11=16000;
fu11=14000;
fo11=18000;
f11=fu11:fo11;
H11=1+ gain11.*(cos((f11-fm11)./((fo11-fu11)/2) .*(pi/2)));
H(fu11:fo11)=H(fu11:fo11).* H11;
%===========================
% Band 11 Ende





subplot(3,1,1);
semilogx(f,H);
axis([0,20000, 0, 15]);
xlabel('Frequenz [Hz]');
ylabel('Verstärkung');
legend('Spektrum mit Cos-Fenster',[2]);

omega=f./(fs/2);
[B,A]=yulewalk(20,omega,H);
[H_yule,f_yule]=freqz(B,A,20000);
subplot(3,1,2);
plot(abs(f_yule.*(fs/(2*pi))),20*log10(abs(H_yule)));
axis([-5,20000, -5, 15]);
xlabel('Frequenz [Hz]');
ylabel('Amplitude [dB]');
legend('Spektrum mit Yule-Walker',[3]);

subplot(3,1,3);
semilogx(abs(f_yule.*(fs/(2*pi))),20*log10(abs(H_yule)));

xlabel('Frequenz [Hz]');
ylabel('Amplitude [dB]');
legend('Spektrum mit Yule-Walker. X-Achse logarithmisch',[3]);

% abtasten und abspielen
% =====================================================
x=wavread('sound3.wav');
z=filter(B,A,x); 
sound(z,22000); 