适合参数

Question

我是Matlab的新手，我正在尝试使用我在网上找到的代码。 我试图拟合HydrodynamicSpectrum描述的图表。但是在输入fvA和fmA之后，我没有考虑它，我也试图获得这个值的拟合参数。

我尝试删除它们，更改它们。但没有一个工作。我想知道这里是否有人能够指出我正确的方向来解决这个问题。

specFunc = @(f, para)HydrodynamicSpectrum(f, [para fvA fmA]);
[fit.AXfc, fit.AXD] = NonLinearFit(fit.f(indXY), fit.AXSpec(indXY), specFunc, [iguess_AXfc iguess_AXD]);
[fit.AYfc, fit.AYD] = NonLinearFit(fit.f(indXY), fit.AYSpec(indXY), specFunc, [iguess_AYfc iguess_AYD]);
[fit.ASumfc, fit.ASumD] = NonLinearFit(fit.f(indSum), fit.ASumSpec(indSum), specFunc, [iguess_ASumfc iguess_ASumD]);
predictedAX = HydrodynamicSpectrum(fit.f, [fit.AXfc fit.AXD fvA fmA]);
predictedAY = HydrodynamicSpectrum(fit.f, [fit.AYfc fit.AYD fvA fmA]);
predictedASum = HydrodynamicSpectrum(fit.f, [fit.ASumfc fit.ASumD fvA fmA]);

function spec = HydrodynamicSpectrum(f, para);

fc = para(1);
D = para(2);
fv = para(3);
fm = para(4);
f = abs(f);  %Kludge!
spec = D/pi^2*(1+sqrt(f/fv))./((fc - f.*sqrt(f./fv) - (f.^2)/fm).^2 + (f + f.*sqrt(f./fv)).^2);


function [fc, D, sfc, sD] = NonLinearFit(f, spec, specFunc, init);

func = @(para, f)spec./specFunc(f, para);

[paraFit, resid, J] = nlinfit(f, ones(1, length(spec)), func, init);

fc = paraFit(1);
D = paraFit(2);

ci = nlparci(real(paraFit), real(resid), real(J));  % Kludge!!

sfc = (ci(1,2) - ci(1,1))/4;
sD = (ci(2,2) - ci(2,1))/4;

Answer 1

[paraFit, resid, J] = nlinfit(f, ones(1, length(spec)), func, init);

看起来您使用此行获得了拟合参数。而且你正在进一步处理它们以从函数中获取其他东西。你也可以修改你的第二个功能来解决它​​们。

由于评论很少，而且问题似乎是针对特定应用的，因此我对您提供的内容没有多大帮助。