1
00:00:10,840 --> 00:00:15,730
Let's recap all that we have learned in this course about common filters, starting with the linnear.

2
00:00:15,730 --> 00:00:22,060
Com and filter, the linnear common filter solves the estimation problem of estimating the current state

3
00:00:22,060 --> 00:00:28,240
of the system by minimizing the expected error over all measurements it receives subject to the system

4
00:00:28,240 --> 00:00:29,600
dynamics sharing here.

5
00:00:29,980 --> 00:00:31,770
So here is the estimated error.

6
00:00:31,930 --> 00:00:38,890
We want to minimize the estimated error for the expected error across all the measurements subject to

7
00:00:38,890 --> 00:00:41,220
the system module dynamics here.

8
00:00:42,100 --> 00:00:46,510
So you can see that these models are linear, which is why we're using the linear Matilda.

9
00:00:47,170 --> 00:00:52,660
The filter assumes that all the noise in the system or any of the error in the system follows these

10
00:00:52,660 --> 00:00:53,800
sets of assumptions here.

11
00:00:54,220 --> 00:00:59,800
So it assumes that each of them are uncorrelated and they can be modeled as a Gaussian distribution

12
00:01:00,010 --> 00:01:03,130
with zero meaning and a given covariance.

13
00:01:06,350 --> 00:01:12,260
The actual common of equations used to solve the filter are shown here, so inside the prediction step,

14
00:01:12,470 --> 00:01:16,580
we take the current state, we predict it forward using our process model.

15
00:01:16,850 --> 00:01:20,660
We add on the effects of the controls to predict the step forward in time.

16
00:01:21,050 --> 00:01:24,610
When we do this, we also want to predict the covariance forward in time.

17
00:01:24,620 --> 00:01:30,740
So we use again the process model and the previous experience and predict it forward in time and add

18
00:01:30,740 --> 00:01:33,170
on the uncertainty inside that prediction step.

19
00:01:33,920 --> 00:01:37,980
When we get the measurement, we apply the measurement set of equations here.

20
00:01:38,240 --> 00:01:42,230
So using the current state and the measurement model, we predict what the measurement should be.

21
00:01:42,410 --> 00:01:46,730
We look at the innovation, the difference between the actual measurement and our predicted measurement

22
00:01:46,880 --> 00:01:53,540
to calculate the innovation and then using this innovation, value and outcome and gain value here we

23
00:01:53,540 --> 00:01:59,990
can correct for the estimated state based on the error as we do the correction here, we also want to

24
00:01:59,990 --> 00:02:03,590
correct and update the covariance matrix using this last equation.

25
00:02:05,580 --> 00:02:07,110
So we're not running the common field.

26
00:02:07,410 --> 00:02:12,690
We also want to make sure the filters tuned correctly and one way is to validate it, and that is to

27
00:02:12,690 --> 00:02:18,300
look at the minimum mean squared error, which is the error between the common filter state estimates

28
00:02:18,300 --> 00:02:19,440
and the true states.

29
00:02:19,680 --> 00:02:24,450
And we want to make this as small as possible for the whole run in the filter like we've done in the

30
00:02:24,450 --> 00:02:25,200
simulation.

31
00:02:25,450 --> 00:02:29,190
We want to minimize the error between the estimates and the three states.

32
00:02:29,640 --> 00:02:31,500
So this is how we can validate the filter.

33
00:02:32,610 --> 00:02:37,470
We're not using the filter to make sure our main squared error is as small as possible.

34
00:02:37,770 --> 00:02:41,480
We want to make sure that the measurement innovation is always close to zero.

35
00:02:41,490 --> 00:02:42,260
So it has zero.

36
00:02:43,020 --> 00:02:46,830
So this is one of the operating conditions that the filter assumes.

37
00:02:46,830 --> 00:02:52,380
So we can calculate the innovation for a series of measurements and we want to make sure that is very

38
00:02:52,380 --> 00:02:53,050
close to zero.

39
00:02:53,050 --> 00:02:58,620
I mean, we also want to ensure that the measurement innovation variance is close to the actual innovation

40
00:02:58,620 --> 00:02:59,160
variance.

41
00:02:59,460 --> 00:03:04,980
So inside the common filter, we calculate our innovation and what we want to do is to make sure the

42
00:03:04,980 --> 00:03:11,460
actual innovation covariance matches what the measured innovation covariance is from the series of innovation

43
00:03:11,460 --> 00:03:13,140
measurements over the time period.

44
00:03:14,590 --> 00:03:19,960
So once these two conditions are met, we have a well working and true in common field which should

45
00:03:19,960 --> 00:03:23,530
minimize the main square area over the whole profile.

46
00:03:26,570 --> 00:03:32,330
Now for the extended calm and filter, the extended calm, I feel that solves the same problem as linear

47
00:03:32,330 --> 00:03:37,820
coming theater, except now we're using a non-linear system model, dynamic equations showing here.

48
00:03:38,000 --> 00:03:41,960
So instead of the linear model, we now have these general nonlinear models.

49
00:03:42,440 --> 00:03:47,050
We still use all the same assumptions for the noise and error variances.

50
00:03:47,510 --> 00:03:51,710
However, we've only change these equations over here for the system model.

51
00:03:52,760 --> 00:03:56,480
The actual equations used to implement the common filter are very similar.

52
00:03:57,230 --> 00:04:00,520
Now we use the full nonlinear system to predict forward.

53
00:04:00,770 --> 00:04:06,890
However, when we do the covariance prediction forward, we're using the Jacobean of these nonlinear

54
00:04:06,890 --> 00:04:07,330
functions.

55
00:04:07,330 --> 00:04:14,360
So the linear approximation of the nonlinear function, instead of using the linear system model for

56
00:04:14,360 --> 00:04:18,680
the innovation, we use the full measurement model to calculate the innovation.

57
00:04:18,920 --> 00:04:24,320
But then for the rest of the equations here is very similar as the linear computer swapping out the

58
00:04:24,320 --> 00:04:30,320
linear measurement model with the linear approximation of the Caribbean for the measurement model.

59
00:04:30,770 --> 00:04:35,780
So the equation shown on this slide here, talk about the Jacobean of the different process model and

60
00:04:35,780 --> 00:04:36,560
measurement model.

61
00:04:36,860 --> 00:04:43,580
So the Jacobean is simply just a partial derivative of the output of the process model with respect

62
00:04:43,580 --> 00:04:44,430
to the states.

63
00:04:44,450 --> 00:04:47,090
And this is for the state process model Jacobean.

64
00:04:47,510 --> 00:04:51,700
You can also calculate the Jacobean of the process model with respect to the noise input.

65
00:04:51,740 --> 00:04:57,290
So this term over here, by taking the partial derivative with respect to the noise for the measurement

66
00:04:57,290 --> 00:05:02,900
model, we have the measurement model Jacobean, which is the Jacobean of this function here with respect

67
00:05:02,900 --> 00:05:09,560
to the actual state inputs to this vector over here, the measurement noise Jakobsson model is the Jacobean

68
00:05:09,560 --> 00:05:11,480
of this measurement model.

69
00:05:11,660 --> 00:05:15,080
But respect to the input noise, the over here.

70
00:05:20,620 --> 00:05:26,110
And then lastly, for the unsent to come and filter, we want to take our state prediction augmented

71
00:05:26,140 --> 00:05:32,860
with the process, modernize, to form this augmented vector, and we do the same thing for the covariance

72
00:05:32,860 --> 00:05:33,310
matrix.

73
00:05:33,310 --> 00:05:36,370
We augment the current covariance with the noise covariance.

74
00:05:36,370 --> 00:05:37,840
The forms are going to covariance.

75
00:05:38,350 --> 00:05:45,730
Then we use these two augmented matrixes to calculate the signal points according to this set of equations

76
00:05:45,730 --> 00:05:46,120
here.

77
00:05:46,390 --> 00:05:52,570
So we use the square root of the covariance matrix to work out the deltas of the similar points between

78
00:05:52,570 --> 00:05:53,560
the main vector.

79
00:05:54,220 --> 00:05:59,170
Once we have the signal points, we then use these weightings and we can use them inside the prediction

80
00:05:59,170 --> 00:05:59,470
step.

81
00:05:59,660 --> 00:06:06,610
So we predict forward each sigma point through the process model and then we can and then we can recover

82
00:06:06,610 --> 00:06:10,900
our predicted state by taking the weighted estimate of all the single points.

83
00:06:10,900 --> 00:06:13,420
And we can do the same thing for the covariance matrix.

84
00:06:14,920 --> 00:06:20,140
When it comes to the measurement step, we do the same thing is that we augment the vector now with

85
00:06:20,140 --> 00:06:22,870
the measurement, noise and measurement covariance.

86
00:06:23,110 --> 00:06:28,300
We form the signal points again using our square root of the covariance matrix.

87
00:06:28,660 --> 00:06:33,520
And then we take these augmented signal points and we put them through the nonlinear measurement model.

88
00:06:34,150 --> 00:06:39,970
Once we do this, we can use the weights again to recover the predicted measurement, calculate the

89
00:06:39,970 --> 00:06:45,130
innovation, and then we can form the measurement covariance matrix and across covariance matrix by

90
00:06:45,130 --> 00:06:51,910
taking the weights and the signal points, again, taking the variances between the state and the measurement

91
00:06:51,910 --> 00:06:53,320
or the measurement and the measurement.

92
00:06:55,140 --> 00:06:59,460
Once we have all those information, then we can apply the uncertainty coming through the update step,

93
00:06:59,700 --> 00:07:06,810
which again is just these series of equations here, which are just a slightly different form of the

94
00:07:06,810 --> 00:07:10,290
extended of filter and linear common filter update step equations.

95
00:07:13,160 --> 00:07:19,580
So in summary, the Linnear midfielder is the best linnear estimate of a linear systems, it is stable

96
00:07:19,580 --> 00:07:23,240
for any initial conditions and it's stable for any error probation's.

97
00:07:23,780 --> 00:07:25,820
So this is one advantage of the linnear coming.

98
00:07:26,360 --> 00:07:29,480
If the system is linear, it has some very nice properties.

99
00:07:30,080 --> 00:07:34,480
Now, moving on to non-linear systems, we start to look at the extended midfielder.

100
00:07:34,970 --> 00:07:40,250
Now, unfortunately, this is not the best estimate of a non-linear systems because there is no best

101
00:07:40,250 --> 00:07:42,050
estimate of a non-linear systems.

102
00:07:42,830 --> 00:07:47,510
We can't guarantee stability and the filter can diverge because of the non-linear effects.

103
00:07:47,990 --> 00:07:53,810
The state estimates must always be close to the true state or the field or diverge because we're using

104
00:07:53,810 --> 00:07:56,530
linear approximations of the nonlinear system.

105
00:07:57,260 --> 00:08:03,740
So overall, this is a first order approximate accuracy for the nonlinear state and covariance transformations

106
00:08:03,980 --> 00:08:05,690
due to the Jacobins.

107
00:08:07,040 --> 00:08:11,780
So now looking at a different form for non-linear filters, we're looking at the unscented coming filter.

108
00:08:12,260 --> 00:08:17,930
So the unscented coming has very similar properties as the extended economy of Europe, except now it

109
00:08:17,930 --> 00:08:23,120
is third order approximation accuracy for the nonlinear state and covariance transforms.

110
00:08:23,400 --> 00:08:28,550
So we have a higher order of accuracy here, which is where the advantages of the uncertainty come a

111
00:08:28,550 --> 00:08:29,380
few to come from.

112
00:08:29,810 --> 00:08:35,540
And when we're using the unscented computer, it does not require the calculation of the Jacobean.

113
00:08:35,690 --> 00:08:41,510
We can use a process model and measurement model directly, which can sometimes aid in the computational

114
00:08:41,510 --> 00:08:45,080
complexity or at least aid in the implementation complexity.

115
00:08:47,100 --> 00:08:52,410
Congratulations on completing the course of all the advance common Philadelphia since the Fusion.

116
00:08:52,590 --> 00:08:56,040
I hope you have learned a lot and you have got quite a bit out of this cause.

117
00:08:57,730 --> 00:09:00,010
I hope to see you in one of my other courses.
