46 #define ABS(x) (((x) > 0) ? (x) : (-(x)))
54 #define M_PI 3.14159265358979323846
76 std::vector<double> code;
80 cerr <<
"file not found" << endl;
81 std::ifstream input (name);
84 while(getline(input, lineData))
86 code.push_back(::atof(lineData.c_str()));
116 double re_abs_hat_alpha=0;
117 double im_abs_hat_alpha=0;
120 for (
unsigned k=0; k<code.size(); k++)
122 im_abs_hat_alpha =im_abs_hat_alpha+code[k]*sin(-xi*(k+0.5));
123 re_abs_hat_alpha =re_abs_hat_alpha+code[k]*cos(-xi*(k+0.5));
128 im_abs_hat_alpha =im_abs_hat_alpha*sin(xi/2)
130 re_abs_hat_alpha =re_abs_hat_alpha*sin(xi/2)
134 return(pow(im_abs_hat_alpha*im_abs_hat_alpha+
135 re_abs_hat_alpha*re_abs_hat_alpha,0.5));
167 pow(sin((xi+epsilon)/2)/((xi+epsilon)/2),2)-
169 pow(sin(xi/2)/(xi/2),2))/epsilon);}
double abs_hat_alpha(std::vector< double > code, double xi)
std::vector< double > fileToVector(const char *name)
Given the file name oof a text file containing a flutter code store the values in a vector of doubles...
double w_prime_estimator(double xi, std::vector< double > code, double epsilon)