blob: ed0048742b99c2c12513716d4757209283472de6 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
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')
|
