我是一名新的MATLAB用户,我需要找到一个方程式的四个系数来计算电动机的磁芯损耗。
我已将测量数据绘制在图表上,需要使用这些结果来定义此图表的等式。
核心损失的等式:
根据测量结果,我为 f 的不同值绘制了 Pfe (Bm)的值,但我需要 a , b , e 和 x 。
使用等式:
我可以通过非线性回归分析计算误差的最小值 E ,其中 Pfei 是我的测量值, Pfei *是我的估计值(在这种情况下,我可能会'猜测'系数的初始值)。
如何使用fminsearch
函数计算最小误差值,从而计算出上述系数的值?
答案 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