From f35ea80659b6a49d1bb2ce1d7d002583f3f40947 Mon Sep 17 00:00:00 2001 From: prashantsinalkar Date: Tue, 10 Oct 2017 12:38:01 +0530 Subject: updated the code --- 22/CH4/EX4.27/ch4ex27.sce | 21 +++++++++++---------- 1 file changed, 11 insertions(+), 10 deletions(-) (limited to '22/CH4/EX4.27/ch4ex27.sce') diff --git a/22/CH4/EX4.27/ch4ex27.sce b/22/CH4/EX4.27/ch4ex27.sce index 0389caa79..85993f0a8 100755 --- a/22/CH4/EX4.27/ch4ex27.sce +++ b/22/CH4/EX4.27/ch4ex27.sce @@ -1,10 +1,11 @@ -omega_0=2*%pi*60; theta = [60 80 87]*(%pi/180); -omega = (0:0.5:1000)'; mag = zeros(3,length(omega)); -s=poly(0,'s') -for m =1:length(theta) - H=syslin('c',((s^2+omega_0^2)/(s^2+2*omega_0*cos(theta(m))*s +omega_0^2))); - bode(H,10,100); -end -f=omega/((2*%pi))plot(f,mag(1,:),'k-',f mag(2,:),'k--',f,mag(3,:),'k-.'); -xlabel('f[hz]'); ylabel('|H(j2/pi f)|'); -legend('\theta=60^\circ','\theta = 80^\circ','\theta = 87^\circ',0) \ No newline at end of file +omega_0=2*%pi*60; theta = [60 80 87]*(%pi/180); +omega = (0:0.5:1000)'; mag = zeros(3,length(omega)); +s=poly(0,'s') +for m =1:length(theta) + H=syslin('c',((s^2+omega_0^2)/(s^2+2*omega_0*cos(theta(m))*s +omega_0^2))); + bode(H,10,100); +end +f=omega/((2*%pi)); +plot(f,mag(1,:)','k-',f,mag(2,:)','k--',f,mag(3,:)','k-.'); +xlabel('f[hz]'); ylabel('|H(j2/pi f)|'); +legend('\theta=60^\circ','\theta = 80^\circ','\theta = 87^\circ',1) \ No newline at end of file -- cgit