Der Spannungsteiler

Mathi-Homework oder "Alleingelassen in der bösen Welt des Matlab"

Hier mein .m-File:

%Spannungsteiler - Berechnung
%Hausaufgabe Mathi TSU 28.3.2000
%Martin Hess
%27.3.2000

%Konstanten
clear;
Us=10; %Speisspannung
Rp=10000; %Potentiometer-Widerstand
a=linspace(1,270,270); %Drehbereich Poti

%Belastungen
Rl(1)=1;
Rl(2)=10*Rp; %Belastungen aus Aufgabenstellung
Rl(3)=5*Rp;
Rl(4)=Rp;
Rl(5)=0.5*Rp;
Rl(6)=0.1*Rp;

%Berechnungen
Ra=(Rp*a)./270; %Ra in Abhängigkeit des Drehwinkels
Va=Us.*(((Rp.*a)./270)./Rp); %Unbelasteter Spannungs-Teiler
for i=2:6,
Ua(i,:)=Us.*(1./(1./Ra+1./Rl(i)))./((Rp-Ra)+1./(1./Ra+1./Rl(i)));
Uabw(i,:)=Va-Ua(i,:);
xmax(i)=max(Uabw(i,:));
ymax(i)=find(max(Uabw(i,:)) == Uabw(i,:));
end;

%Plot
figure; %Neues Fenster
clf reset; %reset Fenster
subplot(2,1,1); %erster Subplot
grid; %Raster ein
hold on; %hold
plot(a,Va,'r',a,Ua(2,:),'b',a,Ua(3,:),'g',a,Ua(4,:),'y',a,Ua(5,:),'m',a,Ua(6,:),'c')
xlabel('Winkel [°]'); %Bezeichnung der Achsen
ylabel('Spannung [V]');
legend('Unbelastet','Rl=10*Rp','Rl=5*Rp','Rl=Rp','Rl=0.5*Rp','Rl=0.1*Rp',0);
subplot(2,1,2); %zweiter subplot, analog wie oben
grid;
hold on;
plot(a,Uabw(2,:),'b',a,Uabw(3,:),'g',a,Uabw(4,:),'y',a,Uabw(5,:),'m',a,Uabw(6,:),'c')
plot(ymax(2),xmax(2),'+k'); %plotten der Maximas, funktioniert nicht in einem Argument
plot(ymax(3),xmax(3),'+k');
plot(ymax(4),xmax(4),'+k');
plot(ymax(5),xmax(5),'+k');
plot(ymax(6),xmax(6),'+k');
xlabel('Winkel [°]');
ylabel('Spannungsabweichung [V]');

(Sorry für die hässliche Darstellung, aber irgendwie will Billy's Gratis-HTML-Editor nicht so wie ich will...)