modelica - Modelicaで平均値と標準偏差を計算する

Question

Modelica で連続信号の平均値 (mean) と標準偏差 (StdDev) を計算する最良の方法は何ですか。平均値と StdDev は一定期間 T について計算する必要があります。すなわち、tT から t まで。

Answer 1

以下は、問題に対する個別の解決策です。Modelica では、1 つの連続入力信号と 2 つの連続出力信号を持つブロックとしてコーディングされています。離散化は、Modelica 組み込み関数サンプルを使用して実行されます。

block MeanStdDevDiscr 
"Determines the mean value and standard deviation of a signal for a fixed time interval T."
  extends Modelica.Blocks.Interfaces.BlockIcon;
  import SI = Modelica.SIunits;

  parameter SI.Time T = 0.1 
    "Time interval used for calculating mean value and standard deviation";
  parameter Integer n = 10 "number of increments in T";

  Modelica.Blocks.Interfaces.RealInput u "signal input"
    annotation (Placement(transformation(extent={{-140,-20},{-100,20}})));
  Modelica.Blocks.Interfaces.RealOutput[2] y 
    "y[1] = average value; y[2] = standard deviation"
    annotation (Placement(transformation(extent={{100,-10},{120,10}})));

protected 
  parameter SI.Time dt = T/n "Precision of monitor";
  Real[n] uArray;

initial equation 
  uArray = ones(n)*u;

equation 
  when sample(0, dt) then
    uArray[1] = u;
    for j in 2:n loop
      uArray[j] = pre(uArray[j-1]);
    end for;
  end when;
  y[1] = sum(uArray)/n; // mean value
  y[2] = sqrt(sum((uArray .- y[1]).^2)/n); // standard deviation

  annotation (Diagram(graphics));
end MeanStdDevDiscr;