ex5_5 - pause clf perform time domain analysis use the...

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

View Full Document Right Arrow Icon
% Example 5.5 RC = input('Input RC '); % perform Frequency Domain analysis w = -40:.3:40; X = 2*sin(w/2)./w; H = (1/RC)./(j*w+1/RC); Y = X.*H; magY = abs(Y); angY = 180*angle(Y)/pi; % plot X(w) clf subplot(211),plot(w,abs(X)); ylabel('|X|') xlabel('Frequency (rad/sec)') title(['Example 5.5, RC = ' num2str(RC)]); subplot(212),plot(w,180*angle(X)/pi); ylabel('Angle(X), degrees') xlabel('Frequency (rad/sec)') pause % plot H(w) subplot(211),plot(w,abs(H)); ylabel('|H|') xlabel('Frequency (rad/sec)') title(['Example 5.5, RC = ' num2str(RC)]); subplot(212),plot(w,180*angle(H)/pi); ylabel('Angle(H), degrees') xlabel('Frequency (rad/sec)') pause % plot Y(w) subplot(211),plot(w,magY); ylabel('|Y|') xlabel('Frequency (rad/sec)') title(['Example 5.5, RC = ' num2str(RC)]); subplot(212),plot(w,angY); ylabel('Angle(Y), degrees') xlabel('Frequency (rad/sec)')
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
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: pause clf % % perform time domain analysis % use the numerical convolution described in Example 3.10 % for accuracy, choose T <= RC/10 if RC >= .1, T = .01; kh = 0:400; % defines indices for h, % corresponds to t=0 to t=4 h = exp(-1/RC*kh*T)/RC; kx = -60:100; % defines indices for x, % corresponds to t=.6 to t=1 x = [zeros(1,10) ones(1,101) zeros(1,50)]; elseif RC >=0.01, T = 0.001; kh = 0:400; h = exp(-1/RC*kh*T)/RC; kx = -600:1000; x = [zeros(1:100) ones(1,1001) zeros(1,500)]; else, error('RC is too small for accurate results using the numerical convolution method') end y = conv(x,h*T); ky = kx(1)+kh(1):kx(length(kx))+kh(length(kh)); clf subplot(211), plot(ky*T,y) title(['Pulse Response, RC = ' num2str(RC)]) ylabel('y(t)') xlabel('Time (sec)') subplot(111)...
View Full Document

This note was uploaded on 03/24/2010 for the course CENG 4331 taught by Professor Maryrandolph-gips during the Fall '09 term at UH Clear Lake.

Page1 / 2

ex5_5 - pause clf perform time domain analysis use the...

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

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