Записать выражения в python

@fyria30 · Регистрация: 11.07.2015

Студворк — интернет-сервис помощи студентам

Добрый день! Никак не удается правильно записать выражение в python, не могу найти ошибку, при том что в Matlab оно записывается верно. Само выражение (+прикрепил на листке, чтобы было понятнее):

g(t)=sin(2*pi*(f1*t+(f2-f1))*t^2/(2*T)

В MatlabРеализации:

Matlab M
1
2
3
4
for i  = 1 : RecordSec*Rate
    t(i) = 1/Rate*i;
    data(i)=sin(2*pi*(f1*t(i)+(f2-f1)*t(i)^2/(2*RecordSec)))
    end

в Python:

Python
1
2
3
    for i ,_ in enumerate(SweptSine):
         t = i * 1. / rate
         SweptSine[i] = math.sin(2*math.pi*t * ((f1*t)+(((f2-f1)*math.pow(t,2))/(2*RecordSec))))

Смысл всего этого, запись этого сигнала в wav формат. На выходе получаю разные сигналы. Спасибо!

@dondublon · 17.05.2018, 16:14

Насколько я помню, в матлабе поддерживаются векторные операции.
В python+numpy тоже. Разве что sin лучше взять из numpy, а не math. И в квадрат возводить без pow, а если с ним - тоже из np.
Так с чем проблема?

@fyria30 · 17.05.2018, 16:20 **[ТС]**

Сообщение от dondublon

Насколько я помню, в матлабе поддерживаются векторные операции.
В python+numpy тоже. Разве что sin лучше взять из numpy, а не math. И в квадрат возводить без pow, а если с ним - тоже из np.
Так с чем проблема?

То что вы описали, может повлиять на результат вычисления ?

Проблема в том, что получив таким образом массив данных и записав его в wav, я получаю разные результаты если сравнивать matlab и python. Запись данных в файл вроде бы не портачит, поэтому я грешу на массив данных. Вот ищу из-за чего получаются разные данные.

@dondublon · 17.05.2018, 16:29

fyria30, не должно, по идее.
Насчёт ошибки - надо аккуратно и пошагово всё проверять, тут знание питона не поможет. Пока я вижу, что у вас, вроде, лишнее t там затесалось, ну и скобки проверить желательно.

@Garry Galler · 17.05.2018, 16:31

Python
1
2
3
4
# ваша транспиляция выражения
sin(2*  pi * t * (( f1 * t) + (((f2 - f1) * t ** 2)) / (2 * RecordSec))))
# оригинальное выражение из matlab 
sin(2 * pi *       (f1 * t  +   (f2 - f1) * t ** 2   / (2 * RecordSec)))

Это все, чем могу помочь.
P.S. Две звездочки - это оператор возведения в степень.

@fyria30 · 17.05.2018, 16:37 **[ТС]**

Сообщение от dondublon

fyria30, не должно, по идее.
Насчёт ошибки - надо аккуратно и пошагово всё проверять, тут знание питона не поможет. Пока я вижу, что у вас, вроде, лишнее t там затесалось, ну и скобки проверить желательно.

Сообщение от Garry Galler

Это все, чем могу помочь.
P.S. Две звездочки - это оператор возведения в степень.

О боже, какая глупость получилась, глаз замылился уже, я все на скобки и порядок вычисления грешил, несколько часов просидел. Спасибо огромное!

Новые блоги и статьи Все статьи Все блоги /
Знаешь почему 90% людей редко бывают счастливыми? kumehtar 14.04.2026 Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.	Фиксация колонок в отчете СКД Maks 14.04.2026 Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .	Настройки VS Code Loafer 13.04.2026 { "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .	Оптимизация кода на разграничение прав доступа к элементам формы Maks 13.04.2026 Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .
Контроль заполнения и очистка дат в зависимости от значения перечислений Maks 12.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .	Архитектура слоя интернета для сервера-слоя. Hrethgir 11.04.2026 В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .	Подстановка значения реквизита справочника в табличную часть документа Maks 10.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача №1: при указании работ (справочник РаботыПоРемонтуСпецтехники),. . .	Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	17.05.2018, 16:14
	Насколько я помню, в матлабе поддерживаются векторные операции. В python+numpy тоже. Разве что sin лучше взять из numpy, а не math. И в квадрат возводить без pow, а если с ним - тоже из np. Так с чем проблема? 0

@dondublon 4653 / 2073 / 366 Регистрация: 17.03.2012 Сообщений: 10,183 Записей в блоге: 6
	17.05.2018, 16:29
	fyria30, не должно, по идее. Насчёт ошибки - надо аккуратно и пошагово всё проверять, тут знание питона не поможет. Пока я вижу, что у вас, вроде, лишнее t там затесалось, ну и скобки проверить желательно. 1

Записать выражения в python

Решение