ps3._soln

ps3._soln - y = 1 - exp(-zeta*x) .* ( cos_term +...

Info iconThis preview shows pages 1–23. 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 DocumentRight 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 DocumentRight Arrow Icon
Background image of page 4
Background image of page 5

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

View Full DocumentRight Arrow Icon
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 DocumentRight 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 DocumentRight 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 DocumentRight 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 DocumentRight Arrow Icon
Background image of page 14
EE 428 Problem Set 3 Solutions -15- Problem 13: (part 2) >> zeta_for_wntr_1p8 = 5.7879e-001
Background image of page 15

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

View Full DocumentRight Arrow Icon
EE 428 Problem Set 3 Solutions -16- Problem 13: (part 2) % m-file for Problem Set 3, Problem 13 Part 2 % clear the workspace and close all figures clear close all % specify values of zeta to consider zeta_array = [0: .01 : 0.99]; % generate a vector for storing the value of wn*tr in wn_tr = zeros(1,length(zeta_array)); % for each value of zeta, determine the rise-time for n = 1 : length(zeta_array) % strip off value of zeta zeta = zeta_array(n); % calculate y(x) over the range 0 < x < wn*tp = wn*pi/wd x_max = pi / sqrt(1 -zeta^2); x = linspace(0,x_max,1000); % evaluate y(x) at given value of zeta cos_term = cos( sqrt(1 - zeta^2) * x); sin_term = sin( sqrt(1 - zeta^2) * x);
Background image of page 16
Background image of page 17

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

View Full DocumentRight Arrow Icon
Background image of page 18
Background image of page 19

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

View Full DocumentRight Arrow Icon
Background image of page 20
Background image of page 21

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

View Full DocumentRight Arrow Icon
Background image of page 22
Background image of page 23
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: y = 1 - exp(-zeta*x) .* ( cos_term + (zeta/sqrt(1-zeta^2)) * sin_term); % compute wn*tr I = find( (y &gt;= 0.1) &amp; (y &lt;= 0.9) ); wn_tr(n) = x( max(I) ) - x ( min(I) ); end % use linear interpolation to identify the value of zeta % for which wntr = 1.8 zeta_for_wntr_1p8 = interp1(wn_tr ,zeta_array, 1.8, 'linear' ) % plot wn*tr as function of zeta figure(1) handle = plot(zeta_array, wn_tr); set(handle, 'LineWidth' , 2, 'MarkerSize' , 10); set(gca, 'FontSize' , 14, 'FontName' , 'times new roman' ); title( 'Determination of Rise-Time as a Function of \zeta' , ... 'FontSize' , 14, 'FontName' , 'Times New Roman' ) ylabel( 'wn * tr' , 'FontSize' , 14, 'FontName' , 'Times New Roman' ) xlabel( '\zeta' , 'FontSize' , 14, 'FontName' , 'Times New Roman' )...
View Full Document

Page1 / 23

ps3._soln - y = 1 - exp(-zeta*x) .* ( cos_term +...

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

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