Filter für Preemphase

Moin Moin!

Ich bin wieder bei der Implementierung des Preemphase-Filters angekommen, nachdem ich jetzt alle Klausuren überstanden habe und den "Rest" meines Programmes halbwegs auf Vordermann gebracht habe.

Und ja, leider habe ich immer noch Fragen (Hilfe - ich fühle mich so dumm).

Zunächst hat sich was geändert: Die Abtastrate beträgt nicht mehr die hohen 300kHz sondern 44,1kHz. Das liegt daran, dass der Filter jetzt an einer anderen Stelle im Programm implementiert wird: Direkt nach dem Einlesen der .wav-Dateien.

Also zu den Fragen:
1. @Lucki: noch einmal wegen der Normierung: War das ein Ausprobieren nach dem Motto: "da muss zu Anfang eine 1 stehen, tuts aber nicht, also teile ich durch das Eingangssignal!"? Und wie kommt man darauf? Wäre diese Normierung dann auch im späteren Programm mit den Musikdateien notwendig?

2. Vergleich mit MatLab: Ich habe versucht (mit Hilfe) mir in MatLab die Koeffizienten für einen Filter berechnen zu lassen. Dies wurde per frequency-sampling durchgeführt (wenn ich nicht vollkommen daneben liege). Hier brauche ich aber mindestens einen Filter der Ordnung 30, damit das Ergebnis brauchbar aussieht. Diese Koeffizienten habe ich mit ausgeben lassen und in LabVIEW dann als Koeffizienten eingefügt.
Mit diesem Versuch lande ich aber ziemlich auf der Klappe, wie man im angehängten VI erkennen kann. Wieso brauche ich in MatLab so viele Koeffizienten und wieso sieht das Ergebnis in MatLab sehr gut, in LabVIEW dann aber wieder total daneben aus?


Immerhin glaube ich, mir Eurer Hilfe langsam den Horizont erkennen zu können

@RoLi: Danke für die Hinweise auf den Tief- und Hochpass. Ich denke, dass ich die vergessen hätte und dadurch später ordentlich geflucht hätte.

mfg
Sönke

PS: Da ich anscheinend keine MatLab-Dateien hochladen darf, kopier ich kurz hier den Quelltext rein, vielleicht/hoffentlich kann jemand damit etwas anfangen:

f0=44100; % Abtastrate: da direkt die .wav Datei bearbeitet wird: 44,1kHz
tau=50; % Zeitkonstante in Microsekunden (für Europa 50, für Amerika 75)
N=30; % Filterordnung

% gewünschte Übertragungsfunktion
fm=0:0.001:16; % Frequenzpunkte: von 0Hz bis 16kHz in 1Hz-Schritten
dp=10*log10(1+10^(-5)*(2*tau*fm).^2); % corresponding magnitudes [dB]

% define frequency points
fs=0:.01:f0/2000; %Frequenzpunkte zum Arbeiten
As=1+10^(-5)*(2*tau*fs).^2; %gewünschte Amplitudenanhebung

% FIR Filter entwerfen (linearer Phasengang)
b=fir2(N,(fs/(f0/2000)),As);
b=b/sum(b); % Die Verstärkung soll 1 betragen

[H,f]=freqz(b,1,1024,f0); % compute transfer function
figure(1)
plot(f/1000,10*log10(abs(H)),'b',fs,10*log10(As),'r'%,fm,dp,'g'
grid on
axis([0.1 11 0 18])
legend('entworfen', 'erwünscht'
set(gca, 'XScale', 'log'