Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
40
Добавлен:
10.02.2015
Размер:
30.29 Кб
Скачать

p=[1 5 6];

r=roots(p);

r1=-2; r2=-3;

dp=polyder(p);

a1=polyval(dp,r1);

a2=polyval(dp,r2);

b0=4;

c1=b0./a1;

c2=b0./a2;

t=[0:0.01:5];

x=c1.*exp(r1.*t)+c2.*exp(r2.*t);

plot(t,x),grid on, xlabel('Time(sec)'), ylabel('x(t)')

c3=c1./(r1);

c4=c2./(r2);

a0=6;

c0=b0./a0;

x1=c0+c3.*exp(r1.*t)+c4.*exp(r2.*t);

plot(t,x1),grid on, xlabel('Time(sec)'), ylabel('x1(t)')

p=[1 2 17];

r=roots(p);

r1=-1+4.*i;

dp=polyder(p);

a1=polyval(dp, r1);

p1=[1 -1];

b=1;

dp1=polyder(p1);

b0=polyval(dp1, b);

c1=b0./a1;

t =[0:0.01:5];

x= exp(-1.*t).*sin(2.*t);

plot(t,x), grid on, xlabel('Time(sec)'), ylabel('x(t)')

c3=c1./r1;

a0=17;

c0=b0./a0;

x1=c0-0.4.*exp(-1.*t).*cos(2.*t)-0.2.*exp(-1.*t).*sin(2.*t);

plot(t,x1), grid on, xlabel('Time(sec)'), ylabel('x1(t)')

Соседние файлы в папке Новая папка
  • #
    10.02.201525.45 Кб39DZ.mdl
  • #
    10.02.2015441 б38DZ1.m
  • #
    10.02.2015491 б38DZ1_2.m
  • #
    10.02.201567.53 Кб39dz_1.docx
  • #
    10.02.201560.09 Кб37L1.mdl
  • #
    10.02.201530.29 Кб40L2.docx
  • #
    10.02.2015231 б37L2_1.m
  • #
    10.02.2015141 б38L2_2.m
  • #
    10.02.2015243 б38L2_3.m
  • #
    10.02.2015151 б38L2_4.m
  • #
    10.02.2015449.29 Кб39L4.docx