summaryrefslogtreecommitdiff
path: root/22/CH4
diff options
context:
space:
mode:
Diffstat (limited to '22/CH4')
-rwxr-xr-x22/CH4/EX4.27/ch4ex27.sce21
1 files changed, 11 insertions, 10 deletions
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