MATLAB Answers

fitting a scatter plot to power law

조회 수: 1(최근 30일)
Joseph
Joseph 24 Sep 2019
Edited: Star Strider 24 Sep 2019
I have an array A(5,2). how can i fit it to a power law y= c1*x^n1+c2*x^n2 in wich C1, C2, n1 and n2 are real numbers and yet to be determined from the fit?
Thank you

  댓글 수: 0

로그인 to comment.

채택된 답변

Star Strider
Star Strider 24 Sep 2019
Star Strider 님이 편집함. 24 Sep 2019
One approach:
objfcn = @(b,x) b(1).*x.^b(3) + b(2).*x.^b(4);
x = 1:10; % Create Data
y = rand(1, 10); % Create Data
B0 = ones(4,1);
[B,rsdnrm] = fminsearch(@(b) norm(y - objfcn(b,x)), B0);
fprintf(1, 'c_1 = %12.6f\nc_2 = %12.6f\nn_1 = %12.6f\nn_2 = %12.6f\n', B)
Experiment to get the result you want.
EDIT — (24 Sep 2019 at 15:19)
I have no idea how ‘A’ fits into this, since it was not described in any detail.
A guess:
[B,rsdnrm] = fminsearch(@(b) norm(A(:,2) - objfcn(b,A(:,1))), B0);
This fits the second column of ‘A’ to the first column, using the function. Reverse the column designations if necessary.

  댓글 수: 0

로그인 to comment.

More Answers (0)

이 질문에 답변하려면 로그인을(를) 수행하십시오.

태그


Translated by