Parameter optimization for a model that relates two datasets

Question

Andrew Johnson 2014-10-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/157038-parameter-optimization-for-a-model-that-relates-two-datasets

评论： Andrew Johnson 2014-10-8

Hello! This is a seemingly simple problem I have been stuck on for some time now. I would greatly appreciate any ideas or help!

I have a model of the form:

Y = A*X^-1 + B*X^4*F(X)

-----------

A and B are unknown parameters.

Y is a measured frequency distribution dataset of the form (Y(i), f(Y(i))).

X is another measured frequency distribution dataset (X(i), g(X(i))).

F is a known (albeit complicated) function of X.

-----------

I would like to find values for A and B such that the right hand side of the function transforms the frequency distribution (X,g(X)) into something that looks very nearly like the measured (Y,f(Y)).

I have tried a nested for-loop "brute-force" method which tests all combinations of A and B within reasonable bounds, but it has proven too sensitive to the discretization scale chosen for A and B. Is there any matlab solution to this sort of problem? Thank you in advance!

Andrew

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

请先登录，再进行评论。

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

Answer 1

Jeremy Kemmerer 2014-10-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/157038-parameter-optimization-for-a-model-that-relates-two-datasets#answer_153696

Do you know for certain that the parameters A and B exist? It sounds like you may have an underdetermined system.

One quick thing to try is to form a matrix Z with three columns: Y, X^-1, and X^4*F(X) and check the rank of the matrix using the function “ rank ”.

If Z is full rank (rank is equal to 3 in this case), this tells you that the columns are not linearly dependent, and hence the parameters you are looking for won’t exist.

For more information on the command “ rank ”, please see the documentation pages: http://www.mathworks.com/help/matlab/ref/rank.html