動作的合成.doc

Size: px

Start display at page:

Download "動作的合成.doc"

院桐
7 years ago
Views:

1 MATLAB MATLAB MATLAB xyz - MATLAB - MATLAB MATLAB function demo_mot1(v) %demo_mot1.m -1-

2 x = -pi:pi/20:pi; t=0.1; h=plot(x,cos(x),'linewidth',5); delx=v*pi/360*t; while 1 x=x+delx; set(h,'xdata',x,'ydata',cos(x)); drawnow pause(t); %if x>2*pi, x=0;; 1. X cosx plot x,y plot xdata ydaya h (handle) MATLAB plot set -2-

3 set(h,'xdata',x,'ydata',cos(x)); x cos(x) x=x+delx; demo_mot1 (v) D (Lorenz) y(t) A y (Euler's method) %demo_mot2.m A = [ -8/3 0 0; ; ]; y = [ ]'; h = 0.01; grid on; p = plot3(y(1),y(2),y(3),'.',... 'EraseMode','none','MarkerSize',5); % Set EraseMode to none axis([ ]) hold on -3-

4 for i=1:4000 A(1,3) = y(2); A(3,1) = -y(2); ydot = A*y; y = y + h*ydot; set(p,'xdata',y(1),'ydata',y(2),'zdata',y(3)) % Change coordinates drawnow p = plot3(y(1),y(2),y(3),'.',... 'EraseMode','none','MarkerSize',5); % Set EraseMode to none p y xyz y(1)y(2)y(3) 2(!) MATLAB MATLAB MATLAB ( -4-

5 %demo_mot3.m p = plot3(y(1),y(2),y(3),'square',... 'EraseMode','background','MarkerSize',10,... 'MarkerEdgeColor',[1.7.7],'MarkerFaceColor',[1.7.7]); for i=1:4000 A(1,3) = y(2); A(3,1) = -y(2); ydot = A*y; y = y + h*ydot; set(p,'xdata',y(1),'ydata',y(2),'zdata',y(3)) drawnow hold off; "pink eraser" ( ) (Xor Erase Mode) ( '') n = 20 s =.02 x = rand(n,1)-0.5; y = rand(n,1)-0.5; -5-

6 h = plot(x,y,'.'); axis([ ]) axis square grid off set(h,'erasemode','xor','markersize',18) while 1 drawnow x = x + s*randn(n,1); y = y + s*randn(n,1); set(h,'xdata',x,'ydata',y) demo_mot4.m %Brownian motion clf; n=20; s=0.02; x = rand(n,1)-0.5; y = rand(n,1)-0.5; h = plot(x,y,'.'); axis([ ]) axis square grid off set(h,'erasemode','xor','markersize',18) grid on; title('press Ctl-C to stop'); while 1 drawnow x = x + s*randn(n,1); y = y + s*randn(n,1); set(h,'xdata',x,'ydata',y) -6-

7 Drawnow -7-

8 h = image(a(:,:,:,1),'erasemode','none'); for i = 2:x set(h,'cdata',a(:,:,:,i)) drawnow X = [ ; ; ]; -8-

9 h = image(x); colormap(colorcube(12)) xlabel x; ylabel y get(h,'xdata') ans = 1 4 get(h,'ydata') ans = 1 3 X = [ ; ; ]; image(x,'xdata',[-1 2],'YData',[2 4]); colormap(colorcube(12)) xlabel x; ylabel y -9-

10 nframes = 50; x = rand(n,1)-0.5; y = rand(n,1)-0.5; h = plot(x,y,'.'); set(h,'markersize',18); axis([ ]) axis square grid off for k = 1:nframes x = x + s*randn(n,1); y = y + s*randn(n,1); set(h,'xdata',x,'ydata',y) M(k) = getframe; movie(m,30) movie(m) movie(m,n) movie(m,n,fps) movie(h,...) movie(h,m,n,fps,loc) -10-

11 Z = peaks; surf(z); axis tight set(gca,'nextplot','replacechildren'); % Record the movie for j = 1:20 surf(sin(2*pi*j/20)*z,z) F(j) = getframe; % Play the movie twenty times movie(f,20) M = moviein(n) M = moviein(n,h) M = moviein(n,h,rect) -11-

12 FFT MATLAB axis equal M = moviein(16); set(gca,'nextplot','replacechildren') for j = 1:16 plot(fft(eye(j+16))) M(:,j) = getframe; set(gca,'nextplot','replacechildren') -12-

13 MATLAB ( QuickTime) 30 movie(m,30) ( uicontrols ) j axis equal M = moviein(16,gcf); set(gca,'nextplot','replacechildren') h = uicontrol('style','slider','position',... [ ],'Min',1,'Max',16) for j = 1:16 plot(fft(eye(j+16))) set(h,'value',j) M(:,j) = getframe(gcf); clf; axes('position',[ ]); movie(m,30) -13-

一课程基本情况课程名称工程应用数学 ( 计算机类 ) 编码 51611026 所属部门工业中心课程所属专业课程所属模块数学计算机类任课教师情况 ( 人数 ) 教授副教授讲师助教 3

一课程基本情况课程名称工程应用数学 ( 计算机类 ) 编码 51611026 所属部门工业中心课程所属专业课程所属模块数学计算机类任课教师情况 ( 人数 ) 教授副教授讲师助教 3 附表深圳职业技术学院文化育人示范课程建设项目申请书课程名称工程应用数学 ( 计算机类 ) 课程性质课程负责人所属专业所属部门基础课郑红数学工业中心填表日期 2015-12-28 深圳职业技术学院制二一五年十二月 1 一课程基本