ps4_soln - y_est_R51 = step(b0[1 a1_R51 a0 tsim y_est_R1000...

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

View Full Document Right Arrow Icon
Background image of page 1

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

View Full Document Right Arrow Icon
Background image of page 2
Background image of page 3

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

View Full Document Right Arrow Icon
Background image of page 4
EE 350 Problem Set 4 Solutions -5- Problem 17: (part 4)
Background image of page 5

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

View Full Document Right Arrow Icon
EE 350 Problem Set 4 Solutions -6- Problem 17: m-file % EE 350 Fall 2007 % Problem Set 4, Problem 17 % measured resistance values R51 = 50.2; R1000 = 983; R2200 = 2180; RL = 41.4; % load unit-step response of RC circuit load RC_resp % find rise-time tr (in ms) yf = y(end); ind = find( (y >= 0.1*yf) & (y<=0.9*yf) ); tr_RC = 1e3 * ( max( t(ind) ) - min( t(ind) ) ); % estimate the value of C C = tr / R2200 /log(9) % load unit-step response of RL circuit load RL_resp % find rise-time tr (in ms) yf = y(end); ind = find( (y >= 0.1*yf) & (y<=0.9*yf) ); tr_LC = 1e3 * ( max( t(ind) ) - min( t(ind) ) ); % estimate the value of L L = tr * (RL + R51) /log(9) % estimate RLC unit-step response: b0 = 1 / L / C; a0 = b0; a1_R51 = (R51 + RL)/L; a1_R1000 = (R1000 + RL)/L; a1_R2200 = (R2200 + RL)/L; tsim = linspace(0,0.14, length(t));
Background image of page 6
Background image of page 7

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

View Full Document Right Arrow Icon
Background image of page 8
Background image of page 9

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

View Full Document Right Arrow Icon
Background image of page 10
Background image of page 11

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

View Full Document Right Arrow Icon
Background image of page 12
Background image of page 13

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

View Full Document Right Arrow Icon
Background image of page 14
Background image of page 15

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

View Full Document Right Arrow Icon
Background image of page 16
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: y_est_R51 = step(b0, [1, a1_R51, a0], tsim); y_est_R1000 = step(b0, [1, a1_R1000, a0], tsim); y_est_R2200 = step(b0, [1, a1_R2200, a0], tsim); % load the measured responses load RLC_R51_resp ; y_meas_R51 = y; load RLC_R1000_resp ; y_meas_R1000 = y; load RLC_R2200_resp ; y_meas_R2200 = y; % compare the simulated and measured responses t = t*1e3; tsim = tsim*1e3; % show time in ms figure(1) subplot(311) plot(t, y_meas_R51, '-r' , tsim, y_est_R51, '--g' ); title( 'RLC Zero-State Unit-Step Response' ) ylabel( 'y(t) [V]' ) legend( 'measured' , 'R = 51 \Omega' ) subplot(312) plot(t, y_meas_R1000, '-r' , tsim, y_est_R1000, '--g' ); ylabel( 'y(t) [V]' ) legend( 'measured' , 'R = 1000 \Omega' ) subplot(313) plot(t, y_meas_R2200, '-r' , tsim, y_est_R2200, '--g' ); ylabel( 'y(t) [V]' ) legend( 'measured' , 'R = 2200 \Omega' ) xlabel( 'time [ms]' )...
View Full Document

{[ snackBarMessage ]}

Page1 / 16

ps4_soln - y_est_R51 = step(b0[1 a1_R51 a0 tsim y_est_R1000...

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

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