How to find parameters of the following coefficient based on optimization tool?

Question

Parham Babakhani Dehkordi 2016-4-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/278434-how-to-find-parameters-of-the-following-coefficient-based-on-optimization-tool

回答： Parham Babakhani Dehkordi 2016-4-12

Hi everyone, I have the following input and function, my purpose is to calculate the value of a,m and n which are scalars using curve fitting in matlab. I have:

ew=[0.679 0.613 0.613 0.501 0.497 0.50 0.389 0.386 0.386 0.391 0.390 0.288 0.288 0.288 0.289 0.289]; Hw=[0.757 0.669 0.71 0.577 0.57 0.60 0.50 0.55 0.52 0.50 0.476 0.353 0.388 0.379 0.431 0.371]; k=[0.217 0.217 0.217 0.163 0.163 0.163 0.131 0.131 0.131 0.131 0.131 0.109 0.109 0.109 0.109 0.109];

Hw=[ew*(1+0.36(1-ew))]*C; C=exp[-a*(k^m)*(1-ew)^n];

any help would be appreciated.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Alan Weiss 2016-4-11

Please format your question using the {} Code button so that we might read it more easily.

Alan Weiss

MATLAB mathematical toolbox documentation

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Roger Stafford 2016-4-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/278434-how-to-find-parameters-of-the-following-coefficient-based-on-optimization-tool#answer_217496

在 MATLAB Online 中打开

Define the following sum of squares of differences:

   ew=[0.679 0.613 0.613 0.501 0.497 0.50  0.389 0.386 0.386 0.391 0.390 0.288 0.288 0.288 0.289 0.289]; 
   Hw=[0.757 0.669 0.71  0.577 0.57  0.60  0.50  0.55  0.52  0.50  0.476 0.353 0.388 0.379 0.431 0.371]; 
   k= [0.217 0.217 0.217 0.163 0.163 0.163 0.131 0.131 0.131 0.131 0.131 0.109 0.109 0.109 0.109 0.109];
   function s = myfun(a,m,n)
     s = sum((Hw-ew.*(1+0.36(1-ew)).*exp[-a*(k.^m).*(1-ew).^n)).^2);
   return

and then use that subfunction for optimization of a, m, and n using 'fminsearch' or 'fminunc'. Their descriptions can be found at:

   http://www.mathworks.com/help/matlab/ref/fminsearch.html
   http://www.mathworks.com/help/optim/ug/fminunc.html

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Parham Babakhani Dehkordi 2016-4-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/278434-how-to-find-parameters-of-the-following-coefficient-based-on-optimization-tool#answer_217600

在 MATLAB Online 中打开

Thank you Roger for your help,but I still have a problem to get a solution, I arranged my code with initial guesses for a=0.31,m=1 and n=0.7 and also I considered a=x(1), m=x(2) and n=x(3) as follows:

for i=1:16
x =fminsearch(@(x) (Hw(i)-(ew(i)*(1+0.36*(1-ew(i)))))*exp(-x(1)*k(i)^x(2)*(1-ew(i))^x(3)), [0.31 1.0 0.7])
end
there is an error which says Error: File: fminsearch.m Line: 31 Column: 5
Function definitions are not permitted in this context.

Actually I could not properly use handle function directly to pass it through function function (in this case fminsearch). Sorry Alan, I am new user to Matlab.