fminsearch for multi-input and single output
조회 수: 8 (최근 30일)
이전 댓글 표시
I have unknown function with two inputs ('x1' and 'x2') and one output ('Y') and a constant parameter 'p'. Each input is a vector of order [3000 x 3] and a output vector with order [3000 x 3]. Constant is [1 x 1]
best guess of unknown function is Y = (1-p)* 0.003* x1 + 275*p* x2.^4
i want to find a optimum value of parameter 'p' for which the error between the (Ymes - ) is minimum.
Ymes , X1 and X2 are known.
i tried fminsearch but it is not working
Ymes = rand(3000,3);
x1= rand(3000,3);
x2 = rand(3000,3);
fun = @(p) sum((Ymes - (1-p)*0.003* x1 + 275*p*x2.^4)).^2;
pguess = 1500;
[p,fminres] = fminsearch(fun,pguess)
댓글 수: 0
답변 (1개)
Torsten
2018년 2월 28일
편집: Torsten
2018년 2월 28일
Ymes = rand(3000,3);
x1= rand(3000,3);
x2 = rand(3000,3);
Ymes = Ymes();
x1 = x1();
x2 = x2();
fun = @(p) sum((Ymes - ((1-p)*0.003* x1 + 275*p*x2.^4)).^2);
%starting guess
pguess = [0.15];
%optimise
[p,fminres] = fminsearch(fun,pguess)
or simply:
Ymes = rand(3000,3);
x1= rand(3000,3);
x2 = rand(3000,3);
Ymes = Ymes();
x1 = x1();
x2 = x2();
p = (-0.003*x1+275*x2.^4)\(Ymes-0.003*x1)
댓글 수: 7
Torsten
2018년 3월 1일
My guess was that you want to fit
Ymes(i,j)
against
(1-p)*0.003* x1(i,j) + 275*p*x2(i,j).^4
for 1<=i<=3000 and 1<=j<=3.
This is what the code I suggested does.
If you want to do something else, you will have to clarify.
Best wishes
Torsten.
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!