diff options
Diffstat (limited to '3812/CH4/EX4.14.a/4_14_a.sce')
-rw-r--r-- | 3812/CH4/EX4.14.a/4_14_a.sce | 22 |
1 files changed, 22 insertions, 0 deletions
diff --git a/3812/CH4/EX4.14.a/4_14_a.sce b/3812/CH4/EX4.14.a/4_14_a.sce new file mode 100644 index 000000000..8f810210f --- /dev/null +++ b/3812/CH4/EX4.14.a/4_14_a.sce @@ -0,0 +1,22 @@ +//Example 4_14_a
+//determine the fourier series coeffient of x(n)
+clc;
+clear;
+N=6;
+n=0:0.01:N;
+Wo=2*%pi/N;
+xn=1*ones(1,length(n))+1*cos(Wo*n);
+for k=0:N-2
+C(k+1,:)=exp(-sqrt(-1)*Wo*n.*k);
+a(k+1)=xn*C(k+1,:)'/length(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:$)];
+Mag_ak=abs(ak);
+k=-(N-2):(N-2);
+plot2d3('gnn',k,Mag_ak,5)
+xtitle('abs(ak)','k','ak')
|