Guest User

Untitled

a guest
Jul 22nd, 2018
101
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.58 KB | None | 0 0
  1. for i=1:10
  2. for j=1:3
  3. ......
  4. r(1)=i-1;
  5. r(2)=j-1;
  6. r(3)=q-1;
  7. ....
  8.  
  9. fwrap = @(ab) wrap_linear_dynamic(r,ab);
  10. [ab] = fminsearch(fwrap,zeros(1,columns(ab)),optimset('MaxIter',2000,'Display','off'));
  11.  
  12.  
  13. function [result] = wrap_linear_dynamic(r,ab)
  14. global ayw shifty;
  15. b = ab(1, columns(ab)-r(columns(r)):columns(ab));
  16. a = ab(1, 1: columns(ayw)- columns(b));
  17. result=criterion_linear_dynamic(a,b);
  18. end
  19.  
  20. function [result]=criterion_linear_dynamic(a,b,Dy,Dx)
  21. ....
  22. global ayw shifty;
  23. .....
  24. u = uf(a,b);
  25. ...
  26. w = wf(a,b,Dx,Dy);
  27. result = u*(w^-1);
  28. clear w u Dx Dy;
  29. end
  30.  
  31. function [resultw]=wf(a,b,Dx,Dy)
  32. d = 1:columns(a);
  33. if max(Dx) > 0
  34. Do = (Dx.^2)/(Dy.^2);
  35. else
  36. Do = Dx;
  37. end
  38. resultw = Dy*(1+dot(b',b)+sum(Do(:,d)*dot(a(:,d)',a(:,d))));
  39. end
  40.  
  41. function [resultu]=uf(a,b)
  42. global ayw shifty;
  43. ab = vertcat(a',b');
  44. resultu = sum((shifty-ayw*(ab)).^2);
  45. end
  46.  
  47. ??? Error using ==> max
  48. Out of memory. Type HELP MEMORY for your options.
  49. Error in ==> criterion_linear_dynamic>wf at 68
  50. if max(Dx) > 0
  51. Error in ==> criterion_linear_dynamic at 35
  52. w = wf(a,b,Dx,Dy);
  53. Error in ==> wrap_linear_dynamic at 8
  54. result=criterion_linear_dynamic(a,b);
  55. Error in ==>
  56. optim_linear_dynamic>@(ab)wrap_linear_dynamic(r,ab) at 28
  57. fwrap = @(ab) wrap_linear_dynamic(r,ab);
  58. Error in ==> fminsearch at 191
  59. fv(:,1) = funfcn(x,varargin{:});
  60. Error in ==> optim_linear_dynamic at 29
  61. [ab] = fminsearch(fwrap,zeros(1,columns(ab)),optimset('MaxIter',2000,'Display','off'));%,'PlotFcns',@optimplotfval));
  62. Error in ==> optim_values_test7d at 21
  63. [curcrit,a,b,curcrit_test_data,std_error]
  64. = optim_linear_dynamic(x,y,rz);
Add Comment
Please, Sign In to add comment