Skip to content

Commit

Permalink
Modify gyro_data.m and README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
xinli committed Jun 22, 2017
1 parent 354631a commit 249cf8a
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 4 deletions.
2 changes: 1 addition & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@

![](https://github.com/XinLiGitHub/GyroAllan/raw/master/Software/untitled.bmp) <br>

量化噪声 X轴:0.458971 Y轴:0.564662 Z轴:1.118367 单位:deg/h <br>
量化噪声 X轴:0.458971 Y轴:0.564662 Z轴:1.118367 单位:urad <br>
角度随机游走 X轴:0.231291 Y轴:0.273737 Z轴:0.532305 单位:deg/h^0.5 <br>
零偏不稳定性 X轴:8.264598 Y轴:8.849770 Z轴:7.367757 单位:deg/h <br>
角速率游走 X轴:42.156012 Y轴:40.870664 Z轴:32.268721 单位:deg/h/h^0.5 <br>
Expand Down
6 changes: 3 additions & 3 deletions Software/gyro_data.m
Original file line number Diff line number Diff line change
Expand Up @@ -16,19 +16,19 @@
C(2, :) = nihe(A', (B(:,2)').^2, 2)';
C(3, :) = nihe(A', (B(:,3)').^2, 2)';

Q = sqrt(abs(C(:, 1) / 3)); %量化噪声,单位:deg/h
Q = sqrt(abs(C(:, 1) / 3)); %量化噪声,单位:urad
N = sqrt(abs(C(:, 2) / 1)) / 60; %角度随机游走,单位:deg/h^0.5
Bs = sqrt(abs(C(:, 3))) / 0.6643; %零偏不稳定性,单位:deg/h
K = sqrt(abs(C(:, 4) * 3)) * 60; %角速率游走,单位:deg/h/h^0.5
R = sqrt(abs(C(:, 5) * 2)) * 3600; %速率斜坡,单位:deg/h/h

fprintf('量化噪声 X轴:%f Y轴:%f Z轴:%f 单位:deg/h\n', Q(1), Q(2), Q(3));
fprintf('量化噪声 X轴:%f Y轴:%f Z轴:%f 单位:urad\n', Q(1), Q(2), Q(3));
fprintf('角度随机游走 X轴:%f Y轴:%f Z轴:%f 单位:deg/h^0.5\n', N(1), N(2), N(3));
fprintf('零偏不稳定性 X轴:%f Y轴:%f Z轴:%f 单位:deg/h\n', Bs(1), Bs(2), Bs(3));
fprintf('角速率游走 X轴:%f Y轴:%f Z轴:%f 单位:deg/h/h^0.5\n', K(1), K(2), K(3));
fprintf('速率斜坡 X轴:%f Y轴:%f Z轴:%f 单位:deg/h/h\n', R(1), R(2), R(3));

D(:, 1) = sqrt(C(1,1)*A.^(-2) + C(1,2)*A.^(-1) + C(1,3)*A.^(0) + C(1,4)*A.^(1) + C(1,5)*A.^(2)); %生成拟合函数
D(:, 1) = sqrt(C(1,1)*A.^(-2) + C(1,2)*A.^(-1) + C(1,3)*A.^(0) + C(1,4)*A.^(1) + C(1,5)*A.^(2)); %生成拟合函数
D(:, 2) = sqrt(C(2,1)*A.^(-2) + C(2,2)*A.^(-1) + C(2,3)*A.^(0) + C(2,4)*A.^(1) + C(2,5)*A.^(2));
D(:, 3) = sqrt(C(3,1)*A.^(-2) + C(3,2)*A.^(-1) + C(3,3)*A.^(0) + C(3,4)*A.^(1) + C(3,5)*A.^(2));

Expand Down

0 comments on commit 249cf8a

Please sign in to comment.