信号滤波matlab代码
时间: 2023-11-25 17:44:00 浏览: 36
这里提供一些常见的信号滤波matlab代码:
1. 均值滤波
```
function y = meanFilter(x, windowSize)
% x为输入信号,windowSize为滑动窗口大小
n = length(x);
y = zeros(n,1);
halfWindowSize = floor(windowSize/2);
for i = 1+halfWindowSize:n-halfWindowSize
y(i) = mean(x(i-halfWindowSize:i+halfWindowSize));
end
```
2. 中值滤波
```
function y = medianFilter(x, windowSize)
% x为输入信号,windowSize为滑动窗口大小
n = length(x);
y = zeros(n,1);
halfWindowSize = floor(windowSize/2);
for i = 1+halfWindowSize:n-halfWindowSize
y(i) = median(x(i-halfWindowSize:i+halfWindowSize));
end
```
3. 高斯滤波
```
function y = gaussianFilter(x, sigma)
% x为输入信号,sigma为高斯核的标准差
n = length(x);
y = zeros(n,1);
kerSize = 2*ceil(3*sigma)+1; % 根据3sigma原则计算核的大小
ker = fspecial('gaussian', kerSize, sigma); % 计算高斯核
y = conv(x, ker, 'same'); % 用卷积运算实现滤波
```
4. 带通滤波
```
function y = bandpassFilter(x, f1, f2, fs)
% x为输入信号,f1和f2为带通滤波器的频率范围,fs为采样频率
n = length(x);
y = zeros(n,1);
[b,a] = butter(2, [f1 f2]/(fs/2), 'bandpass'); % 用butter函数计算滤波器系数
y = filtfilt(b, a, x); % filtfilt函数可以保证滤波后的信号相位不变
```
以上是常用的几种滤波方法,可以根据具体需求选择合适的滤波方法。