コード置換ライブラリ非対応関数をコード置換する方法

Question

Shoumei 2018-2-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/384122-

编辑： Atsushi Matsumoto 2018-3-23

Cortex-M4向けにコード置換ライブラリを使用してターゲットプロセッサ向けに最適化されたCコードを生成したいのですが、MATLABファイル内で使用しているfiltfilt関数が対応していないようです。filtfilt関数の処理は普通の畳み込み演算＋アルファの処理を行っているので、コード置換で最適化できそうなものですが、何か良い方法はありませんか？

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Atsushi Matsumoto 2018-3-23

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/384122-#answer_311483

编辑：Atsushi Matsumoto 2018-3-23

在 MATLAB Online 中打开

コード置換に対応していない関数は、内部処理をコード置換に対応した関数で書き換えることで対応させることができる場合があります。 filtfiltやsgolayfiltなどのIIR/FIRフィルタベースの処理関数はソースコードをカスタマイズすることでコード置換に対応させることができます。

&nbsp

filtfiltのソースコードをIIRフィルタのSystem Object: dsp.BiquadFilterを使って書き換えた例を添付します。

この例では2次のIIRフィルタを使ってフィルタリングする例となっていますので、異なるフィルタを使用している場合は、一部の修正が必要になると思います。

function y = filtfiltSysO(x, sos, g)
%FILTFILT Zero-phase forward and reverse digital IIR filtering.
% filter, reverse data, filter again, and reverse data again
persistent Hd1 Hd2
if isempty(Hd1)
  Hd1 = dsp.BiquadFilter('Structure','Direct form II transposed', ...
      'SOSMatrix',sos,'ScaleValues',g);
  Hd2 = dsp.BiquadFilter('Structure','Direct form II transposed', ...
      'SOSMatrix',sos,'ScaleValues',g);
end
nfact = single(6);
y = [2*x(1)-x(nfact+1:-1:2); x; 2*x(end)-x(end-1:-1:end-nfact)]; 
% y = filter(b(:,ii),a(:,ii),y,zi(:,ii)*y(1));
y = step(Hd1, y);
y = y(end:-1:1);
% y = filter(b(:,ii),a(:,ii),y,zi(:,ii)*y(1));
y = step(Hd2, y);
% retain reversed central section of y
y = y(end-nfact:-1:nfact+1);