%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                                                            %
%                Løsningsforslag til oppgave (b) i øving 2,                  %
%                Bølgefysikk, høst 2009                                      %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Vi skal plotte ladningsamplituden q0, strømamplituden i0 = omega*q0, og
% fasevinkelen alfa, alle tre som funksjoner av vinkelfrekvensen omega, omkring
% resonansfrekvensen omega0 = 1/sqrt(L*C) = 1E6 pr sek. Halvverdibredden
% Delta omega = R/L har her verdien 1E4 pr sek, slik at et passende intervall
% for omega er fra 9.5E5 til 10.5E5. Det skulle holde med 100 datapunkter på
% dette intervallet:
omega = linspace(9.5E5,10.5E5,100);
% La oss for oversiktens skyld bruke symboler for de ulike størrelsene som
% inngår, og gi dem de aktuelle tallverdiene, i SI-enheter:
V0 = 10;
R = 1.0;
L = 0.1E-3;
C = 10E-9;
% Kretsens reaktans er X(omega) = omega*L-1/omega*C:
X = omega.*L - 1./(omega.*C);
% Pass på punktumene...!
% Strømamplituden:
i0 = V0./(sqrt(R^2+X.^2));
% Ladningsamplituden:
q0 = i0./omega;
% Fasekonstanten:
alfa = atan(X./R);
% Kommandoen figure ber om ny figur
% Først q0:
figure;
plot(omega,q0);
xlabel('vinkelfrekvens (1/s)');
ylabel('ladningsamplitude (C)');
% Deretter i0:
figure;
plot(omega,i0);
xlabel('vinkelfrekvens (1/s)');
ylabel('stromamplitude (A)');
% Deretter alfa:
figure;
plot(omega,alfa);
xlabel('vinkelfrekvens (1/s)');
ylabel('fasekonstant (rad)');