mixpp: work/git/mixpp/pmsm/sim_profiles.h Source File

00001 /*
00002   \file
00003   \brief Common simulation profiles for PMSM experiments
00004   \author Vaclav Smidl.
00005 
00006   -----------------------------------
00007   BDM++ - C++ library for Bayesian Decision Making under Uncertainty
00008 
00009   Using IT++ for numerical operations
00010   -----------------------------------
00011 */
00012 
00013 void sim_profile_steps1(double &Ww,bool load=false) {
00014 
00015         if (t>0.2) x[8]=1.2;    // 1A //0.2ZP
00016         if (t>0.4) x[8]=10.8;   // 9A
00017         if (t>0.6) x[8]=25.2;  // 21A
00018 
00019         if (t>0.7) Ww=2.*M_PI*10.;
00020         if (t>1.0) x[8]=1.2;    // 1A
00021         if (t>1.2) x[8]=10.8;   // 9A
00022         if (t>1.4) x[8]=25.2;  // 21A
00023 
00024         if (t>1.6) Ww=2.*M_PI*50.;
00025         if (t>1.9) x[8]=1.2;    // 1A
00026         if (t>2.1) x[8]=10.8;   // 9A
00027         if (t>2.3) x[8]=25.2;  // 21A
00028 
00029         if (t>2.5) Ww=2.*M_PI*100;
00030         if (t>2.8) x[8]=1.2;    // 1A
00031         if (t>3.0) x[8]=10.8;   // 9A
00032         if (t>3.2) x[8]=25.2;  // 21A
00033 
00034         if (t>3.4) Ww=2.*M_PI*150;
00035         if (t>3.7) x[8]=1.2;    // 1A
00036         if (t>3.9) x[8]=10.8;   // 9A
00037         if (t>4.1) x[8]=25.2;  // 21A
00038 
00039         if (t>4.3) Ww=2.*M_PI*0;
00040         if (t>4.8) x[8]=-1.2;    // 1A
00041         if (t>5.0) x[8]=-10.8;   // 9A
00042         if (t>5.2) x[8]=-25.2;  // 21A
00043 
00044         if (t>5.4) Ww=2.*M_PI*(-10.);
00045         if (t>5.7) x[8]=-1.2;    // 1A
00046         if (t>5.9) x[8]=-10.8;   // 9A
00047         if (t>6.1) x[8]=-25.2;  // 21A
00048 
00049         if (t>6.3) Ww=2.*M_PI*(-50.);
00050         if (t>6.7) x[8]=-1.2;    // 1A
00051         if (t>6.9) x[8]=-10.8;   // 9A
00052         if (t>7.1) x[8]=-25.2;  // 21A
00053 
00054         if (t>7.3) Ww=2.*M_PI*(-100.);
00055         if (t>7.7) x[8]=-1.2;    // 1A
00056         if (t>7.9) x[8]=-10.8;   // 9A
00057         if (t>8.1) x[8]=-25.2;  // 21A
00058         if (t>8.3) x[8]=10.8;   // 9A
00059         if (t>8.5) x[8]=25.2;  // 21A
00060         
00061         if (t>9) Ww=2.*M_PI*0;
00062         //remove load if load=false
00063         if (!load) x[8]=0.0;
00064 }