Перевод Matlab @ в код Python (RuntimeWarning: недопустимое значение, встречающееся в true_divide)

Question

Перевод Matlab @ в код Python (RuntimeWarning: недопустимое значение, встречающееся в true_divide)

1

Я перевожу код Matlab на python и застрял здесь. Я думаю, что с @это создает локальную переменную, верно? Я попытался передать код с помощью вложенной функции, но, похоже, возникла проблема с обработкой одной из переменных.

#Code variable:

#LL0 = first tentative value
#Sad and f = returns of psd estimation via welch method
#Sad = is the psd of discrete time-series
#f = is an array of the sampled frequencies
#u_mean = is a scalar

% Matlab Code
f = @(LL, Sad, n, U_media)sum((Sad - 4 * n * LL / U_media .* (1 + 70.8 * (n * LL / U_media).^2).^(-5/6)).^2);
fun = @(LL)f(LL, Sad, n, U_media);
LL = fminsearch(fun, LL0);

def f1(LL, Sad, n, u_mean):
    a = sum((Sad - 4 * n * LL / u_mean * (1 + 70.8 * (n * LL / u_mean)**2)**(-5/6))**2)
    return(a)
f2 = lambda LL, Sad, u_mean, f: f1(LL, Sad, n, u_mean)
fun = lambda LL: f2(LL, Sad, n, u_mean)
LL = scipy.optimize.fmin(func=fun, x0=LL0, maxfun=100000, xtol=1e-6, maxiter=10000, disp=True)

Код работает, но итерации дают мне этот вывод

RuntimeWarning: invalid value encountered in true_divide
  sum((S_adim - 4 * f * LL / u_mean * (1 + 70.8 * (f * LL / u_mean**2)**(-5/6)))**2)
RuntimeWarning: invalid value encountered in reduce
  return ufunc.reduce(obj, axis, dtype, out, **passkwargs)
Warning: Maximum number of iterations has been exceeded.

Значение, которое я получаю взамен, совпадает с тем, которое я использовал в качестве приблизительного значения. Кто-нибудь может мне помочь? заранее спасибо

Giovanni Frison 18 янв. 2019, в 09:19

Источник

0

пожалуйста, укажите, является ли приведенный выше код matlab или python. Это похоже на Matlab, хотя.
Jean-François Fabre♦ 18 янв. 2019, в 07:53
0

Поиск в Google "Matlab Arrobas означает сайт: stackoverflow.com" дал мне 2 отличных ответа. Используйте это в следующий раз, прежде чем задавать тот же вопрос. Спасибо
Jean-François Fabre♦ 18 янв. 2019, в 07:57
0

Привет, спасибо за быстрый ответ. Кстати, я вполне понимаю, что @ делает в matlab, хотя это и казалось хорошим заголовком ... В любом случае, моя проблема заключается в том, как перевести его на python, и ни один из предложенных вами ответов не решит мою проблему.
Giovanni Frison 18 янв. 2019, в 08:06
0

.* кажется умножением матриц. Не уверен, что вы правильно перевели это на python. и в вашем случае @ создает локальную функцию (но вы правильно поняли эту часть)
Jean-François Fabre♦ 18 янв. 2019, в 08:25
0

также «это не похоже на работу» не является правильным вопросом. Пожалуйста, включите минимальный воспроизводимый пример, показывающий, что не работает, чтобы другие могли воспроизвести. Выход из закрытого банкомата.
Jean-François Fabre♦ 18 янв. 2019, в 08:27
0

Спасибо за вашу помощь. Извините, но я здесь новенький. Надеюсь теперь все нормально
Giovanni Frison 18 янв. 2019, в 09:11
0

@Jean-FrançoisFabre Жан-Франсуа Фабр, не могли бы вы хотя бы снять отметку как дубликат?
Giovanni Frison 18 янв. 2019, в 09:35
0

да, потому что после изменений кажется, что это не
Jean-François Fabre♦ 18 янв. 2019, в 09:36
0

@Jean-FrançoisFabre Jean-FrançoisFabre .* - это поэлементное умножение в MATLAB
Ander Biguri 18 янв. 2019, в 11:16
0

Это численная проблема, которую мы не сможем решить без действительных значений переменных. Вы u_mean , на NaN или на ноль, u_mean становится этим значением
Ander Biguri 18 янв. 2019, в 11:18

Показать ещё 8 комментариев

Теги:

python

matlab

scipy

fminsearch

1 ответ

Ещё вопросы

пожалуйста, укажите, является ли приведенный выше код matlab или python. Это похоже на Matlab, хотя.
Поиск в Google "Matlab Arrobas означает сайт: stackoverflow.com" дал мне 2 отличных ответа. Используйте это в следующий раз, прежде чем задавать тот же вопрос. Спасибо
Привет, спасибо за быстрый ответ. Кстати, я вполне понимаю, что @ делает в matlab, хотя это и казалось хорошим заголовком ... В любом случае, моя проблема заключается в том, как перевести его на python, и ни один из предложенных вами ответов не решит мою проблему.
.* кажется умножением матриц. Не уверен, что вы правильно перевели это на python. и в вашем случае @ создает локальную функцию (но вы правильно поняли эту часть)
также «это не похоже на работу» не является правильным вопросом. Пожалуйста, включите минимальный воспроизводимый пример, показывающий, что не работает, чтобы другие могли воспроизвести. Выход из закрытого банкомата.
Спасибо за вашу помощь. Извините, но я здесь новенький. Надеюсь теперь все нормально
@Jean-FrançoisFabre Жан-Франсуа Фабр, не могли бы вы хотя бы снять отметку как дубликат?
да, потому что после изменений кажется, что это не
@Jean-FrançoisFabre Jean-FrançoisFabre .* - это поэлементное умножение в MATLAB
Это численная проблема, которую мы не сможем решить без действительных значений переменных. Вы u_mean , на NaN или на ноль, u_mean становится этим значением

Giovanni Frison · Answer 1 · 2019-01-19T14-52-00.000Z

наконец-то нашел решение; я не очень уверен, почему теперь это работает, думаю, что я назначал одну и ту же переменную дважды. В любом случае, это решение, к которому я пришел:

    def f(LL, Sad, n, u_mean):
        f = sum((Sad - (4 * n * LL / u_mean) *
                 ((1 + 70.8 * ((n * LL / u_mean) ** (2)))**(-5 / 6)))**2)
        return(f)

    fun = lambda LL: f(LL, Sad, n, u_mean)
    res = scipy.optimize.minimize(fun=fun, x0=LL0)
    Lux = res.x

Кстати, spyder бездействует постоянно говорит, что я не должен использовать лямбда-функцию, а вместо этого определить функцию... может кто-нибудь сказать мне, почему? Я смог перевести "веселье" в функцию; Я пробовал с:

def fun(LL):
   f(LL, Sad, n, u_mean)
   return(f)

res = scipy.optimize.minimize(fun=fun, x0=LL0)

но это не работает Было бы хорошо, если кто-нибудь может сказать мне, почему.

Спасибо за вашу помощь