нужно написать прогу который по временным трендам сможет идентифицировать объект управления по методу наименьших квадрантов

@Tiger92 · Регистрация: 26.05.2013

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

нужно написать прогу который по временным трендам сможет идентифицировать объект управления по методу наименьших квадратов желательно с подписью каждой действии (или любой другой метод , но только с описанием )

вот мой вариант:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
function [s2] = mnk(num,den)
s1=tf(num,den);
T_end=60; % интервал измерений
dt=0.2; % шаг дискретизации
t=0:dt:T_end; % массив дискретного времени
N=length(t); % размер выборки
u=ones(N,1); % массив значений единичного входного воздействия
y=lsim(s1,u,t); %массив значений выходного воздействия
n = max([length(num)-1,length(den)-1]);
I=diag([1 1 1]);
i=1; % начальный шаг
P=1000*I; % начальное приближение
betta=[0;0;0];
bet(i,:)=betta; % массив оценок параметров
% формирование расширенной матрицы
for i=n+1:N-1
R=[y(i+n-2:-1:i-1);u(i+n-2:-1:i)]'; 
gamma=P*R'/(R*P*R'+1);
betta=betta+gamma*(y(i+1)-R*betta);
P=(I-gamma*R)*P;
bet(i,:)=betta;
end;
% расчет параметров непрерывной модели
T1=dt^2/(1-betta(1)-betta(2)); 
T2=(betta(2)*T1+T1+dt^2)/dt;
K=betta(3)*T1/dt^2;
s2=tf([K],[T1 T2 1]);%передаточная функция идентифицированного объекта
y2=lsim(s2,u,t);
% сравнение переходных характеристик объекта и модели
plot(t,y,t,y2,':'); 
grid;

, но преподу не понравилось, так как тут заранее указывается объект, а ему нужно чтоб это прога работала по неким временным трендам ( то есть по переходной характеристике -некая кривая )...

Суть таков : есть некая переходная характеристика —некая кривая , из этой кривой через равные шаги -шаг дискретизации , сняты точки ( например :y ( 0 1 2 3 4 5 6 7 8 7 6 5 4 5 5 5 5 5 5 5) и время t (1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 и тд)
рисунок ( переходная характеристика) примерно будет так выглядит:

нужно написать прогу который по временным трендам сможет идентифицировать объект управления по методу наименьших квадрантов

так вот задача заключается в том что зная значения этих точек определить объект, то есть передаточную функцию
есть математическое описание метода наименьших квадратов (страница 79)

diligenskaya_tutorial_2009.pdf

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

Зосима · 03.06.2013, 17:51

Ну смотри, сделать МНК не очень трудно, можно даже использовать встроенную ф-цию polyfit, также можно взять Фурье преобразование/проинтегрировать и получать различные виды выражений описания цепи, но вопрос в другом - как имея формулу идентифицировать звено?

@Tiger92 · 03.06.2013, 20:14 **[ТС]**

Зосима, под формулы вы имеете ввиду трансфер функцию? или я не понял... начальные данные это не формула а уже некий временной тренд ...

Зосима · 04.06.2013, 20:19

Странно, вчера вроде отвечал...

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

под формулы вы имеете ввиду трансфер функцию?

нет!

по заданным отсчетам (тренду) мы получаем только коэфф-ты апроксимирующего многочлена.

Новые блоги и статьи Все статьи Все блоги /
Воспроизведение звукового файла с помощью SDL3_mixer при касании экрана Android 8Observer8 26.01.2026 Содержание блога SDL3_mixer - это библиотека я для воспроизведения аудио. В отличие от инструкции по добавлению текста код по проигрыванию звука уже содержится в шаблоне примера. Нужно только. . .	Установка Android SDK, NDK, JDK, CMake и т.д. 8Observer8 25.01.2026 Содержание блога Перейдите по ссылке: https:/ / developer. android. com/ studio и в самом низу страницы кликните по архиву "commandlinetools-win-xxxxxx_latest. zip" Извлеките архив и вы увидите. . .	Вывод текста со шрифтом TTF на Android с помощью библиотеки SDL3_ttf 8Observer8 25.01.2026 Содержание блога Если у вас не установлены Android SDK, NDK, JDK, и т. д. то сделайте это по следующей инструкции: Установка Android SDK, NDK, JDK, CMake и т. д. Сборка примера Скачайте. . .	Использование SDL3-callbacks вместо функции main() на Android, Desktop и WebAssembly 8Observer8 24.01.2026 Содержание блога Если вы откроете примеры для начинающих на официальном репозитории SDL3 в папке: examples, то вы увидите, что все примеры используют следующие четыре обязательные функции, а. . .
моя боль iceja 24.01.2026 Выложила интерполяцию кубическими сплайнами www. iceja. net REST сервисы временно не работают, только через Web. Написала за 56 рабочих часов этот сайт с нуля. При помощи perplexity. ai PRO , при. . .	Модель сукцессии микоризы anaschu 24.01.2026 Решили писать научную статью с неким РОманом	http://iceja.net/ математические сервисы iceja 20.01.2026 Обновила свой сайт http:/ / iceja. net/ , приделала Fast Fourier Transform экстраполяцию сигналов. Однако предсказывает далеко не каждый сигнал (см ограничения http:/ / iceja. net/ fourier/ docs ). Также. . .	http://iceja.net/ сервер решения полиномов iceja 18.01.2026 Выкатила http:/ / iceja. net/ сервер решения полиномов (находит действительные корни полиномов методом Штурма). На сайте документация по API, но скажу прямо VPS слабенький и 200 000 полиномов. . .

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	03.06.2013, 17:51
	Ну смотри, сделать МНК не очень трудно, можно даже использовать встроенную ф-цию polyfit, также можно взять Фурье преобразование/проинтегрировать и получать различные виды выражений описания цепи, но вопрос в другом - как имея формулу идентифицировать звено? 0

@Tiger92 0 / 0 / 0 Регистрация: 26.05.2013 Сообщений: 10
	03.06.2013, 20:14 [ТС]
	Зосима, под формулы вы имеете ввиду трансфер функцию? или я не понял... начальные данные это не формула а уже некий временной тренд ... 0