设计一50Hz工频干扰数字陷波器 matlab代码

以下是一个50Hz工频干扰数字陷波器的Matlab代码示例：

% 设计数字陷波器
fs = 1000; % 采样率
f0 = 50; % 工频干扰频率
w0 = f0 / (fs/2); % 归一化的工频干扰频率

% 模拟陷波器的设计
Q = 30; % 陷波器的质量因数
bw = w0 / Q; % 陷波器的带宽
[num,den] = iirnotch(w0,bw); % IIR数字陷波器设计

% 对一个信号进行陷波
t = 0:1/fs:1; % 时间向量
x = sin(2*pi*50*t) + sin(2*pi*150*t); % 信号包含50Hz和150Hz两个频率分量

y = filter(num,den,x); % 陷波后的信号

% 绘制信号和陷波后的信号
subplot(2,1,1);
plot(t,x);
xlabel('时间 (s)');
ylabel('幅值');
title('原始信号');

subplot(2,1,2);
plot(t,y);
xlabel('时间 (s)');
ylabel('幅值');
title('陷波后的信号');

这段代码使用了Matlab内置的`iirnotch`函数来设计数字陷波器，然后使用`filter`函数对一个包含50Hz和150Hz两个频率分量的信号进行陷波。最后使用`subplot`函数将原始信号和陷波后的信号绘制在同一个图像中进行比较。