Matlab: ノイズ除去のための簡単な一次遅れフィルタ(mファイル)

【必要なToolbox】

Control System Toolbox

ソースコード

前提条件: 

data: フィルタ処理前のデータがワークスペースに読み込み済み

%% 初期化

clc;                   % コマンドウィンドウの初期化

close all;          % 全てのグラフをclose

clear;                % 全変数クリア

 

%% 変数定義

s = tf('s');                   % ラプラス演算子
T = 10;                       % 時定数
F = 1/(1+T*s);            % 一次遅れフィルタ

 

%% フィルタ処理

data_Filter = lsim(F, data, 0:length(data)-1);          % フィルター処理後のデータ

 

%% グラフ化

figure;
plot(data);
hold on;
plot(data_Filter , 'm:')

legend('data', 'data-Filter')