График в MATLAB, как нарисовать? Опции

[Stefan1]

Доброго времени суток.
У меня есть зависимость x от у, которая изменяется с течением времени. Как мне построить график зависимости х от у в последний момент времени.

Сообщение отредактировал Stefan1 - Nov 20 2012, 09:56

[V_G]

Ну, так пишите plot(x[последние моменты],y[последние моменты]),
где последние моменты - диапазон индексов массивов x,y, соответствующих по-вашему последним моментам времени.

[Stefan1]

Ну, так пишите plot(x[последние моменты],y[последние моменты]),
где последние моменты - диапазон индексов массивов x,y, соответствующих по-вашему последним моментам времени.

Спасибо, заработало!

Сообщение отредактировал Stefan1 - Nov 21 2012, 17:04

[Stefan1]

По ходу решения задачи опять возник вопрос: а можно ли из синусоидальной функции с наличием высших гармоник u в зависимости от времени выделить первую гармонику и построить ее на графике также в зависимости от времени?
Я раскладываю в ряд Фурье функцию и затем выделяю из нее первую гармонику:

uf=fft(u); u1=uf(2)/256*2; uf=[];

Просто uf(2) - не строится. Как мне построить первую гармонику напряжения u в зависимости от времени (отбросив высшие гармоники)?

[Xenia]

По ходу решения задачи опять возник вопрос: а можно ли из синусоидальной функции с наличием высших гармоник u в зависимости от времени выделить первую гармонику и построить ее на графике также в зависимости от времени?
Я раскладываю в ряд Фурье функцию и затем выделяю из нее первую гармонику:
uf=fft(u); u1=uf(2)/256*2; uf=[];
Просто uf(2) - не строится. Как мне построить первую гармонику напряжения u в зависимости от времени (отбросив высшие гармоники)?

Наверное, вычислить обратное преобразование Фурье - ifft.
Например, так:
Стираем (зануляем) в массиве uf все частоты, кроме нужной вам
u1(1)=0
u1(3:length(u))=0
или в одну строку, если она вас не устрашит
u1([1 3:length(u)])=0
написала length(u1), т.к. не знаю точно, сколько у вас частот, вы же можете подставить конкретную константу,
потом
y=ifft(u1);
а дальше строим график обычным образом.

[Stefan1]

Наверное, вычислить обратное преобразование Фурье - ifft.
Например, так:
Стираем (зануляем) в массиве uf все частоты, кроме нужной вам
u1(1)=0
u1(3:length(u))=0
или в одну строку, если она вас не устрашит
u1([1 3:length(u)])=0
написала length(u1), т.к. не знаю точно, сколько у вас частот, вы же можете подставить конкретную константу,
потом
y=ifft(u1);
а дальше строим график обычным образом.

Все работает, только вот что-то не похоже на 1-ю гармонику (т.е. синусоиду). Может быть из-за того что я не то число ставлю вместо length(u1), т.к. у меня сигнал модулирован низкой частотой и первая гармоника здесь соответствует модулирующиму сигналу, а не несущей. Как бы теперь выделить первую гармонику несущей - наверно в ручную подбирать? т.к. частот модулирующего сигнала я не знаю.

Сообщение отредактировал Stefan1 - Nov 22 2012, 10:42