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!