MATLAB fminsearch方程使用四个匿名参数

时间:2013-04-26 15:28:20

标签: matlab fminsearch

我是一名新的MATLAB用户,我需要找到一个方程式的四个系数来计算电动机的磁芯损耗。

我已将测量数据绘制在图表上,需要使用这些结果来定义此图表的等式。

核心损失的等式:

enter image description here

根据测量结果,我为 f 的不同值绘制了 Pfe (Bm)的值,但我需要 a , b e x

使用等式:

enter image description here

我可以通过非线性回归分析计算误差的最小值 E ,其中 Pfei 是我的测量值, Pfei *是我的估计值(在这种情况下,我可能会'猜测'系数的初始值)。

如何使用fminsearch函数计算最小误差值,从而计算出上述系数的值?

1 个答案:

答案 0 :(得分:0)

您可以执行以下代码之类的操作。向量fitParams将具有a,b,e和x的值。

function [fitParams, fval] = callMinEx()
Bm = 1:100;

% fake data.  Use your actual data for pfei
x0 = [2,3,4,5];
pfei = coreLoss(Bm, x0) + 5e10*rand(1,length(Bm));

% call fminsearch with some initial parameters
startParams = [3,3,3,3];
[fitParams,fval] = fminsearch(@func2Minimize, startParams);

% plot data and fit
plot(Bm,pfei,'r*');
hold on
plot(Bm, coreLoss(Bm, x0));
legend('data','fit')

     % The equation to be minimized by fminsearch
    function epsValue = func2Minimize(params)
        pfeStar = coreLoss(Bm,params);
        epsValue = sum(((pfei - pfeStar) ./ pfei).^ 2);
    end

    % core loss function
    function pfei = coreLoss(Bm, x0)
        a = x0(1);
        b = x0(2);
        e = x0(3);
        x = x0(4);

        f = 100;

        pfei =  a * f * Bm .^ x + b * (f * Bm).^2 + e * (f * Bm).^1.5;
    end

end