MATLAB2009_1_1自行车的运动.docVIP

自行车轮装饰物的运动轨迹 {问题[数学模型] 假设路面曲线是y = f(x)，自行车在运动过程中，车轮在路面上作无滑动的滚动，与路面只有一个接触点。车轮可当作一个半径为R的圆，运动时与地面相切。饰物P离圆心的距离为r。如A1图所示，CP与圆相交于A，A是初始接触点。当圆滚过θ角之后，接触点移到A处，圆从C移到C处，饰物从P移到P处。 曲线弧元为 设A点的坐标为(x0，y0)，则有弧长关系 (1) 设圆心C的坐标为(xC，yC)，C在曲线的法线上，因此联立方程 解得圆心坐标 (2) (3) 设P的坐标为(x，y)，则 x = xC – rsin(θ - φ)，y = yC + rcos(θ - φ) 其中θ是圆心角，φ是法线与竖直方向的夹角，也是切线与水平方向的夹角，称为仰角。由(1)式可得圆心角 (4) 根据导数的定义可得仰角 φ = arctan[f (x0)] (5) 饰物的轨迹方程为 (6a) (6b) [解析](1)设y = 0，路面是一条直线，则f (x) = 0，φ = 0，θ = x0/R，饰物的轨迹方程为 x = x0 – rsin(x0/R)，y = R – rcos(x0/R) (7a) 或 x = Rθ – rsinθ，y = R – rcosθ (7b) 当r R时，上式是短摆线方程；如果r = R，上式就是典型的摆线方程。 [算法](1)取R为长度单位，则短摆线方程可化为 x* = θ – r*sinθ，y* = 1 – r*cosθ (7*) 其中x* = x/R，y* = y/R，r* = r/R。θ是自变量，r*是可调节的参数。 [图示](1)如M1_1_1a图所示，当r R时，短摆线处处光滑。如M1_1_1b图所示，当r = R时，摆线在x/R = 2nπ处出现尖点。如M1_1_1c图所示，当r R时，即：饰物在车轮之外(例如火车车轮边缘的情况)，摆线在x/R = 2nπ附近出现回旋。 M1_1_1a图 M1_1_1b图 [程序]zqy1_1_1bicycle.m %摆线的轨迹 clear %清除变量 r=input(请输入半径r/R:); %键盘输入半径比 th=linspace(0,2*2*pi,1000); %参数向量 x=th-r*sin(th); %横坐标向量 y=1-r*cos(th); %纵坐标向量 figure %建立图形窗口 plot([0,th(end)],[0,0]) %画横线 hold on %保持图像 plot(x,y) %画摆线 axis equal %使坐标间隔相等 grid on %加网格 fs=16; %字体大小 title(摆线的轨迹,FontSize,fs) %显示标题 xlabel(\itx/R,FontSize,fs) %x标签 ylabel(\ity/R,FontSize,fs) %y标签 text(0,0,[\itr/R\rm=,num2str(r)],FontSize,fs)%显示字符串 M1_1_1c图 M1_1_2图 [程序]zqy1_1_2bicycle.m %摆线的轨迹族 clear %清除变量 r=[0.5,1,1.5]; %半径比向量 th=linspace(0,2*2*pi,1000); %参