From 2647f7e4d4636157d9d5aab7bc0fcb24cb17b630 Mon Sep 17 00:00:00 2001 From: "Robin. P. Clark" Date: Tue, 10 Sep 2019 09:42:57 +0100 Subject: [PATCH] first attempt at zg. Did not make much difference although increased the gain. Perhaps the zeros needs to be closed in frequency to the dc line These ones were at 7/8 which is 29 degrees approx. Probably not enough to make a real difference. --- tp.gpt | 2 +- two_pole_7_8.c | 5 ++++- 2 files changed, 5 insertions(+), 2 deletions(-) diff --git a/tp.gpt b/tp.gpt index a0bd1ef..8c98e46 100644 --- a/tp.gpt +++ b/tp.gpt @@ -3,7 +3,7 @@ #plot "tp.dat" using 1:3 title "input", "tp.dat" using 1:5 title "two pole", "tp.dat" using 1:6 title "34LAG", "tp.dat" using 1:7 title "78LAG" #plot "tp.dat" using 1:3 title "input" with linespoints, "tp.dat" using 1:5 title "two pole" with linespoints, "tp.dat" using 1:7 title "78LAG" with linespoints, "tp.dat" using 1:9 title "78LAG TWICE" with linespoints, "tp.dat" using 1:11 title "two pole 15 16" plot "tp.dat" using 1:3 title "input" with linespoints, "tp.dat" using 1:5 title "two pole" with linespoints, "tp.dat" using 1:11 title "two pole 15 16" with linespoints, "tp.dat" using 1:13 title "zero only" with linespoints -plot "tp.dat" using 1:3 title "input" with linespoints, "tp.dat" using 1:5 title "two pole" with linespoints, "tp.dat" using 1:15 title "zero only^2" with linespoints, "tp.dat" using 1:17 title "zero only^3" with linespoints,"tp.dat" using 1:19 title "zero only^4" with linespoints +plot "tp.dat" using 1:3 title "input" with linespoints, "tp.dat" using 1:5 title "two pole" with linespoints, "tp.dat" using 1:15 title "7 8 two pole ZG" with linespoints, "tp.dat" using 1:17 title "zero only^3" with linespoints,"tp.dat" using 1:19 title "zero only^4" with linespoints !sleep 10 !sleep 10 diff --git a/two_pole_7_8.c b/two_pole_7_8.c index 94323ca..91abe26 100644 --- a/two_pole_7_8.c +++ b/two_pole_7_8.c @@ -80,6 +80,9 @@ two_pole_7_8 ( int16_t input ) { } +// Added zero guard. tried a pair of zeros near unit circle +// 7/8 in ADJ, i.e hyp = 1 so on unit circle. +// int16_t /* squared version of LAG_7_8 */ two_pole_7_8_zg ( int16_t input ) { // with zero guard @@ -287,7 +290,7 @@ int main () { - zo2 = zero_only_2 (val); + zo2 = two_pole_7_8_zg(val); // zero_only_2 (val); zo3 = zero_only_3 (val); zo4 = zero_only_4 (val);