HW6Problem2

HW6Problem2 - clc; clear all; close all; num=[0 0.6042 0...

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

View Full Document Right Arrow Icon
clc; clear all ; close all ; num=[0 0.6042 0 4.1572 0 6.1591]; den=[1 3.4630 7.1704 10.2787 8.6926 6.1591]; %% Part a lambda2=roots(den) % S domain roots lam=roots(num) % S domain zeros %% Part b figure() trf=tf(num,den); bode(trf); %% part c figure() %LT=2.*(z-1).*z./(3.*z-1); for r=0.8:0.05:1 p=1; for th=0:pi/36:2*pi z=r*exp(sqrt(-1)*th); LT(p)=(2.*(z-1).*z)./(3.*z-1); p=p+1; end hold all ; plot(real(LT),imag(LT)) end for th=0:pi/36:2*pi p=1; for r=0.8:0.05:1 z=r*exp(sqrt(-1)*th); LT(p)=2.*(z-1).*z./(3.*z-1); p=p+1;
Background image of page 1

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

View Full DocumentRight Arrow Icon
end hold all ; plot(real(LT),imag(LT)) end title( 'LambdaT plot' ) Xlabel( 'real axis' ) ylabel( 'Imaginary axis' ) %% Part d T1=0.7; % absolutely unstable T2=0.1; % relatively stable and absolutely stable T3=0.42; % relatively unstable and absolutely stable T=[T1,T2,T3] for j=1:1:3 figure() for r=0.8:0.1:1 p=1; for th=0:pi/36:2*pi z=r*exp(sqrt(-1)*th); LT(p)=2.*(z-1).*z./(3.*z-1); p=p+1; end hold all ; plot(real(LT),imag(LT)) end for th=0:pi/36:2*pi p=1; for r=0.8:0.05:1 z=r*exp(sqrt(-1)*th); LT(p)=2.*(z-1).*z./(3.*z-1); p=p+1; end hold all ;
Background image of page 2
Image of page 3
This is the end of the preview. Sign up to access the rest of the document.

This note was uploaded on 02/21/2011 for the course EE 456 taught by Professor Dr.yang during the Fall '09 term at SUNY Buffalo.

Page1 / 4

HW6Problem2 - clc; clear all; close all; num=[0 0.6042 0...

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

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