%Spannungsteiler.m
%Belasteter Spannungsteiler
%Christian Widmer
%16.04.2002

figure (clf);

% Variablen
Us = 10; %Speisespannung
Rp = 10000; % [Ohm]Widerstand des Potis
Rl = 10000; % Lastwiderstand
min_Winkel=0;%Minimale Winkeleinstellung des Potis
max_Winkel=270;%Maximale Winkeleinstellung de Potis
Raster_Winkel=1;% Ratser des Potis
%Zusatz für Aufgabe 2.
Rlast0 = inf;
Rlast1 = 10 * Rp;
Rlast2 = 5 * Rp;
Rlast3 = Rp;
Rlast4 = 0.5 * Rp;
Rlast5 = 0.1 * Rp;
Rlast6 = 0.001 * Rp;

%Allgemeine Berechnungen
Winkel = min_Winkel:Raster_Winkel:max_Winkel;
Rv = (Rp.*(max_Winkel-Winkel))./max_Winkel;% Rv berechnen

% Ra Berechnen
Ra = Rp .* Winkel./ max_Winkel; % Ra im zusammenhang mit Potiposition (Winkelabhängig)

% Rp Berechnen

Ra_parallel_Rlast0 = 1./((1./Ra) + (1/Rlast0)); % Ra Parralel Rlast0 
Ra_parallel_Rlast1 = 1./((1./Ra) + (1/Rlast1)); % Ra Parralel Rlast1 
Ra_parallel_Rlast2 = 1./((1./Ra) + (1/Rlast2)); % Ra Parralel Rlast2 
Ra_parallel_Rlast3 = 1./((1./Ra) + (1/Rlast3)); % Ra Parralel Rlast3 
Ra_parallel_Rlast4 = 1./((1./Ra) + (1/Rlast4)); % Ra Parralel Rlast4 
Ra_parallel_Rlast5 = 1./((1./Ra) + (1/Rlast5)); % Ra Parralel Rlast5 
Ra_parallel_Rlast6 = 1./((1./Ra) + (1/Rlast6)); % Ra Parralel Rlast6 

% Ua Berechnen

Ua1 = Us.*Ra_parallel_Rlast0./(Rv+Ra_parallel_Rlast0); % Berechnung von Ua1
Ua2 = Us.*Ra_parallel_Rlast1./(Rv+Ra_parallel_Rlast1); % Berechnung von Ua2
Ua3 = Us.*Ra_parallel_Rlast2./(Rv+Ra_parallel_Rlast2); % Berechnung von Ua3
Ua4 = Us.*Ra_parallel_Rlast3./(Rv+Ra_parallel_Rlast3); % Berechnung von Ua4
Ua5 = Us.*Ra_parallel_Rlast4./(Rv+Ra_parallel_Rlast4); % Berechnung von Ua5
Ua6 = Us.*Ra_parallel_Rlast5./(Rv+Ra_parallel_Rlast5); % Berechnung von Ua6
Ua7 = Us.*Ra_parallel_Rlast6./(Rv+Ra_parallel_Rlast6); % Berechnung von Ua7

% Abweichung zum Unbelasteten Poti ( Rl= Unendlich)

Abweichung1 = (Ua1 - Ua2); % U Unbelastet - U Belastet = Abweichung
Abweichung2 = (Ua1 - Ua3); % U Unbelastet - U Belastet = Abweichung
Abweichung3 = (Ua1 - Ua4); % U Unbelastet - U Belastet = Abweichung
Abweichung4 = (Ua1 - Ua5); % U Unbelastet - U Belastet = Abweichung
Abweichung5 = (Ua1 - Ua6); % U Unbelastet - U Belastet = Abweichung
Abweichung6 = (Ua1 - Ua7); % U Unbelastet - U Belastet = Abweichung

% Maximale Stellen suchen

[max_Abweichung1, Indiz1] = max(Abweichung1);
[max_Abweichung2, Indiz2] = max(Abweichung2);
[max_Abweichung3, Indiz3] = max(Abweichung3);
[max_Abweichung4, Indiz4] = max(Abweichung4);
[max_Abweichung5, Indiz5] = max(Abweichung5);
[max_Abweichung6, Indiz6] = max(Abweichung6);



% erste Grafik erstellen
clf; % Grafikfenster löschen
subplot (2,1,1);
plot(Winkel, Ua1,Winkel,Ua2,Winkel,Ua3,Winkel,Ua4,Winkel,Ua5,Winkel,Ua6,Winkel,Ua7);
grid;
legend('Rl=Unendlich','Rl=10*Rp','Rl=5*Rp','Rl=Rp','Rl=0.5*Rp','Rl=0.1*Rp','Rl=0.001*Rp');
xlabel('Potiwinkel [°]');
ylabel('Ausgangsspannung [V]');
title ('Belastetes Potentimeter');

%zweite Grafik erstellen
subplot (2,1,2);
plot (Winkel,Abweichung1,Winkel,Abweichung2,Winkel,Abweichung3,Winkel,Abweichung4,Winkel,Abweichung5,Winkel,Abweichung6,Indiz1,max_Abweichung1,'rx',Indiz2,max_Abweichung2,'rx',Indiz3,max_Abweichung3,'rx',Indiz4,max_Abweichung4,'rx',Indiz5,max_Abweichung5,'rx',Indiz6,max_Abweichung6,'rx',[Indiz1,Indiz2],[max_Abweichung1,max_Abweichung2],'k',[Indiz2,Indiz3],[max_Abweichung2,max_Abweichung3],'k',[Indiz3,Indiz4],[max_Abweichung3,max_Abweichung4],'k',[Indiz4,Indiz5],[max_Abweichung4,max_Abweichung5],'k',[Indiz5,Indiz6],[max_Abweichung5,max_Abweichung6],'k');
grid;
legend('Abw. zu 10*Rp','Abw. zu 5*Rp','Abw. zu Rp','Abw. zu 0.5*Rp','Abw. zu 0.1*Rp','Abw. zu 0.001*Rp',2);
xlabel('Winkel [°]');
ylabel('Abweichung zu unbel. Poti [V]');
