Нелинейные припадки в Октаве - PullRequest
Новая
       34

Нелинейные припадки в Октаве

0 голосов
/ 09 июня 2018

В настоящее время я использую функцию nonlin_curvefit из пакета 'Optim' GNU Octave, чтобы соответствовать данным.Но в этот раз мне также понадобилась неопределенность возвращаемых параметров, чтобы определить качество подгонки.После прочтения документации я связал с помощью функции curvefit_stat.Однако всякий раз, когда я получаю ошибки, используя эту функцию, я не могу понять смысл сообщения об ошибке.Я использую Octave 4.2.2 из репозитория Ubuntu 18.04 по умолчанию.

Автономный минимальный пример и сообщения об ошибках ниже.Стартовые параметры init_cvg обычно дают хороший результат, в то время как использование init_dvg обычно приводит к плохим подгонкам: 

1;

x = linspace(-2*pi, 2*pi, 600);
ydata = 3.4*sin(1.6*x) + rand(size(x))*1.3;

f = @(p, x) p(1)*sin(p(2).*x);
function y  = testfun(p ,x)
  y = p(1).*sin(p(2).*x);
endfunction

init_cvg = [1; 1.1];
init_dvg = [1; 1.0];

[pc, mod_valc, cvgc, outpc] = nonlin_curvefit(f, init_cvg, x, ydata);
[pd, mod_vald, cvgd, outpd] = nonlin_curvefit(f, init_dvg, x, ydata);


hold off
plot(x, yd, "b.");
hold on;
plot(x, mod_valc, "r-");
plot(x, mod_vald, "color", [0.9 0.4 0.3]);

settings = optimset("ret_covp", true);
covpc = curvefit_stat(f, pc, x, ydata, settings);
covpd = curvefit_stat(f, pd, x, ydata, settings);

puts("sqrt(diag(covpc))")
sqrt(diag(covpc))

puts("sqrt(diag(covpd))")
sqrt(diag(covpd))

Первое сообщение об ошибке появляется, когда я использую f в качестве функции модели, второе - когда яиспользуйте testfun вместо: 

>> curvefit_stat_TEST
error: label of objective function must be specified
error: called from
    __residmin_stat__ at line 566 column 7
    curvefit_stat at line 56 column 7
    curvefit_stat_TEST at line 25 column 7
>> curvefit_stat_TEST
error: 'p' undefined near line 8 column 7
error: called from
    testfun at line 8 column 5
    curvefit_stat_TEST at line 25 column 7
>>

Может ли кто-нибудь подтвердить эту ошибку?

Буду признателен за любую помощь.

1 Ответ

0 голосов
/ 11 июня 2018

Я нашел проблему.Мне нужно было добавить "abjf_type", "wls" в качестве аргументов для optimset.

...