Teilentladungsmessung Ladungstrennung

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Wattläufer

Forum-Newbie


	Beiträge: 1

	Anmeldedatum: 28.02.21

	Wohnort: ---

	Version: ---

Verfasst am: 28.02.2021, 17:47 Titel: Teilentladungsmessung Ladungstrennung

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

Funktion ohne Link?

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


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 |

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.