How do I create a series of variables (V0, V1, ...Vn) which act as independent variables in conjunction with a series of sums?

2 次查看（过去 30 天）
Domantas Laurinavicius
编辑： Domantas Laurinavicius 2021-2-4
Basically, I have a symsum which I intend on feeding into a fittype() function.
The fittype is of the form (V0/2) + (V1/2)*cos(1*x) + (V2/2)*cos(2*x) +... (Vn/2)*cos(n*x) where each Vn is an independent variable in a non-linear least squares.
Those Vn variables are going into a matrix so the values need to be extractable, right now I'm using coeffvalues(fitresult) or manually getting them out with fitresult.V0, V1 etc (that could be a question on its own).
I have the sum of cos funcs
syms x n
Func = symsum((1/2)*cos(n*x),n,[1 5])
but I can't figure out how to (1) create an array of variables V0:Vn and outputting them as variables and (2) how to combine those with a symsum output into the fittype given above.

请先登录，再回答此问题。

采纳的回答

Chaitanya Mallela
编辑：Chaitanya Mallela 2021-2-4
fittype function accepts character array as input argument but the symsum function gives symbolic variable. To apply fittype to this function you need to split the symsum expression into terms and convert them to character array and generate independent variable Vn as coefficients to fittype function.
syms x n
Func = cell2sym(children(symsum((1/2)*cos(n*x),n,[1 5])));
g = fittype(arrayfun(@char,Func,'UniformOutput',false),'coefficients',arrayfun(@char,sym('V',[1,5]),'UniformOutput',false))
  1 个评论
Domantas Laurinavicius
Domantas Laurinavicius 2021-2-4
编辑：Domantas Laurinavicius 2021-2-4
That works wonderfully, thank you very much! Now, Why is the output all out of order?
g =
Linear model:
g(V1,V2,V3,V4,V5,x) = V1*cos(2*x)/2 + V2*cos(3*x)/2 + V3*cos(4*x)/2 + V4*cos(x)/2 + V5*1/2
I would've expected V1 to be paired with cos(1*x) and V2 with cos(2*x) but matlab seems to preferentially place the lowest index term last.
Regardless, when I put this equation as the fittype(), I get the message:
Warning: The given fit options specify a NonlinearLeastSquares fit but the model specifies a LinearLeastSquares fit.
Any hints there?
Thanks again!
Edit: thought the full context would help. I'm inputting that linear model to this set of cftool code:
ft = fittype( '(V0/2) + (V1/2)*cos(1*x) + (V2/2)*cos(2*x) + (V3/2)*cos(3*x) + (V4/2)*cos(4*x) + (V5/2)*cos(5*x);', 'independent', 'x', 'dependent', 'y' )
opts = fitoptions( 'Method', 'NonlinearLeastSquares' );
opts.Display = 'Off';
opts.StartPoint = [0.421761282626275 0.915735525189067 0.792207329559554 0.959492426392903 0.655740699156587 0.0357116785741896];
% Fit model to data
[fitresult, gof] = fit( xData, yData, ft, opts );

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Linear and Nonlinear Regression 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by