低通数字滤波器实现
本节的MATLAB源代码覆盖低通数字滤波器matlab代码。它涵盖了滤波器设计使用firrcos函数和matlab FDA的工具。
使用FIRRCOS MATLAB函数
FFT_size = 1024;
t = 0:0.001:(FFT_size-1) / 1000;
x =罪(2 *π* 50 * t) +罪(2 *π* 200 * t);
多项式系数= firrcos (50100, 1000);
y = conv (x,多项式系数);
日元= fft (y, 1024);
图;情节(abs (fft (x, 1024)));
人物,情节(abs (y1));
我们以输入信号为正弦波形的两座山峰50赫兹和200赫兹。因此,如果我们通过这个信号通过低通滤波器的频率高于截止频率应该被消除。这里FFT的长度是1024。在这个实现中,我们使用MATLAB预定义函数生成系数滤波器设计。
多项式系数= firrcos (n, fc, bw, fs);
在那里,
' n ' -过滤器的顺序。n是一个正整数。
“俱乐部”——指定滤波器的截止频率。俱乐部必须大于0小于fs / 2。
“bw”——指定的过渡带宽。bw必须大于0。
[fc-bw / 2, fc + bw / 2)必须下降范围(0,fs / 2)。
“fs”——指定采样频率。fs必须大于0。
我们也使用“conv”预定义函数卷积系数得到的输入信号时域输出。
y = conv (x,多项式系数);
在那里,
“x”——输入信号。
多项式系数,系数。
卷积后得到的输出(y),那么情节人物FFT输出(y)与FFT输入(x)。FFT是将时域信号转换为频域信号。
图1:输入信号在频域(x) 50赫兹和200赫兹。
右边的山峰在上面的图是镜像的原始输入信号。
下面是输出信号在频域(y)。因此,从这两个数据我们可以得出结论,峰值在200赫兹正在消除,因为它是大于截止频率(100赫兹)。
注意:-
假设如果我们给一个输入信号110赫兹和截止频率为100赫兹,然后输入信号输出将不会完全消除,因为它需要一段时间才能完全腐烂。这个衰减时间称为过渡时期。
图2:输出信号在频域(y)。
在MATLAB实现低通数字滤波器使用FDA工具
关闭所有;
clc;
FFT_size = 1024;
t = 0:0.001:(FFT_size-1) / 1000;
x1 =罪(2 *π* 15 * t) +罪(2 *π* 100 * t);
x1 = x1 /;
% % fs = 1000 Hz,成就= 50 Hz, fstop = 100 Hz
b = (-0.00075204 -0.0034964 -0.0025769 -0.0029542 -0.0019023 7.5785 0.0029216 e - 005
0.0060432 0.0086036 0.0096315 0.0082863 0.0041324 -0.0026099 -0.010933 -0.019081 -0.024793 -0.025704 -0.019846 -0.0061564 0.015154 0.042456 0.072834 0.1025 0.12742 0.14402 0.14984 0.14402 0.12742 0.1025 0.072834 0.042456 0.015154 -0.0061564 -0.019846 -0.025704 -0.024793 -0.019081 -0.010933 -0.0026099 0.0041324 0.0082863 0.0096315 0.0086036 0.0060432 0.0029216 7.5785 e - 005 -0.0019023 -0.0029542 -0.0025769 -0.0034964 -0.00075204);
b_dec = b * (32768);
b_deca = b_dec (:);
x_dec = x1 * (32768);
x_decaa = x_dec (:);
x = [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 x1);k = 55:1077
y (k) = (x (k) * b (1)) +…
b (x (k - 1) * (2)) +…
(x (k-2) * b (3)) +…
b (x (k3) * (4) +…
b (x (k-4) * (5)) +…
(x(6年级)* b (6)) +…
b (x (k-6) * (7)) +…
b (x (k-7) * (8) +…
b (x(学习)* (9))+…
b (x (k - 9) * (10)) +…
b (x (k-10) * (11)) +…
b (x (k-11) * (12)) +…
b (x (k - 12) * (13)) +…
b (x (k-13) * (14)) +…
b (x (k-14) * (15)) +…
b (x (k-15) * (16)) +…
b (x (k-16) * (17)) +…
b (x (k-17) * (18)) +…
b (x (k-18) * (19)) +…
b (x (k-19) * (20)) +…
b (x (k-20) * (21)) +…
b (x (k-21) * (22)) +…
b (x (k-22) * (23)) +…
b (x (k-23) * (24)) +…
b (x (k-24) * (25) +…
(x (k-25) * b (26)) +…
b (x (k-26) * (27)) +…
b (x (k-27) * (28) +…
b (x (k-28) * (29)) +…
b (x (k-29) * (30) +…
b (x (k 30) * (31)) +…
b (x (k-31) * (32) +…
b (x (k-32) * (33)) +…
b (x (k-33) * (34) +…
b (x (k-34) * (35)) +…
b (x (k-35) * (36)) +…
b (x (k-36) * (37)) +…
b (x (k-37) * (38)) +…
b (x (k-38) * (39)) +…
b (x (k-39) * (40)) +…
b (x (k-40) * (41)) +…
b (x (k-41) * (42)) +…
b (x (k-42) * (43)) +…
b (x (k-43) * (44)) +…
b (x (k-44) * (45)) +…
b (x (k-45) * (46)) +…
b (x (k-46) * (47)) +…
b (x (k-47) * (48)) +…
b (x (k-48) * (49)) +…
b (x (k-49) * (50));
y_dec = y * 32768;
y_deca = y_dec (:);
结束
z = y(51:结束);
日元= fft (z, 1024);
图;情节(abs (fft (x1, 1024)));
人物,情节(abs (y1));
y1_dec = 32768日元*;
y1_deca = y1_dec (:);
图3:输入信号在频域15赫兹和100赫兹。
图4:输出信号在频域(y) 15赫兹。从上面两个图我们可以得出结论,在100赫兹频率的输入信号是要切断输出。
在这个实现中,我们使用的是食品及药物管理局(过滤器设计分析)工具生成滤波器系数。这个过滤器的规格是:
Fs = 1000 Hz。
成就= 50赫兹。
Fstop = 100 Hz。
公布:FDA工具来生成低通数字滤波器的系数。
图6:低通数字滤波器的结构。
在这个实现FFT的长度是1024。这里,罪恶的输入信号由两个峰波15赫兹和100赫兹。
x1 =罪(2 *π* 15 * t) +罪(2 *π* 100 * t);
我们扩展输入信号,因为当输入(x)范围(我不在Q-15格式。e, 1到0.9999999)。因此,我们即输入x除以2.,x1 /。
现在,系数分配给“b”。现在我们将做补零的x1因为在前面输入信号滤波器方程也乘以系数。过滤方程是专为“50”低通数字滤波器的顺序。“k”值从55,因为原始输入信号从55位置在MATLAB和时光结束1077位置。
实现过滤方程如下:
Y (n) = x (n) * b (1) + x (n - 1) * b (2) + ..... + x (n n) * b (n)
如果1号输入然后乘以系数和剩余的1日系数乘以以前输入的“0”值。
如果第二输入然后第二输入系数乘以1日和第一输入然后乘以第二系数和剩余的系数的增加与先前输入的“0”值。
它会在相似的方式,直到输入过来的系数得到定期重复。然后我们输入和输出值的FFT的阴谋。
MATLAB的有用链接代码
参考以下以及链接提到左侧面板上有用的MATLAB代码。
OFDM序言代时间估计corr频率下估计corr信道估计11 WLAN通道PN序列生成OFDMA Tx的处方AES DES载波聚合CCDF冷杉过滤器IIR滤波器低通冷杉维特比译码器CRC8 CRC32
射频和无线教程
无线局域网802.11交流802.11广告wimax无线个域网z - waveGSMLTEUMTS蓝牙超宽频物联网卫星天线雷达