clc
clear
%Part 1B
f1b= [5 10 20 40 60 70 80 90 100 200];
f1b= transpose(f1b);
f1b= f1b.*1000;

vin1b= [];
vin1b= transpose(vin1b);

vout1b= [];
vout1b= transpose(vout1b);

p1b= [];
p1b= transpose(p1b);

mag1b= 20.*log10(abs(vout1b./vin1b));

figure(1)
plot(f1b, mag1b)
title('Magnitude vs. Frequency Bode Plot for Passive Low Pass Filter')
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')

figure(2)
plot(f1b, p1b)
title('Phase vs. Frequency Bode Plot for Passive Low Pass Filter')
xlabel('Frequency (Hz)')
ylabel('Phase (deg)')
%%
%Part 1D
f1d= [5 10 20 40 60 70 80 90 100 200];
f1d= transpose(f1d);
f1d= f1d.*1000;

vin1d= [];
vin1d= transpose(vin1d);

vout1d= [];
vout1d= transpose(vout1d);

p1d= [];
p1d= transpose(p1d);

mag1d= 20.*log10(abs(vout1d./vin1d));

figure(3)
plot(f1d, mag1d)
title('Magnitude vs. Frequency Bode Plot for Passive High Pass Filter')
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')

figure(4)
plot(f1d, p1d)
title('Phase vs. Frequency Bode Plot for Passive High Pass Filter')
xlabel('Frequency (Hz)')
ylabel('Phase (deg)')
%%
%Part 2B
f2b= [100 200 400 600 800 1e3 2e3 4e3 6e3 8e3 10e3 20e3 40e3 60e3];
f2b= transpose(f2b);

vin2b= [];
vin2b= transpose(vin2b);

vout2b= [];
vout2b= transpose(vout2b);

p2b= [];
p2b= transpose(p2b);

mag2b= 20.*log10(abs(vout2b./vin2b));

figure(5)
plot(f2b, mag2b)
title('Magnitude vs. Frequency Bode Plot for Active Low Pass Filter')
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')

figure(6)
plot(f2b, p2b)
title('Phase vs. Frequency Bode Plot for Active Low Pass Filter')
xlabel('Frequency (Hz)')
ylabel('Phase (deg)')
%%
%Part 2C
%300 mV
f2c300= [400 600 800 1e3 2e3 4e3 6e3 8e3 10e3 20e3 40e3 60e3];
f2c300= transpose(f2c300);

vin2c300= [];
vin2c300= transpose(vin2c300);

vout2c300= [];
vout2c300= transpose(vout2c300);

p2c300= [];
p2c300= transpose(p2c300);

mag2c300= 20.*log10(abs(vout2c300./vin2c300));

figure(7)
plot(f2c300, mag2c300)
title('Magnitude vs. Frequency Bode Plot for Active High Pass Filter at 300 mV')
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')

figure(8)
plot(f2c300, p2c300)
title('Phase vs. Frequency Bode Plot for Active High Pass Filter at 300 mV')
xlabel('Frequency (Hz)')
ylabel('Phase (deg)')
%%
%Part 2C
%1 V
f2c1= [400 600 800 1e3 2e3 4e3 6e3 8e3 10e3 20e3 40e3 60e3];
f2c1= transpose(f2c1);

vin2c1= [];
vin2c1= transpose(vin2c1);

vout2c1= [];
vout2c1= transpose(vout2c1);

p2c1= [];
p2c1= transpose(p2c1);

mag2c1= 20.*log10(abs(vout2c1./vin2c1));

figure(9)
plot(f2c1, mag2c1)
title('Magnitude vs. Frequency Bode Plot for Active High Pass Filter at 1 V')
xlabel('Frequency (Hz)')
ylabel('Magnitude (dB)')

figure(10)
plot(f2c1, p2c1)
title('Phase vs. Frequency Bode Plot for Active High Pass Filter at 1 V')
xlabel('Frequency (Hz)')
ylabel('Phase (deg)')