matlab代码运行的问题

原来是小飞

运行以下求解最大李雅普诺夫指数的程序是，总是提示：   ??? Undefined function or variable 'P'.请问如何解决这个问题。程序代码如下：

1. function lambda_1=largest_lyapunov_exponent(data,N,m,tau,P)
   
2. %the function is used to calcultate largest lyapunov exponent with the
   
3. %mended algorithm,which put forward by lv jing hu.
   
4. %data:the time series
   
5. %N:the length of data
   
6. %m:enbedding dimention
   
7. %tau:time delay
   
8. %P:the mean period of the time series,calculated with FFT
   
9. %lambda_1:return the largest lyapunov exponent
   
10. %skyhawk
    
11. delt_t=1;
    
12. Y=reconstitution(data,N,m,tau );%reconstitute state space
    
13. M=N-(m-1)*tau;%M is the number of embedded points in m-dimensional space
    
14. for j=1:M
    
15.     d_max=1e+100;
    
16.     for jj=1:M                                              %寻找相空间中每个点的最近距离点，并记下
    
17.         d_s=0;                                              %该点下标
    
18.         if abs(j-jj)>P                                      %限制短暂分离
    
19.             for i=1:m
    
20.                 d_s=d_s+(Y(i,j)-Y(i,jj))*(Y(i,j)-Y(i,jj));
    
21.                 d_min=d_max;
    
22.                 if d_s<d_min
    
23.                    d_min=d_s;
    
24.                    idx_j=jj;
    
25.                end
    
26.             end
    
27.         end
    
28.     end
    
29. %     index(j)=idx_j;
    
30.     max_i=min((M-j),(M-idx_j));%计算点j的最大演化时间步长i
    
31.     for k=1:max_i              %计算点j与其最近邻点在i个离散步后的距离
    
32.         d_j_i=0;
    
33.         for kk=1:m
    
34.             d_j_i=d_j_i+(Y(kk,j+k)-Y(kk,idx_j+k))*(Y(kk,j+k)-Y(kk,idx_j+k));
    
35.             d(k,j)=d_j_i;
    
36.         end
    
37.     end
    
38. end
    
39. 
    
40. %对每个演化时间步长i，求所有的j的lnd(i,j)平均
    
41. [l_i,l_j]=size(d);
    
42. for i=1:l_i
    
43.     q=0;
    
44.     y_s=0;
    
45.     for j=1:l_j
    
46.         if d(i,j)~=0
    
47.             q=q+1;
    
48.             y_s=y_s+log(d(i,j));
    
49.         end
    
50.     end
    
51.     y(i)=y_s/(q*delt_t);
    
52. end
    
53. x=1:length(y);
    
54. pp=polyfit(x,y,1);
    
55. lambda_1=pp(1);
    
56. yp=polyval(pp,x);
    
57. plot(x,y,'-o',x,yp,'--')

复制代码

本人的版本为64位win7，matlab7.6.0(R2008a)

321forever

P 是函数的输入之一，用的时候没有赋值

wanxiaoliu

木有看出来

szl1988226

321forever 发表于 2013-10-30 14:57

                               
登录/注册后可看大图

P 是函数的输入之一，用的时候没有赋值

强！好眼力

ChaChing

??? Undefined function or variable 'P'

Ref:
2.常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
From http://home.vibunion.com/blog-63979-18250.html

原来是小飞

321forever 发表于 2013-10-30 14:57

                               
登录/注册后可看大图

P 是函数的输入之一，用的时候没有赋值

赋值后，还是不行，求指导。

1. function lambda_1=largest_lyapunov_exponent(data,N,m,tau,P)
   
2. %the function is used to calcultate largest lyapunov exponent with the
   
3. %mended algorithm,which put forward by lv jing hu.
   
4. %data:the time series
   
5. %N:the length of data
   
6. %m:enbedding dimention
   
7. %tau:time delay
   
8. %P:the mean period of the time series,calculated with FFT
   
9. %lambda_1:return the largest lyapunov exponent
   
10. %skyhawk
    
11. data=load('b.txt');
    
12. N=length(data);
    
13. m=2;
    
14. tau=6;
    
15. p=6;
    
16. delt_t=1;
    
17. Y=reconstitution(data,N,m,tau );%reconstitute state space
    
18. M=N-(m-1)*tau;%M is the number of embedded points in m-dimensional space
    
19. for j=1:M
    
20.     d_max=1e+100;
    
21.     for jj=1:M                                              %寻找相空间中每个点的最近距离点，并记下
    
22.         d_s=0;                                              %该点下标
    
23.         if abs(j-jj)>P                                      %限制短暂分离
    
24.             for i=1:m
    
25.                 d_s=d_s+(Y(i,j)-Y(i,jj))*(Y(i,j)-Y(i,jj));
    
26.                 d_min=d_max;
    
27.                 if d_s<d_min
    
28.                    d_min=d_s;
    
29.                    idx_j=jj;
    
30.                end
    
31.             end
    
32.         end
    
33.     end
    
34. %     index(j)=idx_j;
    
35.     max_i=min((M-j),(M-idx_j));%计算点j的最大演化时间步长i
    
36.     for k=1:max_i              %计算点j与其最近邻点在i个离散步后的距离
    
37.         d_j_i=0;
    
38.         for kk=1:m
    
39.             d_j_i=d_j_i+(Y(kk,j+k)-Y(kk,idx_j+k))*(Y(kk,j+k)-Y(kk,idx_j+k));
    
40.             d(k,j)=d_j_i;
    
41.         end
    
42.     end
    
43. end
    
44. 
    
45. %对每个演化时间步长i，求所有的j的lnd(i,j)平均
    
46. [l_i,l_j]=size(d);
    
47. for i=1:l_i
    
48.     q=0;
    
49.     y_s=0;
    
50.     for j=1:l_j
    
51.         if d(i,j)~=0
    
52.             q=q+1;
    
53.             y_s=y_s+log(d(i,j));
    
54.         end
    
55.     end
    
56.     y(i)=y_s/(q*delt_t);
    
57. end
    
58. x=1:length(y);
    
59. pp=polyfit(x,y,1);
    
60. lambda_1=pp(1);
    
61. yp=polyval(pp,x);
    
62. plot(x,y,'-o',x,yp,'--')

复制代码

请帮忙看看，谢谢了。

ChaChing

1.求助完整格式：出错代码和出错提示
2.直觉楼主不是很清楚m-function与m-script的差异, 建议看下帮助文献

321forever

原来是小飞 发表于 2013-11-6 16:43

                               
登录/注册后可看大图

赋值后，还是不行，求指导。请帮忙看看，谢谢了。

lz 是怎么调用的程序

牛小贱

ChaChing 发表于 2013-11-2 09:19

                               
登录/注册后可看大图

Ref:
2.常见的程序出错问题整理 http://www.chinavib.com/thread-46001-1-1.html
From http://home.chi ...

It is wonderful!!谢谢分享