Matlab信号处理1：模拟去除信号噪声-Toy模板网

这篇具有很好参考价值的文章主要介绍了Matlab信号处理1：模拟去除信号噪声。希望对大家有所帮助。如果存在错误或未考虑完全的地方，请大家不吝赐教，您也可以点击"举报违法"按钮提交疑问。

一个简单的示例：

设计Matlab程序，去掉左图信号中大于20Hz的频率，使其变为右图。

Matlab信号处理1：模拟去除信号噪声,Matlab信号处理,matlab,信号处理,信号系统,数据信号处理

% 1.生成原信号并绘图
% 生成10s的信号，每隔0.0001s生成1次
t1 = 0:0.0001:10; 
% 原信号
y1 = 3*sin(2*pi*10*t1)+sin(2*pi*40*t1)+sin(2*pi*200*t1); 

% 原信号时域图 & 设置横纵轴范围 & 设置图名、横纵轴标签
subplot(321);
plot(t1,y1);
xlim([0 1]);
ylim([-6 6]);
title('原信号时域');
xlabel('时间（s）');
ylabel('幅值');

% 设置横纵轴分度值
set(gca,'XTick',0:0.25:1); 
set(gca,'YTick',-6:3:6);

% 2.模拟采样
% 采样时长1s，采样频率100Hz
t2 = 0:0.01:1;  
y2 = 3*sin(2*pi*10*t2)+sin(2*pi*40*t2)+sin(2*pi*200*t2); 

% 采样信号图 & 设置横纵轴范围 & 设置图名、横纵轴标签 & 设置横纵轴分度值
subplot(322);
plot(t2,y2);
xlim([0 1]);
ylim([-6 6]);
title('采样信号波形');
xlabel('时间（s）');
ylabel('幅值');
set(gca,'XTick',0:0.25:1); 
set(gca,'YTick',-6:3:6);

% 3.fft
% 采样后信号的数据长度
dataLength = length(y2);

% 对采样后的信号进行fft
y2FFT = fft(y2,dataLength);

% 求幅值
mag = abs(y2FFT);

% 幅值归一化
mag = mag*2/dataLength;
pha = angle(y2FFT)*180/pi;

% 无效相位置0
for i = 1:dataLength
    if (mag(1,i)<0.3)
        pha(1,i) = 0;
    end
end

% fft后的序列坐标
n = 0:dataLength-1;

% 采样频率
fs = 1/0.01;

% 序列频率
f = (0:dataLength-1)*fs/dataLength;

% fft频域图
% 绘制幅度谱，stem：绘制离散信号图
subplot(323);
stem(f(1:dataLength/2),mag(1:dataLength/2));
xlim([0 50]);
ylim([0 4]);
set(gca,"XTick",0:10:50);
xlabel("频率（Hz）");
title('幅度谱');

% 绘制相位谱
subplot(324);
stem(f(1:dataLength/2),pha(1:dataLength/2));
xlim([0 50]);
set(gca,'XTick',0:10:50);
xlabel("频率（Hz）");
ylabel("相位");
title('相位谱');

% 低通滤波器滤波  详见下面描述
% 载入设计的滤波器
load("my_filter.mat");  % 另存为的名字

% 滤波
% 滤波器系数存放于filtercoe数组中
filtercoe = my_filter;  % 设计完滤波器，导出的名字
y_Filtered = filter(filtercoe,1,y2);

% 滤波后的时域图
subplot(326);
plot(t2,y_Filtered);
xlim([0 1]);
ylim([-6 6]);
set(gca,'XTick',0:0.25:1); 
set(gca,'YTick',-6:3:6);
title('滤波后时域');
xlabel('时间（s）');
ylabel('幅值');

注：

1. 第85行低通滤波器的设计需要在命令行窗口输入filterDesigner以打开滤波器设计窗口，如下图：

Matlab信号处理1：模拟去除信号噪声,Matlab信号处理,matlab,信号处理,信号系统,数据信号处理