moler写的人行走的程序

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1. function walker(g)
   
2. % WALKER Human gait.
   
3. % This model, developed by Nikolaus Troje, is a five-term Fourier series
   
4. % with vector-valued coefficients that are the principal components for
   
5. % data obtained in motion capture experiments involving subjects wearing
   
6. % reflective markers walking on a treadmill. The components, which are
   
7. % also known as "postures" or "eigenwalkers", correspond to the static
   
8. % position, forward motion, sideways sway, and two hopping/bouncing
   
9. % movements that differ in the phase relationship between the upper and
   
10. % lower portions of the body. The postures are also classified by gender.
    
11. % Sliders allow you to vary the amount that each component contributes to
    
12. % the overall motion. A slider setting greater than 1.0 overemphasizes
    
13. % the characteristic. Can you see whether positive values of the gender
    
14. % coefficient correspond to male or female subjects?
    
15. %
    
16. % References:
    
17. % http://www.bml.psy.ruhr-uni-bochum.de/Demos
    
18. % http://www.biomotionlab.de/Text/WDP2002_Troje.pdf
    
19. % http://journalofvision.org/2/5/2
    
20. 
    
21. clf
    
22. shg
    
23. set(gcf,'doublebuf','on','color','w','name','Walker','numbertitle','off')
    
24. set(gca,'pos',get(gca,'pos')+[0 .07 0 0])
    
25. 
    
26. % The body is represented by 15 points in three space, i.e. a vector of
    
27. % length 45. The data consists of F, five vectors describing the average
    
28. % female and M, five vectors describing the average male. Four linked
    
29. % segments, indexed by L, are the head, torso, arms, and legs.
    
30. 
    
31. % Initial view
    
32. 
    
33. load walkers
    
34. X = reshape((F(:,1)+M(:,1))/2,15,3);
    
35. L = {[1 5],[5 12],[2 3 4 5 6 7 8],[9 10 11 12 13 14 15]};
    
36. for k = 1:4
    
37. p(k) = line(X(L{k},1),X(L{k},2),X(L{k},3),'marker','o', ...
    
38. 'markersize',10,'linestyle','-','erasemode','background');
    
39. end
    
40. set(p(1),'tag','head','userdata',zeros(1,3));
    
41. axis([-750 750 -750 750 0 1500])
    
42. set(gca,'xtick',[],'ytick',[],'ztick',[])
    
43. view(160,10)
    
44. 
    
45. % Sliders and controls
    
46. 
    
47. labels = {'speed','stride','sway','hop','bounce','gender'};
    
48. for j = 1:6
    
49. switch j
    
50. case 1, smin = 0; start = 1; smax = 3;
    
51. case 6, if nargin == 0, g = 0; end
    
52. smin = -3; start = g; smax = 3;
    
53. otherwise, smin = -2; start = 1; smax = 2;
    
54. end
    
55. txt = uicontrol('style','text','string',sprintf('%4.2f',start), ...
    
56. 'back','w','units','norm','pos',[.16*j-.10 .11 .08 .03]);
    
57. sliders(j) = uicontrol('style','slider','units','norm','back','w', ...
    
58. 'pos',[.16*j-.13 .07 .14 .03],'min',smin,'max',smax,'val',start, ...
    
59. 'sliderstep',[1/(4*smax),1/(10*smax)],'userdata',txt,'callback',...
    
60. 'set(get(gco,''userd''),''str'',sprintf(''%4.2f'',get(gco,''val'')))');
    
61. uicontrol('style','text','string',labels{j},'back','w', ...
    
62. 'units','norm','pos',[.16*j-.12 .02 .10 .04])
    
63. end
    
64. stop = uicontrol('style','toggle','units','norm','pos',[.91 .94 .08 .05], ...
    
65. 'backgr','w','fontw','bold','string','stop');
    
66. uicontrol('style','radio','units','norm','pos',[.015 .96 .03 .03], ...
    
67. 'userdata',H,'background','white', ...
    
68. 'callback',['p1 = findobj(''tag'',''head''); if get(gco,''val''),' ...
    
69. 'set(p1,''userd'',get(gco,''userd''),''marker'',''none''),' ...
    
70. 'else, set(p1,''userd'',zeros(1,3),''marker'',''o''), end']);
    
71. 
    
72. % uicontrol('style','text','units','norm','pos',[.72 .25 .25 .08], ...
    
73. % 'backgr','white','fontangle','italic', ...
    
74. % 'fontsize',get(0,'defaultuicontrolfontsize')-2, ...
    
75. % 'string',{'Click on the figure','to change the view'})
    
76. cameratoolbar setmode orbit
    
77. 
    
78. % Start walkin'...
    
79. 
    
80. period = 151.5751;
    
81. omega = 2*pi/period;
    
82. t = 0;
    
83. while get(stop,'value') == 0
    
84. s = cell2mat(get(sliders,'value'));
    
85. t = t + s(1);
    
86. c = [sin(omega*t); cos(omega*t); sin(2*omega*t); cos(2*omega*t)];
    
87. X = (F+M)/2 + s(6)*(F-M)/2;
    
88. w = [1; s(2:5).*c];
    
89. X = reshape(X*w,15,3);
    
90. H = get(p(1),'userdata');
    
91. e = ones(size(H,1),1);
    
92. XH = [H+X(e,:); X(5,:)];
    
93. set(p(1),'xdata',XH(:,1),'ydata',XH(:,2),'zdata',XH(:,3))
    
94. for k = 2:4
    
95. set(p(k),'xdata',X(L{k},1),'ydata',X(L{k},2),'zdata',X(L{k},3));
    
96. end
    
97. pause(.0001)
    
98. end;
    
99. cameratoolbar close
    
100. set(stop,'val',0,'str','close','fontw','bold','callb','close(gcf)')

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：Moler Mol OLE relationship coefficients 程序 Moler

高手真的是很多哦 见识一下

是啊，很强。
是从看他为其他人写“动态随机优化“的程序来找他的帖子的。
论坛里强人多啊。

牛人 厉害厉害