I have constructed a curve fit "func" with the following line.
func=fit(x,y,'pchipinterp')
Now I want to use "func" inside an anonymous function as e.g
func2= @ (x,c) x-(c*x+x^2)/func(x)
Then I want to find the root of func2=0 using e.g FZERO as
c=2;
func3=@ (x) func2(x,c)
sol=fzero(func3,<interval>)
Could it be causing any trouble that i am using a cfit object inside the function "func2" ? My solution output is making me suspicious. Note that func2 is not the actual function i'm using but of a similar form.
Any help is much appreciated!
Best regards

댓글 수: 1

Adam Danz
Adam Danz 2020년 4월 29일
"I am not sure this is a correct way of doing it ..."
Doing what?
"...because the solution is not what at all what you would expect?"
What is the solution and what do you expect?
We obviously cannot run the code nor do we know what the solution is or what you expect the solution to be.

댓글을 달려면 로그인하십시오.

 채택된 답변

Adam Danz
Adam Danz 2020년 4월 29일

0 개 추천

My blind guess is that x or c is a vector in which case,
func2= @ (x,c) x-(c.*x+x.^2)./func(x)
% ^ ^ ^

댓글 수: 4

Philip Croné
Philip Croné 2020년 4월 29일
I am wondering whether or not you can use a cfit object the way I have tried to illustrate.
Best regards
Adam Danz
Adam Danz 2020년 4월 29일
As long as the cfit obj is defined prior to the anonymous function, there are no problems with that. You can think of the cfit obj as a function. You're providing it with inputs and it produces outputs.
Whatever you're seeing that makes you suspicious, it's not because you're using a cfit obj within an anonymous function.
Adam Danz
Adam Danz 2020년 4월 29일
편집: Adam Danz 2020년 4월 29일
Whenever you have a suspicion, simplify the problem and test it with known inputs and outputs.
% Create a cfit object
f = fittype('a*x');
c = cfit(f,2); % This just multiplies x by 2
% Create anon function that uses the cfit obj
fun = @(a,x) c(x) + a; % multiply x by 2 and add a
% Test it with known inputs
fun(10,[2 4 6])
ans =
14 % 2*2 + 10
18 % 4*2 + 10
22 % 6*2 + 10
Philip Croné
Philip Croné 2020년 4월 29일
Thank you!
Best regards

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 Fit Postprocessing에 대해 자세히 알아보기

제품

릴리스

R2017b

질문:

2020년 4월 29일

편집:

2020년 4월 29일

Community Treasure Hunt

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

Start Hunting!

Translated by