updated the code

1 files changed, 31 insertions, 31 deletions

diff --git a/40/CH10/EX10.5b/Exa_10_5b.sce b/40/CH10/EX10.5b/Exa_10_5b.sce
index 7a3f04ee9..b98473816 100755
--- a/40/CH10/EX10.5b/Exa_10_5b.sce
+++ b/40/CH10/EX10.5b/Exa_10_5b.sce
@@ -1,32 +1,32 @@
-//Fir high pass filter design by frequency sampling

-z=%z;

-N=10;

-magHk=[0 0 0 1 1 1 1 1 0 0];

-k=[0:5 -4:-1:-1];

-fik=(-%pi*k'*(N-1)/N)+(0.5*%pi);

-for i=1:length(fik)

-    H1k(i)=magHk(i)*exp(%i*fik(i));

-end

-H1n=(dft(H1k,1));

-H2k=H1k;

-H2k(3)=0.5*%e^(-%i*1.3*%pi);

-H2k(9)=0.5*%e^(%i*1.3*%pi);

-H2n=(dft(H2k,1));

-H1Z=0;H2Z=0;

-for i=1:length(H1n)

-    H1Z=H1Z+H1n(i)*z^(-i); 

-end

-for i=1:length(H2n)

-    H2Z=H2Z+H2n(i)*z^(-i); 

-end

-F=0:0.01:1;

-F1=0:0.1:0.9;

-H1F=abs(horner(H1Z,exp(%i*2*%pi*F')));

-H2F=abs(horner(H2Z,exp(%i*2*%pi*F')));

-a=gca();

-plot2d(F1,magHk);

-plot2d(F,H2F);

-plot2d(F,H1F);

-xlabel('Digital Frequency F');

-ylabel('magnitude');

+//For high pass filter design by frequency sampling
+z=%z;
+N=10;
+magHk=[0 0 0 1 1 1 1 1 0 0];
+k=[0:5 -4:-1:-1];
+fik=(-%pi*k'*(N-1)/N)+(0.5*%pi);
+for i=1:length(fik)
+    H1k(i)=magHk(i)*exp(%i*fik(i));
+end
+H1n=(fft(H1k,1));
+H2k=H1k;
+H2k(3)=0.5*%e^(-%i*1.3*%pi);
+H2k(9)=0.5*%e^(%i*1.3*%pi);
+H2n=(fft(H2k,1));
+H1Z=0;H2Z=0;
+for i=1:length(H1n)
+    H1Z=H1Z+H1n(i)*z^(-i); 
+end
+for i=1:length(H2n)
+    H2Z=H2Z+H2n(i)*z^(-i); 
+end
+F=0:0.01:1;
+F1=0:0.1:0.9;
+H1F=abs(horner(H1Z,exp(%i*2*%pi*F')));
+H2F=abs(horner(H2Z,exp(%i*2*%pi*F')));
+a=gca();
+plot2d(F1,magHk);
+plot2d(F,H2F);
+plot2d(F,H1F);
+xlabel('Digital Frequency F');
+ylabel('magnitude');
 xtitle('Low pass filter using frequency sampling');
\ No newline at end of file