HAPSRG / Mbed 2 deprecated HAPStail

Dependencies:   mbed MatrixMath LPS25HB_I2C LSM9DS1 Matrix2 PIDcontroller LoopTicker SBUS_without_mainfile UsaPack solaESKF_wind Vector3 CalibrateMagneto FastPWM

Committer:
osaka
Date:
Mon Mar 28 09:11:55 2022 +0000
Revision:
116:248f2d8cc11e
Parent:
115:ceb9ef04bb26
Child:
117:8cedc59454c4
always I control

Who changed what in which revision?

UserRevisionLine numberNew contents of line
cocorlow 56:888379912f81 1 #include "global.hpp"
cocorlow 56:888379912f81 2
cocorlow 56:888379912f81 3 // 割り込まれた時点での出力(computeの結果)を返す関数
cocorlow 56:888379912f81 4 void calcServoOut()
cocorlow 56:888379912f81 5 {
cocorlow 56:888379912f81 6 if(sbus.failSafe == false)
cocorlow 56:888379912f81 7 {
cocorlow 56:888379912f81 8 // sbusデータの読み込み
cocorlow 56:888379912f81 9 for (int i = 0; i < 16;i ++)
cocorlow 56:888379912f81 10 {
cocorlow 56:888379912f81 11 rc[i] = 0.65f * mapfloat(float(sbus.getData(i)),368,1680,-1,1) + (1.0f - 0.65f) * rc[i]; // mapped input
cocorlow 56:888379912f81 12 }
cocorlow 56:888379912f81 13 }
NaotoMorita 62:ef10fd919f7b 14 //pc.printf("rc:%f\r\n",rc[1]);
NaotoMorita 57:923e3df16159 15
osaka 85:0b14a2a600fc 16 //姿勢角の取得
osaka 85:0b14a2a600fc 17 euler = eskf.computeAngles();
NaotoMorita 107:78e6f7bee68e 18 Matrix gyroBias = eskf.getGyroBias();
osaka 85:0b14a2a600fc 19 rpy.x = euler(1,1);
osaka 85:0b14a2a600fc 20 rpy.y = euler(2,1);
osaka 85:0b14a2a600fc 21 rpy.z = euler(3,1);
osaka 85:0b14a2a600fc 22
osaka 88:0fc5df2fddcb 23 if (updateValues.calibrationFlag == 1111)
osaka 88:0fc5df2fddcb 24 {
osaka 88:0fc5df2fddcb 25 alignment();
osaka 88:0fc5df2fddcb 26 }
osaka 88:0fc5df2fddcb 27 else
osaka 88:0fc5df2fddcb 28 {
osaka 88:0fc5df2fddcb 29 rpy.x -= roll_offset;
osaka 88:0fc5df2fddcb 30 rpy.y -= pitch_offset;
osaka 88:0fc5df2fddcb 31 }
osaka 88:0fc5df2fddcb 32
NaotoMorita 57:923e3df16159 33 // 自身の位置に応じてエレベータ舵角を決定する
osaka 115:ceb9ef04bb26 34 float derc = -rc[1];
osaka 115:ceb9ef04bb26 35 //pc.printf("derc: %f\r\n", derc);
osaka 115:ceb9ef04bb26 36 float darc = -rc[0] + rc[1];
osaka 115:ceb9ef04bb26 37 //pc.printf("darc: %f\r\n", darc);
cocorlow 63:851e96f54a86 38 deobj = 0.0f;
NaotoMorita 57:923e3df16159 39 switch(pos_tail)
NaotoMorita 57:923e3df16159 40 {
NaotoMorita 57:923e3df16159 41 case 0: //left
osaka 87:981895e1d4f2 42 deobj = derc+darc;
NaotoMorita 57:923e3df16159 43 break;
NaotoMorita 57:923e3df16159 44 case 1: //center
NaotoMorita 71:62eb45ecffe9 45 deobj = derc;
NaotoMorita 57:923e3df16159 46 break;
NaotoMorita 57:923e3df16159 47 case 2: //right
osaka 87:981895e1d4f2 48 deobj = derc-darc;
NaotoMorita 57:923e3df16159 49 break;
NaotoMorita 57:923e3df16159 50 default: // error situation
NaotoMorita 57:923e3df16159 51 deobj = 0.0f;
NaotoMorita 57:923e3df16159 52 break;
NaotoMorita 57:923e3df16159 53 }
NaotoMorita 74:af3b7cedff56 54
NaotoMorita 107:78e6f7bee68e 55 drobj = rc[3];
NaotoMorita 107:78e6f7bee68e 56
NaotoMorita 74:af3b7cedff56 57
NaotoMorita 74:af3b7cedff56 58 float objgain = 15.0f*M_PI/180.0f;
osaka 115:ceb9ef04bb26 59 float pitchobj = objgain * (deobj + updateValues.de_command);
osaka 115:ceb9ef04bb26 60 if (pitchobj > objgain)
osaka 115:ceb9ef04bb26 61 {
osaka 115:ceb9ef04bb26 62 pitchobj = objgain;
osaka 115:ceb9ef04bb26 63 }
osaka 115:ceb9ef04bb26 64 if (pitchobj < -objgain)
osaka 115:ceb9ef04bb26 65 {
osaka 115:ceb9ef04bb26 66 pitchobj = -objgain;
osaka 115:ceb9ef04bb26 67 }
osaka 116:248f2d8cc11e 68 //pc.printf("pitchobj: %f\r\n", pitchobj);
osaka 116:248f2d8cc11e 69
osaka 116:248f2d8cc11e 70 Matrix vihat = eskf.getVihat();
osaka 116:248f2d8cc11e 71 float vihat_norm = vihat(1, 1) * vihat(1, 1) + vihat(2, 1) * vihat(2, 1) + vihat(3, 1) * vihat(3, 1);
osaka 116:248f2d8cc11e 72 //pc.printf("rc[2]: %f\r\n", rc[2]);
osaka 116:248f2d8cc11e 73 pitchPID.setGain(6.36f, 10.6f*0.5f,0.0);
osaka 116:248f2d8cc11e 74 pitchratePID.setGain(0.9540f,0.0f,0.0f);
osaka 116:248f2d8cc11e 75
NaotoMorita 109:27ae949bc38e 76 yawratePID.setGain(2.0f,0.0f,0.0f);
NaotoMorita 108:e582f8bd4729 77 pitchPID.setProcessValue(rpy.y);
NaotoMorita 108:e582f8bd4729 78 pitchratePID.setProcessValue((gyro.y-gyroBias(2,1)));
NaotoMorita 107:78e6f7bee68e 79 yawratePID.setProcessValue(gyro.z-gyroBias(3,1));
NaotoMorita 57:923e3df16159 80 //目標値のセット
NaotoMorita 57:923e3df16159 81 pitchPID.setSetPoint(pitchobj); //目標の設定
NaotoMorita 71:62eb45ecffe9 82
cocorlow 63:851e96f54a86 83 de = pitchPID.compute()+pitchratePID.compute();
NaotoMorita 74:af3b7cedff56 84 if(de<-1.0f){de = -1.0f;}
NaotoMorita 74:af3b7cedff56 85 if(de>1.0f){de = 1.0f;}
osaka 116:248f2d8cc11e 86 //pc.printf("%f %f %f %f %f\r\n",pitchobj,deobj,de,drobj,dr);
NaotoMorita 107:78e6f7bee68e 87
NaotoMorita 107:78e6f7bee68e 88 dr = yawratePID.compute()+drobj;
NaotoMorita 107:78e6f7bee68e 89 if(dr<-1.0f){dr = -1.0f;}
NaotoMorita 107:78e6f7bee68e 90 if(dr>1.0f){dr = 1.0f;}
osaka 82:3daf57031f65 91
NaotoMorita 108:e582f8bd4729 92 scaledServoOut[0]=-de;
NaotoMorita 107:78e6f7bee68e 93 scaledServoOut[1]=dr;
osaka 112:5a82800e1485 94 scaledMotorOut[0]=rc[2]+drobj*0.2f+updateValues.dT_command;
osaka 112:5a82800e1485 95 scaledMotorOut[1]=rc[2]-drobj*0.2f+updateValues.dT_command;
cocorlow 56:888379912f81 96
cocorlow 56:888379912f81 97 float LP_servo = 0.2;
cocorlow 56:888379912f81 98 for(int i = 0; i < sizeof(servoOut)/sizeof(servoOut[0]); i++)
cocorlow 56:888379912f81 99 {
cocorlow 56:888379912f81 100 servoOut[i] = LP_servo*(mapfloat(scaledServoOut[i],-1,1,servoPwmMin,servoPwmMax))+(1.0-LP_servo)*servoOut[i];
cocorlow 56:888379912f81 101 if(servoOut[i]<servoPwmMin)
cocorlow 56:888379912f81 102 {
cocorlow 56:888379912f81 103 servoOut[i] = servoPwmMin;
cocorlow 56:888379912f81 104 }
cocorlow 56:888379912f81 105 if(servoOut[i]>servoPwmMax)
cocorlow 56:888379912f81 106 {
cocorlow 56:888379912f81 107 servoOut[i] = servoPwmMax;
cocorlow 56:888379912f81 108 }
cocorlow 56:888379912f81 109 }
cocorlow 56:888379912f81 110
NaotoMorita 107:78e6f7bee68e 111 float LP_motor = 0.2;
NaotoMorita 107:78e6f7bee68e 112 for(int i = 0; i < sizeof(motorOut)/sizeof(motorOut[0]); i++)
NaotoMorita 107:78e6f7bee68e 113 {
NaotoMorita 107:78e6f7bee68e 114 motorOut[i] = LP_motor*(mapfloat(scaledMotorOut[i],-1,1,motorPwmMin,motorPwmMax))+(1.0-LP_motor)*motorOut[i];
NaotoMorita 107:78e6f7bee68e 115 if(motorOut[i]<motorPwmMin)
NaotoMorita 107:78e6f7bee68e 116 {
NaotoMorita 107:78e6f7bee68e 117 motorOut[i] = motorPwmMin;
NaotoMorita 107:78e6f7bee68e 118 }
NaotoMorita 107:78e6f7bee68e 119 if(motorOut[i]>motorPwmMax)
NaotoMorita 107:78e6f7bee68e 120 {
NaotoMorita 107:78e6f7bee68e 121 motorOut[i] = motorPwmMax;
NaotoMorita 107:78e6f7bee68e 122 }
NaotoMorita 107:78e6f7bee68e 123 }
NaotoMorita 107:78e6f7bee68e 124
cocorlow 56:888379912f81 125 servo.pulsewidth_us(servoOut[0]);
NaotoMorita 107:78e6f7bee68e 126 rudServo.pulsewidth_us(servoOut[1]);
NaotoMorita 107:78e6f7bee68e 127 motor1.pulsewidth_us(motorOut[0]);
NaotoMorita 107:78e6f7bee68e 128 motor2.pulsewidth_us(motorOut[1]);
NaotoMorita 73:be7a8b8188de 129 send2center();
cocorlow 56:888379912f81 130 }