Построение графика функции по точкам

@Singer5278 · Регистрация: 10.12.2016

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

Привет всем!
И снова проблема с графиком, но на этот раз серьезнее...
Объясняю суть задачи.
У меня есть уравнение и все дано, кроме t и U (написал в коде (1)). Из данного уравнения мне нужно получить одно решение из нескольких, которое ближе всего к предыдущей точке. Начальная (предыдущая) точка задается n=0, от нее мы начинаем искать ближнее решение.
Моя цель в данном случае осуществить кусок программы, который построил бы зависимость U от t по точкам (U;t), которые формируются так:
t увеличивается шагом с каждым разом (в цикле);
U каждый раз ищется так: с помощью функции fzero в окрестности точки 0 ищем первое ближайшее решение, а вот во второй и последующий разы необходимо найти U уже от получившегося значения n(i) в предыдущем случае. В примере (2) я написал, что мне нужно, но это должно быть циклом, а я не могу понять, как это сделать. То есть получившееся n нужно как-то сохранить и в следующий раз подставить в fzero(...,n). И вот из всех таких получившихся точек мне надо нарисовать график. Как-то так.

Matlab M

(1)
Imax = 0.0037;
Umax = 0.1;
Imin = 0.0005;
Umin = 0.2;
U0 = 2;
A = 2;
f = 5;
R = 1000;
t0 = 0;
t1 = 1;
Uvd = animatedline(handles.Uvd_t);
tdt = t0;
n = 0;
while tdt<t1
    f = @(Uvd1)(U0+A.*cos(2.*pi.*f.*tdt)-Uvd1)/R-(Imax.*(Uvd1/(Umax)).*exp(1-(Uvd1/(Umax)))+Imin.*exp(2.*(Uvd1-(Umin)))+(10^(-10)).*exp(40.*Uvd1));
    addpoints(Uvd,tdt,fzero(f,n));
    tdt=tdt+0.0001;
    ...???
end
 
(2)
addpoints(Uvd,tdt,fzero(f,n));
tdt=tdt+0.0001;
n1 = fzero(f,n);
addpoints(Uvd,tdt,fzero(f,n1));
tdt=tdt+0.0001;
n2 = fzero(f,n1);
addpoints(Uvd,tdt,fzero(f,n2));
tdt=tdt+0.0001;
n3 = fzero(f,n2);
%и так далее... но сделать нужно все в цикле

Подскажите, как можно такое провернуть. Ну никак не могу додуматься. Большое спасибо заранее за помощь!

Добавлено через 2 часа 20 минут
Хотел сделать такой цикл, но он не работает. Просчитывает n до конца и все, а дальше в цикл по кругу почему-то получившееся значение не идет(

Matlab M

Uvd = animatedline(handles.Uvd_t);
tdt = t0;
n = 0;
while tdt<t1
    f = @(Uvd1)(U0+A.*cos(2.*pi.*f.*tdt)-Uvd1)/R-(Imax.*(Uvd1/(Umax)).*exp(1-(Uvd1/(Umax)))+Imin.*exp(2.*(Uvd1-(Umin)))+(10^(-10)).*exp(40.*Uvd1));
    addpoints(Uvd,tdt,fzero(f,n));
    n = fzero(f,n);
    tdt=tdt+0.0001;
end

@АлександрКом · 25.12.2018, 11:08

Singer5278, может так:

Matlab M

Uvd = animatedline(handles.Uvd_t);
tdt = t0;
n = 0;
i = 1;
while tdt<t1
    f = @(Uvd1)(U0+A.*cos(2.*pi.*f.*tdt)-Uvd1)/R-(Imax.*(Uvd1/(Umax)).*exp(1-(Uvd1/(Umax)))+Imin.*exp(2.*(Uvd1-(Umin)))+(10^(-10)).*exp(40.*Uvd1));
    addpoints(Uvd,tdt,fzero(f,n(i)));
    n(i+1) = fzero(f,n(i));
    tdt=tdt+0.0001;
    i = i + 1;
end

	25.12.2018, 11:08

Построение графика функции по точкам

Решение