diff options
| author | leochanj105 <leochanj@live.unc.edu> | 2020-10-19 23:09:30 -0400 |
|---|---|---|
| committer | leochanj105 <leochanj@live.unc.edu> | 2020-10-20 02:40:39 -0400 |
| commit | f618466c25d43f3bae9e40920273bf77de1e1149 (patch) | |
| tree | 460e739e2165b8a9c37a9c7ab1b60f5874903543 /SD-VBS/benchmarks/localization/src/matlab/mclWhole.m | |
| parent | 47ced4e96bbb782b9e780e8f2cfc637b2c21ff44 (diff) | |
initial sd-vbs
initial sd-vbs
add sd-vbs
sd-vbs
Diffstat (limited to 'SD-VBS/benchmarks/localization/src/matlab/mclWhole.m')
| -rw-r--r-- | SD-VBS/benchmarks/localization/src/matlab/mclWhole.m | 18 |
1 files changed, 18 insertions, 0 deletions
diff --git a/SD-VBS/benchmarks/localization/src/matlab/mclWhole.m b/SD-VBS/benchmarks/localization/src/matlab/mclWhole.m new file mode 100644 index 0000000..4e78b8f --- /dev/null +++ b/SD-VBS/benchmarks/localization/src/matlab/mclWhole.m | |||
| @@ -0,0 +1,18 @@ | |||
| 1 | function retW=mclWhole(quat, Ovel, STDDEV_GPSVel, accl, STDDEV_ACCL) | ||
| 2 | n=size(quat,1); | ||
| 3 | OvelNorm=norm(Ovel); | ||
| 4 | if (OvelNorm>0.5) | ||
| 5 | Ovel=Ovel/norm(Ovel); | ||
| 6 | %trick | ||
| 7 | %quat=addEulNoise(quat, STDDEV_GPSVel); | ||
| 8 | orgWorld=quatRot([1 0 0],quatConj(quat)); | ||
| 9 | p1=get3DGaussianProb(orgWorld, Ovel, eye(3)./STDDEV_GPSVel); | ||
| 10 | else | ||
| 11 | p1=zeros(n,1); | ||
| 12 | end | ||
| 13 | |||
| 14 | gravity=quatRot(ones(n,1)*[0 0 -9.8], quat); | ||
| 15 | p2=get3DGaussianProb(gravity, accl, eye(3)./(STDDEV_ACCL)); | ||
| 16 | retW=p1+p2; | ||
| 17 | retW=retW/sum(retW); | ||
| 18 | \ No newline at end of file | ||
