MATLAB: Export Fourier series coefficient into workspace variables

cftoolcurve fittingCurve Fitting Toolboxfourier seriestext editing

Hi guys, i need to export the results from a Fourier series a0,a1,b1 into variables of the workspace, in order to further elaborate them. If i use the CFTOOL manually there is no problem, since i can save this data output easily. The problem comes when i try to do it in textcode from the editor, since i have to postprocess a huge amount of data files and was trying to generate a code to do this process automatically. For the Fourier series i am using the following code:

[fitobject]=fit(var1,var2,'fourier1') [fitresult,gof]=fitobject(a0,a1,b1)

with a0,a1,b1 as global var. Thx to those willing to help a disperate newby 😀

Best Answer

 a0 = fitobject.a0;
 a1 = fitobject.a1;
 b1 = fitobject.b1;
 w  = fitobject.w;

Related Solutions

MATLAB: How to export a fourier series in its function form after using the curve-fitting method

save theta_fit(t) fitdata; load fitdata;

MATLAB: Zero values of a fitted curve

There's a paper by Prof. Boyd from the University of Michigan that appears to address this:

"Computing the zeros, maxima and inflection points of Chebyshev, Legendre and Fourier series: solving transcendental equations by spectral interpolation and polynomial rootfinding" http://link.springer.com/article/10.1007%2Fs10665-006-9087-5#page-1

A summary seems to be given here: http://math.stackexchange.com/questions/370996/roots-of-a-finite-fourier-series

I think this does what you want (not that the roots of the fourier series are given in the variable t):

% fourier coefficients
% (dummy data)
a = [1, 1, 1, 1];
b = [0, 0, 1, 0];
% Note:
% a(1) corresponds to a_0 in the reference
% a(2:end) corresponds to a(1:N) in the reference;
N = numel(a) - 1;
h = [a(end:-1:2) + 1i*b(end:-1:2), 2*a(1),  (a(2:end) - 1i*b(2:end))];
B = diag(ones(2*N - 1, 1), 1);
B(end, :) = -h(1:end - 1)/(a(end) - 1i*b(end));
z = eig(B);
t = angle(z);
% note: duplicate t's are imgaginary roots
% eliminate duplicates within tolerance
tol = 10*eps;
t = sort(t);
ind = find(diff(t)./abs(t(2:end)) < tol);
ind = [ind; ind + 1];
t(ind(:)) = [];
% plot results
% f = ugly inline fourier series anonymous function
f = @(x, a, b) sum(bsxfun(@times, a(:).', cos(x(:)*(0:numel(a) - 1))) + bsxfun(@times, b(:).', sin(x(:)*(0:numel(b)-1))), 2);
x = linspace(-pi, pi, 128 + 1);
hf = figure; 
ha = axes;
axis([-pi, pi, -sum(a+b)-0.5, sum(a+b)+0.5]);
hold(ha, 'on');
plot(x, f(x, a, b));
plot(t, f(t, a, b), 'o');
plot([-pi, pi], [0, 0], '--c');

I probably won't be around much in the next two days, but feel free to leave a comment, and I'll try to get back to you.

Best Answer

Related Solutions

MATLAB: How to export a fourier series in its function form after using the curve-fitting method

MATLAB: Zero values of a fitted curve

Related Question