2022年完整word版,matlab课程设计 .pdf

《2022年完整word版,matlab课程设计 .pdf》由会员分享，可在线阅读，更多相关《2022年完整word版,matlab课程设计 .pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、物理与电子科学学院摆线方程 MATLAB仿真实验报告实验名称：摆线方程 MATLB 仿真试验日期： 2013年 12 月 2 日专业：电子信息工程姓名：刘斌班级：物电 1105 班学号： 2011112030560 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 8 页 - - - - - - - - - 一、实验目的熟悉掌握数学方程仿真的技能， 即通过对生活实例摆线进行数学建模，然后用计算机来模拟方程的大概形态，并进行仿真和分析。仿真后针对模型分析、 算法程序和数学方程

2、的核心问题做出评价，对实现实际生活中的复杂度和技术难度较大的摆线类似问题的解决积累经验。二、实验内容摆线的绘制。如下图所示：当圆轮在平面上滚动时， 轮上任一点所画出的轨迹称为摆线。如果这一点不在圆周上而在圆内，则生成内摆线 ; 如果该点在圆外， 即离圆心距离大于半径，则生成外摆线。对于后一种情况，可由火车轮来想象。其接触轨道的部分，并不是直径最大处，其内侧的直径还要大一些，以防止车轮左右出轨。在这部分边缘上的点就形成外摆线。通过建模，设定参变方程，完成直线摆线，内摆线和外摆线的绘制和仿真，得出仿真结果图样，加深对摆线的数学形式的深入了解。名师资料总结 - - -精品资料欢迎下载 - - - -

3、 - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 8 页 - - - - - - - - - 三、实验步骤建 模：概括几种情况，其普遍方程可表示为: 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、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*thet

5、a(i);-sin(n*theta(i) cos(n*theta(i)*x0;y0; xr=AA(1,:); yr=AA(2,:); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 8 页 - - - - - - - - - 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,

6、(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=

7、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(epicyclo

8、id.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=cos(ra+rb)*theta(i)/rb) -sin(ra+rb)*theta(i)/rb);sin(ra+rb)*theta(i)/rb) cos(ra+rb)*theta(i)/rb)*xb;yb; xbb=AA(1,:); ybb=AA(2,:); xbb=xbb+(ra+rb)*cos(theta(i); ybb=ybb+(ra+rb)

9、*sin(theta(i); plot(xbb,ybb,b,linewidth,2); plot(xbb(1),ybb(1),o,linewidth,2,markeredgecolor,k,markerfacecolor,g); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 8 页 - - - - - - - - - plot(xbb(1) xbb(m-1)/2),ybb(1) ybb(m-1)/2),k,xbb(m-1)/4) xbb(3*(m-1)/4),ybb(m

10、-1)/4) ybb(3*(m-1)/4),k) axis(-1.2*(ra+rb),1.2*(ra+rb),-1.2*(ra+rb),1.2*(ra+rb); axis equal axis off F=getframe(gcf); mov=addframe(mov,F); hold off end mov = close(mov); 【3】内圆摆线function hypocycloid(ra,rb) fig=figure; set(fig,Color,1,1,1); theta=0:pi/100:2*pi; xa=ra*cos(theta); % 定圆 x 坐标ya=ra*sin(the

11、ta); % 定圆 y 坐标xb=rb*cos(theta); % 动圆 x 坐标yb=rb*sin(theta); % 动圆 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(hypocycloid.avi); for i=1:m plot(xax,0,0,k,0 0,yax,

12、k); hold on plot(xa,ya,k,linewidth,2); plot(xc(1,1:i),yc(1,1:i),r,linewidth,2); AA=cos(ra-rb)*theta(i)/rb) sin(ra-rb)*theta(i)/rb);-sin(ra-rb)*theta(i)/rb) cos(ra-rb)*theta(i)/rb)*xb;yb; xbb=AA(1,:); ybb=AA(2,:); xbb=xbb+(ra-rb)*cos(theta(i); ybb=ybb+(ra-rb)*sin(theta(i); plot(xbb,ybb,b,linewidth,2)

13、; plot(xbb(1),ybb(1),o,linewidth,2,markeredgecolor,k,markerfacecolor,g); plot(xbb(1) xbb(m-1)/2),ybb(1) ybb(m-1)/2),k,xbb(m-1)/4) xbb(3*(m-1)/4),ybb(m-1)/4) ybb(3*(m-1)/4),k) axis equal axis off 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 8 页 - - - - - - - -

14、 - F=getframe(gcf); mov=addframe(mov,F); hold off end mov = close(mov); 四、程序运行结果：1、直线摆线令 r=2，n=5 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 8 页 - - - - - - - - - 2、外圆摆线令 ra=5,rb=1 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - -

15、 第 7 页，共 8 页 - - - - - - - - - 3、内圆摆线令 ra=8,rb=2 五、结果分析：1. 仿真出来的摆线曲线与理论曲线一样，证明计算理论值与仿真值是非常匹配的，误差很少。2仿真过程中 r 和 rb 分别表示转动的圆的半径，其数值越小，所得仿真图样的“花瓣”越多。 Ra是大圆半径，其数值的大小，决定了动圆的运动轨迹的弯曲程度。3. 此次试验最大的乐趣在于，与以前matlab 仿真实验相比较，系统会自动生成动画文件，可以打开播放，直接观察仿真图样的生成过程。 4.再一次切身的感受到matlab 功能的强大，使用之方便。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 8 页 - - - - - - - - -