static void _Main(string[] args) { const int N = 1024; double[] x = new double[N]; IFilter filter = new ShelvingEqualizer(Common.Normalize(1000), 1); // インパルス応答取得 filter.Clear(); /* filter.GetValue(1); for(int i=0; i<16; ++i) filter.GetValue(0); */ filter.GetValue(1); for(int i=1; i<N; ++i) x[i] = filter.GetValue(0); #if TIME_SEQUENCE int n = N; double[] t = new double[n]; for(int i=0; i<n; ++i) t[i] = i; Graph.GraphForm gf = new Graph.GraphForm(); gf.Graph.AddEntry(t, x, new Pen(Color.Crimson)); gf.Graph.SetXAxis(0, n, 5, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Graph.SetYAxis(0, 2, 5, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Size = new Size(640, 480); gf.ShowDialog(); #elif PHASE Spectrum spec = Spectrum.FromTimeSequence(x, 1); int n = spec.Count; double[] t = new double[n]; for(int i=0; i<n; ++i) t[i] = i; double[] y = spec.GetPhase(); Spectrum.Unwrap(y); Graph.GraphForm gf = new Graph.GraphForm(); gf.Graph.AddEntry(t, y, new Pen(Color.Crimson)); gf.Graph.SetXAxis(0, n, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Graph.SetYAxis(-50, 0, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Size = new Size(640, 480); gf.ShowDialog(); #elif GAIN Spectrum spec = Spectrum.FromTimeSequence(x, 1); int n = spec.Count; double[] t = new double[n]; for(int i=0; i<n; ++i) t[i] = i; double[] y = spec.GetAmplitude(); Graph.GraphForm gf = new Graph.GraphForm(); gf.Graph.AddEntry(t, y, new Pen(Color.Crimson)); gf.Graph.SetXAxis(0, n, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Graph.SetYAxis(0, 4, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Size = new Size(640, 480); gf.ShowDialog(); #endif #if false const int N = 512; double[] x = new double[N]; double[] x0 = new double[N]; double[] y0 = new double[N]; /* Random rnd = new Random(); x0[0] = 1; x0[1] = 1; y0[0] = 0; y0[1] = 0; for(int i=1; i<N/2; ++i) { x0[2*i] = rnd.NextDouble(); x0[2*i+1] = rnd.NextDouble(); y0[2*i] = -x0[2*i+1]; y0[2*i+1] = x0[2*i]; } Fft fft = new Fft(N); fft.Invert(x0); fft.Invert(y0); for(int i=0; i<N; ++i) x[i] = x0[i]; //*/ //* for(int i=0; i<N; ++i) { x[i] = 100; for(int k=1; k<200; ++k) x[i] += short.MaxValue / 4 / k * Math.Cos(Math.PI / 512 * 32 * i * k); x0[i] = x[i]; y0[i] = 0; for(int k=1; k<200; ++k) y0[i] += short.MaxValue / 4 / k * Math.Sin(Math.PI / 512 * 32 * i * k); } //*/ double[] y = Spectrum.HilbertTransform(x); using(StreamWriter writer = new StreamWriter("z.csv", false, Encoding.Default)) { /* for(int i=0; i<x.Length; ++i) writer.Write("{0},{1},{2},{3},{4}\n", i, x0[i], x[i], y0[i], y[i]); //*/ //* Spectrum a0 = Spectrum.FromTimeSequence(x0); Spectrum a = Spectrum.FromTimeSequence(x); Spectrum b0 = Spectrum.FromTimeSequence(y0); Spectrum b = Spectrum.FromTimeSequence(y); for(int i=0; i<a.Count; ++i) { writer.Write("{0},{1},{2},{3},{4}\n", i, a0[i].Abs, a[i].Abs, b0[i].Abs, b[i].Abs); // writer.Write("{0},{1},{2},{3},{4}\n", i, a0[i].Arg, a[i].Arg, b0[i].Arg, b[i].Arg); } //*/ } #elif false const int N = 1024; double[] lp = new double[N]; IFilter lpf = new LowPassFir(16, Common.Normalize(2000)); lp[0] = lpf.GetValue(1); for(int i=1; i<N; ++i) lp[i] = lpf.GetValue(0); double[] bp = new double[N]; IFilter bpf = new BandPassFir(16, Common.Normalize(2000), Common.Normalize(8000)); lp[0] = bpf.GetValue(1); for(int i=1; i<N; ++i) bp[i] = bpf.GetValue(0); double[] hp = new double[N]; IFilter hpf = new HighPassFir(16, Common.Normalize(8000)); lp[0] = hpf.GetValue(1); for(int i=1; i<N; ++i) hp[i] = hpf.GetValue(0); Fft fft = new Fft(N); fft.Transform(1, lp); fft.Transform(1, bp); fft.Transform(1, hp); using(StreamWriter writer = new StreamWriter("z.csv", false, Encoding.Default)) { double la = Common.Amp(lp[0]); double ba = Common.Amp(bp[0]); double ha = Common.Amp(hp[0]); writer.Write("{0},{1},{2},{3}\n", 0, la, ba, ha); for(int i=1; i<N/2; ++i) { la = Common.Amp(lp[2*i] * lp[2*i] + lp[2*i+1] * lp[2*i+1]); ba = Common.Amp(bp[2*i] * bp[2*i] + bp[2*i+1] * bp[2*i+1]); ha = Common.Amp(hp[2*i] * hp[2*i] + hp[2*i+1] * hp[2*i+1]); writer.Write("{0},{1},{2},{3}\n", i, la, ba, ha); } la = Common.Amp(lp[1]); ba = Common.Amp(bp[1]); ha = Common.Amp(hp[1]); writer.Write("{0},{1},{2},{3}\n", N/2, la, ba, ha); } #elif false uint length; double[] l; double[] r; FormatHeader header; using(WaveReader reader = new WaveReader(@"TD100V00H030.wav")) { header = reader.Header; length = reader.Length; int tmp = reader.Read(length, out l, out r); } IFilter lpf = new LowPassFir(16, Common.Normalize(4500)); IFilter bpf = new BandPassFir(16, Common.Normalize(4500), Common.Normalize(10000)); IFilter hpf = new HighPassFir(16, Common.Normalize(10000)); for(int i=0; i<l.Length; ++i) { double lp = lpf.GetValue(l[i]); double bp = bpf.GetValue(l[i]); double hp = hpf.GetValue(l[i]); l[i] = lp + bp + hp; } if(r != null) for(int i=0; i<l.Length; ++i) { double lp = lpf.GetValue(r[i]); double bp = bpf.GetValue(r[i]); double hp = hpf.GetValue(r[i]); r[i] = lp + bp + hp; } using(WaveWriter writer = new WaveWriter("z.wav", header)) { writer.Write(l, r); } #elif false using(StreamWriter writer = new StreamWriter("z.csv", false, Encoding.Default)) { for(int i=0; i<l.Length; ++i) { writer.Write("{0},{1},{2}\n", i, l[i], r[i]); } } #endif }
static void _Main(string[] args) { const int N = 1024; double[] x = new double[N]; IFilter filter = new ShelvingEqualizer(Common.Normalize(1000), 1); // インパルス応答取得 filter.Clear(); /* * filter.GetValue(1); * for(int i=0; i<16; ++i) * filter.GetValue(0); */ filter.GetValue(1); for (int i = 1; i < N; ++i) { x[i] = filter.GetValue(0); } #if TIME_SEQUENCE int n = N; double[] t = new double[n]; for (int i = 0; i < n; ++i) { t[i] = i; } Graph.GraphForm gf = new Graph.GraphForm(); gf.Graph.AddEntry(t, x, new Pen(Color.Crimson)); gf.Graph.SetXAxis(0, n, 5, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Graph.SetYAxis(0, 2, 5, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Size = new Size(640, 480); gf.ShowDialog(); #elif PHASE Spectrum spec = Spectrum.FromTimeSequence(x, 1); int n = spec.Count; double[] t = new double[n]; for (int i = 0; i < n; ++i) { t[i] = i; } double[] y = spec.GetPhase(); Spectrum.Unwrap(y); Graph.GraphForm gf = new Graph.GraphForm(); gf.Graph.AddEntry(t, y, new Pen(Color.Crimson)); gf.Graph.SetXAxis(0, n, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Graph.SetYAxis(-50, 0, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Size = new Size(640, 480); gf.ShowDialog(); #elif GAIN Spectrum spec = Spectrum.FromTimeSequence(x, 1); int n = spec.Count; double[] t = new double[n]; for (int i = 0; i < n; ++i) { t[i] = i; } double[] y = spec.GetAmplitude(); Graph.GraphForm gf = new Graph.GraphForm(); gf.Graph.AddEntry(t, y, new Pen(Color.Crimson)); gf.Graph.SetXAxis(0, n, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Graph.SetYAxis(0, 4, 4, new Font(FontFamily.GenericSerif, 12), new SolidBrush(Color.Black)); gf.Size = new Size(640, 480); gf.ShowDialog(); #endif #if false const int N = 512; double[] x = new double[N]; double[] x0 = new double[N]; double[] y0 = new double[N]; /* * Random rnd = new Random(); * x0[0] = 1; * x0[1] = 1; * y0[0] = 0; * y0[1] = 0; * for(int i=1; i<N/2; ++i) * { * x0[2*i] = rnd.NextDouble(); * x0[2*i+1] = rnd.NextDouble(); * y0[2*i] = -x0[2*i+1]; * y0[2*i+1] = x0[2*i]; * } * Fft fft = new Fft(N); * fft.Invert(x0); * fft.Invert(y0); * for(int i=0; i<N; ++i) * x[i] = x0[i]; * //*/ //* for (int i = 0; i < N; ++i) { x[i] = 100; for (int k = 1; k < 200; ++k) { x[i] += short.MaxValue / 4 / k * Math.Cos(Math.PI / 512 * 32 * i * k); } x0[i] = x[i]; y0[i] = 0; for (int k = 1; k < 200; ++k) { y0[i] += short.MaxValue / 4 / k * Math.Sin(Math.PI / 512 * 32 * i * k); } } //*/ double[] y = Spectrum.HilbertTransform(x); using (StreamWriter writer = new StreamWriter("z.csv", false, Encoding.Default)) { /* * for(int i=0; i<x.Length; ++i) * writer.Write("{0},{1},{2},{3},{4}\n", i, x0[i], x[i], y0[i], y[i]); * //*/ //* Spectrum a0 = Spectrum.FromTimeSequence(x0); Spectrum a = Spectrum.FromTimeSequence(x); Spectrum b0 = Spectrum.FromTimeSequence(y0); Spectrum b = Spectrum.FromTimeSequence(y); for (int i = 0; i < a.Count; ++i) { writer.Write("{0},{1},{2},{3},{4}\n", i, a0[i].Abs, a[i].Abs, b0[i].Abs, b[i].Abs); // writer.Write("{0},{1},{2},{3},{4}\n", i, a0[i].Arg, a[i].Arg, b0[i].Arg, b[i].Arg); } //*/ } #elif false const int N = 1024; double[] lp = new double[N]; IFilter lpf = new LowPassFir(16, Common.Normalize(2000)); lp[0] = lpf.GetValue(1); for (int i = 1; i < N; ++i) { lp[i] = lpf.GetValue(0); } double[] bp = new double[N]; IFilter bpf = new BandPassFir(16, Common.Normalize(2000), Common.Normalize(8000)); lp[0] = bpf.GetValue(1); for (int i = 1; i < N; ++i) { bp[i] = bpf.GetValue(0); } double[] hp = new double[N]; IFilter hpf = new HighPassFir(16, Common.Normalize(8000)); lp[0] = hpf.GetValue(1); for (int i = 1; i < N; ++i) { hp[i] = hpf.GetValue(0); } Fft fft = new Fft(N); fft.Transform(1, lp); fft.Transform(1, bp); fft.Transform(1, hp); using (StreamWriter writer = new StreamWriter("z.csv", false, Encoding.Default)) { double la = Common.Amp(lp[0]); double ba = Common.Amp(bp[0]); double ha = Common.Amp(hp[0]); writer.Write("{0},{1},{2},{3}\n", 0, la, ba, ha); for (int i = 1; i < N / 2; ++i) { la = Common.Amp(lp[2 * i] * lp[2 * i] + lp[2 * i + 1] * lp[2 * i + 1]); ba = Common.Amp(bp[2 * i] * bp[2 * i] + bp[2 * i + 1] * bp[2 * i + 1]); ha = Common.Amp(hp[2 * i] * hp[2 * i] + hp[2 * i + 1] * hp[2 * i + 1]); writer.Write("{0},{1},{2},{3}\n", i, la, ba, ha); } la = Common.Amp(lp[1]); ba = Common.Amp(bp[1]); ha = Common.Amp(hp[1]); writer.Write("{0},{1},{2},{3}\n", N / 2, la, ba, ha); } #elif false uint length; double[] l; double[] r; FormatHeader header; using (WaveReader reader = new WaveReader(@"TD100V00H030.wav")) { header = reader.Header; length = reader.Length; int tmp = reader.Read(length, out l, out r); } IFilter lpf = new LowPassFir(16, Common.Normalize(4500)); IFilter bpf = new BandPassFir(16, Common.Normalize(4500), Common.Normalize(10000)); IFilter hpf = new HighPassFir(16, Common.Normalize(10000)); for (int i = 0; i < l.Length; ++i) { double lp = lpf.GetValue(l[i]); double bp = bpf.GetValue(l[i]); double hp = hpf.GetValue(l[i]); l[i] = lp + bp + hp; } if (r != null) { for (int i = 0; i < l.Length; ++i) { double lp = lpf.GetValue(r[i]); double bp = bpf.GetValue(r[i]); double hp = hpf.GetValue(r[i]); r[i] = lp + bp + hp; } } using (WaveWriter writer = new WaveWriter("z.wav", header)) { writer.Write(l, r); } #elif false using (StreamWriter writer = new StreamWriter("z.csv", false, Encoding.Default)) { for (int i = 0; i < l.Length; ++i) { writer.Write("{0},{1},{2}\n", i, l[i], r[i]); } } #endif }