HW2sol - 1(c clear all t1 = 0 0.01 1.99 t2 = 2 0.01 10 y1 = 0.5*exp-t1*cos(2*t1(pi/2 y2 = 0.5*exp-t2*cos(2*t2(pi/2 2.5*exp(t2-2*cos(2(t2-2(pi/2 t

Info iconThis preview shows pages 1–21. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 1(c) clear all t1 = 0 : 0.01 : 1.99; t2 = 2 : 0.01 : 10; y1 = 0.5*exp(-t1).*cos(2*t1+(pi/2)); y2 = 0.5*exp(-t2).*cos(2*t2+(pi/2)) + 2.5*exp(-(t2-2)).*cos(2*(t2-2)-(pi/2)); t = [t1, t2]; y = [y1, y2]; plot(t,y), xlabel('t'), ylabel('y(t)') 1(d) 1(e) Initial conditions set as derived in part (d). Max step size (in "simulation, configuration parameters") set to 50 ms to get smooth looking result. 1(f) If we double R, we get: Notice we have more damping. If we halve R, we get: Notice we have less damping. 2 (b) clear all k1 = 0 : 479; y1 = (250/0.003)*(1.005.^k1 - 1.002.^k1); k2 = 480 : 542; y2 = (250/0.003)*(1.005.^k2 - 1.002.^k2) - (5000*(1.002^480)/0.003)*(1.005.^(k2-480) 1.002.^(k2-480)); k = [ k1, k2 ]; y = [ y1, y2 ]; plot(k,y), xlabel('k'), ylabel('y(k)'), title('"plot" used instead of "stem" due to large time interval') 2 (c) 2 (d) 3 (b) clear all t1 = 0 : 0.01 : 0.99; t2 = 1 : 0.01 : 10; y1 = 1+3.1623*exp(-t1).*cos(t1+1.8925); y2 = 1+3.1623*exp(-t2).*cos(t2+1.8925)-(1+1.4142*exp(-(t2-1)).*cos((t2-1)+2.3562)); t = [t1 t2]; y = [y1 y2]; plot(t,y), xlabel('t'), ylabel('y(t)'), title('prob3b') 3 (c) Note: initial conditions are set up on integrator. 3 (d) 3 (e) Note: A,B,C,D and x(0) set up in state-space block. 3 (f) 4 (b) clear all k1 = 0; k2 = 1 : 10; y1 = -4 + 7/2 + 3/2; y2 = (7/2)*(0.5).^k2 + (3/2)*(-0.5).^k2; k = [k1 k2]; y = [y1 y2]; stem(k,y) xlabel('k'), ylabel('y(k)'), title('prob4b') 4 (c) Note: x(0) set up on delay element. 4 (d) 4 (e) Note: A,B,C,D and x(0) set up in state-space block. 4 (f) ...
View Full Document

This note was uploaded on 10/28/2008 for the course ECE 3085 taught by Professor Taylor during the Spring '08 term at Georgia Institute of Technology.

Page1 / 21

HW2sol - 1(c clear all t1 = 0 0.01 1.99 t2 = 2 0.01 10 y1 = 0.5*exp-t1*cos(2*t1(pi/2 y2 = 0.5*exp-t2*cos(2*t2(pi/2 2.5*exp(t2-2*cos(2(t2-2(pi/2 t

This preview shows document pages 1 - 21. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online