fourierseriesapproach.cpp

来自「何光渝所著《计算机常用数值算法与程序(C++版)》的随书源代码。 是想学习C+」· C++ 代码 · 共 34 行

//FourierSeriesApproach.cpp	傅里叶级数逼近

#include <iostream>			//模板类输入输出流标准头文件
#include <Transform.h>		//数学变换头文件
using namespace std;		//名字空间

void main(void)
{	
	valarray<long double> a(51);
	valarray<long double> b(51);
	valarray<long double> f(101);
	
	long double h = 6.283185306 / 101.0;

	for(int i=0; i<101; i++)
    {
		long double ld = (i + 0.5) * h;
		f[i] = ld * ld;
	}
	
  	cout.setf(ios::fixed);	//输出数据为定点法
	cout.precision(8);		//精度8位

    FourierSeriesApproach(f, a, b);

	cout << endl;

	for (i=0; i<51; i++)
		cout << "a(" << i << ") = " << a[i]
						<< "\t b(" << i << ") = " << b[i] << endl;
}