http://www.gomatlab.de/index.php goMatlab.de - Deutschsprachiges Forum für die Benutzung der Software MATLAB / Simulink. Hilfe bei Bedienungsfragen oder bei der Umsetzung technischer Probleme. de http://www.gomatlab.de/index.php http://www.gomatlab.de/templates/vSpiralFull/images/_header_logo.gif 0 http://www.gomatlab.de/viewtopic.php?p=200560#200560 Wed, 02 Jun 2021 22:18:54 +0200 Hallo ich weiß nicht genau ob ich hier im richtigen Forum bin, aber ich bin am verzweifeln kann mir jemand vielleicht den Code schicken der hinter der hough Transformations funktion steht da ich gerne den Code für ein anderes Projekt abwandeln möchte aber meine selbst programmierte bough transformation funktioniert nicht 0 http://www.gomatlab.de/viewtopic.php?p=192700#192700 Wed, 13 Nov 2019 08:27:44 +0100 Hallo zusammen, ich befasse mich seit neuen erst mit matlab und versuche aktuell die Phasenverschiebung zweier Signale zu berechnen. Mein Ansatz war folgender: - Ich habe die Nullstellen beider Signale ermittelt - Anschließend berechne ich mir die Differenz zwischen den Nullstellen und lasse mir die Signalwerte der beiden Signale für die Nullstellen eines Signals ausgeben. - Jetzt dachte ich ich könnte mir die Differenz der Werte berechnen und mir darauf mit Hilfe der Abtastrate meine Phasenverschiebung ermitteln. Aber irgendwie kommt da nichts sinnvolles zustande. Hat jemand möglicherweise einen besseren Ansatz für mein Problem oder hatte schonmal so ein Problemr? vielen Dank :) 18 http://www.gomatlab.de/viewtopic.php?p=192241#192241 Thu, 03 Oct 2019 16:19:19 +0200 Guten Tag zusammen, Ich möchte gern Daten einer Tabelle (1.Spalte: Zeit, 2.Spalte: Geschwindigkeit) über eine FFT in den Frequenzbereich bringen und anschliessend die Geschwindigkeit integrieren, um auf die Wegamplitude (Strecke) zu kommen und wieder über ifft rücktransformieren - ich habe dazu folgenden Code erarbeitet: clear all close all clc dataset=xlsread("PK_F00_MatLab.xlsx","PK_F01","A2:B31251"); t=dataset(:,1) %Definition der ersten Spalte für Zeit v_01=dataset(:,2) %Definition der zweiten Spalte für Geschwindigkeit L_v=20*log10((v_01/1000)/5e-08); subplot(2,2,1) plot(t,L_v); xlabel("t in s") %x-Achsenbeschriftung ylabel("L_v in dB") %y-Achsenbeschriftung title("Zeitbereich - Schallschnellepegel") %Diagrammtitel f=fft(t); %Fast-Fouriertransformation von Zeit in Frequenz subplot(2,2,2) %siehe oben plot(abs(f),L_v); xlim([0 60]); xlabel("f in Hz") ylabel("L_v in dB") title("Frequenzbereich - Schallschnellepegel") dt=t(2)-t(1); s_01=dt*cumtrapz(v_01); %Integral für Wertetabellen ( subplot(2,2,3) %siehe oben plot(abs(f),s_01); xlim([0 50]); xlabel("f in Hz") ylabel("s in mm") title("Frequenzbereich - Strecke") T=ifft(f); %inverse Fast-Fouriertransformation von Frequenz zurück in Zeit subplot(2,2,4) %siehe oben plot(T,s_01); xlim([0 0.5]) xlabel("t in s") ylabel("s in mm") title("Zeitbereich - Strecke") ... 0 http://www.gomatlab.de/viewtopic.php?p=190461#190461 Sat, 25 May 2019 13:18:07 +0200 ich brauche Hilfe ich muss diese VBA skript in Matlab schreiben: anbei ist die VBA skipt und meine Matlab code danke sehr im voraus If anz < 0 Then Ti0_mK = "-" Exit Function End If If anz > 8 Then Ti0_mK = "n. zul." Exit Function End If Pi = 3.14159265358979 'Zahl pi (Konstante) l = 0.1 'fiktive Matrizenlänge 1 [m], kürzt sich später raus alpha_i = 1E+16 'Wärmeübergangskoeffizient Pressteil (Stahl) unendlich hoch [W/(m²*K)] a40_lambda = 0 'a0-Wert für Kern hier angegeben, da sonst zu viele Elemente in Klammer oben (Variablenangabe) Ai0 = Pi * di0 / 1000 * l 'Mantelfläche Kern-Innen-Durchmesser [m²] Aa0 = Pi * da0 / 1000 * l 'Mantelfläche Kern-Außen-Durchmesser [m²] Aa1 = Pi * da1 / 1000 * l Aa2 = Pi * da2 / 1000 * l Aa3 = Pi * da3 / 1000 * l Aa4 = Pi * da4 / 1000 * l Aa5 = Pi * da5 / 1000 * l Aa6 = Pi * da6 / 1000 * l Aa9 = Pi * da9 / 1000 * l 'fiktive Mantelfläche Außen Werkzeuaufnahme [m²] T0_Start = 1 T1_Start = 1 T2_Start = 1 T3_Start = 1 T4_Start = 1 T5_Start = 1 T6_Start = 1 T7_Start = 1 T8_Start = 1 T9_Start = 1 Ti0_ = 0 Ta0_ = 0 Ta1_ = 0 Ta2_ = 0 Ta3_ = 0 Ta4_ = 0 Ta5_ = 0 Ta6_ = 0 Ta7_ = 0 Ta8_ = 0 Ta9_ = 0 T0 = Ta0_ - (Ta0_ - Ti0_) / 2 T1 = Ta1_ - (Ta1_ - Ta0_) / 2 T2 = Ta2_ - (Ta2_ - Ta1_) / 2 T3 = Ta3_ - (Ta3_ - Ta2_) / 2 T4 = Ta4_ - (Ta4_ - Ta3_) / 2 T5 = Ta5_ - (Ta5_ - Ta4_) / 2 T6 = Ta6_ - (Ta6_ - Ta5_) / 2 T7 = Ta7_ - (Ta7_ - Ta6_) / 2 T8 = Ta8_ - (Ta8_ - Ta7_) / 2 T9 = Ta9_ - (Ta9_ - Ta8_) / 2 While T0_Start > T0 + 0.5 Or T0_Start < T0 - 0.5 T0_Start = T0 'Wärmeleitfähigkeit [W/(m*K)] Kern bei vorliegender Temperatur im Kern lambda_0 = a40_lambda * T0 ^ 4 + a30_lambda * T0 ^ 3 + a20_lambda * T0 ^ 2 + a10_lambda * T0 + a00_lambda 'Wärmeleitfähigkeit [W/(m*K)] 1. Arm.-Ring bei vorliegender Temperatur im 1. Arm ... 1 http://www.gomatlab.de/viewtopic.php?p=187565#187565 Wed, 07 Nov 2018 09:47:40 +0100 Hallo! Ich habe mir ein Minimalbeispiel der Filterung mit dem window-sinc_filter zusammengebaut und dabei festgestellt, dass das gefilterte Signal mehr Samples enthält als das Ausgangssignal. Das liegt wahrscheinlich an der hohen Filterordnung des Bandpasses vermute ich. Das Problem der ganzen Geschichte liegt nun darin, dass es sich um ein Zeitsignal handelt, welches gefiltert werden soll. Wie kann ich nun also die Werte wieder richtig "sortieren" sodass ich im gefilterten Signal die korrekten Zeitpunkte erhalte? Gruß Christian load sensorData fs = Fs; order = 500; fc1 = 700; fc2 = 800; filter_type = 'bandpass'; window_type = 'hamming'; analyse_plot = 'y'; % Filterfrequenzantwort darstellen % Berechnung der Impulsantwort = Filterkoeffizienten % Bandpass im Bereich 10...100 Hz h = window_sinc_filter(order, fc1, fc2, Fs, filter_type, window_type, analyse_plot); % Signal mit dem Bandpass filtern gefiltertes_Signal = FFT_Faltung(s1,h); % äquivalent zur Faltung im Zeitbereich gefiltertes_Signal2 = conv(s1,h); % Faltung im Frequenzbereich ist jedoch bei höherer Filterordnung schneller als conv() figure hold on plot(s1) plot(gefiltertes_Signal) 0 http://www.gomatlab.de/viewtopic.php?p=181647#181647 Wed, 17 Jan 2018 17:19:01 +0100 // 3d-Analyse 6 // WAV Dateien mit Abtastfrequenzen von 22005 und 96000 Hz können genutzt werden // Vermutlich können auch Dateien mit anderen Abtastfrequenzen genutzt werden. // Abtastfrequenzen von 22000 bis 96000 sind möglich. Ich bevorzuge 48000. // Es werden FFT Frequenzen von 1 bis 250 Hz berechnet // Es werden Zeitabschnitte von 60 Sekunden berechnet // Das Signal kann von 1 bis 10 verstärkt werden // Es kann zwischen Kanal 1 und Kanal 2 gewechselt werden // Es wird im Textfeld Position Sek. die Laufzeit der WAV Datei angezeigt // Die Start-Position des zu berechnenden Abschnittes wird mit dem Zähler ausgewählt // Es kann zwischen FFT Darstellung und Power Darstellung gewählt werden // Es kann die Grafik Anzeige auf der Frequenz Achse gezoomt werden // Die Zoomfrequenz kann mit dem Zähler eingestellt werden // Mit der Taste Zoom Grafik wird die eingestellte Zoom Frequenz aktiviert // Mit der Taste Save Grafik wird die angezeigte Grafik in einer gif Datei gespeichert // Mit der Taste Start werden alle oben eingestellten Voreinstellungen aktiviert und neu berechnet // mit der Playposition wird ein Startpunkt ausgewählt und eine Sequenz // von 10 Sekunden Länge wird von diesem Startpunkt an abgespielt // Mit der Taste Neu wird eine neue WAV Datei ausgewählt // Nit der Taste Exit wird die Sitzung beendet ///////////////////////////////////////////////////////////////////////////// mclose('all'); clear; xdel(); z=struct('pfad','','datei','','pfaddatei','','info','','laufzeit',0,'samplerate',0,'ueberdeckung',0,'gleich',0,'posi',1,'start1',0,'start2',0,'wa',0,'ff',0,'power',0,'aaff',0,'aapp',0,'fo1',0,'fo2',500,'zeit',0,'a1',0,'a2',0,'a3',0,'b1',0,'b2',0,'b3',0,'h1',0,'graph1',0,'f',0,'cm',0,'kanal',1,'verstaerkung',1,'wa',0,'darst',1,'spiel',0,'spiel1',0); ///////////////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////// ... 0 http://www.gomatlab.de/viewtopic.php?p=177267#177267 Mon, 10 Jul 2017 01:59:23 +0200 SCILAB ist klasse!!! Der Skript läuft unter Scilab 6 Eine Audiodatei, Kanal 1, wird je Sekunde analysiert. Die gesamte Audio-Datei kann von Sekunde zu Sekunde durchgeblättert werden. Die entsprechende Sekunde wird in folgender Weise analysiert. Grafisch angezeigt werden Wav, Gleichspannungsanteil, FFT, Leistungsspektrum, und mit der FFT wird das ursprüngliche Wav Signal wieder berechnet. Mit einstellbarer oberer Grenzfrequenz. // Signalanalyse 6 mclose('all'); clear; xdel(); z=struct('wa',0,'ft',0,'ftl',0,'ftr',0,'st1',0,'st2',0,'lz',0,'g',0,'da',0,'pf',0,'re',0,'sf',0,'pfda',0,'info',0,'vs',1,'h1',0,'f1',0,'h2',0,'f2',0,'h3',0,'f3',0,'h4',0,'f4',0,'h5',0,'f5',0,'po',0,'fi',0); ///////////////////////////////////////////////////////////////////////////// function [z]=neudatei(z) ////////////////////////////////////////////////////////////////////////////// // WAV Datei auswählen D:\Mdaten\Testsignale/100Hz-03A_g.wav [z.da,z.pf]=uigetfile('*.wav','D:\Mdaten'); //WAV Datei wählen if isequal(z.da,''); //Ist eine Date gewählt? else z.pfda=z.pf+'\'+z.da; //Vollständigen Pfad erstellen z.info = wavread(z.pfda,'info'); //WAV Info`s laden disp('//////////////////////////////////////////////////////////') disp('WAV-Datei: '+z.pfda); //Ausgabe vollständiger Pfad disp('Kanäle: '+string(z.info(2))); //Ausgabe Anzahl der Kanäle disp('Samplerate: '+string(z.info(3))); //Ausgabe Samplerate disp('Dateilänge : '+string(z.info(8))); //Ausgabe Dateilänge disp('Laufzeit : '+string(z.info(8)/z.info(3))+' Sekunden'); //Laufzeit in Sek. disp('//////////////////////////////////////////////////////////') z.lz=floor(z.info(8)/z.info(3)); // lauzeit in sekunden z.sf=z.info&# ... 7 http://www.gomatlab.de/viewtopic.php?p=172345#172345 Fri, 20 Jan 2017 09:04:41 +0100 Hallo DSP, der Thread ist ja leider schon etwas älter, jedoch hoffe ich doch auf eine Antwort. Meine Frage betrifft den Korrekturfaktor (N / sum(window)). Dieser erscheint mir sinnvoll ich würde aber gerne wissen wie man darauf kommt. Ich habe schon in vielen Foren nach dem Korrekturfaktor gesucht, aber oft wird der einfach weggelassen oder was ich auch schon gesehen habe ist einfach ein konstanter Faktor von z.B. 2. Literatur dazu habe ich noch nicht durchstöbert dafür bin ich im Alltag zu weit weg von dem Thema Signalverarbeitung, wäre aber für Referenzen sehr dankbar. Vielen Dank im Voraus, [/quote]Gruß Brian 2 http://www.gomatlab.de/viewtopic.php?p=168164#168164 Wed, 31 Aug 2016 09:59:10 +0200 So, nun habe ich den weiteren Fehler erkannt. Die Samplerate und die FFT-Buffergröße spielen eine ganz wichtige Rolle und sollten bei verschiedenen Fragestellungen richtig gewählt werden. Ich habe nun folgendes festgestellt, in meinem Beispiel ist die Samplerate kleiner als die FFT-Buffergröße. Damit bekommt man eine Frequenzauflösung von 0.15625Hz. D.h. 5 Hz lassen sich gut auf der x-achse darstellen, denn 5 Hz ist das 32-fache von der Frequenzauflösung. Das Einzige, die Amplitude stimmt nicht und die ist scheinbar auch von der Samplerate und FFT-Buffergröße abhängig. Nehme ich z.B. eine Samplerate von 160 Hz und eine FFT-Buffergröße von ebenfalls 160, dann habe ich eine Frequenzauflösung von 1 Hz. Hier stimmt die x-Achse haargenau und die y-achse zeigt fast exakt die gleiche Amplitude im Amplitudenspektrum, wie im Zeitbereich an. Falls jemand noch ergänzende Anmerkungen hat, sind diese gerne willkommen! 16 http://www.gomatlab.de/viewtopic.php?p=143477#143477 Sun, 21 Dec 2014 19:06:19 +0100 Hat alles wunderbar funktioniert!:-)