本文整理汇总了C++中rates函数的典型用法代码示例。如果您正苦于以下问题：C++ rates函数的具体用法？C++ rates怎么用？C++ rates使用的例子？那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了rates函数的15个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: rates
/*
* Attitude rates controller.
* Input: '_rates_sp' vector, '_thrust_sp'
* Output: '_att_control' vector
*/
void
MulticopterAttitudeControl::control_attitude_rates(float dt)
{
/* reset integral if disarmed */
if (!_armed.armed || !_vehicle_status.is_rotary_wing) {
_rates_int.zero();
}
/* current body angular rates */
math::Vector<3> rates;
rates(0) = _ctrl_state.roll_rate;
rates(1) = _ctrl_state.pitch_rate;
rates(2) = _ctrl_state.yaw_rate;
/* angular rates error */
math::Vector<3> rates_err = _rates_sp - rates;
_att_control = _params.rate_p.emult(rates_err) + _params.rate_d.emult(_rates_prev - rates) / dt + _rates_int +
_params.rate_ff.emult(_rates_sp - _rates_sp_prev) / dt;
_rates_sp_prev = _rates_sp;
_rates_prev = rates;
/* update integral only if not saturated on low limit and if motor commands are not saturated */
if (_thrust_sp > MIN_TAKEOFF_THRUST && !_motor_limits.lower_limit && !_motor_limits.upper_limit) {
for (int i = 0; i < 3; i++) {
if (fabsf(_att_control(i)) < _thrust_sp) {
float rate_i = _rates_int(i) + _params.rate_i(i) * rates_err(i) * dt;
if (PX4_ISFINITE(rate_i) && rate_i > -RATES_I_LIMIT && rate_i < RATES_I_LIMIT &&
_att_control(i) > -RATES_I_LIMIT && _att_control(i) < RATES_I_LIMIT) {
_rates_int(i) = rate_i;
}
}
}
}
}示例2: rates
/*
* Attitude rates controller.
* Input: '_rates_sp' vector, '_thrust_sp'
* Output: '_att_control' vector
*/
void
MulticopterAttitudeControl::control_attitude_rates(float dt)
{
/* reset integral if disarmed */
if (!_armed.armed) {
_rates_int.zero();
}
/* current body angular rates */
math::Vector<3> rates;
rates(0) = _v_att.rollspeed;
rates(1) = _v_att.pitchspeed;
rates(2) = _v_att.yawspeed;
/* angular rates error */
math::Vector<3> rates_err = _rates_sp - rates;
_att_control = _params.rate_p.emult(rates_err) + _params.rate_d.emult(_rates_prev - rates) / dt + _rates_int;
_rates_prev = rates;
/* update integral only if not saturated on low limit */
if (_thrust_sp > MIN_TAKEOFF_THRUST) {
for (int i = 0; i < 3; i++) {
if (fabsf(_att_control(i)) < _thrust_sp) {
float rate_i = _rates_int(i) + _params.rate_i(i) * rates_err(i) * dt;
if (isfinite(rate_i) && rate_i > -RATES_I_LIMIT && rate_i < RATES_I_LIMIT &&
_att_control(i) > -RATES_I_LIMIT && _att_control(i) < RATES_I_LIMIT) {
_rates_int(i) = rate_i;
}
}
}
}
}示例3: initmodel
static void initmodel(double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {
int _i; double _save;{
h_q = h_q0;
m_q = m_q0;
{
ena = 67.0 ;
temperature = celsius + 273.15 ;
rates ( _threadargs_ ) ;
rates ( _threadargs_ ) ;
m_q = m_inf ;
h_q = h_inf ;
}
}
}示例4: _ode_spec1
/*CVODE*/
static int _ode_spec1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {int _reset = 0; {
rates ( _threadargs_ ) ;
Dm_q = rate_m_q ;
Dh_q = rate_h_q ;
}
return _reset;
}示例5: main
int main()
{
double rate = 0.0;
double vol = 0.1;
int N = 3;
int M = 1;
int b = 700;
double maturity = 1.0;
double initialSpot = 100.0;
double timeStep = (maturity / (double)N);
double m = timeStep*(rate - vol*vol / 2.0);
//std::ofstream fichier("C:\\Users\\Lucien\\Desktop\\e.txt", std::ios::out | std::ios::trunc);
for (int i = 0; i < 40; i++)
{
double initialSpots[1]{80 + (double)i};
//False ditribution in the estimator
lognormal f(m, vol*sqrt(timeStep));
lognormal g(0.0, 5.0*vol);
normalGen gGenerator(ENG(time(0)), normalDist(0.0, 5.0*vol));
RandomGenerationTool kit(f, gGenerator, g);
BlackScholesPathGenerator paths(rate, vol, N, b, M, maturity, kit);
BlackScholesRate rates(N, maturity, rate);
AmericanPutPayoff payoff(100.0, maturity, N, b, M);
Estimator estimator(paths, payoff, rates);
std::cout << estimator.computePrice(initialSpots) << std::endl;
}
//fichier.close();
system("pause");
}示例6: states
/*END CVODE*/
static int states (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { {
rates ( _threadargscomma_ v ) ;
n = n + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / taun)))*(- ( ( ( ninf ) ) / taun ) / ( ( ( ( - 1.0) ) ) / taun ) - n) ;
l = l + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / taul)))*(- ( ( ( linf ) ) / taul ) / ( ( ( ( - 1.0) ) ) / taul ) - l) ;
}
return 0;
}示例7: _ode_matsol1
static int _ode_matsol1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {
rates ( _threadargscomma_ v , cai ) ;
Dn = Dn / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_n )) ;
Dh = Dh / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_h )) ;
Dc = Dc / (1. - dt*( ( ( ( - 1.0 ) ) ) / tau_c )) ;
return 0;
}示例8: states
/*END CVODE*/
static int states (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { {
rates ( _threadargs_ ) ;
m_q = m_q - dt*(- ( rate_m_q ) ) ;
h_q = h_q - dt*(- ( rate_h_q ) ) ;
}
return 0;
}示例9: _ode_spec1
/*CVODE*/
static int _ode_spec1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {int _reset = 0; {
rates ( _threadargscomma_ v ) ;
DX = ( Xinf - X ) / Xtau ;
DY = ( Yinf - Y ) / Ytau ;
}
return _reset;
}示例10: _ode_matsol1
static int _ode_matsol1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {
rates ( _threadargscomma_ v ) ;
Dm = Dm / (1. - dt*( ( ( ( - 1.0 ) ) ) / mtau )) ;
Dh = Dh / (1. - dt*( ( ( ( - 1.0 ) ) ) / htau )) ;
Dn = Dn / (1. - dt*( ( ( ( - 1.0 ) ) ) / ntau )) ;
return 0;
}示例11: _ode_spec1
/*CVODE*/
static int _ode_spec1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {int _reset = 0; {
rates ( _threadargscomma_ v ) ;
Dm = ( minf - m ) / mtau ;
Dh = ( hinf - h ) / htau ;
}
return _reset;
}示例12: _ode_spec1
/*CVODE*/
static int _ode_spec1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {int _reset = 0; {
rates ( _threadargs_ ) ;
Dm = ( mInf - m ) / mTau ;
Dh = ( hInf - h ) / hTau ;
}
return _reset;
}示例13: states
/*END CVODE*/
static int states (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) { {
rates ( _threadargs_ ) ;
m = m + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / mTau)))*(- ( ( ( mInf ) ) / mTau ) / ( ( ( ( - 1.0) ) ) / mTau ) - m) ;
h = h + (1. - exp(dt*(( ( ( - 1.0 ) ) ) / hTau)))*(- ( ( ( hInf ) ) / hTau ) / ( ( ( ( - 1.0) ) ) / hTau ) - h) ;
}
return 0;
}示例14: _ode_matsol1
static int _ode_matsol1 (double* _p, Datum* _ppvar, Datum* _thread, _NrnThread* _nt) {
rates ( _threadargscomma_ v ) ;
DX = DX / (1. - dt*( ( ( ( - 1.0 ) ) ) / Xtau )) ;
DY = DY / (1. - dt*( ( ( ( - 1.0 ) ) ) / Ytau )) ;
DZ = DZ / (1. - dt*( ( ( ( - 1.0 ) ) ) / Ztau )) ;
return 0;
}示例15: _ode_spec1
/*CVODE*/
static int _ode_spec1 () {_reset=0;
{
rates ( _threadargscomma_ v ) ;
Dn = ( ninf - n ) / taun ;
}
return _reset;
}