回 帖 发 新 帖 刷新版面

主题:请教关于最小二乘的一些问题

各位高手,问题的描述如下:

定义函数文件 function y=Qa(np)
请注意,返回的y是一个数组,即np函数的数组。数组中元素的个数大于1,现在我想求Qa最小时的np,但是,返回的np应该是和Qa同样大小的数组,即对应一个y中的函数,求一次最小,返回一个np。请问应该用哪个函数,fminsearch不行,它不能运算函数矩阵。请各位高手赐教。


注:不能使用符号变量,因为我在Qa函数中已经使用了符号变量x,然后利用了inline函数,如果把np定义为符号变量,那么Qa函数中的inline后获得的函数,有两个变量x和np,无法进一步运算。

谢谢各位!

回复列表 (共1个回复)

沙发

条件有点含糊,np是数组,我猜想你的意思应该是每一行代表一个自变量,每一列是一个取值序列。
经过你Qa函数映射后的解集也应该是同样大小的数组。问题是现在是在解集中取出每一行的最小的数,所对应自变量的值。
这样描述完不就是很简单的问题吗?
[minValue minIndex]=min(y');  %注意min在处理矩阵时,是按列处理的,所以要倒置
for i=1:size(np,1)
     minX(i)=np(i,minIndex(i));    %minIndex是联系解集和自变量集合的下标组合
end


我来回复

您尚未登录,请登录后再回复。点此登录或注册