Tải bản đầy đủ (.docx) (8 trang)

Project 5: Zero padding and decimation properties5.1. Pad with Zeros -BÁO CÁO BÀI TẠP DÀI XỬ LÝ THÔNG TIN

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (148.98 KB, 8 trang )

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)





×