Project 5: Zero padding and decimation properties
5.1. Pad with Zeros
- Tạo tín hiệu rời rạc x=sin(w0n) với 16 điểm.
>> w=2*pi/sqrt(17);
>> N=16;
>> n=0:N-1;
>> X=sin(n*w);
>> stem(X)
Đồ thị:
Tính DFT của tín hiệu rời rạc trên sử dụng thuật toán FFT
Tính DFT 16 điểm
>> X1=fft(X)
X1 =
Columns 1 through 4
0.4817
0.4939 + 0.1488i 0.5441 + 0.3652i 0.7404 + 0.9381i
Columns 5 through 8
-2.5203 - 7.2727i 0.0287 - 0.7333i 0.1631 - 0.3198i 0.2028 - 0.1344i
Columns 9 through 12
0.2128
0.2028 + 0.1344i 0.1631 + 0.3198i 0.0287 + 0.7333i
Columns 13 through 16
-2.5203 + 7.2727i 0.7404 - 0.9381i 0.5441 - 0.3652i 0.4939 - 0.1488i
Phổ:
>> stem(X1)
FFT với 32 điểm.
Để tính FFT với 32 điểm, ta thêm vào tín hiệu mẫu 16 điểm 0 vào sau tín hiệu gốc
Để thêm ta sử dụng lệnh
%tạo tín hiệu có độ dài 32 (thêm 16 điểm 0 vào cuối tín hiệu mẫu)%
>> XX1=[X,zeros(1,16)];
>> stem(XX1)
Đồ thị:
Tính FFT của 32 điểm
>> XX11=fft(XX1)
XX11 =
Columns 1 through 4
0.4817
0.5850 - 0.0712i 0.4939 + 0.1488i 0.7612 - 0.2414i
Columns 5 through 8
0.5441 + 0.3652i 1.3590 - 0.5573i 0.7404 + 0.9381i 5.4938 - 2.2564i
Columns 9 through 12
-2.5203 - 7.2727i -3.9147 + 1.3820i 0.0287 - 0.7333i -1.7764 + 0.4706i
Columns 13 through 16
0.1631 - 0.3198i -1.3255 + 0.2157i 0.2028 - 0.1344i -1.1824 + 0.0647i
Columns 17 through 20
0.2128
-1.1824 - 0.0647i 0.2028 + 0.1344i -1.3255 - 0.2157i
Columns 21 through 24
0.1631 + 0.3198i -1.7764 - 0.4706i 0.0287 + 0.7333i -3.9147 - 1.3820i
Columns 25 through 28
-2.5203 + 7.2727i 5.4938 + 2.2564i 0.7404 - 0.9381i 1.3590 + 0.5573i
Columns 29 through 32
0.5441 - 0.3652i 0.7612 + 0.2414i 0.4939 - 0.1488i 0.5850 + 0.0712i
>> stem(XX11)
FFT 64 điểm
>> XX2=[X,zeros(1,48)];
>> stem(XX2)
>> XX21=fft(XX2)
>> stem(XX21)
FFT 256 điểm.
>> XX3=[X,zeros(1,240)];
>> stem(XX3)
>> XX31=fft(XX3);
>> stem(XX31)