1 files changed, 23 insertions, 0 deletions
diff --git a/3812/CH4/EX4.2/4_2.sce b/3812/CH4/EX4.2/4_2.sce
new file mode 100644
index 000000000..ed0048742
--- /dev/null
+++ b/3812/CH4/EX4.2/4_2.sce
@@ -0,0 +1,23 @@
+//Example 4_2
+//Finad DTFS coefficients of periodic square wave
+clc;
+clear;
+N=10;
+N1=2;
+Wo=2*%pi/N;
+xn=ones(1,length(N));
+n=-(2*N1+1):(2*N1+1);
+a(1)=(2*N1+1)/N;
+for k=1:2*N1
+a(k+1)=sin((2*%pi*k*(N1+0.5))/N)/sin(%pi*k/N);
+a(k+1)=a(k+1)/N;
+if(abs(a(k+1))<=0.1)
+a(k+1)=0;
+end
+end
+a=a';
+a_conj=conj(a);
+ak=[a_conj($:-1:1),a(2:$)];
+k=-2*N1:2*N1;
+plot2d3('gnn',k,abs(ak))
+xtitle('Magnitude spectrum','k','jakj')