1 files changed, 43 insertions, 0 deletions
diff --git a/135/CH10/EX10.1/EX1.sce b/135/CH10/EX10.1/EX1.sce
new file mode 100755
index 000000000..8c59bb156
--- /dev/null
+++ b/135/CH10/EX10.1/EX1.sce
@@ -0,0 +1,43 @@
+// Example 10.1: Asymptotic magnitude and phase response curves
+clc, clear
+w=[0:70];
+// Asymptotic magnitude response curve
+for i=1:length(w)
+    a(i)=32;
+    if w(i)<10 then
+        b(i)=0;
+        c(i)=0;
+    elseif w(i)<50
+        b(i)=14*(w(i)-10)/40;
+        c(i)=0;
+    else
+        b(i)=20*log10(w(i)/10);
+        c(i)=-20*log10(w(i)/50);
+    end
+end
+A=a+b+c;
+plot2d(w,A,rect=[0,0,70,50]);
+xtitle("Asymptotic magnitude response curve","ω(rad/sec)","20 log |A(jω)| in dB");
+// Asymptotic phase response curve
+scf(1);
+w=[1:600];
+for i=1:length(w)
+    if w(i)<1 then
+        theta1(i)=0;
+    elseif w(i)<5
+        theta1(i)=31.45*(w(i)-1)/4;
+        theta2(i)=0;
+    elseif w(i)<100
+        theta1(i)=45*log10(w(i)/10);
+        theta2(i)=-45*log10(w(i)/50);
+    elseif w(i)<500
+        theta1(i)=90;
+        theta2(i)=-58.55-31.45*(w(i)-100)/400;
+    else
+        theta1(i)=90;
+        theta2(i)=-90;
+    end
+end
+theta=theta1+theta2;
+plot(w,theta);
+xtitle("Asymptotic phase curve response","ω(rad/sec)","θ(ω)")
+\ No newline at end of file