Моделирование скорости гармонической системы генератора

Question

Я являюсь написать моделирования для получения истинной скорости гармонической системы генератора, как

где P = [p1 p2, p2 p3] можно найти с помощью Rung-Кутта Метод интеграции с P (0) = [1 0; 0 1]

Теперь, я хочу написать код matlab, чтобы получить истинное положение z системы и оценить x inwhich x = [x1 x2] '. Это мой код для поиска x и z. Однако результат неверен. Не могли бы вы помочь мне изменить его. Это мой код MATLAB

function z=getPos(x) 
    %% Function to get velocity 
    v=0+1*randn; %r=1; 
    z=[0 1]*x+v; 
return 

Это мой результат

Посмотрим помочь мне правильно реализовать для решения. Спасибо

Answer 1

Ваш код может возвращать отрицательное значение, так как randn является нормально распределенным случайным числом со стандартным отклонением 1. Теоретически это может вернуть любое число, и вероятность того, что будет ~ 15%, будет меньше -1.

Вы также верны, чем следует обновлять x на каждом шаге. Я бы принял x как вход в getPos и ​​xdot в качестве вывода. Затем добавьте еще один шаг для обновления позиции.