|
|
|
Ergebnis von xcorr interpretieren |
|
| RamonaH |
Forum-Newbie
|
 |
Beiträge: 1
|
|
|
|
Anmeldedatum: 24.02.14
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
 |
|
Verfasst am: 26.02.2014, 17:31
Titel: Ergebnis von xcorr interpretieren
|
 |
| |
 |
 |
Hallo,
ich bin noch eher dabei mich mit den basic functions von matlab vertraucht zu machen, muss aber dringend meine Daten analysieren.
Meine Daten (Druckmessungen) bestehen deutlich aus 2 Phasen, und mein Problem liegt darin, diese von einander zu trennen.
Ein Ansatz war nun um ein spaeteres Stueck des Signals (in der Grafik rot) zu definieren und ueber andere Sequenzen zu legen. Hierbei wuerde ich gerne die Korrelation bestimmen (mit xcorr).
Das schaut dann so aus wie in der Datei.
Das dazugehoerige Skript schaut aus wie folgt:
load('20130626EP3 agonist GP1 2nd resp.mat');
%define pressure and time
p=data;
tijd=1:length(data);
t=(tijd./1000)+(start_sample/1000);
%choose selected area
pseg3=p(end-370000:end-100000);
tseg3=t(end-370000:end-100000);
%choose areas for cross corrleation
pseg4=p(end-640000:end-370000);
tseg4=t(end-640000:end-370000);
pseg5=p(1:1+270000);
tseg5=t(1:1+270000);
%plot data set, selected area and areas for cross correlation
figure (1)
plot(t,p);
hold on
plot(tseg3,pseg3,'r');
plot(tseg4,pseg4,'g')
plot(tseg5,pseg5,'m')
title('selected part');
%xcorr: kruiscorrelatie tussen signaal p en pseg, MAXLAG mee te geven
% data to compare from 2nd phase
x=p(end-640000:end-370000);
F=xcorr(x,pseg3,50000);
z=(-50000:1:50000);
figure(2)
hold on
plot(z,F,'g')
kk=max(F)
lkk=find(F==kk)-500001;
%data to compare from initial phase
x1=p(1:1+270000)
G=xcorr(x1,pseg3,50000);
z=(-50000:1:50000);
hold on
plot(z,G,'m')
km=max(G)
lkm=find(G==km)-500001;
Mein Problem besteht nun darin, dass ich eingentlich das Resultat nicht wirklich interpretieren kann (ich wuerde einen niedriegeren Wert fuer das erste Datenstueck erwarten als fuer das zweite).
Desweiteren muessten sich die Werte die auf der y-aches aufgetragen warden zwischen 0 und 1 befinden, oder irre ich mich da?
Ganz liebe Gruesse und vielen Dank fuer eure Hilfe,
Ramona
| Beschreibung: |
|
Download |
| Dateiname: |
ScreenHunter_370 Feb. 26 16.27.jpg |
| Dateigröße: |
286.75 KB |
| Heruntergeladen: |
343 mal |
_________________
If you can't explain it simply, you don't understand it well enough.
-Albert Einstein-
|
|
|
|
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
|
|
Impressum
| Nutzungsbedingungen
| Datenschutz
| FAQ
|  RSS
Hosted by:
Copyright © 2007 - 2025
goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
|
|
