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!