预测代码的优化问题

原来是小飞

以下是用支持向量机预测混沌时间序列的一个代码，可是预测时间要好长，上午10:40运行，到下午16:55还是没有出结果，大家看下，可否帮忙优化下？代码如下：

1. clc;
   
2. clear all;
   
3. 
   
4. load data20.txt;
   
5. [nrow,ncolumn]=size(data20);
   
6. x=zeros(nrow,3);
   
7. tao=24*6;
   
8. for i=1:nrow
   
9.     x(i,1)=data20(i);   
   
10.     if mod(i,tao)==0
    
11.         x(i,2)=fix(i/tao);
    
12.         x(i,3)=24*6;
    
13.     else
    
14.         x(i,3)=mod(i,tao);
    
15.         x(i,2)=fix(i/tao)+1;
    
16.     end
    
17. end
    
18. tau=12;
    
19. m_demension=14;
    
20. % cur_day=30;
    
21. % datanum_day=144;
    
22. %preday_index=2;
    
23. p=3;
    
24. cur_index=2736;%x(t)中的t
    
25. n_used_datset=cur_index +p;
    
26. n_record=cur_index - (m_demension -1)*tau;
    
27. 
    
28. 
    
29. for i=1:n_record
    
30.     for j=1:m_demension
    
31.         row_t_index=cur_index -(i -1);
    
32.         jindex=row_t_index-(m_demension -j)*tau;
    
33.         x_jihe(2580-(i-1),j)=x(jindex,1);
    
34.     end
    
35.     kk=row_t_index + p;
    
36.     y_jihe(2580-(i-1),1)=x(kk,1);
    
37. end
    
38. x_train=x_jihe(1:(n_record-1),:);
    
39. y_train=y_jihe(1:(n_record-1),:);
    
40. x_fina= x_jihe(n_record,:);
    
41. y_fina= y_jihe(n_record,:);
    
42. %%%%yuce
    
43. y_pre_test=0;
    
44. e=0;
    
45. 
    
46. gam=10;
    
47. sig2=2;
    
48. [gam0,sig20]=bay_initlssvm({x_train,y_train,'f',gam,sig2,'RBF_kernel'});
    
49. [model,gamopt]=bay_optimize({x_train,y_train,'f',gam0,sig20,'RBF_kernel','preprocess'},2);
    
50. [cost_l3,sig2opt]=bay_optimize({x_train,y_train,'f',gamopt,sig20,'RBF_kernel','preprocess'},3);
    
51. gamopt
    
52. sig2opt
    
53. 
    
54. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
55. 
    
56. x_test=x_fina;
    
57. y_test=y_fina;
    
58. %%%%%%%%%%%%%%%%%%%%
    
59. type = 'function approximation';%这个引号里面的内容说明了用于回归运算，还有另一种查手册把。
    
60.    [alpha,b] = trainlssvm({x_train,y_train,type,gamopt,sig2opt,'RBF_kernel'});
    
61.     x_single_test=x_fina;
    
62.     y_pre_test= simlssvm({x_train,y_train,type,gamopt,sig2opt,'RBF_kernel','preprocess'},{alpha,b},x_single_test);
    
63.     actal=y_pre_test;
    
64.     %fe(i)=e(i)/y_test_fgyh*100;
    
65. % FCaculate=max(e);
    
66. % Femain=sum(abs(fe))/31
    
67. %%%%%1—————输出求得的极值
    
68. % plot(1:max_iter,actal,'-')
    
69. % hold on
    
70. % plot(1:max_iter,predict_load,'-.d')
    
71.    
    
72.    
    

复制代码

PS：x_jihe为2580x14的矩阵；y_jihe为2580x1的矩阵。谢谢大家了！急求