ps3._soln

ps3._soln - y = 1 exp-zeta*x cos_term(zeta/sqrt(1-zeta^2...

This preview shows pages 1–23. Sign up to view the full content.

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

View Full Document

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

View Full Document

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

View Full Document

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

View Full Document

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

View Full Document

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

View Full Document

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

View Full Document
EE 428 Problem Set 3 Solutions -15- Problem 13: (part 2) >> zeta_for_wntr_1p8 = 5.7879e-001

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

View Full Document
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);

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

View Full Document

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

View Full Document

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

View Full Document
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 >= 0.1) & (y <= 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

{[ snackBarMessage ]}

Page1 / 23

ps3._soln - y = 1 exp-zeta*x cos_term(zeta/sqrt(1-zeta^2...

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

View Full Document
Ask a homework question - tutors are online