Построение графика производной для каждого графика из 3 режимов

@8-bit · Регистрация: 02.01.2013

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

Дорогие друзья, прошу помочь доделать три программы. В каждой программе выводится по одному графику и нужно чтобы каждая программа выводила не один, а два графика в разных окнах . Второй график должен иметь вид графика производной от первого графика.
Заранее очень благодарен.

Режим 1,2 и 3 это отдельные программы.

1) Свободные незатухающие колебания:

Matlab M

dt=0.0001;
t=0:dt:20;
y=zeros(1,length(t));
y(1)=1;
y1=zeros(1,length(t));
y1(1)=0;
for i=1:length(t)-1
  y1(i+1)=y1(i)+dt*f1(t(i),y1(i),y(i));
  y(i+1)=y(i)+dt*y1(i);
end
x=0:20;
 f=cos(x)
 figure(1);
 hold on
plot(t,y,'LineWidth',4);
plot(x,f,'Marker','o','LineStyle','none','Color','red','LineWidth',3);
 hold off

2) Затухающее колебательное движение:

Matlab M

dt=0.0001;
t=0:dt:20;
y=zeros(1,length(t));
y(1)=1;
y1=zeros(1,length(t));
y1(1)=0;
for i=1:length(t)-1
  y1(i+1)=y1(i)+dt*f1(t(i),y1(i),y(i));
  y(i+1)=y(i)+dt*y1(i);
end
x=0:20;
 f=exp(-x).*(cos(3*x)-(1/3)*sin(3*x));
 figure(1);
 hold on
plot(t,y,'LineWidth',4);
plot(x,f,'Marker','o','LineStyle','none','Color','red','LineWidth',3);
 hold off

3) Вынужденное колебание без сопротивления

Matlab M

dt=0.0001;
t=0:dt:20;
y=zeros(1,length(t));
y(1)=1;
y1=zeros(1,length(t));
y1(1)=0;
for i=1:length(t)-1
  y1(i+1)=y1(i)+dt*f1(t(i),y1(i),y(i));
  y(i+1)=y(i)+dt*y1(i);
end
x=0:20;
 f=cos(2*x)+(sin(2*x))/8-(x.*cos(2*x))/4;
 figure(1);
 hold on
plot(t,y,'LineWidth',4);
plot(x,f,'Marker','o','LineStyle','none','Color','red','LineWidth',3);
 hold off

Зосима · 02.01.2013, 13:30

Чтобы вывести два графика в разных окнах нужно в конце каждой программки написать:

Matlab M

figure
plot(t,y,'LineWidth',4);
figure
plot(x,f,':or','LineWidth',3);

Кстать, любопытно взглянуть на ф-цию f1(t,y1,y), можешь выложить?

@8-bit · 02.01.2013, 17:30 **[ТС]**

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

1) Свободные незатухающие колебания:

function f=f1(t,y1,y)
f=-y;

2) Затухающее колебательное движение:

function f=f1(t,y1,y)
f=-10*y-2*y1;

3) Вынужденное колебание без сопротивления

function f=f1(t,y1,y)
f=sin(2*t)-4*y;

@8-bit 0 / 0 / 0 Регистрация: 02.01.2013 Сообщений: 7
	02.01.2013, 17:30 [ТС]	3
	Спасибо, что откликнулись и за вашу помощь. Вот три функции к каждой программе. Как я и говорил уже ,нужно построить к каждому выдаваемому графику его производную в виде графика. На бумаге это не так уж и сложно, даже, могу выложить фотографию. Это к тому, что в крайнем случае 2-ой график можно просто подогнать по фотографии. Если моя просьба возможна к осуществлению, прошу помочь. 1) Свободные незатухающие колебания: function f=f1(t,y1,y) f=-y; 2) Затухающее колебательное движение: function f=f1(t,y1,y) f=-10y-2y1; 3) Вынужденное колебание без сопротивления function f=f1(t,y1,y) f=sin(2t)-4y; 0

Опции темы