完整word版,matlab课程设计.docVIP

物理与电子科学学院 摆线方程 MATLAB 仿真 实验报告 实验名称： 摆线方程 MATLB 仿真 试验日期： 2013 年12 月2 日 专 业： 电子信息工程 姓 名： 刘 斌 班 级： 物电 1105 班 学 号： 2011112030560 一、 实验目的 熟悉掌握数学方程仿真的技能， 即通过对生活实例摆线进行数学 建模，然后用计算机来模拟方程的大概形态，并进行仿真和分析。仿 真后针对模型分析、 算法程序和数学方程的核心问题做出评价， 对实 现实际生活中的复杂度和技术难度较大的摆线类似问题的解决积累 经验。 二、 实验内容 摆线的绘制。如下图所示： 当圆轮在平面上滚动时， 轮上任一点所画出的轨迹称为摆线。 如果这一点不在圆周上而在圆内， 则生成内摆线 ; 如果该点在圆外， 即离圆心距离大于半径，则生成外摆线。对于后一种情况，可由火车轮来想象。其接触轨道的部分，并不是直径最大处，其内侧的直径还要大一些，以防止车轮左右出轨。在这部分边缘上的点就形成外摆线。 通过建模，设定参变方程，完成直线摆线，内摆线和外摆线的绘制和仿真，得出仿真结果图样，加深对摆线的数学形式的深入了解。 三、 实验步骤 建 模： 概括几种情况，其普遍方程可表示为 : xA = rt-Rsin t yA = r-Rcos t 可由这组以 t 为参数的方程分析其轨迹 。 即用 matlab 以动画的方式绘制出摆线 x=a(t-sint) y=a(1-cost) （a 自己赋值）的渐屈线 MATLAB程序： 【1】直线摆线 function cycloid(r,n) fig=figure; set(fig,Color,[1,1,1]); r=100; n=4; theta=0:pi/100:2*pi; x=r*(n*theta-sin(n*theta));% 摆线 x 坐标 y=r*(1-cos(n*theta)); % 摆线 y 坐标 x0=r*sin(theta-pi);% 圆 x 坐标 y0=r*cos(theta-pi); % 圆 y 坐标 xa=[-r,8*n*r];% 坐标轴 x ya=[-2.5*r,5*r];% 坐标轴 y m=length(theta); mov=avifile(cycloid.avi); for i=1:m plot(xa,[0,0],k,[0,0],ya,k); hold on AA=[cos(n*theta(i)) sin(n*theta(i));-sin(n*theta(i)) cos(n*theta(i))]*[x0;y0]; xr=AA(1,:); yr=AA(2,:); xr=xr+r*n*theta(i); yr=yr+r; plot(xr,yr); plot(xr(1,[1,(m-1)/2]),yr(1,[1,(m-1)/2]),k); plot(xr(1,[(m-1)/4,(m-1)*3/4]),yr(1,[(m-1)/4,(m-1)*3/4]),k); plot(xr(1),yr(1),og,linewidth,2,markeredgecolor,k,markerfacecolor,g); plot(x(1:i),y(1:i),r,linewidth,2); axis equal axis off F=getframe(gcf); mov=addframe(mov,F); hold off end mov = close(mov); 【2】外圆摆线 function epicycloid(ra,rb) fig=figure; set(fig,Color,[1,1,1]); theta=0:pi/50:2*pi; xa=ra*cos(theta); % 定圆 x 坐标 ya=ra*sin(theta); % 定圆 y 坐标 xb=rb*cos(theta-pi); % 动圆 x 坐标 yb=rb*sin(theta-pi); % 动圆 y 坐标 xc=(ra+rb)*cos(theta)-rb*cos((ra+rb)*theta/rb); % 摆线 x 坐标 yc=(ra+rb)*sin(theta)-rb*sin((ra+rb)*theta/rb); % 摆线 y 坐标 xax=[-1.2*ra,1.2*ra]; yax=[-1.2*ra,1.2*ra]; m=length(theta); mov=avifile(epicycloid.avi); for i=1:m plot(xax,[0,0],k,[0 0],yax,k); hold on plot(xa,ya,k,linewidth,2); plot(xc(1,[1:i]),yc(1,[1:i]),r,linewidth,2); AA=[c