iPython notebook как альтернатива Matlab
#1

Бывает что нуджен инструмент по типу матлаба, там с фильтрами поковырятся, можетчисла там какие-то поперереставлять. Ну или просто нужен калкулятор работающий с комплексными числами :)
Помимо безплатного Octave есть еще достаточно популярный и все так же безплатный iPython.

Для установки всех этих дел на диск можно пойти одним из путей описаным здесь http://ipython.org/install.html
Я предпочитаю вот этот пакет, он все ставит сам https://store.enthought.com/downloads/ более чем достаточно для среднестатистического аудиофила :)
Далее, можно запускать как прилодение, а можно запускать как "записную книжку" и работать прямо из браузера. Для этого надо выполнить из желаемой директории ipython notebook --pylab=inline

В результате видим в браузере что-то типа этого:


Файлы вложений Эскизы(ов)
   

"Найкраще сало то ковбаса." (с)
The following 1 user says Thank You to БендеровецЪ for this post:
  • EDWARD (11-22-2014)
Ответ
#2

Как первый пример что можно сделать - затянуть звуковой файл и провести пару операций с ним. Для жтого кликаем на "новую записную книжку" и вводим такой код:


import numpy as np
import scipy as sp
import matplotlib.pyplot as plt
from scipy import signal
from scipy.io import wavfile

fs, sig = wavfile.read('EXIT.wav')
print fs
print sig.shape

plt.figure(figsize=(20,5))
plot( sig[0: 100000])
plt.xlabel('Sample')
plt.ylabel('Mag')

n_steps = (np.power(2, 16))
offset = 1000
freq = linspace(0, fs, n_steps)
window = np.zeros((n_steps, 2))
window[:, 0] = (np.blackman(n_steps))
window[:, 1] = (np.blackman(n_steps))
spectr = np.zeros((n_steps, 2))
spectr[:, :] = np.float64(np.multiply(window[:, :], sig[offset: offset + n_steps ,:] ))

spectrum = np.zeros((n_steps, 2))
print spectrum.shape

for i in range(2):
spectrum[:, i] = (20*np.log10(np.absolute(np.fft.fft(spectr[:, i])) + 1e-15) )
spectrum[:, i] -= np.amax(spectrum)

figure()
plt.figure(figsize=(15,7))
plt.plot(freq, spectrum[:, :]);
plt.xscale('log')
plt.xlim(1e2,fs/2)
plt.ylim(-150, 0)
plt.grid(which='both')
xlabel('Freq, Hz');
ylabel('dB');

print 'ch0 max sample = {} '.format(np.amax(sig[: ,0]))
print 'ch0 min sample = {} '.format(np.amin(sig[: ,0]))
print 'ch1 max sample = {} '.format(np.amax(sig[: ,1]))
print 'ch1 min sample = {} '.format(np.amin(sig[: ,1]))


Для удобства не обязательно вводить все это в одну ячейку, а можно разбить на несколько. Ячейки исполняются последовательно по одной, с нажатием шифт+энтер

На выходе это даст пару графиков, сам вид вейвформы (точнее ее части) и fft куска фонограммы


Файлы вложений Эскизы(ов)
       

"Найкраще сало то ковбаса." (с)
Ответ
#3

Если у вас есть кое какя измерительная апаратура подключаемая к компу и поддерживающая VESA, то с помощью пайтона тоже моджно ею управлять.
Общее описание - http://pyvisa.readthedocs.org/en/latest/
Для начало надо будет установить что-то с нейшенал инструмент VESA. Я ставил Agilent IO Libraries - http://www.keysight.com/en/pd-1985909-pn...330.977662&cc=AU&lc=eng
Когда все установлено и настрено это работает приблизительно вот так (пример с мостом E4980A).


Файлы вложений Эскизы(ов)
   

"Найкраще сало то ковбаса." (с)
The following 2 users say Thank You to БендеровецЪ for this post:
  • begemot (11-22-2014), Сергей Кор (11-08-2015)
Ответ
#4

Если непример хочется что-то пофильтровать - без проблем. В библиотеках уже есть синтез разных фильтров. Вот пример наложения Батеворта второго порядка на использованый ранее файл. Я правда int16 в дабл флоат не переводил.


Файлы вложений Эскизы(ов)
   

"Найкраще сало то ковбаса." (с)
Ответ
#5

123


Файлы вложений Эскизы(ов)
       

"Найкраще сало то ковбаса." (с)
The following 1 user says Thank You to БендеровецЪ for this post:
  • naroznyv (05-12-2015)
Ответ
#6

БендеровецЪ Написал:123
"иголоки" при переключении на выходе boost dc-dc LT3467 конвертера + их FFT. На выходе регулятора было 0805 4.7мкф и 22мкф +1ом послед для демфирования.

художник творит для тех, кто имеет похожую/близкую конструкцию мозга, позволяющую воспринимать созданное
Ответ
#7

Шота на форуме кроме как про бухло ваще ничего не пишут. Запощу пару картинок-иллюстраций к срачу на одном не безисвесном форуме :) На первой картинке можно с легкостью узнать синусоиду с частотой 22035Гц просемплированую с частотой 44100Гц :)


Файлы вложений Эскизы(ов)
                       

"Найкраще сало то ковбаса." (с)
The following 2 users say Thank You to БендеровецЪ for this post:
  • naroznyv (11-09-2015), EDWARD (11-10-2015)
Ответ
#8

БендеровецЪ Написал:Шота на форуме кроме как про бухло ваще ничего не пишут. Запощу пару картинок-иллюстраций к срачу на одном не безисвесном форуме :) На первой картинке можно с легкостью узнать синусоиду с частотой 22035Гц просемплированую с частотой 44100Гц :)
с первой понятно, далее ?

художник творит для тех, кто имеет похожую/близкую конструкцию мозга, позволяющую воспринимать созданное
Ответ
#9

Далее класика - зеро стафинг. Рассово верный метод :)

"Найкраще сало то ковбаса." (с)
Ответ
#10

Чуть более наглядная картика


Файлы вложений Эскизы(ов)
   

"Найкраще сало то ковбаса." (с)
The following 2 users say Thank You to БендеровецЪ for this post:
  • naroznyv (11-09-2015), EDWARD (11-10-2015)
Ответ
#11

Картинка спектра сигнала ~15kHz при fs 44.1kHz и "прозеростафленого" до 352.8kHz. Отрицательные частоты не убирал для наглядности.


Файлы вложений Эскизы(ов)
   

"Найкраще сало то ковбаса." (с)
The following 2 users say Thank You to БендеровецЪ for this post:
  • naroznyv (11-09-2015), EDWARD (11-10-2015)
Ответ
#12

Особо упоротые теперь говорят что нельзя востановить точно фазу близко к fs/2... Ну вот просемплированый, проаспсемпленый и отфильтрованый синус и косинус. Выглядят достаточно ортогонально :)


Файлы вложений Эскизы(ов)
   

"Найкраще сало то ковбаса." (с)
The following 1 user says Thank You to БендеровецЪ for this post:
  • naroznyv (11-10-2015)
Ответ
#13

БендеровецЪ Написал:Особо упоротые теперь говорят что нельзя востановить точно фазу близко к fs/2
интересно с чего такие предположения взяли

"The universal aptitude for ineptitude makes any human accomplishment an incredible miracle." John Stapp
The following 1 user says Thank You to nazar for this post:
  • naroznyv (11-10-2015)
Ответ
#14

Я думаю что лучше спросить прямо там в теме. Веселее будет :)

"Найкраще сало то ковбаса." (с)
Ответ
#15

У кого вы там спрашивать сбрались?
Я там вообще уже писать что-то боюсь, тем более что все равно не в коня корм....
Ответ
#16

Я согласен что их репертуар поднадоел. Не генерят ничего нового.

"Найкраще сало то ковбаса." (с)
Ответ
#17

Сергей, я пытаюсь "следить" за темой, но не в курсе проблемы. Проясните для невкурсячих, к чему эти все кривульки?
Ответ
#18

В двух словах: https://en.wikipedia.org/wiki/Noise_shaping и https://en.wikipedia.org/wiki/Delta-sigma_modulation
Я, на пальцах, лучше не обьясню.

"Найкраще сало то ковбаса." (с)
Ответ
#19

Покупаешь или собираешь на трёх ОУ (обычно) нотчницы и режешь нойс, чо там непонятного.
Ответ
#20

Какие оу, какой нотч! ты хоть понял о чем речь? Собсна можно не отвечать, и так видно :)

"Найкраще сало то ковбаса." (с)
Ответ


Возможно похожие темы ...
Тема / Автор Ответы Просмотры Последний пост

Перейти к форуму:


Пользователи, просматривающие эту тему: 2 Гость(ей)