% A-Bewertungsfilter im Laplace-Bereich
s = tf('s');
ka = 7.39705*10^9;
% Quelle: Wiki
A_filter = (ka * s^4) / ((s + 129.4)^2 * (s + 676.7) * (s + 4636) * (s + 76655)^2);

fa = 30000; % Hz
df = 2^15/fa;
f11 = df:df:15000; % Hz
w = 2*pi .* f11; % Umwandlung in rad/s
% A-Berwertungsfilter bisher
Ra = ( (12200^2 .* f11.^4) ./ ((f11.^2 + 20.6^2) .* (f11.^2 + 12200^2) .* (sqrt(f11.^2 + 107.7^2)) .* (sqrt(f11.^2 + 737.9^2))) );

figure(1)
% Bodediagramm mit Matlabbefehl 
bode(A_filter);
grid on;

% Bodediagramm selber erstellen und A_filter und Ra vergleichen
figure(2) 
[mag phase] = bode(A_filter,w); 
% Betrag 
subplot(211) ; 
semilogx([w./(2*pi)],20*log10(mag(1,:))) ;
hold on;
% hier muss 2.00 dB addiert werden um das gleiche Spektrum wie A_filter zu
% haben -> siehe Quelle:
% http://www.sengpielaudio.com/BerechnungDerBewertungsfilter.pdf
semilogx(f11,20*log10(Ra) + 2.00,'r--') ;
hold off;
title('Bode Diagramm') ; 
xlabel('Frequenz in [Hz]') ; 
ylabel('Magnitude in [dB]') ; 
grid on; 
hold on; 
subplot(212) ; 
semilogx([w./(2*pi)],unwrap(angledim(phase(1,:),'radians','degrees'))) ; 
xlabel('Frequenz in [Hz]') ; 
ylabel('Phase in [Grad]') ; 
grid on; 
hold off;