diff options
Diffstat (limited to '3812/CH4/EX4.12/4_12_b.sce')
| -rw-r--r-- | 3812/CH4/EX4.12/4_12_b.sce | 21 |
1 files changed, 21 insertions, 0 deletions
diff --git a/3812/CH4/EX4.12/4_12_b.sce b/3812/CH4/EX4.12/4_12_b.sce new file mode 100644 index 000000000..84a6c32e5 --- /dev/null +++ b/3812/CH4/EX4.12/4_12_b.sce @@ -0,0 +1,21 @@ +//Example 4_12_b
+clc;
+clear;
+N=8;
+n=0:0.01:N;
+Wo=3*%pi/N;
+xn=1*ones(1,length(n))+1*sin(Wo*n+%pi/4);
+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')
|