माडेलिका (SystemModeler) में सफेद शोर

Question

मैं एक अनुकरण करने के लिए माप शोर जोड़ने के लिए कोशिश कर रहा हूँ उत्पन्न करें। उदाहरण के लिए सिमुलिंक करना संभव है लेकिन मॉडलिका और सिस्टममोडेलर में अधिक कठिन लगता है।

यह कैसे करें इस पर कोई विचार है?

Answer 1

एक वैकल्पिक Modelica.Blocks.Noise उपयोग करने के लिए बाहरी कोड अपने आप को (माडेलिका स्टैंडर्ड लाइब्रेरी 3.2.2 3 अप्रैल, वर्ष 2016 में जारी की में जोड़ा लेखन से बचने के लिए है, यानी जब मूल प्रश्न पूछा गया तो इससे मदद नहीं मिली होगी)। Modelica.Blocks.Noise की

एक लाभ यह है कि नमूना, कई बीज, आदि के साथ मुश्किल सवाल हल कर रहे हैं।

Answer 2

आप बाहरी सी-कोड के माध्यम से Wolfram SystemModeler में सफेद शोर जोड़ सकते हैं।

माडेलिका कोड (मैं कोड से आरेख एनोटेशन हटा दिया है, ताकि इसे पढ़ने के लिए आसान हो सकता है):

package WhiteNoise "Package for generating white noise" 
    extends Modelica.Icons.Library; 

    block NoiseNormal "Normally distributed random noise" 
    parameter Real mean=0 "Mean value of random noise"; 
    parameter Real stdev=1 "Standard deviation of random noise"; 
    parameter Real tSample=0.01 "Noise sample time"; 
    Modelica.Blocks.Interfaces.RealOutput y; 
    equation 
    when initial() then 
     WhiteNoise.initRandomNormal(); 
    end when; 
    when sample(0, tSample) then 
     y=mean + stdev*WhiteNoise.RandomNormal(time); 
    end when; 
    end NoiseNormal; 

    function initRandomNormal 
    external "C" ext_initRandomNormal() annotation(Include="#include \"ext_initRandNormal.c\""); 
    end initRandomNormal; 

    function RandomNormal 
    output Real y; 
    input Real u; 
    external "C" y=ext_RandomNormal(u) annotation(Include="#include \"ext_RandNormal.c\""); 
    end RandomNormal; 

end WhiteNoise; 

बाहरी कोड:

ext_intRandNormal.c

#include <math.h> 
#include <limits.h> 

void ext_initRandomNormal() 
{ 
    srand(time(NULL)); 
} 

ext_RandNormal.c

#include <math.h> 
#include <limits.h> 
double ext_RandomNormal(double timein) 

{ 
    unsigned int seed = 0; 
    double v1, v2, r; 

    timein /= 100; 
    seed = (timein - floor(timein)) * UINT_MAX; 

    do 
    { 
     v1 = 2 * ((double) rand()) /((double) RAND_MAX) - 1; 
     v2 = 2 * ((double) rand()) /((double) RAND_MAX) - 1; 
     r = v1 * v1 + v2 * v2; 
    } while((r >= 1.0) || (r == 0.0)); 

    return v1 * sqrt(- 2.0 * log(r)/r); 
}