Path: csiph.com!usenet.pasdenom.info!gegeweb.org!de-l.enfer-du-nord.net!feeder2.enfer-du-nord.net!feeds.phibee-telecom.net!newsfeed.xs4all.nl!newsfeed3.news.xs4all.nl!xs4all!post.news.xs4all.nl!not-for-mail Return-Path: X-Original-To: python-list@python.org Delivered-To: python-list@mail.python.org X-Spam-Status: OK 0.077 X-Spam-Evidence: '*H*': 0.85; '*S*': 0.00; 'argument': 0.05; 'responding': 0.07; 'latter': 0.09; 'replied': 0.09; 'question.': 0.14; 'thread': 0.14; 'compute': 0.16; 'element:': 0.16; 'reedy': 0.16; 'index': 0.16; 'wrote:': 0.18; '>>>': 0.22; 'header:User- Agent:1': 0.23; 'header:In-Reply-To:1': 0.27; 'function': 0.29; 'array': 0.29; 'included': 0.31; 'easier': 0.31; 'context.': 0.31; 'call.': 0.33; 'but': 0.35; 'representing': 0.36; "i'll": 0.36; 'to:addr:python-list': 0.38; 'pm,': 0.38; 'track': 0.38; 'to:addr:python.org': 0.39; 'signal': 0.60; 'numbers': 0.61; "you're": 0.61; 'here:': 0.62; 'email addr:gmail.com': 0.63; 'real': 0.63; 'maximum': 0.63; 'skip:n 10': 0.64; 'more': 0.64; 'received:74.208': 0.68; 'said:': 0.68; 'to,': 0.72; "it'd": 0.84; 'received:74.208.4.194': 0.84 Date: Sun, 19 May 2013 21:11:57 -0400 From: Dave Angel User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:17.0) Gecko/20130510 Thunderbird/17.0.6 MIME-Version: 1.0 To: python-list@python.org Subject: Re: Harmonic distortion of a input signal References: In-Reply-To: Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit X-Provags-ID: V02:K0:lGLCkBVknVmYCqMKQrgnuVhw14bZM33UfW4Bz7z3MLX 7icdlq1VBsZ/ll9C4P+i7bTW/Mu2iZofrAqvYXfHH/dfOmavNg GcakIBeAIaUwnOfDidBdDh1nr2tgZLjSMsEqHIPJL5Xrk/+3O7 SxWn16l5fFULP7LSyAk9XWGfIk+pWX9g/aFG2bi1ecv6wU8Ok4 edvPlidMtLqphHYnbCDth3IlE8QCETkZZBFVQy5eP/XtU3JauB znFF5aNdhP/L/8Z7FxSjFHMxaXXTdA1cPqiOqIrmRMari0GfCB XXtubiYdOWwDrkjyKN2SNH8nYWBDRb55h23wZxLykgSEj5DPQ= = X-BeenThere: python-list@python.org X-Mailman-Version: 2.1.15 Precedence: list List-Id: General discussion list for the Python programming language List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Newsgroups: comp.lang.python Message-ID: Lines: 33 NNTP-Posting-Host: 2001:888:2000:d::a6 X-Trace: 1369012338 news.xs4all.nl 16006 [2001:888:2000:d::a6]:32999 X-Complaints-To: abuse@xs4all.nl Xref: csiph.com comp.lang.python:45576 On 05/19/2013 07:36 PM, killybeard91@gmail.com wrote: > One more question. Function np.argmax returns max of non-complex numbers ? > Because FFT array of my signal is complex. > It'd be easier to track the thread if you actually replied to the message you're responding to, and also if you included some context. But I'll paste the latter in here: Terry Reedy said: > Compute the magnitude spectrum: >>> spect = abs(np.fft.fft(x)[:25]) >>> spect > array([ 0. , 31.85194222, 0. , 10.67342282, > 0. , 6.47213595, 0. , 4.69726931, > 0. , 3.73254943, 0. , 3.13762901, > 0. , 2.7436023 , 0. , 2.47213595, > 0. , 2.28230601, 0. , 2.15105461, > 0. , 2.06487174, 0. , 2.01589594, > 0. ]) > Find the index of the maximum element: >>> np.argmax(spect) > 1 Notice that argmax's argument is the result of an abs() call. It's got real numbers representing the magnitude of the various complex numbers. -- DaveA