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

не очень понятно что имееться ввиду

t=1:0.1:4*pi;
y=sin(t);
plot(t(100:116),y(100:116))

если я правильно понял вопрос то так

Имею ввиду вот что: у меня есть некая двухмерная фигура, описываемая функцией f=(x,y), которая зависит еще и от времени по обоим координатам (x и у).

Если я пишу: plot(x, y) - то матлаб строит мне эту фигуру во всем заданном интервале времени. А мне нужно только в конечный момент времени.

Не понял, вам 3-мерный график нужен?
Есть функции surf, mesh.

Нет, не 3-х мерный. У меня есть две разные синусоиды (синусоида тока и синусоида напряжения), которые изменяются со временем. Но мне надо вывести на графике ток относительно напряжения в конечный момент времени, а не за все периоды.

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

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

По ходу решения задачи опять возник вопрос: а можно ли из синусоидальной функции с наличием высших гармоник 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);
а дальше строим график обычным образом.

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

Возьмем 8-точеный пример:

1) Загоним туда косинусоиду (удобна тем, что ее спектр действителен):

2) Делаем fft

Ось симметрии проходит через частоту Найквиста.
Постоянная составляющая своего отражения не имеет.
"Комплексносопряженное" означает, что действительная часть числа сохраняется, а мнимое изменяет свой знак.

3) Обнуляем ненужные гармоники
Кроме u(2) надо оставить незануленной еще и ее отражение u(8).
Когда в массиве не 8 элементов, а больше, то надо рассчитать номер зеркального отражения.
Для первой гармоники он самый последний в ряду:


4) Делаем обратное преобразование Фурье

Готово! Поскольку у нас изначально была только 1-ая гармоника, то она и осталась.

Спасибо за пояснения! Только вот у меня этих элементов - несколько тысяч, как это зеркальное отражение рассчитать то?
И еще один вопрос: если мой сигнал модулирован низкой частотой, то 1-я гармоника будет именно этого модулирующего сигнала, может есть способы его убрать?

Я же сказала - отражение 1-ой гармоники будет в самом конце массива, отражение 2-ой гармоники - 2-ая сзади. В чем проблема? Или вы только с одного бока считать умеете? А может быть за вас программу такую написать, чтобы элементы с заднего конца считала? Тем более что для 1-ой гармоники я ее уже приводила:
u([1 3:length(u)-1])=0
Не заметили, как я единичку вычла? Это чтобы последний элемент не стирать! И эта формула работает при любой длине массива.



Про "модулирован" не знаю. Но думаю, что, как всегда, эту проблему можно решить на модельном примере. Возьмите массив, длиной в целую степень двойки (8, 16, 32, выше будет трудно разбираться), поскольку я не понимаю, как организовано отражение при fft на массивах иной длины. А потом промодулируйте одну частоту (достаточно высокую) другой (более низкой). Только обе частоты выберите с периодом, точно укладывающимся в отрезок. И посмотрите, что будет после fft. Можно для пущей ясности задать им разные амплитуды, чтобы в спектре это выглядело нагляднее. Ну и желательно косинусоиды, чтобы не возиться с мнимой частью.

Гугла, наконец, спросите - http://www.support17.com/component/content....html?task=view - пишет, что при амплитудной модуляции (какая у вас ЧМ или АМ вы так и не сказали) надо искать в спектре не частоту модуляции, а разностные частоты F+f и F-f по обе стороны от несущей F.

Да, так и сделаю. Спасибо за совет.

Есть совет еще лучше - сформулируйте точнее вашу проблему. Скажите прямым текстом, что вам нужна цифровая демодуляция, т.е. вычленение огибающей частоты на фоне несущей цифровыми методами. Тогда вам ответит вместо меня кто-то более сведущий в этом деле. А то ведь вы просите график построить...

Столкнулся опять с проблемой: матлаб почему-то не хочет строить график зависимости импеданса X от частоты f. Где Х - получается в ходе решения нескольких уравнений и в итоге в конечная формула выглядит так:
Х=1/Y+X2; частота f стоит в формулах Y и Х2.
Код графика такой:
plot(f, X, 'b');

Величину Х матлаб считает верно, а вот выводить ее зависимость от частоты - не хочет: на всех частотах приравнивает ее к 0. Кто-нибудь сталкивался с такой проблемой?

А X(f) комплексно-значное? MATLAB не умеет выводить комплексные числа на двумерном графике. Стройте отдельно модуль и аргумент от частоты.

Так вывести получилось, благодарю за подсказку!