rsb.cc

用systemC制作的快速付立叶变换

// rsb.cc
// Pierre Bomel - LESTER/UBS - janvier 2005
//    Formation chez France Telecom R&D
#include <iostream.h>
#include <fstream.h>
#include <math.h>

int main(int argc, char *argv[])
{
	ifstream sig1, sig2;
	sig1.open(argv[1], ios::in);
	sig2.open(argv[2], ios::in);
	int t1, t2;
	double s1, s2;
	double errsum = 0.0;
	double sigsum = 0.0;
	int n = 0;

	while (sig1 >> t1) {
		n++;
		sig2 >> t2;
		sig1 >> s1;
		sig2 >> s2;
		//cout << "s1 = " << s1 << " s2 = " << s2 << endl;

		double delta = s1 - s2;
		delta *= delta;
		errsum += delta;

		double ref = s1 * s1;
		sigsum += ref;
	}
	//cout << "Somme du carr