我有一个名为A 的101x82 大小矩阵。我试图最小化目标函数obj_fun，其值是使用A 间接计算的。
现在为了最小化这个目标函数obj_fun，我需要扰动A 的值。我想检查obj_fun 的值是否下降。如果不是，那么我需要对A 的值进行扰动/更改到某个百分比，以使其最小化obj_fun。继续扰动/改变A 的值，直到我们得到最小值obj_fun。在任何扰动之前，我的 A 平均值约为 1.1529e+003。
有没有人建议我该怎么做？另外，我有点关心速度，即方法/算法不应该太慢。谢谢。
【问题讨论】：
你可以给A添加随机高斯噪声：
A = 0; % seed value for A with something more interesting than 0
best = obj_fun(A);
for iter = 1:max_iter % max_iter should be the maximum number of iterations
newA = A + normrnd(0, 1, size(A));
newObj = obj_fun(A);
if( newObj < best )
best = newObj;
A = newA;
end
end
【讨论】：
normrand 未显示为正确的函数 - normrand not found。另外，我想扰动A 的值，使得结果值的大小也接近A。我A的平均值是~1.1529e+003
normrnd。 normrnd 的前两个参数分别是您想要生成的扰动的平均值和标准差。如果您需要扰动以平均获得该值，则可以将示例中的 0 更改为 1.1529+003。