summaryrefslogtreecommitdiff
path: root/openEMS/Analyse
diff options
context:
space:
mode:
authorRuben Undheim <ruben.undheim@gmail.com>2016-07-05 18:02:38 +0200
committerRuben Undheim <ruben.undheim@gmail.com>2016-07-05 18:02:38 +0200
commitef962f6008f25ab7cbd4ca21bcc72b97a1e2d76f (patch)
tree8149bee93d1a3f91d4503bfb3853adac4af0a85e /openEMS/Analyse
Imported Upstream version 0.0.34
Diffstat (limited to 'openEMS/Analyse')
-rw-r--r--openEMS/Analyse/PlotVoltage.m70
1 files changed, 70 insertions, 0 deletions
diff --git a/openEMS/Analyse/PlotVoltage.m b/openEMS/Analyse/PlotVoltage.m
new file mode 100644
index 0000000..19989c1
--- /dev/null
+++ b/openEMS/Analyse/PlotVoltage.m
@@ -0,0 +1,70 @@
+%close all;
+clear all;
+clc
+
+fmax = 50e6;
+
+figure(1);
+tmpu = load('../tmp/u1');
+tmpi = load('../tmp/i1');
+
+t = tmpu(:,1);
+u = tmpu(:,2);
+
+subplot(2,2,1);
+title('u_1 TD');
+plot(t,u);
+xlabel('t \rightarrow');
+ylabel('ut_1 \rightarrow');
+grid on;
+
+dt=t(2)-t(1);
+u= [u ; zeros(5000,1)];
+L=numel(u);
+t = (1:L)*dt;
+
+f = (0:L-1)/L/dt;
+fu = fft(u)/L;
+subplot(2,2,2);
+title('u_1 FD');
+plot(f(1:L/2),abs(fu(1:L/2)));
+xlabel('f \rightarrow');
+ylabel('|uf_1| \rightarrow');
+grid on;
+
+
+t = tmpi(:,1);
+i = tmpi(:,2);
+
+subplot(2,2,3);
+title('i_1 TD');
+plot(t,i);
+xlabel('t \rightarrow');
+ylabel('it_1 \rightarrow');
+grid on;
+
+dt=t(2)-t(1);
+i = [i; zeros(5000,1)];
+L=numel(i);
+t = (1:L)*dt;
+f = (0:L-1)/L/dt;
+
+fi = fft(i)/L;
+subplot(2,2,4);
+title('i_1 FD');
+plot(f(1:L/2),abs(fi(1:L/2)));
+xlabel('f \rightarrow');
+ylabel('|if_1| \rightarrow');
+grid on;
+
+figure(2);
+subplot(2,1,1);
+plot(f,real(fu./fi));
+xlim([0 fmax]);
+grid on;
+subplot(2,1,2);
+plot(f,imag(fu./fi));
+xlim([0 fmax]);
+grid on;
+
+