Численное дифференцирование

@Zoops · Регистрация: 30.09.2012

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

Доброе время суток, ув. форумчане. У меня есть оцифрованный график перемещения моста при прогибе y(x) (наборы точек). Координату x я заменил временем и получил y(t) (разделил x на скорость ленты устройства, которое записывало колебания моста).

Далее мне нужен график зависимости скорости от времени, т.е. производная от перемещения. Т.к. функцию зависимости y(t) я не знаю - нужно использовать численное дифференцирование. Попробовал diff(x) и, также, нашел в инете, diff(x)/diff(y) - получаются 2 разных графика, какая запись будет правильной? Я склоняюсь к последнему варианту.

И в результате получается кол-во точек на 1 меньше, чем было. Так вот главный вопрос, как мне построить график y'(t), если кол-во точек y' на 1 меньше, чем кол-во точек t.

P.S. Возможна ли вообще такая замена x на t?

Прилагаю архив со значениями x и y

@faraonCar · 09.04.2018, 15:51

Численное дифференцирование

@Zoops · 09.04.2018, 19:18 **[ТС]**

faraonCar, сделал diff(y)./diff(x) и получился график с разрывами и NaN в некоторых значениях... Получается где-то будет деление на 0?

@faraonCar · 10.04.2018, 07:24

Все зависит от данных.

@Zoops · 10.04.2018, 19:36 **[ТС]**

faraonCar, хорошо, тогда вопрос, почему после дифференцирования точек становится на 1 меньше. И я не могу построить график x'(t), т.к. x' = допустим 99 точек, а t = 100 точек.

@faraonCar · 10.04.2018, 20:20

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

почему после дифференцирования точек становится на 1 меньше

Y = diff(X)
Y = [X(2)-X(1) X(3)-X(2) ... X(m)-X(m-1)]

@Zoops · 10.04.2018, 20:27 **[ТС]**

faraonCar, хорошо, тогда как мне построить график x'(t), если точек на 1 меньше у производной? Матлаб пишет, Vectors must be the same lengths.

@faraonCar · 10.04.2018, 20:32

Может это поможет
http://codetown.ru/matlab/chis... cirovanie/

@Zoops · 10.04.2018, 20:39 **[ТС]**

faraonCar, функция неизвестная, есть только набор данных значений x и t.

Добавлено через 1 минуту
faraonCar, polyfit-ом пробовал, очень плохо аппроксимирует.

@faraonCar · 10.04.2018, 20:40

Тогда смотрите ссылку снова, как делал Зосима
Численное дифференцирование

@Zoops · 10.04.2018, 20:42 **[ТС]**

faraonCar, все, я нашел в гугле, как строить: plot(X(1:end-1), dydx);

@faraonCar · 10.04.2018, 20:44

А если

Matlab M
1
diff(Y)./diff(t)

Добавлено через 52 секунды
А при чем тут plot?

@Zoops · 10.04.2018, 20:50 **[ТС]**

faraonCar,

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

faraonCar, хорошо, тогда как мне построить график x'(t), если точек на 1 меньше у производной? Матлаб пишет, Vectors must be the same lengths.

Это был второстепенный вопрос)

Добавлено через 1 минуту
faraonCar, премного благодарен за помощь)

@faraonCar · 10.04.2018, 21:04

Можно было сразу спросить.

Новые блоги и статьи Все статьи Все блоги /
Программная установка даты и запрет ее изменения Maks 02.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .	Вывод данных через динамический список в справочнике Maks 01.04.2026 Реализация из решения ниже выполнена на примере нетипового справочника "Спецтехника" разработанного в конфигурации КА2. Задача: вывести данные из ТЧ нетипового документа. . .	Функция заполнения текстового поля в реквизите формы документа Maks 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .
Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем	10 пpимет, которые всегда сбываются Maks 31.03.2026 1. Чтобы, наконец, пришла маршрутка, надо закурить. Если сигарета последняя, маршрутка придет еще до второй затяжки даже вопреки расписанию. 2. Нaдоели зима и снег? Не надо переезжать. Достаточно. . .	Перемещение выделенных строк ТЧ из одного документа в другой Maks 31.03.2026 Реализация из решения ниже выполнена на примере нетипового документа "ВыдачаОборудованияНаСпецтехнику" с единственной табличной частью "ОборудованиеИКомплектующие" разработанного в конфигурации КА2. . . .

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	09.04.2018, 15:51
	Численное дифференцирование 1

@Zoops 17 / 9 / 2 Регистрация: 30.09.2012 Сообщений: 95
	09.04.2018, 19:18 [ТС]
	faraonCar, сделал diff(y)./diff(x) и получился график с разрывами и NaN в некоторых значениях... Получается где-то будет деление на 0? 0

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	10.04.2018, 07:24
	Все зависит от данных. 0

@Zoops 17 / 9 / 2 Регистрация: 30.09.2012 Сообщений: 95
	10.04.2018, 19:36 [ТС]
	faraonCar, хорошо, тогда вопрос, почему после дифференцирования точек становится на 1 меньше. И я не могу построить график x'(t), т.к. x' = допустим 99 точек, а t = 100 точек. 0

@Zoops 17 / 9 / 2 Регистрация: 30.09.2012 Сообщений: 95
	10.04.2018, 20:27 [ТС]
	faraonCar, хорошо, тогда как мне построить график x'(t), если точек на 1 меньше у производной? Матлаб пишет, Vectors must be the same lengths. 0

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	10.04.2018, 20:32
	Может это поможет http://codetown.ru/matlab/chis... cirovanie/ 1

@Zoops 17 / 9 / 2 Регистрация: 30.09.2012 Сообщений: 95
	10.04.2018, 20:39 [ТС]
	faraonCar, функция неизвестная, есть только набор данных значений x и t. Добавлено через 1 минуту faraonCar, polyfit-ом пробовал, очень плохо аппроксимирует. 0

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	10.04.2018, 20:40
	Тогда смотрите ссылку снова, как делал Зосима Численное дифференцирование 1

@Zoops 17 / 9 / 2 Регистрация: 30.09.2012 Сообщений: 95
	10.04.2018, 20:42 [ТС]
	faraonCar, все, я нашел в гугле, как строить: plot(X(1:end-1), dydx); 0

@faraonCar 164 / 144 / 49 Регистрация: 29.10.2014 Сообщений: 738
	10.04.2018, 21:04
	Можно было сразу спросить. 0