Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Bücher:

Digitale Signalverarbeitung

Fachkräfte:
weitere Angebote

Partner:


Vermarktungspartner


Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Teilentladungsmessung Ladungstrennung

 

Wattläufer
Forum-Newbie

Forum-Newbie


Beiträge: 1
Anmeldedatum: 28.02.21
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 28.02.2021, 18:47     Titel: Teilentladungsmessung Ladungstrennung
  Antworten mit Zitat      
Hallo zusammen,

im Rahmen einer Projektarbeit im Fachbereich Elektrotechnik beschäftige ich mich aktuell mit Teilentladungsmessungen. Diese sollten in MATLAB geplottet werden.

Zur Darstellung des Grundrauschen haben wir folgenden Coce verwendet:

Code:
%
clear all;
close all;
clc

load fisheriris

%Phasenlage der Spannung
file1=fopen('unit1.2.PH', 'r');
Phase=fread(file1, 'double');
fclose(file1);
Phase=Phase*360;

%Entladungen
file2=fopen('unit1.2.Q', 'r');
Charge=fread(file2, 'float', 8);
fclose(file2);

%array
file3=fopen('unit1.2.Q', 'r');
help=fread(file3, 1, 'float');

%charge value
time=fread(file2, 'double', 4);
fclose(file2);
Charge=(Charge*1*10^12);

%Spannungen
file4=fopen('unit1.2.V','r');
Voltage=fread(file4, 'float');
fclose(file4);
VoltageTime=0:48*10^-6:((length(Voltage)-1)*(48*10^-6));

%Diagramm
figure
yyaxis left
semilogy(Phase, abs(Charge),'.k','LineWidth',2,'MarkerSize',10) % plot discharge
ylabel('Entladungen Q / pC ');
hold on
plot(Phase,Charge,'.g') %I assume that Phase is in [0,360}
plot(0:360,max(abs(Charge))*sin((0:360)*pi/180)) %you must create a sinusoidal for references
set(gca,'xtick',[0 90 180 270 360]) %puts adequate tick marks
grid

%some suggested formatting (second y-axis can be formatted accordingly)
yyaxis left

hAx=gca; %get a handle to the current axis
hAx.YColor='black' % y axis color set to black
hAx.XLim=[0 360] % limits of the x-axis
hAx.YLim=[0.1 10^5] % limits of the left y-axis
hAx.FontSize=14;
hAx.XTick=[0 60 120 180 240 300 360];
hAx.XTickLabel=[0 60 120 180 240 300 360];
yyaxis right
hAx=gca;
hAx.YTick=[-200 0 200];
hAx.YLim=[-800 200];

% offset=30;
%
% index= knnsearch(Voltage- mean( Voltage) + offset,0);
% VoltageOnePeriodWindow= Voltage([index:index+416],1); %due to the 48 µs
% VoltagePhase=[ 0:360/417:360-(360/417)] ;

yyaxis right % right y axis of the figure
% VOP= smooth( VoltagePhase, VoltageOnePeriodWindow-mean(VoltageOnePeriodWindow) ,0.25, 'rloess');
%
%   plot(VoltagePhase , VOP, 'r','LineWidth',2); %plot voltage

ylabel('Rippel ΔU DC / V');
xlabel('Phasenwinkel φ/°');
%some suggested formatting (second y-axis can be formatted accordingly)
yyaxis left
hAx=gca; %get a handle to the current axis
hAx.YColor='black' % y axis color set to black
hAx.XLim=[0 360] % limits of the x-axis
hAx.YLim=[0.1 10^5] % limits of the left y-axis
hAx.FontSize=14;
hAx.XTick=[0 60 120 180 240 300 360];
hAx.XTickLabel=[0 60 120 180 240 300 360];
yyaxis right
hAx= gca;
hAx.YTick=[-200 0 200];
hAx.YLim=[-800 200];

a=size(Charge);

predictor=zeros(3867375,1);
Charge=abs(Charge);
 for i=1:3867375
 
if Charge (i)>0.4
 predictor(i,1)=1;
end    
 end
 


Das einzige was leider nicht klappt die Ladung in zwei seperate Grafiken darzustellen, d.h. das Grundrauschen trennen.

Ich hoffe ihr könnt gut Ansätze für das Problem nennen.

Viele Grüße

Wattläufer
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
.


goMatlab ist ein Teil des goForen-Labels
goForen.de goMATLAB.de goLaTeX.de


 Impressum  | Nutzungsbedingungen  | Datenschutz  | Werbung/Mediadaten | Studentenversion | FAQ | goMatlab RSS Button RSS


Copyright © 2007 - 2021 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
Partner: LabVIEWforum.de

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.