|
|
|
@@ -3,19 +3,22 @@ S = asymptote(S_smaller,S_larger,f,Fc); |
|
|
|
P = mod(P,pi); |
|
|
|
|
|
|
|
figure; |
|
|
|
subplot(2,1,1); |
|
|
|
ax1 = subplot(2,1,1); |
|
|
|
semilogx(f,A,f,S) |
|
|
|
ylabel('Magnitude [dB]','Fontsize',10); |
|
|
|
xlabel('Frequency [Hz]','Fontsize',10); |
|
|
|
xlim([f_start,f_stop]); |
|
|
|
legend('Theory','Asymptote'); |
|
|
|
grid on |
|
|
|
title('RLC-circuit Bodeplot Theory') |
|
|
|
subplot(2,1,2); |
|
|
|
ax2 = subplot(2,1,2); |
|
|
|
semilogx(f,P) |
|
|
|
ylabel('Phase [rad]','Fontsize',10); |
|
|
|
xlabel('Frequency [Hz]','Fontsize',10); |
|
|
|
xlim([f_start,f_stop]); |
|
|
|
legend('Theory'); |
|
|
|
grid on |
|
|
|
linkaxes([ax1,ax2],'x') |
|
|
|
|
|
|
|
%% Run the measurements |
|
|
|
if theory_only == false |
|
|
|
@@ -28,24 +31,25 @@ if theory_only == false |
|
|
|
|
|
|
|
%% plot average measurement, theory and asymptote. |
|
|
|
figure; |
|
|
|
subplot(2,1,1); |
|
|
|
ax1 = subplot(2,1,1); |
|
|
|
semilogx(data.frequency,20*log10(A_mean),f,A,f,S) |
|
|
|
ylabel('Magnitude [dB]','Fontsize',10); |
|
|
|
xlabel('Frequency [Hz]','Fontsize',10); |
|
|
|
legend('Average Measurement','Theory','Asymptote'); |
|
|
|
grid on |
|
|
|
title('RLC-circuit Bodeplot') |
|
|
|
subplot(2,1,2); |
|
|
|
ax2 = subplot(2,1,2); |
|
|
|
semilogx(data.frequency,P_mean,f,P) |
|
|
|
ylabel('Phase [rad]','Fontsize',10); |
|
|
|
xlabel('Frequency [Hz]','Fontsize',10); |
|
|
|
legend('Average Measurement','Theory'); |
|
|
|
grid on |
|
|
|
linkaxes([ax1,ax2],'x') |
|
|
|
|
|
|
|
%% plot avearage measurement with standard deviation and theory. |
|
|
|
figure; |
|
|
|
ax = subplot(2,1,1); |
|
|
|
ax.XScale = 'log'; |
|
|
|
ax1 = subplot(2,1,1); |
|
|
|
ax1.XScale = 'log'; |
|
|
|
errorbar(data.frequency,20*log10(A_mean),20*log10(A_std)) |
|
|
|
hold on |
|
|
|
semilogx(f,A) |
|
|
|
@@ -55,8 +59,8 @@ if theory_only == false |
|
|
|
legend('Average Measurement','Theory','Asymptote'); |
|
|
|
grid on |
|
|
|
title('RLC-circuit Bodeplot') |
|
|
|
ax = subplot(2,1,2); |
|
|
|
ax.XScale = 'log'; |
|
|
|
ax2 = subplot(2,1,2); |
|
|
|
ax2.XScale = 'log'; |
|
|
|
errorbar(data.frequency,P_mean,P_std); |
|
|
|
hold on |
|
|
|
semilogx(data.frequency,P_mean,f,P) |
|
|
|
@@ -65,21 +69,23 @@ if theory_only == false |
|
|
|
xlabel('Frequency [Hz]','Fontsize',10); |
|
|
|
legend('Average Measurement','Theory'); |
|
|
|
grid on |
|
|
|
linkaxes([ax1,ax2],'x') |
|
|
|
%% plot avearage measurement with standard deviation. |
|
|
|
figure; |
|
|
|
ax = subplot(2,1,1); |
|
|
|
ax.XScale = 'log'; |
|
|
|
ax1 = subplot(2,1,1); |
|
|
|
ax1.XScale = 'log'; |
|
|
|
errorbar(data.frequency,20*log10(A_mean),20*log10(A_std)) |
|
|
|
ylabel('Magnitude [dB]','Fontsize',10); |
|
|
|
xlabel('Frequency [Hz]','Fontsize',10); |
|
|
|
legend('Average Measurement'); |
|
|
|
grid on |
|
|
|
title('RLC-circuit Bodeplot') |
|
|
|
ax = subplot(2,1,2); |
|
|
|
ax.XScale = 'log'; |
|
|
|
ax2 = subplot(2,1,2); |
|
|
|
ax2.XScale = 'log'; |
|
|
|
errorbar(data.frequency,P_mean,P_std); |
|
|
|
ylabel('Phase [rad]','Fontsize',10); |
|
|
|
xlabel('Frequency [Hz]','Fontsize',10); |
|
|
|
legend('Average Measurement'); |
|
|
|
grid on |
|
|
|
linkaxes([ax1,ax2],'x') |
|
|
|
end |