+ /i(Rt^RIt^RIHtR.t/^^.b^^.b^^.b^^.b^^.b^^.b^.ROb^ ^.b^ ^.b^ ^ .b^ .ROb^ .ROb^.ROb^^.b^.R Ob^^.b^^ .b^.R O^^.^^.^^.^^.^.R O^^.^.R O^.R O^^.^^.^.RO^^.^ .RO/CtRRltR#)zTools for MLS generationN)_max_len_seq_inner max_len_seqcB\P!4P^8Xd\PM\PpVfV\ 9d`\P !\\ P444p\RVP4 RVP4 R24h\P !\ V,V4pM\P!\P !W444RRR 1,p\P!V^84'g0\P!W084'gVP^8d \R4h\P !V4p^V,^, pVfTpM\V4pV^8d \R4hVf)\P !V\P"RR7pM:\P !V\$RR7P'\P"4pVP(^8wgVPV8wd \R 4h\P*!V^8H4'd \R 4h\P,!V\P"RR7p\/W1WV4pWq3#) ay Maximum length sequence (MLS) generator. Parameters ---------- nbits : int Number of bits to use. Length of the resulting sequence will be ``(2**nbits) - 1``. Note that generating long sequences (e.g., greater than ``nbits == 16``) can take a long time. state : array_like, optional If array, must be of length ``nbits``, and will be cast to binary (bool) representation. If None, a seed of ones will be used, producing a repeatable representation. If ``state`` is all zeros, an error is raised as this is invalid. Default: None. length : int, optional Number of samples to compute. If None, the entire length ``(2**nbits) - 1`` is computed. taps : array_like, optional Polynomial taps to use (e.g., ``[7, 6, 1]`` for an 8-bit sequence). If None, taps will be automatically selected (for up to ``nbits == 32``). Returns ------- seq : array Resulting MLS sequence of 0's and 1's. state : array The final state of the shift register. Notes ----- The algorithm for MLS generation is generically described in: https://en.wikipedia.org/wiki/Maximum_length_sequence The default values for taps are specifically taken from the first option listed for each value of ``nbits`` in: https://web.archive.org/web/20181001062252/http://www.newwaveinstruments.com/resources/articles/m_sequence_linear_feedback_shift_register_lfsr.htm .. versionadded:: 0.15.0 Examples -------- MLS uses binary convention: >>> from scipy.signal import max_len_seq >>> max_len_seq(4)[0] array([1, 1, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 0, 0], dtype=int8) MLS has a white spectrum (except for DC): >>> import numpy as np >>> import matplotlib.pyplot as plt >>> from numpy.fft import fft, ifft, fftshift, fftfreq >>> seq = max_len_seq(6)[0]*2-1 # +1 and -1 >>> spec = fft(seq) >>> N = len(seq) >>> plt.plot(fftshift(fftfreq(N)), fftshift(np.abs(spec)), '.-') >>> plt.margins(0.1, 0.1) >>> plt.grid(True) >>> plt.show() Circular autocorrelation of MLS is an impulse: >>> acorrcirc = ifft(spec * np.conj(spec)).real >>> plt.figure() >>> plt.plot(np.arange(-N/2+1, N/2+1), fftshift(acorrcirc), '.-') >>> plt.margins(0.1, 0.1) >>> plt.grid(True) >>> plt.show() Linear autocorrelation of MLS is approximately an impulse: >>> acorr = np.correlate(seq, seq, 'full') >>> plt.figure() >>> plt.plot(np.arange(-N+1, N), acorr, '.-') >>> plt.margins(0.1, 0.1) >>> plt.grid(True) >>> plt.show() Nznbits must be between z and z if taps is NonezEtaps must be non-empty with values between zero and nbits (inclusive)z)length must be greater than or equal to 0c)dtypeorderz'state must be a 1-D array of size nbitszstate must not be all zeros)npintpitemsizeint32int64 _mls_tapsarraylistkeys ValueErrorminmaxuniqueanysizeintonesint8boolastypendimallemptyr)nbitsstatelengthtaps taps_dtype known_tapsn_maxseqs&&&& W/var/www/html/photoedit/myenv/lib/python3.14/site-packages/scipy/signal/_max_len_seq.pyrrsf WWY//14"((J |  !$y~~'7"89J5jnn6F5Gu * 011ACD Dxx %(*5yy$34TrT: 66$(  rvvdl33tyy1}:; ;xx~ XNE ~V A:HI I }RWWC8d#6==bggF zzQ%**-BCC vveqj677 ((6 4C tE3 ?E :)) ) r-) r0)r2r/))r6))r:r9)r*)r.)NNN)__doc__numpyr r__all__rrr)r(rFs42 / )Q )Q )Q )Q )Q )Q )Q  )  )! ) 1# )'); )8:K )  ) 2$ )(*K )9;bT )" )<bT2tR" "r<bT2| rB4bT2{ "r<  ) ur)