change kalman
Dependencies: ADXL345 AigamozuControlPackets HMC5843 ITG3200 MBed_Adafruit-GPS-Library XBee agzIDLIST mbed
Fork of Aigamozu_Robot_ver4 by
main.cpp@16:a07350b3eb64, 2015-05-15 (annotated)
- Committer:
- s1200058
- Date:
- Fri May 15 07:04:15 2015 +0000
- Revision:
- 16:a07350b3eb64
- Parent:
- 15:f3d01f37f00d
- Child:
- 17:e563cf8e6bcf
2015/05/15
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
kityann | 0:daab5accfd83 | 1 | /**********************************************/ |
kityann | 0:daab5accfd83 | 2 | // |
kityann | 0:daab5accfd83 | 3 | // |
kityann | 0:daab5accfd83 | 4 | // |
kityann | 0:daab5accfd83 | 5 | // Program name: Aigamozu ROBOT |
kityann | 0:daab5accfd83 | 6 | // Author: Mineta Kizuku |
kityann | 0:daab5accfd83 | 7 | // |
kityann | 0:daab5accfd83 | 8 | // |
kityann | 0:daab5accfd83 | 9 | /**********************************************/ |
kityann | 0:daab5accfd83 | 10 | |
kityann | 0:daab5accfd83 | 11 | /**********************************************/ |
kityann | 0:daab5accfd83 | 12 | //更新情報 |
kityann | 0:daab5accfd83 | 13 | //2015/05/11 |
kityann | 0:daab5accfd83 | 14 | //ロボットプログラムの作成 |
kityann | 0:daab5accfd83 | 15 | // |
kityann | 4:f36986ceb73d | 16 | //2015/05/13 |
kityann | 4:f36986ceb73d | 17 | //カルマンフィルタの共分散の値を0.0001以下にならないようにした |
kityann | 4:f36986ceb73d | 18 | //共分散の値を10進数に変換するようにした |
kityann | 0:daab5accfd83 | 19 | // |
s1200058 | 6:13f212b75e71 | 20 | //2015/05/13 |
s1200058 | 6:13f212b75e71 | 21 | //何回目のGPSでとられたGPSか確認するようにしました。 |
s1200058 | 6:13f212b75e71 | 22 | // |
kityann | 0:daab5accfd83 | 23 | /**********************************************/ |
kityann | 0:daab5accfd83 | 24 | |
kityann | 0:daab5accfd83 | 25 | #include "mbed.h" |
kityann | 0:daab5accfd83 | 26 | #include "XBee.h" |
kityann | 0:daab5accfd83 | 27 | #include "MBed_Adafruit_GPS.h" |
kityann | 0:daab5accfd83 | 28 | #include "AigamozuControlPackets.h" |
kityann | 0:daab5accfd83 | 29 | #include "agzIDLIST.h" |
kityann | 0:daab5accfd83 | 30 | #include "aigamozuSetting.h" |
kityann | 0:daab5accfd83 | 31 | #include "Kalman.h" |
kityann | 0:daab5accfd83 | 32 | |
kityann | 4:f36986ceb73d | 33 | #define SIGMA_MIN 0.0001 |
kityann | 4:f36986ceb73d | 34 | |
kityann | 0:daab5accfd83 | 35 | //************ID Number***************** |
kityann | 0:daab5accfd83 | 36 | //Robot ID: 'A' ~ 'Z' |
kityann | 0:daab5accfd83 | 37 | //Base ID: 'a' ~ 'a' |
kityann | 0:daab5accfd83 | 38 | //manager ID: '0'(数字のゼロ) |
kityann | 0:daab5accfd83 | 39 | // |
s1200058 | 12:48ef69b4f0e4 | 40 | const char MyID = 'C'; |
kityann | 0:daab5accfd83 | 41 | //************ID Number***************** |
kityann | 0:daab5accfd83 | 42 | |
kityann | 0:daab5accfd83 | 43 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 44 | // |
kityann | 0:daab5accfd83 | 45 | //Pin Setting |
kityann | 0:daab5accfd83 | 46 | // |
kityann | 0:daab5accfd83 | 47 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 48 | VNH5019 agz_motorShield(p21,p22,p23,p24,p25,p26); |
kityann | 0:daab5accfd83 | 49 | |
kityann | 0:daab5accfd83 | 50 | |
kityann | 0:daab5accfd83 | 51 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 52 | // |
kityann | 0:daab5accfd83 | 53 | //Connection Setting |
kityann | 0:daab5accfd83 | 54 | // |
kityann | 0:daab5accfd83 | 55 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 56 | |
kityann | 0:daab5accfd83 | 57 | //Serial Connect Setting: PC <--> mbed |
kityann | 0:daab5accfd83 | 58 | Serial pc(USBTX, USBRX); |
kityann | 0:daab5accfd83 | 59 | |
kityann | 0:daab5accfd83 | 60 | //Serial Connect Setting: GPS <--> mbed |
kityann | 0:daab5accfd83 | 61 | Serial * gps_Serial; |
kityann | 0:daab5accfd83 | 62 | |
kityann | 0:daab5accfd83 | 63 | //Serial Connect Setting: XBEE <--> mbed |
kityann | 0:daab5accfd83 | 64 | XBee xbee(p13,p14); |
kityann | 0:daab5accfd83 | 65 | ZBRxResponse zbRx = ZBRxResponse(); |
kityann | 0:daab5accfd83 | 66 | |
kityann | 0:daab5accfd83 | 67 | //set up GPS module |
kityann | 0:daab5accfd83 | 68 | |
kityann | 0:daab5accfd83 | 69 | //set up AigamozuControlPackets library |
kityann | 0:daab5accfd83 | 70 | AigamozuControlPackets agz(agz_motorShield); |
kityann | 0:daab5accfd83 | 71 | |
kityann | 0:daab5accfd83 | 72 | |
kityann | 0:daab5accfd83 | 73 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 74 | // |
kityann | 0:daab5accfd83 | 75 | //For Kalman data |
kityann | 0:daab5accfd83 | 76 | // |
kityann | 0:daab5accfd83 | 77 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 78 | double sigmaGPS[2][2] = {{250,0},{0,250}}; |
kityann | 0:daab5accfd83 | 79 | double K[2][2][2]= {{{1,0},{0,1}},{{1,0},{0,1}}}; |
kityann | 0:daab5accfd83 | 80 | double sigma[2][2][2]= {{{250,0},{0,250}},{{250,0},{0,250}}}; |
kityann | 0:daab5accfd83 | 81 | double y[2],x[2][2]={0}; |
kityann | 0:daab5accfd83 | 82 | void Kalman(double Latitude,double Longitude,Adafruit_GPS *myGPS); |
kityann | 0:daab5accfd83 | 83 | |
kityann | 0:daab5accfd83 | 84 | |
s1200058 | 2:886fac7f4399 | 85 | ///////////////////////////////////////// |
s1200058 | 2:886fac7f4399 | 86 | // |
s1200058 | 2:886fac7f4399 | 87 | //Plus Speed |
s1200058 | 2:886fac7f4399 | 88 | // |
s1200058 | 2:886fac7f4399 | 89 | //MNUAL_MODEの時にスピードを変える |
s1200058 | 2:886fac7f4399 | 90 | ///////////////////////////////////////// |
s1200058 | 2:886fac7f4399 | 91 | void Plus_Speed(uint8_t *packetdata){ |
s1200058 | 2:886fac7f4399 | 92 | |
s1200058 | 2:886fac7f4399 | 93 | if(agz.nowMode == MANUAL_MODE){ |
s1200058 | 2:886fac7f4399 | 94 | agz.changeSpeed(packetdata); |
s1200058 | 2:886fac7f4399 | 95 | } |
s1200058 | 2:886fac7f4399 | 96 | |
s1200058 | 2:886fac7f4399 | 97 | } |
kityann | 0:daab5accfd83 | 98 | |
kityann | 0:daab5accfd83 | 99 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 100 | // |
s1200058 | 2:886fac7f4399 | 101 | //Send Status |
kityann | 0:daab5accfd83 | 102 | // |
kityann | 0:daab5accfd83 | 103 | //リクエストがきたとき、自分の位置情報などを返信する |
kityann | 0:daab5accfd83 | 104 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 105 | void Send_Status(char SenderIDc){ |
kityann | 0:daab5accfd83 | 106 | XBeeAddress64 send_Address; |
kityann | 0:daab5accfd83 | 107 | if(SenderIDc == '0'){ |
kityann | 0:daab5accfd83 | 108 | send_Address = manager_Address; |
kityann | 0:daab5accfd83 | 109 | } |
kityann | 0:daab5accfd83 | 110 | if(SenderIDc >= 'A' && SenderIDc <= 'Z'){ |
kityann | 0:daab5accfd83 | 111 | send_Address = robot_Address[SenderIDc - 'A']; |
kityann | 0:daab5accfd83 | 112 | } |
kityann | 0:daab5accfd83 | 113 | if(SenderIDc >= 'a' && SenderIDc <= 'z'){ |
kityann | 0:daab5accfd83 | 114 | send_Address = base_Address[SenderIDc - 'a']; |
kityann | 0:daab5accfd83 | 115 | } |
kityann | 0:daab5accfd83 | 116 | //send normal data |
kityann | 0:daab5accfd83 | 117 | //Create GPS Infomation Packet |
s1200058 | 15:f3d01f37f00d | 118 | agz.createReceiveStatusCommand(MyID,SenderIDc, (int)agz.gpsAuto(), |
kityann | 0:daab5accfd83 | 119 | agz.get_agzPoint_lati(),agz.get_agzPoint_longi(), |
kityann | 0:daab5accfd83 | 120 | agz.get_agzPointKalman_lati(),agz.get_agzPointKalman_longi(), |
kityann | 0:daab5accfd83 | 121 | agz.get_agzCov_lati(),agz.get_agzCov_longi()); |
kityann | 0:daab5accfd83 | 122 | |
kityann | 0:daab5accfd83 | 123 | //debug*************************************************** |
kityann | 0:daab5accfd83 | 124 | printf("latitude:%f,longitude:%f\nlatitudeK:%f,longitudeK:%f\nCovlat:%f,Covlongi:%f\n", |
kityann | 0:daab5accfd83 | 125 | agz.get_agzPoint_lati(),agz.get_agzPoint_longi(), |
kityann | 0:daab5accfd83 | 126 | agz.get_agzPointKalman_lati(),agz.get_agzPointKalman_longi(), |
kityann | 0:daab5accfd83 | 127 | agz.get_agzCov_lati(),agz.get_agzCov_longi() |
kityann | 0:daab5accfd83 | 128 | ); |
kityann | 0:daab5accfd83 | 129 | for(int i = 0; i < RECEIVE_STATUS_COMMNAD_LENGTH; ++i) printf("%d ",agz.packetData[i]); |
kityann | 0:daab5accfd83 | 130 | printf("\n"); |
kityann | 0:daab5accfd83 | 131 | //debug end*************************************************** |
kityann | 0:daab5accfd83 | 132 | |
kityann | 0:daab5accfd83 | 133 | //Select Destination |
kityann | 0:daab5accfd83 | 134 | ZBTxRequest tx64request(send_Address,agz.packetData,agz.getPacketLength()); |
kityann | 0:daab5accfd83 | 135 | //Send -> Base |
kityann | 0:daab5accfd83 | 136 | xbee.send(tx64request); |
kityann | 0:daab5accfd83 | 137 | } |
kityann | 0:daab5accfd83 | 138 | |
kityann | 0:daab5accfd83 | 139 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 140 | // |
kityann | 0:daab5accfd83 | 141 | //Get GPS function |
kityann | 0:daab5accfd83 | 142 | // |
kityann | 0:daab5accfd83 | 143 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 144 | |
kityann | 0:daab5accfd83 | 145 | void Get_GPS(Adafruit_GPS *myGPS){ |
kityann | 0:daab5accfd83 | 146 | static bool flag = true; |
kityann | 0:daab5accfd83 | 147 | |
kityann | 0:daab5accfd83 | 148 | if (myGPS->fix) { |
kityann | 0:daab5accfd83 | 149 | agz.nowStatus = GPS_AVAIL; |
kityann | 0:daab5accfd83 | 150 | |
kityann | 0:daab5accfd83 | 151 | if(flag){//初期値設定 |
kityann | 0:daab5accfd83 | 152 | if(myGPS->latitudeH>=36 && myGPS->latitudeH<=38 && myGPS->longitudeH>=138 && myGPS->longitudeH<=141){ |
kityann | 0:daab5accfd83 | 153 | flag = false; |
kityann | 0:daab5accfd83 | 154 | x[0][0]=(double)myGPS->latitudeL; |
kityann | 0:daab5accfd83 | 155 | x[0][1]=(double)myGPS->longitudeL; |
kityann | 0:daab5accfd83 | 156 | } |
kityann | 0:daab5accfd83 | 157 | } |
kityann | 0:daab5accfd83 | 158 | |
kityann | 0:daab5accfd83 | 159 | if(myGPS->longitudeH<138 || myGPS->longitudeH>141 || myGPS->latitudeH<36 || myGPS->latitudeH>38){ |
kityann | 0:daab5accfd83 | 160 | return; |
kityann | 0:daab5accfd83 | 161 | } |
kityann | 0:daab5accfd83 | 162 | //Kalman Filter |
kityann | 0:daab5accfd83 | 163 | Kalman(myGPS->latitudeL,myGPS->longitudeL,myGPS); |
kityann | 0:daab5accfd83 | 164 | |
kityann | 0:daab5accfd83 | 165 | agz.reNewRobotPoint(myGPS->latitudeH,myGPS->latitudeL,myGPS->longitudeH,myGPS->longitudeL); |
kityann | 0:daab5accfd83 | 166 | agz.reNewRobotPointKalman(myGPS->latitudeKH,myGPS->latitudeKL,myGPS->longitudeKH,myGPS->longitudeKL); |
kityann | 0:daab5accfd83 | 167 | } |
kityann | 0:daab5accfd83 | 168 | else agz.nowStatus = GPS_UNAVAIL; |
kityann | 0:daab5accfd83 | 169 | |
kityann | 0:daab5accfd83 | 170 | } |
s1200058 | 2:886fac7f4399 | 171 | |
s1200058 | 2:886fac7f4399 | 172 | ///////////////////////////////////////// |
s1200058 | 2:886fac7f4399 | 173 | // |
s1200058 | 2:886fac7f4399 | 174 | //New Mode |
s1200058 | 2:886fac7f4399 | 175 | // |
s1200058 | 2:886fac7f4399 | 176 | ///////////////////////////////////////// |
s1200058 | 2:886fac7f4399 | 177 | |
s1200058 | 2:886fac7f4399 | 178 | void New_Mode(uint8_t *packetdata){ |
s1200058 | 2:886fac7f4399 | 179 | |
s1200058 | 2:886fac7f4399 | 180 | agz.changeMode(packetdata); |
s1200058 | 2:886fac7f4399 | 181 | |
s1200058 | 2:886fac7f4399 | 182 | } |
kityann | 1:b2b950b916ce | 183 | |
kityann | 1:b2b950b916ce | 184 | ///////////////////////////////////////// |
kityann | 1:b2b950b916ce | 185 | // |
kityann | 1:b2b950b916ce | 186 | //Get Status |
kityann | 1:b2b950b916ce | 187 | // |
kityann | 1:b2b950b916ce | 188 | ///////////////////////////////////////// |
kityann | 1:b2b950b916ce | 189 | void Get_Status(char SenderIDc,uint8_t *packetdata){ |
kityann | 1:b2b950b916ce | 190 | |
kityann | 1:b2b950b916ce | 191 | //マネージャからデータが来たとき |
kityann | 1:b2b950b916ce | 192 | if(SenderIDc == '0'){ |
kityann | 1:b2b950b916ce | 193 | printf("get manager Status\n"); |
kityann | 1:b2b950b916ce | 194 | } |
kityann | 1:b2b950b916ce | 195 | //他のロボットからデータが来たとき |
kityann | 1:b2b950b916ce | 196 | if(SenderIDc >= 'A' && SenderIDc <= 'Z'){ |
kityann | 1:b2b950b916ce | 197 | printf("get other robots Status\n"); |
kityann | 1:b2b950b916ce | 198 | } |
kityann | 1:b2b950b916ce | 199 | //基地局からデータが来たとき |
kityann | 1:b2b950b916ce | 200 | if(SenderIDc >= 'a' && SenderIDc <= 'z'){ |
kityann | 1:b2b950b916ce | 201 | printf("Get Base data\n"); |
kityann | 1:b2b950b916ce | 202 | int id = SenderIDc - 'a'; |
kityann | 1:b2b950b916ce | 203 | agz.reNewBasePoint(id,&packetdata[13],&packetdata[21]); |
kityann | 1:b2b950b916ce | 204 | agz.reNewBasePointKalman(id,&packetdata[29],&packetdata[37]); |
kityann | 1:b2b950b916ce | 205 | |
kityann | 1:b2b950b916ce | 206 | //debug |
kityann | 1:b2b950b916ce | 207 | for(int i = 0;i < 4;i++){ |
kityann | 1:b2b950b916ce | 208 | printf("BASE:%d\n",i); |
kityann | 1:b2b950b916ce | 209 | printf("latitude:%f,longitude:%f\nlatitudeK:%f,longitudeK:%f\n", |
kityann | 1:b2b950b916ce | 210 | agz.get_basePoint_lati(i),agz.get_basePoint_longi(i), |
kityann | 1:b2b950b916ce | 211 | agz.get_basePointKalman_lati(i),agz.get_basePointKalman_longi(i) |
kityann | 1:b2b950b916ce | 212 | ); |
kityann | 1:b2b950b916ce | 213 | } |
kityann | 1:b2b950b916ce | 214 | } |
kityann | 1:b2b950b916ce | 215 | } |
kityann | 1:b2b950b916ce | 216 | |
s1200058 | 2:886fac7f4399 | 217 | void Get_Status_Kalman(char SenderIDc,uint8_t *packetdata){ |
s1200058 | 2:886fac7f4399 | 218 | |
s1200058 | 2:886fac7f4399 | 219 | //マネージャからデータが来たとき |
s1200058 | 2:886fac7f4399 | 220 | if(SenderIDc == '0'){ |
s1200058 | 2:886fac7f4399 | 221 | printf("get manager Status Kalman\n"); |
s1200058 | 2:886fac7f4399 | 222 | } |
s1200058 | 2:886fac7f4399 | 223 | //他のロボットからデータが来たとき |
s1200058 | 2:886fac7f4399 | 224 | if(SenderIDc >= 'A' && SenderIDc <= 'Z'){ |
s1200058 | 2:886fac7f4399 | 225 | printf("get other robots Status Kalman\n"); |
s1200058 | 2:886fac7f4399 | 226 | } |
s1200058 | 2:886fac7f4399 | 227 | //基地局からデータが来たとき |
s1200058 | 2:886fac7f4399 | 228 | if(SenderIDc >= 'a' && SenderIDc <= 'z'){ |
s1200058 | 2:886fac7f4399 | 229 | printf("Get Base data Kalman\n"); |
s1200058 | 2:886fac7f4399 | 230 | int id = SenderIDc - 'a'; |
s1200058 | 2:886fac7f4399 | 231 | agz.reNewBasePointKalman(id,&packetdata[29],&packetdata[37]); |
s1200058 | 2:886fac7f4399 | 232 | |
s1200058 | 2:886fac7f4399 | 233 | //debug |
s1200058 | 2:886fac7f4399 | 234 | for(int i = 0;i < 4;i++){ |
s1200058 | 2:886fac7f4399 | 235 | printf("BASE:%d\n",i); |
s1200058 | 2:886fac7f4399 | 236 | printf("latitudeK:%f,longitudeK:%f\n", |
s1200058 | 2:886fac7f4399 | 237 | agz.get_basePointKalman_lati(i),agz.get_basePointKalman_longi(i)); |
s1200058 | 2:886fac7f4399 | 238 | } |
s1200058 | 2:886fac7f4399 | 239 | } |
s1200058 | 2:886fac7f4399 | 240 | } |
s1200058 | 2:886fac7f4399 | 241 | |
kityann | 1:b2b950b916ce | 242 | ///////////////////////////////////////// |
kityann | 1:b2b950b916ce | 243 | // |
kityann | 1:b2b950b916ce | 244 | //Send_Request_to_base |
kityann | 1:b2b950b916ce | 245 | // |
kityann | 1:b2b950b916ce | 246 | ///////////////////////////////////////// |
kityann | 1:b2b950b916ce | 247 | void Send_Request_Base(int basenumber){ |
kityann | 1:b2b950b916ce | 248 | printf("send\n"); |
kityann | 1:b2b950b916ce | 249 | agz.createRequestCommand(MyID, basenumber); |
kityann | 1:b2b950b916ce | 250 | //Select Destination |
kityann | 1:b2b950b916ce | 251 | ZBTxRequest tx64request(base_Address[basenumber],agz.packetData,agz.getPacketLength()); |
kityann | 1:b2b950b916ce | 252 | //Send -> Base |
kityann | 1:b2b950b916ce | 253 | xbee.send(tx64request); |
kityann | 1:b2b950b916ce | 254 | } |
kityann | 0:daab5accfd83 | 255 | |
s1200058 | 3:1ac506a96fd6 | 256 | ///////////////////////////////////////// |
s1200058 | 3:1ac506a96fd6 | 257 | // |
s1200058 | 3:1ac506a96fd6 | 258 | //auto_Move |
s1200058 | 3:1ac506a96fd6 | 259 | // |
s1200058 | 3:1ac506a96fd6 | 260 | //InAreaかOutAreaの判定 |
s1200058 | 3:1ac506a96fd6 | 261 | //Kalmanを通した値を出力(Baseと自分) |
s1200058 | 3:1ac506a96fd6 | 262 | ///////////////////////////////////////// |
s1200058 | 3:1ac506a96fd6 | 263 | |
s1200058 | 3:1ac506a96fd6 | 264 | void auto_Move(){ |
s1200058 | 3:1ac506a96fd6 | 265 | |
s1200058 | 3:1ac506a96fd6 | 266 | bool result; |
s1200058 | 3:1ac506a96fd6 | 267 | int i; |
s1200058 | 3:1ac506a96fd6 | 268 | |
s1200058 | 3:1ac506a96fd6 | 269 | result = agz.gpsAuto(); |
s1200058 | 3:1ac506a96fd6 | 270 | agz.set_agzAutoGPS(); |
s1200058 | 3:1ac506a96fd6 | 271 | agz.set_agzKalmanGPS(); |
s1200058 | 3:1ac506a96fd6 | 272 | |
s1200058 | 3:1ac506a96fd6 | 273 | if(result == true){ |
s1200058 | 3:1ac506a96fd6 | 274 | printf("Out Area\n"); |
s1200058 | 3:1ac506a96fd6 | 275 | } |
s1200058 | 3:1ac506a96fd6 | 276 | else if(result == false){ |
s1200058 | 3:1ac506a96fd6 | 277 | printf("In Area\n"); |
s1200058 | 3:1ac506a96fd6 | 278 | } |
s1200058 | 3:1ac506a96fd6 | 279 | for(i = 0; i < 4; i++){ |
s1200058 | 3:1ac506a96fd6 | 280 | printf("%d: %f, %f\n", i, agz.get_basePointKalman_lati(i), agz.get_basePointKalman_longi(i)); |
s1200058 | 3:1ac506a96fd6 | 281 | } |
s1200058 | 3:1ac506a96fd6 | 282 | printf("robot: %f, %f\n", agz.get_agzPointKalman_lati(), agz.get_agzPointKalman_longi()); |
s1200058 | 3:1ac506a96fd6 | 283 | |
s1200058 | 3:1ac506a96fd6 | 284 | } |
s1200058 | 3:1ac506a96fd6 | 285 | |
s1200058 | 5:522c47c78401 | 286 | void print_gps(int count){ |
s1200058 | 5:522c47c78401 | 287 | |
s1200058 | 5:522c47c78401 | 288 | printf("%d times:\n", count); |
s1200058 | 5:522c47c78401 | 289 | printf("%f, %f\n", agz.get_agzPointKalman_lati(), agz.get_agzPointKalman_longi()); |
s1200058 | 5:522c47c78401 | 290 | |
s1200058 | 5:522c47c78401 | 291 | } |
kityann | 0:daab5accfd83 | 292 | |
kityann | 0:daab5accfd83 | 293 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 294 | // |
kityann | 0:daab5accfd83 | 295 | //Kalman Processing |
kityann | 0:daab5accfd83 | 296 | // |
kityann | 0:daab5accfd83 | 297 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 298 | |
kityann | 0:daab5accfd83 | 299 | void get_K(){ |
kityann | 0:daab5accfd83 | 300 | double temp[2][2]={ |
kityann | 0:daab5accfd83 | 301 | {sigma[0][0][0]+sigmaGPS[0][0],sigma[0][0][1]+sigmaGPS[0][1]}, |
kityann | 0:daab5accfd83 | 302 | {sigma[0][1][0]+sigmaGPS[1][0],sigma[0][1][1]+sigmaGPS[1][1]} |
kityann | 0:daab5accfd83 | 303 | }; |
kityann | 0:daab5accfd83 | 304 | double ad_bc = temp[0][0]*temp[1][1]-temp[1][0]*temp[0][1]; |
kityann | 0:daab5accfd83 | 305 | K[1][0][0] = sigma[0][0][0]*(1/ad_bc)*(temp[1][1]); |
kityann | 0:daab5accfd83 | 306 | K[1][1][1] = sigma[0][1][1]*(1/ad_bc)*(temp[0][0]); |
kityann | 0:daab5accfd83 | 307 | } |
kityann | 0:daab5accfd83 | 308 | |
kityann | 0:daab5accfd83 | 309 | |
kityann | 0:daab5accfd83 | 310 | void get_x(){ |
kityann | 0:daab5accfd83 | 311 | x[1][0] = x[0][0]+K[1][0][0]*(y[0]-x[0][0]); |
kityann | 0:daab5accfd83 | 312 | x[1][1] = x[0][1]+K[1][1][1]*(y[1]-x[0][1]); |
kityann | 0:daab5accfd83 | 313 | } |
kityann | 0:daab5accfd83 | 314 | |
kityann | 0:daab5accfd83 | 315 | |
kityann | 0:daab5accfd83 | 316 | void get_sigma(){ |
kityann | 0:daab5accfd83 | 317 | double temp[2][2]; |
kityann | 0:daab5accfd83 | 318 | for(int i=0;i<2;i++) { |
kityann | 0:daab5accfd83 | 319 | for(int j=0;j<2;j++) { |
kityann | 0:daab5accfd83 | 320 | for(int k=0;k<2;k++) { |
kityann | 0:daab5accfd83 | 321 | temp[i][j]+=K[1][i][k]*sigma[0][k][j]; |
kityann | 0:daab5accfd83 | 322 | } |
kityann | 0:daab5accfd83 | 323 | } |
kityann | 0:daab5accfd83 | 324 | } |
kityann | 0:daab5accfd83 | 325 | for(int i = 0;i < 2;i++){ |
kityann | 0:daab5accfd83 | 326 | for(int j = 0;j < 2;j++){ |
kityann | 0:daab5accfd83 | 327 | sigma[1][i][j] = sigma[0][i][j]-temp[i][j]; |
kityann | 0:daab5accfd83 | 328 | } |
kityann | 0:daab5accfd83 | 329 | } |
kityann | 0:daab5accfd83 | 330 | } |
kityann | 0:daab5accfd83 | 331 | |
kityann | 0:daab5accfd83 | 332 | void Kalman(double Latitude,double Longitude,Adafruit_GPS *myGPS){ |
kityann | 0:daab5accfd83 | 333 | y[0] = Latitude; |
kityann | 0:daab5accfd83 | 334 | y[1] = Longitude; |
kityann | 0:daab5accfd83 | 335 | //K[t+1] = sigma[t]*Inverse[sigma[t]+sigmaGPS]; |
kityann | 0:daab5accfd83 | 336 | get_K(); |
kityann | 0:daab5accfd83 | 337 | //x[t+1] = x[t]+K[t+1]*(y[t*1]-x[t]); |
kityann | 0:daab5accfd83 | 338 | get_x(); |
kityann | 0:daab5accfd83 | 339 | //sigma[t+1] = sigma[t]-K[t+1]*sigma[t]; |
kityann | 0:daab5accfd83 | 340 | get_sigma(); |
kityann | 0:daab5accfd83 | 341 | |
kityann | 0:daab5accfd83 | 342 | |
kityann | 0:daab5accfd83 | 343 | //kousinn |
kityann | 0:daab5accfd83 | 344 | for(int i = 0;i < 2;i++){ |
kityann | 0:daab5accfd83 | 345 | for(int j = 0;j < 2;j++){ |
kityann | 0:daab5accfd83 | 346 | K[0][i][j]=K[1][i][j]; |
kityann | 0:daab5accfd83 | 347 | x[0][i]=x[1][i]; |
kityann | 0:daab5accfd83 | 348 | sigma[0][i][j]=sigma[1][i][j]; |
kityann | 0:daab5accfd83 | 349 | } |
kityann | 0:daab5accfd83 | 350 | } |
kityann | 0:daab5accfd83 | 351 | |
kityann | 4:f36986ceb73d | 352 | if(sigma[0][0][0] < SIGMA_MIN)sigma[0][0][0]=SIGMA_MIN; |
kityann | 4:f36986ceb73d | 353 | if(sigma[0][1][1] < SIGMA_MIN)sigma[0][1][1]=SIGMA_MIN; |
kityann | 4:f36986ceb73d | 354 | |
kityann | 0:daab5accfd83 | 355 | myGPS->latitudeKH=myGPS->latitudeH;//latitude after filtering |
kityann | 0:daab5accfd83 | 356 | myGPS->longitudeKH=myGPS->longitudeH;//longitude after filtering |
kityann | 0:daab5accfd83 | 357 | myGPS->latitudeKL=(long)x[1][0];//latitude after filtering |
kityann | 0:daab5accfd83 | 358 | myGPS->longitudeKL=(long)x[1][1];//longitude after filtering |
kityann | 0:daab5accfd83 | 359 | |
kityann | 0:daab5accfd83 | 360 | agz.set_agzCov(sigma[0][0][0],sigma[0][1][1]); |
kityann | 0:daab5accfd83 | 361 | } |
kityann | 0:daab5accfd83 | 362 | |
kityann | 0:daab5accfd83 | 363 | |
kityann | 0:daab5accfd83 | 364 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 365 | // |
kityann | 0:daab5accfd83 | 366 | //Main Processing |
kityann | 0:daab5accfd83 | 367 | // |
kityann | 0:daab5accfd83 | 368 | ///////////////////////////////////////// |
kityann | 0:daab5accfd83 | 369 | int main() { |
kityann | 0:daab5accfd83 | 370 | //start up time |
kityann | 0:daab5accfd83 | 371 | wait(3); |
kityann | 0:daab5accfd83 | 372 | //set pc frequency to 57600bps |
kityann | 0:daab5accfd83 | 373 | pc.baud(PC_BAUD_RATE); |
kityann | 0:daab5accfd83 | 374 | //set xbee frequency to 57600bps |
kityann | 0:daab5accfd83 | 375 | xbee.begin(XBEE_BAUD_RATE); |
kityann | 0:daab5accfd83 | 376 | |
kityann | 0:daab5accfd83 | 377 | |
kityann | 0:daab5accfd83 | 378 | //GPS setting |
kityann | 0:daab5accfd83 | 379 | gps_Serial = new Serial(p28,p27); |
kityann | 0:daab5accfd83 | 380 | Adafruit_GPS myGPS(gps_Serial); |
kityann | 0:daab5accfd83 | 381 | Timer refresh_Timer; |
s1200058 | 7:a9a1acc7673b | 382 | const int refresh_Time = 1000; //refresh time in ms |
s1200058 | 3:1ac506a96fd6 | 383 | Timer auto_Timer; |
s1200058 | 3:1ac506a96fd6 | 384 | const int auto_Time = 2000; //refresh time in ms |
s1200058 | 5:522c47c78401 | 385 | int count = 0; |
kityann | 0:daab5accfd83 | 386 | myGPS.begin(GPS_BAUD_RATE); |
kityann | 0:daab5accfd83 | 387 | |
s1200058 | 14:5deb7a4f1cd4 | 388 | Timer collect_Timer; |
s1200058 | 14:5deb7a4f1cd4 | 389 | const int collect_Time = 2000; //when we collect 4 GPS point, we use |
s1200058 | 14:5deb7a4f1cd4 | 390 | int collect_flag = 0; |
s1200058 | 14:5deb7a4f1cd4 | 391 | char collect_state = 'a'; |
s1200058 | 14:5deb7a4f1cd4 | 392 | |
kityann | 0:daab5accfd83 | 393 | char SenderIDc; |
kityann | 0:daab5accfd83 | 394 | //GPS Send Command |
kityann | 0:daab5accfd83 | 395 | myGPS.sendCommand(PMTK_SET_NMEA_OUTPUT_RMCONLY); |
kityann | 0:daab5accfd83 | 396 | myGPS.sendCommand(PMTK_SET_NMEA_UPDATE_1HZ); |
kityann | 0:daab5accfd83 | 397 | myGPS.sendCommand(PGCMD_ANTENNA); |
kityann | 0:daab5accfd83 | 398 | |
kityann | 0:daab5accfd83 | 399 | wait(2); |
kityann | 0:daab5accfd83 | 400 | |
kityann | 0:daab5accfd83 | 401 | //interrupt start |
kityann | 0:daab5accfd83 | 402 | refresh_Timer.start(); |
s1200058 | 3:1ac506a96fd6 | 403 | auto_Timer.start(); |
s1200058 | 16:a07350b3eb64 | 404 | collect_Timer.start(); |
kityann | 1:b2b950b916ce | 405 | |
kityann | 0:daab5accfd83 | 406 | printf("start\n"); |
kityann | 0:daab5accfd83 | 407 | |
kityann | 0:daab5accfd83 | 408 | while (true) { |
kityann | 0:daab5accfd83 | 409 | |
kityann | 0:daab5accfd83 | 410 | //Check Xbee Buffer Available |
kityann | 0:daab5accfd83 | 411 | xbee.readPacket(); |
kityann | 0:daab5accfd83 | 412 | |
kityann | 0:daab5accfd83 | 413 | if (xbee.getResponse().isAvailable()) { |
kityann | 0:daab5accfd83 | 414 | xbee.getResponse().getZBRxResponse(zbRx); |
kityann | 0:daab5accfd83 | 415 | uint8_t *buf = zbRx.getFrameData(); |
kityann | 0:daab5accfd83 | 416 | |
kityann | 0:daab5accfd83 | 417 | if (xbee.getResponse().getApiId() == ZB_RX_RESPONSE) { |
kityann | 0:daab5accfd83 | 418 | xbee.getResponse().getZBRxResponse(zbRx); |
kityann | 0:daab5accfd83 | 419 | uint8_t *buf = zbRx.getFrameData();//フレームデータを格納する |
kityann | 0:daab5accfd83 | 420 | uint8_t *buf1 = &buf[11];//データの部分のみを格納する |
kityann | 0:daab5accfd83 | 421 | SenderIDc = buf1[5];//送信元のIDを取得する |
kityann | 0:daab5accfd83 | 422 | char Command_type =agz.checkCommnadType(buf);//コマンドタイプを取得する |
kityann | 0:daab5accfd83 | 423 | |
kityann | 0:daab5accfd83 | 424 | //Check Command Type |
kityann | 0:daab5accfd83 | 425 | switch(Command_type){ |
kityann | 0:daab5accfd83 | 426 | //Get Request command |
s1200058 | 2:886fac7f4399 | 427 | case MANUAL:{ |
s1200058 | 12:48ef69b4f0e4 | 428 | Plus_Speed(buf); |
s1200058 | 2:886fac7f4399 | 429 | break; |
s1200058 | 2:886fac7f4399 | 430 | } |
kityann | 0:daab5accfd83 | 431 | case STATUS_REQUEST:{ |
kityann | 0:daab5accfd83 | 432 | Send_Status(SenderIDc); |
kityann | 0:daab5accfd83 | 433 | break; |
kityann | 0:daab5accfd83 | 434 | } |
s1200058 | 2:886fac7f4399 | 435 | case CHANGE_MODE:{ |
s1200058 | 12:48ef69b4f0e4 | 436 | New_Mode(buf); |
s1200058 | 2:886fac7f4399 | 437 | break; |
s1200058 | 2:886fac7f4399 | 438 | } |
kityann | 1:b2b950b916ce | 439 | case RECEIVE_STATUS:{ |
kityann | 1:b2b950b916ce | 440 | Get_Status(SenderIDc,buf1); |
kityann | 1:b2b950b916ce | 441 | break; |
kityann | 1:b2b950b916ce | 442 | } |
s1200058 | 2:886fac7f4399 | 443 | case RECEIVE_KALMAN:{ |
s1200058 | 2:886fac7f4399 | 444 | Get_Status_Kalman(SenderIDc, buf1); |
s1200058 | 2:886fac7f4399 | 445 | break; |
s1200058 | 2:886fac7f4399 | 446 | } |
kityann | 0:daab5accfd83 | 447 | default:{ |
kityann | 0:daab5accfd83 | 448 | break; |
kityann | 0:daab5accfd83 | 449 | } |
kityann | 0:daab5accfd83 | 450 | }//endswitch |
kityann | 0:daab5accfd83 | 451 | }//endifZB_RX_RESPONSE |
kityann | 0:daab5accfd83 | 452 | }//endifisAvailable |
kityann | 0:daab5accfd83 | 453 | |
kityann | 0:daab5accfd83 | 454 | myGPS.read(); |
kityann | 0:daab5accfd83 | 455 | //recive gps module |
kityann | 0:daab5accfd83 | 456 | //check if we recieved a new message from GPS, if so, attempt to parse it, |
kityann | 0:daab5accfd83 | 457 | if ( myGPS.newNMEAreceived() ) { |
kityann | 0:daab5accfd83 | 458 | if ( !myGPS.parse(myGPS.lastNMEA()) ) { |
s1200058 | 8:994c73f6fad9 | 459 | continue; |
s1200058 | 8:994c73f6fad9 | 460 | } |
s1200058 | 8:994c73f6fad9 | 461 | else{ |
s1200058 | 5:522c47c78401 | 462 | count++; |
s1200058 | 8:994c73f6fad9 | 463 | } |
kityann | 0:daab5accfd83 | 464 | } |
kityann | 1:b2b950b916ce | 465 | //一定時間ごとに自分のGPSデータを取得し、AigamozuControlPacketsないのagzPointとagzPointKalmanに格納する |
kityann | 0:daab5accfd83 | 466 | if (refresh_Timer.read_ms() >= refresh_Time) { |
kityann | 0:daab5accfd83 | 467 | refresh_Timer.reset(); |
s1200058 | 12:48ef69b4f0e4 | 468 | //print_gps(count); |
kityann | 0:daab5accfd83 | 469 | Get_GPS(&myGPS); |
kityann | 0:daab5accfd83 | 470 | |
kityann | 0:daab5accfd83 | 471 | } |
s1200058 | 14:5deb7a4f1cd4 | 472 | |
s1200058 | 14:5deb7a4f1cd4 | 473 | //get base GPS |
s1200058 | 14:5deb7a4f1cd4 | 474 | if( collect_Timer.read_ms() >= collect_Time){ |
s1200058 | 15:f3d01f37f00d | 475 | collect_Timer.reset(); |
s1200058 | 15:f3d01f37f00d | 476 | |
s1200058 | 15:f3d01f37f00d | 477 | Send_Request_Base(collect_flag); |
s1200058 | 14:5deb7a4f1cd4 | 478 | |
s1200058 | 14:5deb7a4f1cd4 | 479 | collect_flag++; |
s1200058 | 14:5deb7a4f1cd4 | 480 | collect_state++; |
s1200058 | 14:5deb7a4f1cd4 | 481 | |
s1200058 | 14:5deb7a4f1cd4 | 482 | if(collect_flag == 4){ |
s1200058 | 14:5deb7a4f1cd4 | 483 | collect_state = 'a'; |
s1200058 | 14:5deb7a4f1cd4 | 484 | collect_flag = 0; |
s1200058 | 14:5deb7a4f1cd4 | 485 | } |
s1200058 | 14:5deb7a4f1cd4 | 486 | } |
s1200058 | 3:1ac506a96fd6 | 487 | |
s1200058 | 14:5deb7a4f1cd4 | 488 | if(agz.nowMode == AUTO_GPS_MODE && auto_Timer.read_ms() >= auto_Time){ |
s1200058 | 3:1ac506a96fd6 | 489 | auto_Timer.reset(); |
s1200058 | 3:1ac506a96fd6 | 490 | auto_Move(); |
kityann | 1:b2b950b916ce | 491 | } |
kityann | 0:daab5accfd83 | 492 | } |
kityann | 0:daab5accfd83 | 493 | } |