ftj<FddlZddlZddlmZmZmZddlZddlm Z ddlm Z m Z m Z mZmZddlmZmZmZmZmZejjZee GddZedd Gd d Zee Gd dZdS)N)assert_array_almost_equalxp_assert_closemake_xp_test_case) cont2discrete)dlsimss2tfss2zpklsimlti)tf2ssimpulsedimpulsestepdstepceZdZdZdZdZdZdZdZdZ dZ d Z e d d Z d Ze d dZe d dZdS)TestC2Dc4|d|j}|dd|j}|ddgddgddgg}|dgdgd gg}d |dz}|dd }d}t ||||f|d \} } } } } t || t || t || t || t j| |dsJdS)N)dtyper????Qտ)a?g6 <S?zohmethod+=abs_tol)eyefloat64fullasarrayc2drmathisclose)selfxpacbcccdcad_truthbd_truth dt_requestedadbdcddddts c/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/scipy/signal/tests/test_cont2discrete.pytest_zohzTestC2D.test_zohs, VVARZV ( ( WWVS W 3 3 ZZ$sCj3+> ? ? ZZ#w/ 0 0$rvvayy0776#455  "b"b!1<NNNBB!(B///!(B///!"b)))!"b)))|B e<<<<<<<<cR|d}|dd}|ddgddgddgg}|dgdgdgg}d |dz}|dd }|}|d gd gd gg} d} t||||f| d\} } } }}t || t || t || t | |t j|| dsJdS)Nrrrrrrrrrgm?g>."?g4pL ?g.Ȫq¿fohr r"r#r%r'r(r)rr*r+r,r-r.r/r0r1r2r3cd_truthdd_truthr4r5r6r7r8r9s r:test_fohzTestC2D.test_foh(sQ VVAYY WWVS ! ! ZZ$sCj3+> ? ? ZZ#w/ 0 0%rvvayy0776#455:: 120112455  "b"b!1<NNNBB!(B///!(B///!(B///!(B///|B e<<<<<<< ? ? ZZ#u- . .%rvvayy0776#455::x#9::  "b"b!1<(1333BB "(B///!(B///!(B///!(B///|B e<<<<<<< ? ? ZZ#w/ 0 0 ?7763''::Sz #Sz #Sz+,,::w #u &x)**!"b"b!1<(-U<<<BB "(B///!(B///!(B///!(B/////r<cX|d}|dd}|ddgddgddgg}|dgdgdgg}d}d |dz}|dd}|ddgddgddgg} |} t||||f|d \} } } }}t || t || t | | t | |t j||d sJdS)Nrrrrrrrr?eulerr r"r#r?r,r-r.r/r0r1r4r2r3rArBr5r6r7r8r9s r: test_eulerzTestC2D.test_eulerosX VVAYY WWVS ! ! ZZ$sCj3+> ? ? ZZ#w/ 0 0 ?7764((::c{ #Sz #T{,-- "b"b!1<(/111BB "(B///!(B///!(B///!(B///|B e<<<<<<< ? ? ZZ#w/ 0 0 ?7763''::Sz #Sz #Sz+,,::w #u %w())!"b"b!1<(7999BB "(B///!(B///!(B///!(B/////r<cB|d}|dd}|ddgddgddgg}|dgdgdgg}d}d |dz}|dd }|dd gd d gd d gg} |d gd gd gg} t||||f|d\} } } }}t || t || t | | t | |t j||dsJd|dz}|dd}|ddgddgddgg} |dgdgdgg} d }t||||f|d\} } } }}t || t || t | | t | |t j||sJdS)Nrrrrrrrrg?rGgUUUUUU?g?g>X%bilinearr r"r#gffffff?rLrIrJrHrKrMr?rVs r: test_bilinearzTestC2D.test_bilinearsz VVAYY WWVS ! ! ZZ$sCj3+> ? ? ZZ#w/ 0 0 *7769--::Y/ )95 )95788:: 12 ){ 23566!"b"b!1<(2444BB "(B///!(B///!(B///!(B///|B e<<<<<<?7763''::SzC:SzBCC::wx899  "b"b!1<(2444BB "(B///!(B///!(B///!(B///|L"-------r<cb|gd}|gd}|gdg}|gd}d}t||f|d\}}} t||t||tj|| sJdS)N)rrr)rrr)rGgNZۿgϱ)_?)rg]Ogoh?rrr )r(r)rr*r+) r,r-numcdencnumddendr4numdenr9s r:test_transferfunctionzTestC2D.test_transferfunctionszz+++,,zz+++,,zzMMMNOOzzFFFGG D$<eDDD S"!$,,,!$,,,|L"-------r<c|ddg}|d|dz d|dz g}d}|ddg}|d d g}d}d}t|||f|d \} } } } t || t || t j|| sJt j|| sJdS) Nrgy?ryrgwRN?g.-. ?y\9O?(?y\9O?(ֿrr )arraysqrtr(r)rr*r+) r,r-zeros_cpoles_ck_czeros_dpolls_dk_dr4zerospoleskr9s r:test_zerospolesgainzTestC2D.test_zerospolesgains((C;''((C"''!**,dRWWQZZ.?@AA**.0ABCC**DDFGG !7GS"9<).000ua "'5111!'5111|C#####|L"-------r<T)np_onlycZd}d}d}d}t||||\}}t||||\}} } d} d} dd| z | z|zzd| | z|zz z } | |zd| | z|zz z }|d| | z|zz z }|| |z|zz}t| |||\}}t||f| d| \}}}t||t||t| |||\}}}t|| | f| d| \}}}}t||t||t||dS) z7Test method='gbt' with alpha=0.25 for tf and zpk cases.grrrrrNrON)rr r)r)r,r-ABCDcnumcdenczcpckhrPAdBdCdDddnumddenc2dnumc2ddenr9dzdpdkc2dzc2dpc2dks r:test_gbt_with_sio_tf_and_zpkz$TestC2D.test_gbt_with_sio_tf_and_zpks    1aA&& dAq!Q'' B 1u9/A%%!eai!m*; < Ua%!)a-' ( !eai!m# $ R 2r2r** d!$qeLLLf%%%f%%%BB++ B #BB<5NNNdD"D!!!D!!!D!!!!!r<cfd}dgg}dgg}dgg}dgg}d}ddd}|d|dz } ||} t||||f|| | \}} } t||||f| d } d | d d| dd zz}|d d}t | |dd||\}}}d | d d| dd zz}t ||dd S)a Test that the solution to the discrete approximation of a continuous system actually approximates the solution to the continuous system. This is an indirect test of the correctness of the implementation of cont2discrete. c4d|zS)Ng@)sin)tr-s r:uz'TestC2D.test_discrete_approx..u"s66#'?? "r<g{GzrrLrg$@er)TUX0r\r rN)rrx0g-C6?rtol)r(linspacer r)rreshaperravel)r,r-rabcdrrr9u1youtxoutdsysu2t2td2yd2xd2ymids ` r:test_discrete_approxzTestC2D.test_discrete_approxs # # # # # JJy ! ! JJw   JJw   JJw    KK4 % % qTAaD[ QqTTaAq\Q"<<< 4Aq!Q<J777BssGbf$ % ssVdbjjQ&7&72"EEE S# d3B3i$qrr(*+ T555555r<cddgddggddgf}t|d\}}}|dksJt|ddgdt|ddgddggddS) Nrr{Gz?gUxegMbP?rrgT(|)r)r)r,tfrcrdr9s r: test_simo_tfzTestC2D.test_simo_tfBs1v1vA '2t}} S"Tzzzza.T::::q"g;'78tDDDDDDr<cd}ddgddggddgf}t||\}}}|dd|df}t||\}}} |dd|df} t| |\} } } || ksJ|| ksJt|tj|| fdt||dt|| ddS)NrrrgvIh%<=r)r)rnpvstack)r,tsrrcrdr9tf1num1den1dt1tf2num2den2dt2s r:test_multioutputzTestC2D.test_multioutputLs 2wA!Q (2r{{ S"!uQxAc2,,dC!uQxAc2,,dCSyyyySyyyy RYd|445AAAA T....T......r<N)__name__ __module__ __qualname__r;rCrErRrWrZr]rerrskip_xp_backendsrrrrr<r:rrs.===&===.===,0004===00002*.*.*.X . . ....(d###("("$#("T&6&6&6Pd###EE$#Ed###//$#///r<rTzlti currently not supported)rsreasonceZdZdZdZdS) TestC2dLtic"tjddgddgg}tjdgdgg}tjddgg}d}d}tjdd gd d gg}tjd gd gg}t||||} | |} t | j|t | j|t | j|t | jtj | jt| |} t | j|t | j|t | j|t | jtj | jdS)Ng333333ӿ皙?rLgffffffrrg?gГ>>>>r<c tddgddg}|d}tjddg}tjddg}t |j|d t |j|d dS) NrrHrg?g{Gzt?gi ߿gY쭛]g{Gz?)atol)r rrrgrrdrc)r,r-sysnum_resden_ress r: test_c2d_tfzTestC2dLti.test_c2d_tfs3*sCj))ooe$$(C!3455(C!3455 t4444t444444r<N)rrrrrrr<r:rres2???6 5 5 5 5 5r<rcveZdZeddggdddfeddggdddfedgd ddfgZeeeej d ed Z ej d ed Z ej d ed Z dS)TestC2dInvariantsr)rrTrr r)rrTrrr)rrrrzsys,sample_time,samples_numberc~tt|j|}|||z}t ||\}}t t ||dt|\}}t|| z|d dS)Nrr r nr) tuplemapr(aranger rr)lenrr r,r-r sample_timesamples_numbertime_ yout_cont yout_discs r:test_impulse_invariantz(TestC2dInvariants.test_impulse_invariantsC C(())yy((;6sd+++ 9CY G G G"%d))--- 9 ioo&7&7719K9K9M9MNNNNNr<cxtt|j|}|||z}t ||\}}t t ||dt|\}}t| |d dS)Nrrr rr) rrr(rrrr)rrrrs r:test_step_invariantz%TestC2dInvariants.test_step_invariantsC C(())yy((;6C4((( 9Sk%@@@CIINNN 9 ))9Q<+=+=+?+?@@@@@r<cXtt|j|}|||z}t |||\}}}t t ||d|\}}}t||dS)N)rrr>r )r) rrr(rr rr)rrrs r:test_linear_invariantz'TestC2dInvariants.test_linear_invariantsC C(())yy((;6sdd3339aCU C C CtLLL9a ))9??+<+<=====r<N)rrrr casesrr rpytestmark parametrizerrrrr<r:rrs! 1v{{{ # #T2. 1v~~~ & &R0 sLLL ! !3+ Ew)) [=uEEOOFE*)O [=uEEAAFEA [=uEE>>FE>>>r<r)r*numpyrscipy._lib._array_apirrrr scipy.signalrr)rrr r r r r rrrrrrrrrr<r:rs   ------88888888888888>>>>>>>>>>>>>> ;/3O/O/O/O/O/O/O/O/b $'DEEE'5'5'5'5'5'5'5FE'5T3%>%>%>%>%>%>%>%>%>%>r<