Contents

% /home/dimarzio/Documents/working/12303/matlab/bridge.m
% Wed Sep 19 11:30:14 2018
%
% Demonstration of bridge circuit

R1=10e3;R2=10e3;Rx=1161;Rv=[1100:1200];
Vs=10;
y=[Vs/R1;Vs/Rx];

Ammeter

Rm=2;

  for n=1:length(Rv);
  M=[1/R1+1/Rm+1/R2,-1/Rm;-1/Rm,1/Rx+1/Rm+1./Rv(n)];
  x=inv(M)*y;
  v1(n)=x(1);
  v2(n)=x(2);
  end;
im=(v2-v1)/Rm;
figure;plot(Rv,im*1e6,'b-',Rx,0,'bo');grid on;
xlabel('Rv, Ohms');ylabel('im, \mu Amps');
title(['Rm = ',num2str(Rm),' \Omega']);
figure;plot(Rv,(v2-v1)*1e6,'b-',Rx,0,'bo');grid on;
xlabel('Rv, Ohms');ylabel('V2-V1, \mu Volts');
title(['Rm = ',num2str(Rm),' \Omega']);

Voltmeter

Rm=2e6;

  for n=1:length(Rv);
  M=[1/R1+1/Rm+1/R2,-1/Rm;-1/Rm,1/Rx+1/Rm+1./Rv(n)];
  x=inv(M)*y;
  v1(n)=x(1);
  v2(n)=x(2);
  end;
im=(v2-v1)/Rm;
figure;plot(Rv,im*1e9,'b-',Rx,0,'bo');grid on;
xlabel('Rv, Ohms');ylabel('im, nAmps');
title(['Rm = ',num2str(Rm/1e6),' M\Omega']);
figure;plot(Rv,(v2-v1)*1e3,'b-',Rx,0,'bo');grid on;
xlabel('Rv, Ohms');ylabel('V2-V1, mVolts');
title(['Rm = ',num2str(Rm/1e6),' M\Omega']);

Published with MATLAB® R2017a