ハイパス FIR フィルタ。
Inheritance: OddLinearFir
Exemplo n.º 1
0
		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
		}
Exemplo n.º 2
0
		static IFilter GetFilter()
		{
			double Fl = Common.Normalize(5000);
			double Fh = Common.Normalize(15000);

			LowPassFir  lpf = new LowPassFir(16, Fl);
			BandPassFir bpf = new BandPassFir(16, Fl, Fh);
			HighPassFir hpf = new HighPassFir(16, Fh);
			PeakingEqualizer peql = new PeakingEqualizer(Common.Normalize(900), 3, 1.5);
			PeakingEqualizer peqm = new PeakingEqualizer(Common.Normalize(8000), 5, 2);
			PeakingEqualizer peqh = new PeakingEqualizer(Common.Normalize(20000), 2, 4);
			Delay dl = new Delay(4);
			Delay dm = new Delay(8);
			Delay dh = new Delay(16);

#if OUTPUT_COEF
			using(StreamWriter writer = new StreamWriter("coef.txt"))
			{
				Write(writer, "LPF", lpf);
				Write(writer, "BPF", bpf);
				Write(writer, "HPF", hpf);
				Write(writer, "PEQ L", peql);
				Write(writer, "PEQ M", peqm);
				Write(writer, "PEQ H", peqh);
				Write(writer, "DLY L", dl);
				Write(writer, "DLY M", dm);
				Write(writer, "DLY H", dh);
			}
#endif

			IFilter low = new SerialConnector(new IFilter[]{lpf, peql, dl});
			IFilter mid = new SerialConnector(new IFilter[]{bpf, peqm, dm});
			IFilter hig = new SerialConnector(new IFilter[]{hpf, peqh, dh});

			IFilter filter = new Mixer(
				new IFilter[]{low, mid, hig},
				new Double[]{2, 1, 0.5});

			return filter;
		}