1 个投票

I want to execute a code to calculate Largest Lyapunov Exponent in time series. There is a "nlinfit" function in my m-file, but the result is this error:
Error using nlinfit (line 239)
No usable observations after removing NaNs in Y and
in the result of evaluating MODELFUN at the initial
value BETA0.
Error in lyaprosen (line 377)
beta =
nlinfit(K(1:Tl),L(1:Tl),@nonlin1,[betar;randn(1,1)]);
I searched this error, but there was not any results. Can anyone help me to solve this issue?
Thanks; Peyman

4 个评论
显示 2更早的评论 隐藏 2更早的评论

the cyclist
the cyclist 2016-2-11
Are you able to post your code and a small subset of data that will allow us to generate this error for ourselves?
Peyman Ghasemi
Peyman Ghasemi 2016-2-12
Yes, of course. please find the codes in the attachment. copy the files in your directory and then run this command:
[LLE lambda tau m]=lyaprosen(R,0,0);
It should calculate largest Lyapunov exponent.
Thank You.
the cyclist
the cyclist 2016-2-14
I was not able to run your code. I get the error
dist is not included in your installed products. These products offer 'dist':
Neural Network Toolbox
Error in lyaprosen (line 106)
Dmm=dist(EEMmm(:,1:k)');
I may still be able to help. If you halt your code just before the call to nlinfit [line 377 of lyaprosen], and save those variables into a *.mat file, then you could upload the file here and I could just load the variables and run the command.
Peyman Ghasemi
Peyman Ghasemi 2016-2-14
I did it. Because the workspace variables was bigger than 5 MB (They are 7.5 MB!). I upload them on my Dropbox. please download it from this link:
https://www.dropbox.com/s/udzrb8kxpdba221/Lyapunov-2.zip?dl=0
The "lyaprosen" was a function and I hadn't access to variables, So I modified it to a m-file and a function. please load the files in your directory, then load the variables and go to line 377 of "lyaprosen.m".
Thank you very much.

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

the cyclist
the cyclist 2016-2-15

0 个投票

The immediate cause of the problem with nlinfit() is that your variables K and L are empty arrays. They are empty because Tl is an empty array as well.
You can trace this back to the fact that Ldiff is all NaNs, which is due to Lm being NaN, and so on. (That's as far as I went.)
I suggest you use the debugger to diagnose the ultimate cause.

更多回答（1 个）

Peyman Ghasemi
Peyman Ghasemi 2016-5-26
编辑：Peyman Ghasemi 2016-5-26

1 个投票

Finally I solved the problem! My signals had a low domain and they were changed to NaN in the process of the Lyapanov code! The solution is just multiplying the signal in a proper number.

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Aaron Best
Aaron Best 2019-6-3
Hi Peyman,
I am running into the same error could you please go into a little more detial about how you solved this problem?

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Matrix Computations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by