MATLAB选修论文一.docxVIP

Matlab课程设计 2014111079 刘青青 测控二班 一、轮对磁爬式检测试验模型仿真 如图所示,为测量系统的示意图,它由两个能相互转动的连杆,角度编码器和滚轮等组成。O1为固定点,O2点为转动点,O3点为滚轮的中心,连杆的有效长度分别为L1和L2。任一位置时,连杆1相对于某基准位置的角度为θ1,两连杆的相对角度为θ2。其中对于θ1，取垂直方向为基准线，在基准线左侧的为正，基准线右侧的为负；对于θ2，以连杆一为基准线，在基准线左侧的为正，在基准线右侧的为负。这样可以计算求出O3点的坐标值在踏面直角坐标系(XOY)中的位置。当连杆1连续朝一个方向旋转时,可以获得O3点的一系列的坐标点所构成的轨迹坐标值(x3,y3) x3=x1-L1sinθ1-L2sin(θ1+θ2)；y3=y1-L1cosθ1-L2cos(θ1+θ2) (1) 当滚轮与踏面紧密贴靠时,滚轮的包络线就是踏面外形的实际轮廓线。滚轮的包络线在踏面坐标系中的坐标值可以用下式求得 x′3=x3-rsinθ；y′3=y3-rcosθ (2) 式中,θ为滚轮轨迹拟合曲线在某一点法线与垂直轴之间的夹角,可用下式求得： tanθ=dy3/dx3 (3) 式(1)中,x1、L1和L2为已知值,θ1和θ2为测量所得值。 其中L1=65mm;L2=75mm,r=8mm. x1=0;y1=0; 说明：angle.mat是角度数据，共7980个数据，其中前面一半（3990）为θ1，后面一半为θ2，且一一对应 课程设计报告要求： （1）用图形画出踏面外形的实际轮廓线。 （2）在实际轮廓线上，距离最高点O（x方向54mm距离处A）两点垂直高度差； （3）在实际轮廓线上，求比A点高10mm曲线上的B点与O点的x方向距离。 备注：*图形中横坐标方向为X轴 ，纵坐标方向为Y轴 程序代码 clear; ang=load(F:\matlab课件\MATLAB选修课\2015课程设计\仿真1附件\angle.mat); %将已知数据赋给变量ang theta1=ang.angle(1:3990); %前面一半数据（3990）为θ1 theta2=ang.angle(3991:7980); %后面一半数据为θ2 L1=60; L2=72;r=8.5; %题目所给原始数据 x1=0;y1=0; %题目所给原始数据 x3=x1-L1*sin(theta1*pi/180)-L2*sin((theta1+theta2)*pi/180); y3=y1-L1*cos(theta1*pi/180)-L2*cos((theta1+theta2)*pi/180); dy=diff(y3);%求偏导数 dx=diff(x3); p=atan(dy );%求角度 q=p*pi/180;%角度转化为弧度公式 x5=x3(2:3990);%给定矩阵维数 y5=y3(2:3990);%矩阵维数相同 x4=x5-r*sin(q);%题目所给计算公式 y4=y5-r*sin(q); plot(x4,y4);grid on; %绘制踏面外形实际轮廓线 title(MATLAB模拟踏面外形的实际轮廓线); [max_y,i]=max(y4);%求y4中的最大值max_y和位置i max_x=x4(i);%求i位置对应的x x_A=max_x+54;%距离y3最大值位置处54mm处横坐标 x6=x4(2000:3989); y6=y4(2000:3989); y_A=interp1(x6,y6,x_A); h_OA=max_y-y_A%求AO距离 实验结果： h_OA = 26.5663 d_BO = 20.9663