Ejercicios de mediciones eléctricas con M AT LAB r &

Simulink r
Versión 2020

Mediciones Eléctricas 1

Dr. Ing. Jorge Luis Strack - Ing. Juan Francisco Martinez

1
 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Índice
Objetivos 3

Ejercicio 1 3
Programación de un script en MATLAB . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3

Ejercicio 2 4
Gráficos paramétricos - Interacción entre SIMULINK y MATLAB . . . . . . . . . . . . . . . . . . 4

Ejercicio 3 7
Medición de potencia monofásica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

Ejercicio 4 14
Respuesta de los instrumentos a señales de distintas formas de onda. . . . . . . . . . . . . . . . . . 14

Ejercicio adicional 1 18
Cálculo de propagación de incertezas por GUM (magnitudes de entrada no correlacionadas) . . . . 18

Ejercicio adicional 2 20
Análisis de una sonda conectada a un osciloscopio . . . . . . . . . . . . . . . . . . . . . . . . . . . 20

Bibliografı́a recomendada 23

Mediciones Eléctricas 1 Página 2 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Objetivos
Crear modelos de simulación y resolver problemas prácticos.

Afianzar conceptos de electrotecnia y mediciones eléctricas.

Identificar posibles errores o limitaciones en las simulaciones.

Ejercicio 1
Programación de un script en MATLAB
Se realizará un script que determine los errores relativos de método corto y largo, y grafique los mismos en
función de la resistencia medida. Para la resolución se proponen dos metodologı́as. La primera consiste en
realizar operaciones con variables simbólicas y luego realizar sustituciones numéricas. La segunda consiste en
operar directamente con variables numéricas. A partir de cualquiera de las dos metodologı́as se puede obtener
el mismo gráfico. La diferencia radica en que al operar con variables simbólicas se pueden determinar las
expresiones analı́ticas, mientras que para operar con variables numéricas, ya se deben conocer previamente
las expresiones analı́ticas o algebraicas.

% Solucion con variables simbolicas

close a l l ; clear a l l ; c l c ;
% declaracion de variables simbolicas
syms Ra Rv Rm R
5 % error de metodo corto
R = 1/(1/Rm - 1/Rv)
emc = (Rm-R)/R
simplify(emc)
pretty(simplify(emc))
10 % error de metodo largo
R = Rm - Ra
eml = (Rm-R)/R
simplify(eml)
pretty(simplify(eml))
15 % Evaluacion numerica de las funciones analiticas
emc = -subs(emc,{Rv,Rm},{10e3,[1:1:300]});
eml = subs(eml,{Ra,Rm},{0.25,[1:1:300]});
Rm = [1:1:300];
%% Solucion con variables numericas
20 close a l l ; clear a l l ; c l c ;
% declaracion de variables numericas
Rv = 10e3;
Ra = 0.25;
Rm = [1:1:300];
25 emc = Rm/Rv;
eml = Ra./Rm;
%% Grafico
figure
hold on
30 plot(Rm,emc,'b');
plot(Rm,eml,'r');
grid on
xlim([0 100])

Ejercicio 1 continúa en próxima página. . . Mediciones Eléctricas 1 Página 3 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

ylim([0 0.02])
35 t i t l e ('Errores relativos de metodo corto y largo')
xlabel('Resistencia medida [Ohms]')
ylabel('Error relativo')
legend('Metodo corto','Metodo largo')
hold off

Sintaxis de la función subs (sustitución simbólica)

subs(y, {x1 x2 x3 . . . xn}, {a b c . . . n})

donde:

y : función matemática previamente definida

{x1 x3 x3 . . . xn} : vector de variables simbólicas de la función matemática “y” que se

desea especializar en valores numéricos

{a b c d . . . n} : vector de valores numéricos en los que se van a especializar a las

variables simbólicas x1 x3 x3 . . . xn de la función matemática “y”

Ejercicio 2
Gráficos paramétricos - Interacción entre SIMULINK y MATLAB
En éste ejercicio se propone una forma de realizar un análisis del comportamiento de un circuito a partir
de la variación de uno de sus parámetros. A fin de ilustrar la interación entre MATLAB y SIMULINK
y la facilidad con la que se pueden obtener curvas paramétricas, se utilizará un ejemplo bien conocido. Se
obtendrá la curva de regulación de un sistema paralelo, obteniendo la tensión de salida del sistema en función
de la posición del cursor n para diversos factores de carga m.
A continuación se presenta el modelo circuital de un sistema de regulación paralelo que proveé de una tensión
continua variable a dos voltı́metros durante un ensayo de contraste:

Figura 1: Sistema de regulación paralelo para contraste de un voltı́metro.

Donde:

Ejercicio 2 continúa en próxima página. . . Mediciones Eléctricas 1 Página 4 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Tabla 1: Datos del sistema de regulación.

E Tensión de la fuente de corriente continua 30V

ri Resistencia interna de la fuente de corriente continua 1mΩ
R Resistencia de regulación 3800Ω
Rp Resistencia interna del voltı́metro patrón en el alcance de 30 V 10kΩ
Rx Resistencia interna del voltı́metro a contrastar 28,4kΩ
Rc Resistencia de carga Rp //Rx
n Parámetro adimensional que indica la porción de R que queda en paralelo con Rc
m Factor de carga. Es el cociente entre R y Rc

A partir de estos datos, se procede a construir el modelo circuital en SIMULINK con el auxilio de bloques
y funciones de la librerı́a SimPowerSystems.

Figura 2: Modelo de un sistema de regulación paralelo para contraste de un voltı́metro.

Se han utilizado los siguientes bloques:

Configuración para bloques de la librerı́a SimPowerSystems: Library/SimPowerSystems/powergui

Fuente de tensión continua: Library/SimPowerSystems/Electrical Sources/DC Voltage Source

Resistencias: Library/SimpowerSistems/Elements/Series RLC Branch

Referencia (GND): Library/SimpowerSistems/Elements/Ground

Voltı́metro: Library/SimPowerSystems/Measurements/Voltage Measurement

Bloque To Workspace: Library/Simulink/Sinks/To Workspace

Se sugiere guardar el archivo del modelo SIMULINK con el nombre regulacion paralelo.mdl.

Como ejemplo de interacción entre SIMULINK Y MATLAB se escribirá un script MATLAB que realice
reiteradas simulaciones del modelo SIMULINK, modificando en cada una de ellas la posición del cursor n
y el factor de carga m en el circuito, y almacenando los datos obtenidos en una variable para graficarlas al
finalizar el proceso.

Ejercicio 2 continúa en próxima página. . . Mediciones Eléctricas 1 Página 5 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

%% Primera etapa: configuracion del sistema de regulacion paralelo para un

% valor de n

close a l l ; clear a l l ; c l c ;
5

E = 30; % [V]
Ri = 1e-3; % [Ohms]
R = 3800; % [Ohms]
Rp = 10000; % [Ohms]
10 Rx = 28000; % [Ohms]
n = 0.999999;

%% Segunda etapa: configuracion del sistema de regulacion paralelo para ...

distintos valores de n

15 close a l l ; clear a l l ; c l c ;

E = 30; % [V]
Ri = 1e-3; % [Ohms]
R = 3800; % [Ohms]
20 Rp = 10000; % [Ohms]
Rx = 28400; % [Ohms]
vector_n = [0:0.05:1]; % Vector con distintas posiciones del cursor

open_system('regulacion_paralelo')
25

for i=1: s i z e (vector_n,2)

n = vector_n(i);
sim('regulacion_paralelo');
salida_v(i)=salida_v_i.signals.values (end) ;
30 clc;
mensaje = ['Iteracion ' num2str(i)];
disp(mensaje);
end

35 % Curva de regulacion
figure
plot(vector_n,salida_v,'-or')
t i t l e ('Curva de regulacion de un sistema paralelo')
xlabel('Posicion del cursor: n')
40 ylabel('Tension de salida')

%% Tercera etapa: configuracion del sistema de regulacion paralelo para ...

distintos valores de n y de m
close a l l ; clear a l l ; c l c ;
E = 30; % [V]
45 Ri = 1e-3; % [Ohms]
%R = 3800; % [Ohms]
Rp = 10000; % [Ohms]
Rx = 28000; % [Ohms]
vector_n = [0:0.1:1]; % Vector con distintas posiciones del cursor
50 vector_m = [0.25 0.5 1:1:10]; % Vector con distintos factores de carga
Rc = Rp*Rx/(Rp+Rx);

Ejercicio 2 continúa en próxima página. . . Mediciones Eléctricas 1 Página 6 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

open_system('regulacion_paralelo')

55 for j = 1: s i z e (vector_m,2)
m = vector_m(j);
R = m*Rc;
for i = 1: s i z e (vector_n,2)
n = vector_n(i);
60 sim('regulacion_paralelo');
salida_v(i,j)=salida_v_i.signals.values (end) ;
clc;
mensaje = ['Iteracion ' num2str(i)];
disp(mensaje);
65 end
end

% Curva de regulacion
colores(1,:)=[1 0 0];
70 colores(2,:)=[0 1 0];
colores(3,:)=[0 0 1];
colores(4,:)=[0.85 0 0];
colores(5,:)=[0 0.85 0];
colores(6,:)=[0 0 0.85];
75 colores(7,:)=[0.7 0 0];
colores(8,:)=[0 0.7 0];
colores(9,:)=[0 0 0.7];
colores(10,:)=[0.55 0 0];
colores(11,:)=[0 0.55 0];
80 colores(12,:)=[0 0 0.55];
figure
hold on
for j = 1: s i z e (vector_m,2)
plot(vector_n,salida_v(:,j),'-o','color',colores(j,:))
85 end
hold off
t i t l e ('Curva de regulacion de un sistema paralelo')
xlabel('Posicion del cursor: n')
ylabel('Tension de salida')
90 legend('m=0.25','m=0.5','m=1','m=2','m=3','m=4','m=5','m=6','m=7','m=8','m=9','m=10')

Ejercicio 3
Medición de potencia monofásica
Se modelarán los tres circuitos (R, RL y RLC) ensayados en la respectiva práctica de laboratorio, utilizando
los valores medidos en el ensayo. Se observarán las señales temporales de tensión, corriente y potencia y se
medirá la potencia activa, reactiva y aparente en cada uno de ellos utilizando el modo de simulación continua.
Luego, utilizando el modo de simulación fasorial, y con el auxilio de código en MATLAB, se procederá a
realizar el diagrama fasorial de cada circuito. Por último, en el circuito RLC se realizará una medición de
impedancia, a fin de observar que el mismo resuena a 50 Hz.

a) Análisis temporal: circuito monofásico resistivo (R)

Ejercicio 3 continúa en próxima página. . . Mediciones Eléctricas 1 Página 7 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Abra MATLAB y a continuación Simulink. Abra un modelo nuevo y guárdelo con un nombre que lo identi-
fique, por ejemplo R.mdl. Agregue al modelo los siguientes bloques a fin de obtener el circuito de la Fig. 3:
Powergui: Library/SimPowerSystems/powergui
Fuente de tensión CA: Library/SimPowerSystems/Electrical Sources/AC Voltage Source
Carga RLC configurable: Library/SimPowerSystems/Elements/Parallel RLC Load
Tierra: Library/SimpowerSistems/Elements/Ground
Voltı́metro: Library/SimPowerSystems/Measurements/Voltage Measurement
Amperı́metro: Library/SimPowerSystems/Measurements/Current Measurement
Gráfico: Library/Simulink/Sinks/Scope
Indicador: Library/Simulink/Sinks/Display
Multiplicador: Library/Simulink/Commonly Used Blocks/Product
Valor medio: Library/SimPowerSystems/Extra Library/Measurements/Mean Value
Valor RMS: Library/SimPowerSystems/Extra Library/Measurements/RMS
Creador de bus: Library/Simulink/Commonly Used Blocks/Bus Creator
Función: Library/Simulink/User-Defined Functions/Fcn

Scope

+ i
−

V1 A

v + AC Voltage Source
− Pins
Mean
Parallel RLC Load
P
RMS

Vrms

RMS
Continuous S
Irms
powergui
P [W]

Vrms [V]

sqrt(u(1)^2 − u(2)^2)
Irms [A]
Bus Fcn
Creator
S [VA]

Q[VAr]

Figura 3: Modelo de circuito monofásico con carga resistiva.

Configure los bloques de acuerdo a las siguientes indicaciones:

Ejercicio 3 continúa en próxima página. . . Mediciones Eléctricas 1 Página 8 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Powergui: En configuración de parámetros se selecciona tipo de simulación continua.

Fuente de tensión CA: Configurar amplitud (Vrms[V]*sqrt(2)) y frecuencia (50 [Hz]).

Carga RLC configurable: Configure el valor de potencia activa y ponga cero en las potencias reactivas.

Gráfico: Configurar el número de ejes (7) y deshabilitar el lı́mite de datos en el historial.

Indicador: Configurar el formato que crea más conveniente. Establecer el factor de diezmado en 1.

Multiplicador: Por defecto el nº de entradas es igual a dos, pero se puede cambiar según las necesidades.

Valor medio: Ingresar el perı́odo de la señal a medir (1/50 [s]).

Valor RMS: Ingresar la frecuencia de la señal a medir (50 [Hz]).

Creador de bus: También es posible cambiar el nº de entradas, por defecto es igual a dos.

Función: Aquı́ se debe cargar la siguiente ecuación: sqrt(u(1)2 − u(2)2 ) a fin de calcular la potencia
reactiva a partir de la potencia aparente, u(1), y de la potencia activa, u(2). Los ı́ndices (1 y 2) se
establecen de acuerdo al orden en que ingresan las señales al bus de datos.

Por último configure los parámetros generales de la simulación:

Simulation time/Start time: 0 (tiempo inicial)

Simulation time/Stop time: 0.08 (tiempo final)

Solver options/Type: Variable-step (paso variable entre cada iteración)

Solver options/Solver: ode23tb (solver adecuado para este tipo de simulación)

Solver options/Max step size: 1/50/156 (paso máximo de tiempo entre cada iteración)

Solver options/Relative tolerance: 1e-4 (error relativo máximo en cada iteración)

Solver options/Solver reset method: Robust (logra mayor estabilidad en la simulación)

Zero-crossing control/Algorithm: Adaptive (logra mayor estabilidad en la simulación)

Ejecute la simulación, observe las señales graficadas y lea los cuadros indicadores del circuito. Complete la
primera fila de la Tabla 2 y obtenga conclusiones.

Tabla 2: Resumen de los valores medidos en la simulación de los circuitos resistivo, resistivo-inductivo y
resistivo-inductivo con corrección del factor de potencia.

Tipo de Valores medidos

circuito Vrms [V] Irms [A] fv(t) [Hz] fp(t) [Hz] Scarga [VA] Pcarga [W] Qcarga [VAr]
R
RL
RLC

b) Análisis temporal: circuito monofásico resistivo inductivo (RL)

Abra el modelo del ejercicio anterior, y guárdelo con otro nombre, por ejemplo RL.mdl. A fin de apren-
der a utilizar nuevas funciones, agregue al modelo los siguientes bloques para obtener el circuito de la Fig. 4.

Ejercicio 3 continúa en próxima página. . . Mediciones Eléctricas 1 Página 9 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Medidor de P y Q: Library/SimPowerSystems/Extra Library/Measurements/Active & Reactive Power

Multı́metro: Library/SimPowerSystems/Measurements/Multimeter

Demux: Library/Simulink/Commonly Used Blocks/Demux

Ganancia: Library/Simulink/User-Defined Functions/Gain

Mux: Library/Simulink/Commonly Used Blocks/Mux

Configure los bloques de acuerdo a las siguientes indicaciones:

Medidor de P y Q: Ingresar la frecuencia de la señal a medir (50 [Hz]).

Fuente de tensión CA: Setear el campo Measurements, seleccionando Voltage.

Carga RLC configurable: Cargar la potencia reactiva inductiva y setear el campo Measurements, selec-
cionando Branch Current.

Multı́metro: Seleccionar la tensión de la fuente y la corriente en la carga.

Ganancia: Configurar un valor adecuado para amplificar la corriente y graficarla junto a la tensión.

Gráfico: Copiar un nuevo gráfico de un eje.

Indicador: Copiar un nuevo display.

Scope

+ i
−

V1 A

v + AC Voltage Source
− Pins
Mean
Parallel RLC Load
P
RMS
Vrms

RMS
Continuous S
Irms
powergui
P [W]

Vrms [V]

sqrt(u(1)^2 − u(2)^2)
2 Irms [A]
100 Bus Fcn
MultimeterDemux Mux Scope1 Creator
Gain
S [VA]
V
PQ
I

Active & Reactive P[W] Q[VAr]

Power Q[VAr]

Figura 4: Modelo de circuito monofásico con carga resistiva.

Ejercicio 3 continúa en próxima página. . . Mediciones Eléctricas 1 Página 10 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Observe las señales graficadas y lea los cuadros indicadores del circuito. Complete la segunda fila de la
Tabla 2 y obtenga conclusiones.

c) Análisis temporal: circuito monofásico resistivo inductivo con corrección del factor de po-
tencia (RLC).

Abra el modelo del ejercicio anterior, y guárdelo con otro nombre, por ejemplo RLC.mdl. Configure los
bloques de acuerdo a las siguientes indicaciones, a fin de modelar la carga capacitiva y obtener el circuito de
la Fig. 5.

Carga RLC configurable: Cargar la potencia reactiva capacitiva.

Ejecute la simulación, observe las señales graficadas y lea los cuadros indicadores del circuito. Complete la
tercera fila de la Tabla 2 y obtenga conclusiones.

Scope

+ i
−

V1 A
Series RLC Branch
v + AC Voltage Source
− Pins
Mean
Parallel RLC Load
P
RMS

Vrms

RMS
Continuous S
Irms
powergui
P [W]

Vrms [V]

sqrt(u(1)^2 − u(2)^2)
2 Irms [A]
100 Bus Fcn
MultimeterDemux Mux Scope1 Creator
Gain
S [VA]
V
PQ
I

Active & Reactive P[W] Q[VAr]

Power Q[VAr]

Figura 5: Modelo de circuito monofásico con carga resistiva.

d) Análisis fasorial.

A fin de obtener el diagrama fasorial de tensión y corriente de cada circuito, en cada uno de ellos, se
hará una copia y se configurará el bloque powergui para realizar una simulación en el modo phasor a la
frecuencia especı́fica de 50[Hz]. Se quitarán todos los bloques que se usaron en el análisis temporal para
graficar las señales e indicar los valores medidos, y se reemplazaran por los siguientes bloques:

Mux: Library/Simulink/Commonly Used Blocks/Mux

To File: Library/Simulink/Sinks/To File

En la Fig. 6 se observa como deberı́a quedar el circuito RL a modo de ejemplo.

Ejercicio 3 continúa en próxima página. . . Mediciones Eléctricas 1 Página 11 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

datos.mat
+ i
−
Mux To File
V1 A

v + AC Voltage Source
−
Parallel RLC Load Phasors

powergui

Figura 6: Análisis fasorial del modelo RL.

Donde en el bloque To File se configura el nombre del archivo de extensión ”.mat”donde se guardarán los
datos de los fasores de tensión y corriente; y se asignará un nombre a la estructura que contendrá estos datos.
A continuación se presenta un código sencillo de MATLAB conocido comunmente como script de extensión
.m, que importa los datos del archivo de extensión .mat , calcula las potencias P, Q y S, y grafica el diagrama
fasorial de tensión y corriente.

%close all; clear all; clc;

disp('comenzando calculos')
%////////////////////////////////////////
% Carga de datos provenientes de Simulink
5 load datos
Urms = datos.Data(end,1)/sqrt(2);
Uphi = datos.Data(end,2);
Irms = datos.Data(end,3)/sqrt(2);
Iphi = datos.Data(end,4);
10 %////////////////////////////////////////
% Calculo de potencias
S = Urms * Irms;
phi = (Uphi - Iphi)* pi/180;
P = S * cos(phi);
15 Q = S * sin(phi);
FP = cos(phi);
%////////////////////////////////////////
% Calculos para diagrama fasorial
U = Urms*exp(1i*Uphi* pi/180);
20 Ure = real(U);
Uim = imag(U);
I = Irms*exp(1i*Iphi* pi/180);
Ire = real(I);
Iim = imag(I);
25 %////////////////////////////////////////
% Diagrama fasorial
figura = figure(1);
Uesc = 1;
Iesc = 100;
30 quiver(0,0,Ure,Uim,Uesc,'color',[0 0 1],'LineWidth',2); hold on
quiver(0,0,Ire,Iim,Iesc,'color',[1 0 0],'LineWidth',2)
axis equal
xmin = -50; xmax = 500; xlim([xmin,xmax]);
ymin = -500; ymax = 50; ylim([ymin,ymax]);

Ejercicio 3 continúa en próxima página. . . Mediciones Eléctricas 1 Página 12 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

35 grid on
t i t l e ('Diagrama fasorial de tension y corriente');
xlabel('Eje real')
ylabel('Eje imaginario')
legend('U[V]','I[A] x 100','Location','Best')
40 text(Uesc*Ure,Uesc*Uim,'U','color',[0 0 1],'FontSize',20);
text(Iesc*Ire,Iesc*Iim,'I','color',[1 0 0],'FontSize',20);
disp('Calculo terminado')

%////////////////////////////////////////
45 % Exportacion del grafico
% saveas(figura,'diagrama_fasorial.pdf');
% saveas(figura,'diagrama_fasorial.eps','epsc2');
% saveas(figura,'diagrama_fasorial.png');

Ejecute la simulación y el script, observe el diagrama fasorial obtenido y saque conclusiones. ¿Será posible
ejecutar el script en forma automática una vez simulado el modelo?

e) Análisis de impedancia del circuito RLC.

A continuación se medirá la impedancia del circuito en función de la frecuencia. Para realizar este análisis,
agregue al modelo utilizado en el análisis fasorial, el siguiente bloque:

Medición de impedancia: Library/SimPowerSystems/Measurements/Impedance Measurement

Luego realice el cableado correspondiente para obtener el circuito de la Fig. 7.

R1

AC

RLC Load Z
Impedance
Measurement
Continuous

powergui

Figura 7: Análisis frecuencial de la impedancia del modelo RLC.

Ingrese al menu del bloque powergui y realice un barrido en frecuencia de la impedancia del circuito. Obtenga
conclusiones respecto al gráfico observado. Exporte los datos del gráfico a MATLAB e intente reproducir el
gráfico a través de un script.

Mediciones Eléctricas 1 Página 13 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Ejercicio 4
Respuesta de los instrumentos a señales de distintas formas de onda.
El propósito de este ejercicio es modelar voltı́metros de hierro móvil, IPBM e IPBM con rectificador de media
onda, a fin de analizar la respuesta de cada uno de ellos ante una señal de entrada senoidal, cuadrada y
trapezoidal, sin rectificador, con rectificador de media onda y con rectificador de onda completa. El objetivo
es verificar lo visto en la teorı́a y en el laboratorio.

a) Medición de señales sin rectificador

En el entorno de SIMULINK, se procede a realizar el modelo de la Fig. 8. Para realizar el mismo se deben
cargar los siguientes bloques, de los cuáles se detalla su ubicación en las respectivas librerı́as:

Powergui: Library/SimPowerSystems/powergui

Generador de señales: Library/Simulink/Sources/Signal Generator

Fuente de tensión controlada: Library/SimPowerSystems/Electrical Sources/Controlled Voltage Source

Resistencia: Library/SimpowerSistems/Elements/Series RLC Branch

Tierra: Library/SimpowerSistems/Elements/Ground

Voltı́metro: Library/SimPowerSystems/Measurements/Voltage Measurement

Amperı́metro: Library/SimPowerSystems/Measurements/Current Measurement

Gráfico: Library/Simulink/Sinks/Scope

Indicador: Library/Simulink/Sinks/Display

Subsistema: Library/Simulink/Commonly Used Blocks/Subsystem

+
− v

Continuous
+ i
−
powergui A
+

Controlled Voltage Source 1 + 1

Scope
Display Display Display
R
−
s

2 − 2
1
Voltimetro HM2 Voltimetro IPBM DC Voltimetro IPBM AC
Display
2
Signal
Generator Voltimetro HM1
Display1 Display2 Display3 Display4

Figura 8: Modelo para medición de señales sin rectificador.

Para modelar los tres tipos de instrumentos se deben crear tres subsistemas, cuyos circuitos internos se
detallan a continuación:

Ejercicio 4 continúa en próxima página. . . Mediciones Eléctricas 1 Página 14 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

a.1) Modelado de voltı́metro de hierro móvil

1
1

+
RMS −K− 1
Ra − v
Display
V RMS Gain

Rm

2
2

Figura 9: Modelo simplificado de un voltı́metro de hierro móvil.

a.2) Modelado de voltı́metro IPBM DC y AC

1
1
+
1

+
Diode
Mean −K− 1 +
Ra − v Mean −K− 1
Display Ra − v
V Mean Value Gain Display
V Mean Value Gain

Rm
Rm

2
2
−
2

Figura 10: Modelo simplificado de un IPBM en función DC (izquierda) y AC (derecha).

Guarde estos tres bloques en una librerı́a personal abriendo el menu file/New/Library.

Luego de armar el modelo, se procede a configurar cada uno de los bloques, teniendo en cuenta las siguientes
indicaciones:

Powergui: En configuración de parámetros se selecciona tipo de simulación continua.

Generador de señales: Configurar amplitud (50*sqrt(2)[V]), frecuencia (50[Hz]) y forma de onda (senoidal,
rectangular, ó triangular).

Fuente de tensión: No seleccionar la opción de inicializar parámetros.

Resistencia: Cargar un valor elevado (por ej. 10e3 [Ohm]) para modelar la resistencia interna de un
voltı́metro.

Gráfico: Configurar el número de ejes (6) y deshabilitar el lı́mite de datos en el historial.

Indicador: Configurar el formato que crea más conveniente. Establecer el factor de diezmado en 1.

Por último se configuran los parámetros generales de la simulación:

Simulation time/Start time: 0 (tiempo inicial)

Simulation time/Stop time: 0.08 (tiempo final)

Ejercicio 4 continúa en próxima página. . . Mediciones Eléctricas 1 Página 15 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Solver options/Type: Variable-step (paso variable entre cada iteración)

Solver options/Solver: ode23tb (solver adecuado para este tipo de simulación)

Solver options/Max step size: 1/50/156 (paso máximo de tiempo entre cada iteración)

Solver options/Relative tolerance: 1e-4 (error relativo máximo en cada iteración)

Solver options/Solver reset method: Robust (logra mayor estabilidad en la simulación)

Zero-crossing control/Algorithm: Adaptive (logra mayor estabilidad en la simulación)

Finalmente, ejecute la simulación, observe las señales graficadas y lea los cuadros indicadores del circuito.
Complete la Tabla 3 y saque conclusiones.

Tabla 3: Resumen de los valores medidos en la simulación.

Tipo Señal sin rectificar

de señal VHM [V] VIPBM−DC [V] VIPBM−AC [V] VPICO [V] VRMS [V] VMEDIO [V]
Senoidal
Rectangular
Triangular

b) Medición de señales con rectificador de media onda

A fin de analizar la respuesta de un rectificador de onda completa, en el modelo del ejercicio anterior,
se debe reemplazar el bloque del diodo por el siguiente bloque, obteniendo el modelo de la Fig. 11.

Puente universal: Library/SimPowerSystems/Power Electronics/Universal Bridge

Luego de armar el modelo, se procede a configurar el bloque recién visto:

Puente universal: Configurar 2 ramas y seleccionar el dispositivo diodo. Modelar un puente ideal y luego
uno real.

Los parámetros generales de la simulación son similares a los de la simulación anterior.

Ejecute la simulación, observe las señales graficadas y lea los cuadros indicadores del circuito. Complete la
Tabla 4 y obtenga conclusiones.

+
− v

Continuous
+ i
−
powergui A
+

Diode
Controlled Voltage Source 1 + 1
Scope
Display Display Display
R
−
s

2 − 2
1
Voltimetro HM2 Voltimetro IPBM DC Voltimetro IPBM AC
Display
2
Signal
Generator Voltimetro HM1
Display1 Display2 Display3 Display4

Figura 11: Modelo de circuito de corriente alterna con rectificador de onda completa.

Ejercicio 4 continúa en próxima página. . . Mediciones Eléctricas 1 Página 16 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Tabla 4: Resumen de los valores medidos en la simulación.

Tipo Señal rectificada en media onda

de señal VHM [V] VIPBM−DC [V] VIPBM−AC [V] VPICO [V] VRMS [V] VMEDIO [V]
Senoidal
Rectangular
Triangular

c) Medición de señales con rectificador de onda completa

A fin de analizar la respuesta de un rectificador de onda completa, en el modelo del ejercicio anterior,
se debe reemplazar el bloque del diodo por el siguiente bloque, obteniendo el modelo de la Fig. 12.

Puente universal: Library/SimPowerSystems/Power Electronics/Universal Bridge

Luego de armar el modelo, se procede a configurar el bloque recién visto:

Puente universal: Configurar 2 ramas y seleccionar el dispositivo diodo. Modelar un puente ideal y luego
uno real.

Los parámetros generales de la simulación son similares a los de la simulación anterior.

Ejecute la simulación, observe las señales graficadas y lea los cuadros indicadores del circuito. Complete la
Tabla 5 y obtenga conclusiones.

+
− v

Continuous A +
+ i
−
powergui A
+

B −
Controlled Voltage Source 1 + 1
Scope
Display Display Display
Universal Bridge R
−
s

2 − 2
1
Voltimetro HM2 Voltimetro IPBM DC Voltimetro IPBM AC
Display
2
Signal
Generator Voltimetro HM1
Display1 Display2 Display3 Display4

Figura 12: Modelo de circuito de corriente alterna con rectificador de onda completa.

Tabla 5: Resumen de los valores medidos en la simulación.

Tipo Señal rectificada en onda completa

de señal VHM [V] VIPBM−DC [V] VIPBM−AC [V] VPICO [V] VRMS [V] VMEDIO [V]
Senoidal
Rectangular
Triangular

Mediciones Eléctricas 1 Página 17 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Ejercicio adicional 1
Cálculo de propagación de incertezas por GUM (magnitudes de entrada no co-
rrelacionadas)
Se ha estudiado que de acuerdo al método propuesto por la GUM (Guide to the Expression of Uncertainty in
Measurements) y la Norma IRAM 35050:2001, cuando no existe correlación entre las magnitudes de entrada,
la incertidumbre asociada a un mensurando y denominada incertidumbre combinada se calcula como:

Donde:

xi : estimación de entrada i-ésima

u(xi ) : incerteza asociada a la estimación de entrada xi previamente calculada

ci = ∂f (X1 ,X 2 ,...,XN )
∂Xi |X1 =x1 ,X2 =x2 ,...,XN =xN derivada parcial i-ésima de la función mensurando respecto
a la variable de entrada i-esima (referida también como coeficiente de sensibilidad )

En casos en que el mensurando sea función de varias variables y el cálculo manual de las derivadas parciales
se vuelva complejo, puede resultar conveniente el empleo de un programa de aplicación matemática como
MATLAB. A continuación se propone un método para resolver dicho cálculo a partir de un script.

Se tiene una función mensurando con la que se estima en forma indirecta el valor de capacitancia de un
capacitor cerámico, a partir de mediciones de tensión (V1 y V2 ), frecuencia (f ), resistencia (Rref ) y desfasaje
angular entre las tensiones (θ):

Conocidos los valores de incerteza asociada a cada estimación de magnitud de entrada y considerando que
no existe correlación en las mediciones de dichas magnitudes, se puede aplicar el siguiente script MATLAB
para la resolución de los cálculos:

% Mediciones Electricas - 2017

% Script para obtencion de incerteza combinada - Funcion de 5 variables
% Y=Cx
%----------------------------------------------------------------------
5 clear a l l % Borra todas las variables del espacio de trabajo
c l c % limpia la ventana de comandos
% Inicializacion de las variables simbolicas
syms x1 x2 x3 x4 x5
% Carga de valores medidos en variables numericas
10 x1m=1.5; %[Volts pico]
x2m=1.35; %[Volts pico]
x3m=-19.8; %[grados]
x4m=1000; %[Hz]
x5m=296; %[ohm]

Mediciones Eléctricas 1 Página 18 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

15 % Carga de incertezas individuales

u1=0.03; %[Volts]
u2=0.03; %[Volts]
u3=0.5; %[grados]
u4=18.25; %[Hz]
20 u5=2.53; %[ohm]
% Conversion de angulos para calculos
x3m=x3m* pi/180; % grados a radianes
u3=u3* pi/180; % grados a radianes
% Carga de funcion
25 y=-(x1ˆ2+x2ˆ2-2*x1*x2* cos(x3))/(x5*2* pi *x4*x1*x2* sin(x3));
% Calculo de mejor estimacion (especializando la funcion simbolica 'y'
ym=vpa(subs(y,[x1 x2 x3 x4 x5],[x1m x2m x3m x4m x5m]),5)
% Calculo de coeficientes de sensibilidad
c1=vpa(subs( d i f f (y,x1),[x1 x2 x3 x4 x5],[x1m x2m x3m x4m x5m]),5);
30 c2=vpa(subs( d i f f (y,x2),[x1 x2 x3 x4 x5],[x1m x2m x3m x4m x5m]),5);
c3=vpa(subs( d i f f (y,x3),[x1 x2 x3 x4 x5],[x1m x2m x3m x4m x5m]),5);
c4=vpa(subs( d i f f (y,x4),[x1 x2 x3 x4 x5],[x1m x2m x3m x4m x5m]),5);
c5=vpa(subs( d i f f (y,x5),[x1 x2 x3 x4 x5],[x1m x2m x3m x4m x5m]),5);
% Calculo de contribuciones a la incerteza combinada
35 c1u1=vpa(c1*u1,5)
c2u2=vpa(c2*u2,5)
c3u3=vpa(c3*u3,5)
c4u4=vpa(c4*u4,5)
c5u5=vpa(c5*u5,5)
40 % Calculo de incerteza combinada
u_y=vpa(sqrt(c1u1ˆ2+c2u2ˆ2+c3u3ˆ2+c4u4ˆ2+c5u5ˆ2),5);
% Calculo de la incerteza relativa porcentual a la mejor estimacion
er_y_porcentual=vpa((u_y/ym)*100,3)

Sintaxis de la función subs (sustitución simbólica)

subs(y, {x1 x2 x3 . . . xn}, {a b c . . . n})

donde:

y : función matemática previamente definida

{x1 x3 x3 . . . xn} : vector de variables simbólicas de la función matemática “y” que se

desea especializar en valores numéricos

{a b c d . . . n} : vector de valores numéricos en los que se van a especializar a las

variables simbólicas x1 x3 x3 . . . xn de la función matemática “y”

Mediciones Eléctricas 1 Página 19 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Sintaxis de la función vpa

b = vpa(a, d)

donde:

b : número a expresado en formato decimal, con la cantidad de dı́gitos especificada por

el argumento d de la función vpa

a : número original

d : cantidad de dı́gitos con la que se expresará el número d

Ejercicio adicional 2
Análisis de una sonda conectada a un osciloscopio
En este ejercicio se propone una forma de realizar un análisis del comportamiento de un circuito a partir de
la variación de uno de sus parámetros (por ejemplo, el comportamiento de una tensión en función del tiempo
a medida que se varı́a el valor de una capacitancia en el circuito). Para ello, mediante el uso de SIMULINK
y un script MATLAB, se trazarán en este ejemplo una serie de curvas, cada una de ellas con un valor de
parámetro diferente. En este caso, las curvas corresponderán a la tensión a la salida de una sonda pasiva
atenuadora de osciloscopio en función del tiempo, y el parámetro será el valor de la capacitancia de punta
regulable, asociada a dicha sonda.
A continuación se presenta el modelo circuital de una sonda pasiva de tensión atenuadora 10X de un osci-
loscopio:

Figura 13: Modelo de una sonda pasiva de tensión atenuadora 10X conectada a un osciloscopio.

Donde:

Mediciones Eléctricas 1 Página 20 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

Tabla 6: Datos de la sonda.

Vi Tensión a la entrada de la sonda

Vo Tensión a la salida de la sonda (entrada del osciloscopio)
Vs Tensión del generador de señales 1V pico
Rs Resistencia de salida del generador de señales 50Ω
Cp Capacitancia ajustable en la punta de la sonda (parámetro) (11.3...15.3) ∆Cp = 1pF
Rp Resistencia en la punta de la sonda 9M Ω
Co Capacitancia de entrada del osciloscopio 20pF
Ro Resistencia de entrada del osciloscopio 1M Ω
Cc Capacitancia del cable coaxil de la sonda 100pF

A partir de estos datos, se procede a construir el modelo circuital en SIMULINK con el auxilio de bloques
y funciones de la librerı́a SimPowerSystems.

Figura 14: Modelo SIMULINK de una sonda pasiva de tensión atenuadora 10X conectada a un osciloscopio,

Los bloques de la librerı́a SimPowerSystems operan con señales que se podrı́an llamar “señales de potencia”
para facilitar la comprensión, mientras que los bloques generales de SIMULINK operan con señales que se
podrı́an llamar “señales matemáticas”. A fin de diferenciar las señales a nivel SimPowerSystems (“señales
de potencia”) respecto las señales a nivel SIMULINK (“señales matemáticas”), es conveniente crear un
subsistema a partir del circuito dibujado.

Figura 15: Simulación del modelo anterior alojado en un subsistema.

Se sugiere nombrar el subsistema creado con el nombre modelo sonda. En cada simulación se exportarán

Mediciones Eléctricas 1 Página 21 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

los datos desde el modelo de SIMULINK hacia el espacio de trabajo (Workspace) de MATLAB mediante el
bloque To Workspace, en el cuál se indica el nombre de la variable de salida, haciéndo doble click sobre el
mismo. Se sugiere nombrar esta variable como salida vo i.
Se han utilizado los siguientes bloques:

Powergui: Library/SimPowerSystems/powergui

Generador de señales: Library/Simulink/Sources/Signal Generator

Fuente de tension controlada: Library/SimPowerSystems/Electrical Sources/Controlled Voltage Source

Resistencia y capacitor: Library/SimpowerSistems/Elements/Series RLC Branch

Referencia (GND): Library/SimpowerSistems/Elements/Ground

Voltı́metro: Library/SimPowerSystems/Measurements/Voltage Measurement

Subsistema: Library/Simulink/Commonly Used Blocks/Subsystem

Bloque To Workspace: Library/Simulink/Sinks/To Workspace

Se sugiere guardar el archivo del modelo SIMULINK con el nombre sonda 10X tiempo.mdl.
Como ejemplo de interacción entre SIMULINK Y MATLAB se escribirá un script MATLAB que realice
reiteradas simulaciones del modelo SIMULINK, modificando en cada una de ellas el valor de capacitancia
CP en el circuito, y almacenando los datos obtenidos en una variable para graficarlas al finalizar el proceso.

%% Script desde el cual se ejecutaran las simulaciones de SIMULINK

close a l l ; % Cierra todas las ventanas graficas
clear a l l ; % Limpia todas las variables del Workspace
c l c ; % Limpia la ventana de comandos
5 %% Parametros fijos cargados en los bloques de SIMULINK
Vs_amp = 1;
Vs_frec = 1e3;
Rs = 50;
Rp = 9e6;
10 Co = 20e-12;
Ro = 1e6;
Cc = 100e-12;
%% Analisis parametrico
N=5; % Se define la cantidad de simulaciones del modelo
15 leyendas = cell(N,1); % Se define un vector con las leyendas de cada curva
tstop=0.0025; % Parametro cargado en SIMULINK que define el tiempo total de ...
simulacion
C = {'b',[0 .8 0],'r','c',[.8 .2 .9],[.9 .4 .2],[.9 .9 .9],[.7 .7 .7]}; % Vector ...
con los colores de las curvas
open_system('sonda_10x_tiempo') % Abre el modelo de SIMULINK
for i=1:1:N % Se inicia un bucle for de N iteraciones
20 Cp=(10.3+i)*10ˆ(-12); % Se define el valor de Cp en cada iteracion
sim('sonda_10x_tiempo'); % Simula el modelo de SIMULINK
salida_vo(i)=salida_vo_i; % Guarda la salida de cada simulacion
leyendas{i} = ['Cp = ' num2str(Cp*10ˆ12) 'pF']; % Se escribe la leyenda de ...
cada curva
c l c ; % Limpia la ventana de comandos
25 mensaje = ['Iteracion ' num2str(i)];
disp(mensaje);

Mediciones Eléctricas 1 Página 22 de 23

 Ejercicios de mediciones eléctricas con M AT LAB r & Simulink r

end % Fin del bucle for

%% Graficos
figure % Crea una nueva figura
30 hold on; % Permite agregar mas curvas al mismo grafico

for i=1:1:N % Se inicia un bucle for para realizar los N graficos

plot(salida_vo(i).time,salida_vo(i).signals.values,'color',C{i},'LineWidth',3) ...
% Traza una curva en cada ciclo
end; % Fin del bucle for
35

grid on; % Agrega una grilla a la figura

t i t l e ('Tension a la salida de sonda vs. Tiempo') % Agrega un titulo a la figura
xlabel('t[s]'); ylabel('Vo[V]'); % Agrega etiquetas en ambos ejes del grafico
legend(leyendas); % Agrega leyendas a la figura
40 clc
disp('Simulaciones y graficos finalizados')

Bibliografı́a recomendada
SimP owerSystemsT M 5 User’s Guide, Hydro-Québec and The MathWorks, Inc, 2010.

SimscapeT M User’s Guide, The MathWorks, Inc, 2016.

Steven T. Karris, Circuit Analysis I with M AT LAB r Computing and Simulink r /SimP owerSystemsr
Modeling, Orchard Publications, 2009.

Steven T. Karris, Circuit Analysis II with M AT LAB r Computing and Simulink r /SimP owerSystemsr
Modeling, Orchard Publications, 2009.

Justo Roberts, Gustavo Bacino, Electrotecnia. Apunte de cátedra: Introducción a MATLAB

Justo Roberts, Gustavo Bacino, Electrotecnia. Apunte de cátedra: Introducción a SIMULINK

Viktor M. Perelmuter, Electrotechnical Systems. Simulation with Simulink r and SimP owerSystemsT M ,
CRC Press, 2013.

Adrian B. Biran, What Every Engineer Should Know About M AT LAB r and Simulink r , CRC Press,
2010

Holly Moore, M AT LAB r para ingenieros, PEARSON Prentice Hall, 2007

Miguel Ataurima Arellano, M AT LAB r & Simulink r para ingenierı́a. Nivel 1., UCH Mathworks,
2013.

Manuel Gil Rodriguez, Introducción rápida a MATLAB y SIMULINK para ciencia e ingenierı́a, Ed.
Diaz de Santos, 2003.

Mediciones Eléctricas 1 Página 23 de 23