Sie muss im Pfad sein (also am einfachsten im aktuellen Ordner)
Betriebssystem Windows 7
>> [v,tv] = reedMIC("evos.wav","csv",1,1500,50)
warning: wavread is obsolete and will be removed from a future version of Octave, please use audioread instead
warning: called from
wavread at line 62 column 5
reedMIC at line 15 column 9
v = [](0x1)
tv = [](0x1)
>>
>>Ich müsste fast am Ziel sein:
>> [v,tv] = reedMIC("ohne_Prellen.wav","Output.csv",0.015,1.448,200)
error: 'tv' undefined near line 68 column 6
error: called from
reedMIC at line 68 column 1
>>
Ok das bedeutet,
if length(tvup)==length(tvdown)
tv = mean([tvup,tvdown],2);
v = mean([vup,vdown],2);
endBin ich zwischenzeitlich ohne deine Rückmeldung auch draufgekommen. Dass es so empfindlich auf den Schwellenwert reagiert ?!
>> reedMIC("ohne_Prellen.wav","Output.csv",0.02,1.500,200)
Da habe ich mich verzweifelt auch schon versucht - aber immer wenn ich ins Fenster
Diese klappt jedoch nicht, da mit s=0.004 die Oben-Unten-Pärchen noch nicht erkannt werden und s=0.005 bereits die Fehlermeldung auf Grund ungleicher Oben-Unten-Ereignisse erfolgt.
function [v,tv] = reedMIC(wavfile,csvfile,s,U,vmax,varargin)
% [v,tv] = reed(wavfile,csvfile,s,U,vmax) gibt den Geschwindigkeitsvektor v und
% den dazugehörigen Zeitvektor tv zurück. wavfile ist der Pfad zur
% .wav-Datei; csvfile ist der Pfad zur Ausgabedatei; s ist der
% Schwellenwert ab dem ein Peak erkannt werden soll; U ist der
% Radumfang; vmax ist die theoretisch maximale erkennbare Geschwindigkeit
% Es findet keine Filterung der Daten statt.
%
% [v,tv] = reed(...,dofilter) ist dofilter > 0 wird ein
% Bandpass mit den Grenzfrequenzen 3.2kHz und 3.5kHz angewendet.
%
% [v,tv] = reed(...,F3dB1,F3dB2) die Variablen F3dB1 und
% F3dB2 sind die untere bzw. die obere Grenzfrequenz für den Bandpass.
[data,fs]=wavread(wavfile);
dataUF = data;
if length(varargin)>0 && varargin{1}>0
F3dB1 = 3.2e3;
F3dB2 = 3.5e3;
if length(varargin)==2
F3dB1 = varargin{1};
F3dB2 = varargin{2};
end
D = fdesign.bandpass('N,F3dB1,F3dB2',2,F3dB1,F3dB2,fs);
Hd= design(D);
[B,A] = sos2tf(Hd.sosMatrix);
data = filtfilt(B,A,data).*Hd.ScaleValues(1);
end
eventlength=U*fs/vmax*3.6;
ww = eventlength/2;
if mod(ww,2) ~= 0
ww = ww -1;
end
diffdata=diff(data);
t=[1:length(diffdata)]./fs;
tevtup=t(diffdata >s);
tevtdown=t(diffdata <-s);
tevtup([0 diff(tevtup)<ww/fs]>0)=[];
tevtdown([0 diff(tevtdown)<ww/fs]>0)=[];
figure();hold on;
plot(t,diffdata,'b-')
plot(tevtup,0,'r*')
plot(tevtdown,0,'cs')
figure();hold on;
vup=[U./diff(tevtup)*3.6]';
vdown=[U./diff(tevtdown)*3.6]';
tvup=zeros(length(vup),1);
tvdown=zeros(length(vdown),1);
for i=1:length(tevtup)-1
tvup(i) = mean(tevtup(i:i+1));
end
for i=1:length(tevtdown)-1
tvdown(i) = mean(tevtdown(i:i+1));
end
plot(tvup,vup,'r-');
plot(tvdown,vdown,'c-');
if length(tvup)==length(tvdown)
tv = mean([tvup,tvdown],2);
v = mean([vup,vdown],2);
else
tv=unique([tvup;tvdown]);
v=mean([interp1(tvup,vup,tv),interp1(tvdown,vdown,tv)],2);
end
plot(tv,v,'b-');
xlabel('Zeit [s]');ylabel('Geschwindigkeit [km/h]');
export=[tv,v];
if ~isempty(csvfile)
save(csvfile,'export','-ascii','-double','-tabs');
end reedMIC('evos.wav','evos.csv',0.02,1.5,200);Ich denke, dass das Smartphone die Aussteuerung automatisch zurücknimmt. Der bzw. die ersten Impulse haben immer einen hohen Pegel.
Werde ich untersuchen. Es ist jedenfalls nicht linear: während ich bei 5-Sekunden-Abschnitten locker auf das Ergebnis in weniger als einer Minute warten kann, hat es für obige 60s mehr als eine Stunde gedauert und ich war schon am Abbrechen als er dann doch noch fertig wurde. Darf gar nicht an die Messstrecke bei Elfershausen denken, dort sind es jeweils ca. 10 Minuten mit ca. 1000 Radumdrehungen - aber wir müssen diese ja nicht unbedingt von Gnu-Plot zeichnen lassen.
Weder - noch, lediglich etliche Anfangswerte weggelassen, bei welchen der Tic auf Grund des Schwellwerts zusätzlich als Toc erkannt wurde. Dort ist der Algorithmus noch ungeeignet.
)
Mag sein, das kann ich nicht wirklich beurteilen.
Im Grunde hätten wir das damals schon schlussfolgern können:
Wir (und andere HerstellerBei höherer Geschwindigkeit würde diese ominöse bremsende Kraft also abnehmen
Beim Milan Mk2 ist es noch extremer.
) waren aber von der Konstanz des Cw-Werts überzeugt.