你好,
检查这个
清晰的所有
关闭所有
调频= 0.1;%的消息信号频率
Fs = 1;%的采样频率
Ts = 1 / f;%采样时间
n = 0: Ts: 999.9;%时间范围
N =大小(N, 2);% n空间的大小
sine_wave =罪(2 *π*调频* n);%产生正弦信号
随机= 2 *轮(兰德(1、长度(n))) 1;%产生噪音
%计算低通数字滤波器参数
%的假设低频率0赫兹;上频率0.12赫兹;
%通带波纹0.01(线性);阻带波动0.01(线性)
%采样频率1 hz
0.12订单= firpmord ([0]、[1 0], [0.01 - 0.01], 1);
coeffs1 = firls(秩序,[0 0.25 0.26 1],[1 1 0 0]);%计算coeficcients冷杉过滤器- 0.25和0.26弧度/秒的频率
sin_bwlimited =过滤器(coeffs1 1 sine_wave);%过滤正弦信号
random_bwlimited =过滤器(coeffs1 1随机);%过滤噪音
次要情节(5、1、1);%分割窗口
情节(n, sine_wave,“- b”);
持有在;
情节(n, sin_bwlimited,“——r”);
包含(“时间”);
ylabel (“信号”);
标题(正弦信号的带宽有限的);
次要情节(5、1、2);
情节(n,随机的,“- b”);
持有在;
情节(n, random_bwlimited,“——r”);
包含(“时间”);
ylabel (“信号”);
标题(随机噪声带宽有限的);
f3 = 0: Fs / N: Fs / 2 - f / N;%的空间频率域
Y = fft (sin_bwlimited) / N;
次要情节(5、1,3);
情节(f3, 2 * abs (Y (1: N / 2)));
包含(“频率”);
ylabel (“级”);
标题(过滤后的正弦信号的频率响应);
R = fft (random_bwlimited) / N;
次要情节(5、1、4);
情节(f3, 2 * abs (R (1: N / 2)));
包含(“频率”);
ylabel (“级”);
标题(过滤后的随机噪声的频率响应);
次要情节(5、1、5);
情节(coeffs1“o”);
包含(“样本”);
ylabel (“系数”);
标题(“数字滤波器系数”);
