changed values and added RIAA curve to plot
authorNils <forssennils@gmail.com>
Fri, 12 Jul 2024 20:10:05 +0000 (22:10 +0200)
committerNils <forssennils@gmail.com>
Fri, 12 Jul 2024 20:10:05 +0000 (22:10 +0200)
matlab/calculations.m

index b13828636e43a831072c9e2334f529002d6ef9fd..a072e312fb522d8740ec6dd140c834f24f3be122 100644 (file)
@@ -2,8 +2,8 @@ s = tf('s');
 Z_R4 = 180000;\r
 Z_R3 = 10000;\r
 Z_R5 = 4700;\r
-Z_R6 = 16000;\r
-Z_R7 = 330;\r
+Z_R6 = 100000;\r
+Z_R7 = 2200;\r
 Z_R13 = 750;\r
 Z_R14 = 100000;\r
 \r
@@ -22,11 +22,16 @@ H2 = ((1 / (Z_R7 + Z_C3)) + (1 / Z_R6)) * Z_R6;
 H3 = (1 / (Z_R13 * Z_C8)) / (s + (1 / (Z_R13 * Z_C8)));\r
 H4 = (s * Z_R14 * Z_C9 / (1 + s * Z_R14 * Z_C9));\r
 \r
-H_tot = H1 * H2 * H3 * H4;\r
+H_tot = (1 / 2) * H1 * H2 * H3 * H4;\r
+\r
+H_RIAA = 10 * ((s / (2 * pi * 500)) + 1) / (((s / (2 * pi * 50)) + 1)*((s / (2 * pi * 2122)) + 1));\r
+\r
+H_RIAA_amp = 100 * H_RIAA;\r
 \r
 p = bodeoptions('cstprefs');\r
 p.FreqUnits = 'Hz';\r
 p.PhaseVisible = 'off';\r
-p.XLim = {[10,30000]};\r
+p.XLim = {[20,20000]};\r
+\r
 \r
-bodeplot(H_tot,p)
\ No newline at end of file
+bodeplot(H_tot, H_RIAA_amp,p)
\ No newline at end of file