blob: a1251bdb93f516d86193562d0d1a92a972af5953 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
//DTFS coefficients of periodic square wave
clear;
close;
clc;
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','|ak|')
|
