HAPSRG / Mbed 2 deprecated HAPStail

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

Committer:
NaotoMorita
Date:
Fri Jul 08 08:25:42 2022 +0000
Revision:
127:839ae718713b
Parent:
122:2e77adcf0c0d
new

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 119:43ac44c68ff0 55 drobj = rc[3] + updateValues.dr_command;
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 116:248f2d8cc11e 60 //pc.printf("pitchobj: %f\r\n", pitchobj);
NaotoMorita 118:97ffe77b6f38 61
osaka 116:248f2d8cc11e 62 Matrix vihat = eskf.getVihat();
osaka 116:248f2d8cc11e 63 float vihat_norm = vihat(1, 1) * vihat(1, 1) + vihat(2, 1) * vihat(2, 1) + vihat(3, 1) * vihat(3, 1);
NaotoMorita 127:839ae718713b 64 if (vihat_norm > 9.0f || rc[2] > -0.8f)
NaotoMorita 127:839ae718713b 65 {
NaotoMorita 127:839ae718713b 66 pitchPID.setGain(6.36f, 10.6f*0.1f,0.0f);
NaotoMorita 127:839ae718713b 67 pitchratePID.setGain(0.9540f,0.0f,0.0f);
NaotoMorita 118:97ffe77b6f38 68 }else{
NaotoMorita 118:97ffe77b6f38 69 pitchPID.setGain(6.36f, 0.0f,0.0);
NaotoMorita 118:97ffe77b6f38 70 pitchPID.resetIntError();
NaotoMorita 118:97ffe77b6f38 71 pitchratePID.setGain(0.9540f,0.0f,0.0f);
NaotoMorita 127:839ae718713b 72
NaotoMorita 118:97ffe77b6f38 73 }
osaka 116:248f2d8cc11e 74
NaotoMorita 109:27ae949bc38e 75 yawratePID.setGain(2.0f,0.0f,0.0f);
NaotoMorita 108:e582f8bd4729 76 pitchPID.setProcessValue(rpy.y);
NaotoMorita 108:e582f8bd4729 77 pitchratePID.setProcessValue((gyro.y-gyroBias(2,1)));
NaotoMorita 107:78e6f7bee68e 78 yawratePID.setProcessValue(gyro.z-gyroBias(3,1));
NaotoMorita 57:923e3df16159 79 //目標値のセット
NaotoMorita 57:923e3df16159 80 pitchPID.setSetPoint(pitchobj); //目標の設定
NaotoMorita 71:62eb45ecffe9 81
cocorlow 63:851e96f54a86 82 de = pitchPID.compute()+pitchratePID.compute();
NaotoMorita 74:af3b7cedff56 83 if(de<-1.0f){de = -1.0f;}
NaotoMorita 74:af3b7cedff56 84 if(de>1.0f){de = 1.0f;}
osaka 116:248f2d8cc11e 85 //pc.printf("%f %f %f %f %f\r\n",pitchobj,deobj,de,drobj,dr);
NaotoMorita 107:78e6f7bee68e 86
NaotoMorita 107:78e6f7bee68e 87 dr = yawratePID.compute()+drobj;
NaotoMorita 107:78e6f7bee68e 88 if(dr<-1.0f){dr = -1.0f;}
NaotoMorita 107:78e6f7bee68e 89 if(dr>1.0f){dr = 1.0f;}
osaka 82:3daf57031f65 90
NaotoMorita 118:97ffe77b6f38 91 dTl = rc[2]+drobj*0.2f+updateValues.dT_command;
NaotoMorita 118:97ffe77b6f38 92 dTr = rc[2]-drobj*0.2f+updateValues.dT_command;
NaotoMorita 118:97ffe77b6f38 93
NaotoMorita 108:e582f8bd4729 94 scaledServoOut[0]=-de;
NaotoMorita 107:78e6f7bee68e 95 scaledServoOut[1]=dr;
NaotoMorita 118:97ffe77b6f38 96 scaledMotorOut[0]= dTl;
NaotoMorita 118:97ffe77b6f38 97 scaledMotorOut[1]= dTr;
cocorlow 56:888379912f81 98
NaotoMorita 127:839ae718713b 99 float LP_servo = 1.0;
cocorlow 56:888379912f81 100 for(int i = 0; i < sizeof(servoOut)/sizeof(servoOut[0]); i++)
cocorlow 56:888379912f81 101 {
cocorlow 56:888379912f81 102 servoOut[i] = LP_servo*(mapfloat(scaledServoOut[i],-1,1,servoPwmMin,servoPwmMax))+(1.0-LP_servo)*servoOut[i];
cocorlow 56:888379912f81 103 if(servoOut[i]<servoPwmMin)
cocorlow 56:888379912f81 104 {
cocorlow 56:888379912f81 105 servoOut[i] = servoPwmMin;
cocorlow 56:888379912f81 106 }
cocorlow 56:888379912f81 107 if(servoOut[i]>servoPwmMax)
cocorlow 56:888379912f81 108 {
cocorlow 56:888379912f81 109 servoOut[i] = servoPwmMax;
cocorlow 56:888379912f81 110 }
cocorlow 56:888379912f81 111 }
cocorlow 56:888379912f81 112
NaotoMorita 107:78e6f7bee68e 113 float LP_motor = 0.2;
NaotoMorita 107:78e6f7bee68e 114 for(int i = 0; i < sizeof(motorOut)/sizeof(motorOut[0]); i++)
NaotoMorita 107:78e6f7bee68e 115 {
NaotoMorita 107:78e6f7bee68e 116 motorOut[i] = LP_motor*(mapfloat(scaledMotorOut[i],-1,1,motorPwmMin,motorPwmMax))+(1.0-LP_motor)*motorOut[i];
NaotoMorita 107:78e6f7bee68e 117 if(motorOut[i]<motorPwmMin)
NaotoMorita 107:78e6f7bee68e 118 {
NaotoMorita 107:78e6f7bee68e 119 motorOut[i] = motorPwmMin;
NaotoMorita 107:78e6f7bee68e 120 }
NaotoMorita 107:78e6f7bee68e 121 if(motorOut[i]>motorPwmMax)
NaotoMorita 107:78e6f7bee68e 122 {
NaotoMorita 107:78e6f7bee68e 123 motorOut[i] = motorPwmMax;
NaotoMorita 107:78e6f7bee68e 124 }
NaotoMorita 107:78e6f7bee68e 125 }
NaotoMorita 107:78e6f7bee68e 126
cocorlow 56:888379912f81 127 servo.pulsewidth_us(servoOut[0]);
NaotoMorita 107:78e6f7bee68e 128 rudServo.pulsewidth_us(servoOut[1]);
NaotoMorita 107:78e6f7bee68e 129 motor1.pulsewidth_us(motorOut[0]);
NaotoMorita 107:78e6f7bee68e 130 motor2.pulsewidth_us(motorOut[1]);
NaotoMorita 73:be7a8b8188de 131 send2center();
NaotoMorita 118:97ffe77b6f38 132 send2pc();
cocorlow 56:888379912f81 133 }