WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Ausgabe von allen Abschnitte einen bestimmten Werte

 

hedro1234

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 19.11.2022, 16:56     Titel: Ausgabe von allen Abschnitte einen bestimmten Werte
  Antworten mit Zitat      
Hallo Zusammen,

ich bin neu in Matlab und ich brauche Hilfe mit paar Sachen.
ich habe ein Drucksignal, es ist ca. 5min lang. ich habe dieses Signal in 20 Abschnitte geteilt, so dass jedes Abschnitt 1,30 min ist. und der Abstand zwischen jedes Abschnitt 10s ist. (1.Abschnitt von 0 s bis 1,30 min; 2.Abschnitt von 10 s bis 1,40 min ; 3.Abschnitt von 20 s bis 1,50 min; .....usw.)
für jedes Abschnitt habe ich auch einen FFT Analyse durchgeführt.

meine Fragen sind:
1. wie kann ich im Code die Verteilung des Signals in 20 Abschnitte und die FFt und das Plotten vereinfachen
2. wie kann ich im Frequenzbereich von jedes Abschnitt zb. das 1. Wert von allen Abschnitte in einem Diagramm zeigen lassen

ich bin sehr Dankbar, falls mir jemand mir helfen kann

Liebe Gruße
Hedro

Code:


clear
clc

load p.mat

n=2^14;

p0= p.p;


p1= p0(1000:10000);
pabs1= abs(fft(p1,n));

p2= p0(2000:11000);
pabs2= abs(fft(p2,n));

p3= p0(3000:12000);
pabs3= abs(fft(p3,n));

p4= p0(4000:13000);
pabs4= abs(fft(p4,n));

p5= p0(5000:14000);
pabs5= abs(fft(p5,n));

p6= p0(6000:15000);
pabs6= abs(fft(p6,n));

p7= p0(7000:16000);
pabs7= abs(fft(p7,n));

p8= p0(8000:17000);
p2abs8= abs(fft(p8,n));

p9= p0(9000:18000);
pabs9= abs(fft(p9,n));

p10= p0(10000:19000);
pabs10= abs(fft(p10,n));

p11= p0(11000:20000);
pabs11= abs(fft(p11,n));

p12= p0(12000:21000);
pabs12= abs(fft(p12,n));

p13= p0(13000:22000);
pabs13= abs(fft(p13,n));

p14= p0(14000:23000);
pabs14= abs(fft(p14,n));

p15= p0(15000:24000);
pabs15= abs(fft(p15,n));

p16= p0(16000:25000);
pabs16= abs(fft(p16,n));

p17= p0(17000:26000);
pabs17= abs(fft(p17,n));

p18= p0(18000:27000);
pabs18= abs(fft(p18,n));

p19= p0(19000:28000);
pabs19= abs(fft(p19,n));

p20= p0(20000:29000);
pabs20= abs(fft(p20,n));

%%%%%

figure(1)

subplot(4,10,1)
plot(p1)
title ('1')
subplot(4,10,11)
plot(pabs1)
title ('1')
%axis([0 2000 0 3.5*10^4 ])

subplot(4,10,2)
plot(p2)
title ('2')
subplot(4,10,12)
plot(pabs2)
title ('2')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,3)
plot(p3)
title ('3')
subplot(4,10,13)
plot(pabs3)
title ('3')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,4)
plot(p4)
title ('4')
subplot(4,10,14)
plot(pabs4)
title ('4')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,5)
plot(p5)
title ('5')
subplot(4,10,15)
plot(pabs5)
title ('5')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,6)
plot(p6)
title ('6')
subplot(4,10,16)
plot(pabs6)
title ('6')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,7)
plot(p7)
title ('7')
subplot(4,10,17)
plot(pabs7)
title ('7')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,8)
plot(p8)
title ('8')
subplot(4,10,18)
plot(p2abs8)
title ('8')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,9)
plot(p9)
title ('9')
subplot(4,10,19)
plot(pabs9)
title ('9')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,10)
plot(p10)
title ('10')
subplot(4,10,20)
plot(pabs10)
title ('10')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,21)
plot(p11)
title ('11')
subplot(4,10,31)
plot(pabs11)
title ('11')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,22)
plot(p12)
title ('12')
subplot(4,10,32)
plot(pabs12)
title ('12')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,23)
plot(p13)
title ('13')
subplot(4,10,33)
plot(pabs13)
title ('13')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,24)
plot(p14)
title ('14')
subplot(4,10,34)
plot(pabs14)
title ('14')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,25)
plot(p15)
title ('15')
subplot(4,10,35)
plot(pabs15)
title ('15')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,26)
plot(p16)
title ('16')
subplot(4,10,36)
plot(pabs16)
title ('16')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,27)
plot(p17)
title ('17')
subplot(4,10,37)
plot(pabs17)
title ('17')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,28)
plot(p18)
title ('18')
subplot(4,10,38)
plot(pabs18)
title ('18')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,29)
plot(p19)
title ('19')
subplot(4,10,39)
plot(pabs19)
title ('19')
axis([0 2000 0 3.5*10^4 ])

subplot(4,10,30)
plot(p20)
title ('20')
subplot(4,10,40)
plot(pabs20)
title ('20')
axis([0 2000 0 3.5*10^4 ])

sgtitle('druck änderung');

 


p.mat
 Beschreibung:

Download
 Dateiname:  p.mat
 Dateigröße:  35.86 KB
 Heruntergeladen:  192 mal


Harald
Forum-Meister

Forum-Meister


Beiträge: 24.495
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 19.11.2022, 18:00     Titel:
  Antworten mit Zitat      
Hallo,

zur ersten Frage: Variablen nicht durchnumerieren, sondern Matrizen oder Cell Arrays in einer for-Schleife verwenden.

Der erste Teil dürfte sein:
Code:
pCell = cell(1,20);
pabsCell = cell(1,20);
for k = 1:20
pCell{k} = p0(k*1000:k*1000+9000);
pabsCell{k} = abs(fft(pCell{k},n));
end


Bei den ersten Schritten kann auch das MATLAB Onramp recht hilfreich sein.
https://www.mathworks.com/learn/tutorials/matlab-onramp.html

Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
hedro1234

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 19.11.2022, 22:40     Titel:
  Antworten mit Zitat      
hallo Harald,

vielen dank für deine Hilfe und den Kurs in dem Link kann mir auf jeden fall helfen.

aber wie kann ich das plotten von alle Signale in einem Figur zeigen. ich meine zb. 1. Abschnitt von pCell und pabsCell untereinander erstellen. wie in dem Bild im Anhang

Liebe Grüße


Code:


clear
clc

load p.mat

n=2^14;
%%tiefpass
tf=ones(1,9)/9;


p= p.p;
p0= conv(p,tf,'same');


pCell = cell(1,20);
pabsCell = cell(1,20);
for k = 1:20
pCell{k} = p0(k*1000:k*1000+9000);
pabsCell{k} = abs(fft(pCell{k},n));
end


figure(1)
 hold on;
 cellfun(@plot,pCell);

figure(2)
 hold on;
 cellfun(@plot,pabsCell);

 


1.PNG
 Beschreibung:

Download
 Dateiname:  1.PNG
 Dateigröße:  146.36 KB
 Heruntergeladen:  165 mal
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.495
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 19.11.2022, 23:03     Titel:
  Antworten mit Zitat      
Hallo,

auch das plotten würde ich in einer for-Schleife machen, der Einfachheit halber eine für jede Zeile im Plot.

Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
hedro123

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 20.11.2022, 01:24     Titel:
  Antworten mit Zitat      
hallo,

ich habe es nur geschafft, bei verschiede Figure zu plotten



Code:



figure (1)
for c = 1:length(pCell)
subplot(2,10,c)
plot(pCell{c})
hold on
end

figure (2)
for c = 1:length(pabsCell)
subplot(2,10,c)
plot(pabsCell{c})
hold on
end


 
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.495
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 20.11.2022, 07:52     Titel:
  Antworten mit Zitat      
Hallo,

Code:
figure
for c = 1:length(pCell)
subplot(4,10,c)
plot(pCell{c})
subplot(4,10,c+20)
plot(pabsCell{c})
end


Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

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 | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2024 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.