iN" , d Z ddlZddlmZ g dZd Zd Zd Zd Z d Z
d
Z d Z d Z
d
Zd Zd Z G d de Z e Zd Z G d de Z e Zd Zd Zd Zd Zd Zd Zd Z d Z G d d e Z e Z G d d e Z e Z!y)!zI Collection of Model instances for use with the odrpack fitting package.
N)Model)r exponential multilinear unilinear quadratic
polynomialc | d | dd }}|j |j d df }|||z j d z S Nr axis)reshapeshapesum)Bxabs f/Volumes/fast/ai/experiments/voice-extract-mac/.venv/lib/python3.12/site-packages/scipy/odr/_models.py_lin_fcnr
sI Q412qA 1771:q/"A !yyay c t j |j d t }t j ||j f }|j
| j d |j d f S N)nponesr float concatenateravelr )r r r ress r _lin_fjbr! sV
U#A
..!QWWY
(C ;; QWWR[1 22r c | dd }t j ||j d f|j d z d }|j |j S )Nr r r r )r repeatr r )r r r s r _lin_fjdr$ sI !"A
!aggbk^AGGBK/a8A 99QWW r c t | j j dk( r| j j d }nd}t j |dz ft
S N r r )lenr r r r r )datams r _lin_estr+ sF
466<<A FFLLO
77AE8U ##r c | d | dd }}|j |j d df }|t j |t j || z d z S r
r r r r power)r r powersr r s r _poly_fcnr0 * sU Q412qA 1771:q/"A rvva"((1f--A6 66r c t j t j |j d t t j
|| j f }|j | j d |j d f S r )r r r r r r. flatr )r r r/ r s r _poly_fjacbr3 1 sb
.."''!''"+u5 ((1f-224 5C ;; QWWR[1 22r c | dd }|j |j d df }||z }t j |t j ||dz
z d S )Nr r r r- )r r r/ r s r _poly_fjacdr5 7 sU !"A 1771:q/"A F
A
66!bhhq&(++! 44r c D | d t j | d |z z S Nr r r expr r s r _exp_fcnr; @ " Q4"&&1" ""r c D | d t j | d |z z S )Nr r8 r: s r _exp_fjdr> D r< r c t j t j |j d t |t j
| d |z z f }|j
d|j d f S )Nr r r' )r r r r r r9 r )r r r s r _exp_fjbr@ H sZ
.."''!''"+u5q266!A$(;K7KL
MC ;;1772;' ((r c 0 t j ddg S )N ?)r arrayr) s r _exp_estrE M s
88RH r c " e Zd ZdZ fdZ xZS )_MultilinearModela
Arbitrary-dimensional linear model
.. deprecated:: 1.17.0
`scipy.odr` is deprecated and will be removed in SciPy 1.19.0. Please use
`pypi.org/project/odrpack/ `_
instead.
This model is defined by :math:`y=\beta_0 + \sum_{i=1}^m \beta_i x_i`
Examples
--------
We can calculate orthogonal distance regression with an arbitrary
dimensional linear model:
>>> from scipy import odr
>>> import numpy as np
>>> x = np.linspace(0.0, 5.0)
>>> y = 10.0 + 5.0 * x
>>> data = odr.Data(x, y)
>>> odr_obj = odr.ODR(data, odr.multilinear)
>>> output = odr_obj.run()
>>> print(output.beta)
[10. 5.]
c
V t | t t t t
dddd y )Nz Arbitrary-dimensional Linearz y = B_0 + Sum[i=1..m, B_i * x_i]z&$y=\beta_0 + \sum_{i=1}^m \beta_i x_i$nameequTeXequ)fjacbfjacdestimatemeta)super__init__r r! r$ r+ self __class__s r rR z_MultilinearModel.__init__o s.
HHx8; EG Hr __name__
__module__ __qualname____doc__rR
__classcell__rU s @r rG rG R s 8H Hr rG c 6 t j | }|j dk( rt j d|dz }|j t | df }t | dz }|fd}t
t t t ||fdd|dz
z d|dz
z d S ) az
Factory function for a general polynomial model.
.. deprecated:: 1.17.0
`scipy.odr` is deprecated and will be removed in SciPy 1.19.0. Please use
`pypi.org/project/odrpack/ `_
instead.
Parameters
----------
order : int or sequence
If an integer, it becomes the order of the polynomial to fit. If
a sequence of numbers, then these are the explicit powers in the
polynomial.
A constant term (power 0) is always included, so don't include 0.
Thus, polynomial(n) is equivalent to polynomial(range(1, n+1)).
Returns
-------
polynomial : Model instance
Model instance.
Examples
--------
We can fit an input data using orthogonal distance regression (ODR) with
a polynomial model:
>>> import numpy as np
>>> import matplotlib.pyplot as plt
>>> from scipy import odr
>>> x = np.linspace(0.0, 5.0)
>>> y = np.sin(x)
>>> poly_model = odr.polynomial(3) # using third order polynomial model
>>> data = odr.Data(x, y)
>>> odr_obj = odr.ODR(data, poly_model)
>>> output = odr_obj.run() # running ODR fitting
>>> poly = np.poly1d(output.beta[::-1])
>>> poly_y = poly(x)
>>> plt.plot(x, y, label="input data")
>>> plt.plot(x, poly_y, label="polynomial ODR")
>>> plt.legend()
>>> plt.show()
r c 8 t j |ft S )N)r r r )r) len_betas r _poly_estz polynomial.._poly_est s ww{E**r zSorta-general Polynomialz$y = B_0 + Sum[i=1..%s, B_i * (x**i)]z)$y=\beta_0 + \sum_{i=1}^{%s} \beta_i x^i$rI )rN rM rO
extra_argsrP )
r asarrayr aranger r( r r0 r5 r3 )orderr/ r` ra s r r r z s \ ZZ
F
||r1fqj)
^^S[! ,
-F6{QH!) + +[# 9>(1*MG!! %& ' 'r c " e Zd ZdZ fdZ xZS )_ExponentialModela
Exponential model
.. deprecated:: 1.17.0
`scipy.odr` is deprecated and will be removed in SciPy 1.19.0. Please use
`pypi.org/project/odrpack/ `_
instead.
This model is defined by :math:`y=\beta_0 + e^{\beta_1 x}`
Examples
--------
We can calculate orthogonal distance regression with an exponential model:
>>> from scipy import odr
>>> import numpy as np
>>> x = np.linspace(0.0, 5.0)
>>> y = -10.0 + np.exp(0.5*x)
>>> data = odr.Data(x, y)
>>> odr_obj = odr.ODR(data, odr.exponential)
>>> output = odr_obj.run()
>>> print(output.beta)
[-10. 0.5]
c
V t | t t t t
dddd y )N Exponentialzy= B_0 + exp(B_1 * x)z$y=\beta_0 + e^{\beta_1 x}$rI rN rM rO rP )rQ rR r; r> r@ rE rS s r rR z_ExponentialModel.__init__ s.
"*'4&=)G I Jr rV r\ s @r rg rg 4J Jr rg c || d z | d z S r7 r^ r: s r _unilinrm s QqT6AaD=r c V t j |j t | d z S )Nr )r r r r r: s r _unilin_fjdro s
77177E "QqT ))r c t j |t j |j t f }|j d|j z S )N)r' r r r r r r r r _rets r _unilin_fjbrt s;
>>1bggaggu 56 7D <<qww ''r c y)N)rB rB r^ rD s r _unilin_estrv s r c 0 ||| d z | d z z | d z S )Nr r r' r^ r: s r
_quadraticrx s& a!fqtm qt ##r c $ d|z | d z | d z S r& r^ r: s r _quad_fjdrz s Q3qt8ad?r c t j ||z |t j |j t f }|j d|j z S )N) rq rr s r _quad_fjbr} sA
>>1Q3277177E#:; `_
instead.
This model is defined by :math:`y = \beta_0 x + \beta_1`
Examples
--------
We can calculate orthogonal distance regression with an unilinear model:
>>> from scipy import odr
>>> import numpy as np
>>> x = np.linspace(0.0, 5.0)
>>> y = 1.0 * x + 2.0
>>> data = odr.Data(x, y)
>>> odr_obj = odr.ODR(data, odr.unilinear)
>>> output = odr_obj.run()
>>> print(output.beta)
[1. 2.]
c
V t | t t t t
dddd y )NzUnivariate Linearzy = B_0 * x + B_1z$y = \beta_0 x + \beta_1$rI rj )rQ rR rm ro rt rv rS s r rR z_UnilinearModel.__init__ s.
;"-':&9)F H Ir rV r\ s @r r r s 4I Ir r c " e Zd ZdZ fdZ xZS )_QuadraticModela
Quadratic model
.. deprecated:: 1.17.0
`scipy.odr` is deprecated and will be removed in SciPy 1.19.0. Please use
`pypi.org/project/odrpack/ `_
instead.
This model is defined by :math:`y = \beta_0 x^2 + \beta_1 x + \beta_2`
Examples
--------
We can calculate orthogonal distance regression with a quadratic model:
>>> from scipy import odr
>>> import numpy as np
>>> x = np.linspace(0.0, 5.0)
>>> y = 1.0 * x ** 2 + 2.0 * x + 3.0
>>> data = odr.Data(x, y)
>>> odr_obj = odr.ODR(data, odr.quadratic)
>>> output = odr_obj.run()
>>> print(output.beta)
[1. 2. 3.]
c
V t | t t t t
dddd y )N Quadraticzy = B_0*x**2 + B_1*x + B_2z&$y = \beta_0 x^2 + \beta_1 x + \beta_2rI rj )rQ rR rx rz r} r rS s r rR z_QuadraticModel.__init__E s.
iy9%5 GI Jr rV r\ s @r r r * rk r r )"rZ numpyr scipy.odr._odrpackr __all__r r! r$ r+ r0 r3 r5 r; r> r@ rE rG r r rg r rm ro rt rv rx rz r} r r r r r r^ r r r s $ !3
$73 5##)
"H "HJ ! ?'D J JF ! *(
$(
Ie IF
Je JF
r