I understand you’re trying to estimate the clock offset and clock skew of Pairwise Broadcast Synchronization using Truncated Levy distributed random delays. However, the mean square error is not decreasing as expected, and so you’re seeking changes to improve the results.
To improve the results and decrease the MSE, you can try the following changes:
- Modify the optimization options: You can experiment with different optimization options to improve the convergence of the ‘lsqnonlin’ function. Try changing the line options = optimoptions(@lsqnonlin,'Algorithm','trust-region-reflective') to options = optimoptions(@lsqnonlin,'Algorithm','levenberg-marquardt') or options = optimoptions(@lsqnonlin,'Algorithm','sqp').
- Increase the number of random data points: You can increase the number of random data points to get a larger sample size and potentially improve the estimation accuracy. Modify the line ‘n = 100’ to a higher value, such as ‘n = 500’ or ‘n = 1000’.
- Adjust the range of random clock offset values: The code generates random clock offset values between -1010^-3 and 1010^-3. You can try narrowing down this range to focus on values closer to the expected offset. Modify the lines ‘a = -10*10^-3’ and ‘b = 10*10^-3‘ to smaller values, such as ‘a = -5*10^-3’ and ‘b = 5*10^-3’.
- Adjust the random clock skew values range: The code generates random clock skew values between 0.99 and 1.01. You can try narrowing down this range to focus on values closer to the expected skew. Modify the lines ‘a1 = 0.99’ and ‘b1 = 1.01’ to smaller values, such as ‘a1 = 0.995’ and ‘b1 = 1.005’.
