机械原理课程设计matlab凸轮代码
机械原理凸轮编程
s=351−10×3Ψ−ππ3+15×3Ψ−ππ4−6×3Ψ−ππ5v=−35ω1π330×3Ψ−ππ2−60×3Ψ−ππ3+30×3Ψ−ππ4a=−35ω12π2960×3Ψ−ππ−180×3Ψ−ππ2+120×3Ψ−ππ3fai0=80/180*pi;fai01=60/180*pi;fais=100/180*pi;fais1=120/180*pi;h=35;t1=0.:0.01:fai0;s1=h/2*(1-cos(pi*t1./fai0));t2=fai0:0.01:fai0+fais;s2=h;t3=fai0+fais:0.01:fai0+fais+fai01;s3=h*(1-10*((t3-fai0-fais)/fai01).^3+15*((t3-fai0-fais)/fai01).^4-6*( (t3-fai0-fais)/fai01).^5);t4=fai0+fais+fai01:0.01:2*pi;s4=0;t1=t1/pi*180;t2=t2/pi*180;t3=t3/pi*180;t4=t4/pi*180;subplot(1,2,1);plot(t1,s1);hold on;plot(t2,s2);hold on;plot(t3,s3);hold on;plot(t4,s4);grid on;%从动件位移，速度，加速度与凸轮转角fai的函数关系function [s,v,a]=wyh(fai)h=35; %行程fai0=80/180*pi; %升程运动角fai01=60/180*pi; %回程运动角fais=100/180*pi; %远休止角fais1=120/180*pi; %近休止角w=1; %令凸轮转动角速度为1rad/sif fai>=0 && fai<=fai0s=h/2*(1-cos(pi*t1/fai0));v=pi*h*w/(2*fai0)*sin(pi*fai/fai0));a=pi^2*h*w^2/(2*fai0^2)*cos(pi*fai/fai0);elseif fai>fai0 && fai<fai0+faiss=h;v=0;a=0;elseif fai>=fai0+fais && fai<=fai0+fais+fai01s=h*(1-10*((t3-fai0-fais)/fai01).^3+15*((t3-fai0-fais)/fai01).^4-6*((t3-fai0-fais)/fai01).^5);v=-h*w/fai01(30*((t3-fai0-fais)/fai01).2-60*((t3-fai0-fais)/fai01).^3+30*((t3-fai0-fais)/fai01).^4);a=-h*w^2/fai01(60*((t3-fai0-fais)/fai01)-180*((t3-fai0-fais)/fai01).^2+120*((t3-fai0-fais) /fai01).^3);elseif fai > fai0+fais+fai01 && fai <=fai0+fai01+fais+fais1s=0;v=0;a=0;endend%作从动件位移，速度，加速度曲线s=zeros(1,36001);v=zeros(1,36001);a=zeros(1,36001);for i=1:36001fai=(i-1)/36000*2*pi;[s(i),v(i),a(i)]=sva(fai);endfai=0:0.01:360;subplot(1,3,1);plot(fai,s);xlabel('凸轮转角φ(°)');ylabel('从动件位移s(mm)');title('位移曲线');grid on;subplot(1,3,2);plot(fai,v);xlabel('凸轮转角φ/°');ylabel('从动件速度v(mm/s)');title('速度曲线');grid on;subplot(1,3,3);plot(fai,a);xlabel('凸轮转角φ/°');ylabel('从动件加速度a(mm/s^2)'); title('加速度曲线');grid on;。
基于matlab的GUI设计——机械原理教学演示系统
机械原理教学演示系统——基于matlab的GUI设计xxx指导老师：20年月日目录一、功能简介 (3)二、总界面 (3)三、凸轮模块 (4)四、齿轮模块 (6)五．参考书目 (6)六．附录（部分程序源代码） (7)一、功能简介本系统能实现机械原理教学过程中凸轮模块与齿轮模块的设计与运动仿真，加深对机械原理课程学习的理解。
二、总界面总界面标题设置：set(gcf,'name','机械原理教学演示系统 made by 翟鲁鑫');背景图片设置：ximg=imread('机械原理课本.jpg');imshow(ximg);背景声音播放：Fs=44100;[ywav,Fs]=wavread('夜的钢琴曲 - 六3.wav');sound(ywav,Fs);到各个模块：到凸轮模块clcclose(gcf);clear allGUItulun到齿轮模块clcclose(gcf);clear allGUIchilun关闭系统：clcquestion='真的要退出吗？';title='确认退出？';button=questdlg(question,title,'是','否','是'); switch buttoncase'是'clear allclosecase'否'returnend三、凸轮模块设计要点：1.背景声音设置方法同主界面2.推程角、远休角、回程角之和不能大于360度的判别条件；基圆半径、滚子半径、行程不能为0的判别条件sr0=get(handles.edit2,'string');r0=str2num(sr0);if isequal(r0,0)errordlg('基圆半径不能为0,请重新输入','出错');returnendsrr=get(handles.edit3,'string');rr=str2num(srr);if isequal(rr,0)errordlg('滚子半径不能为0,请重新输入','出错');returnendsh=get(handles.edit4,'string');h=str2num(sh);if isequal(h,0)errordlg('行程不能为0,请重新输入','出错');returnendn3=phi01+phi02+phi03; %推程角、远休止角与回程角的总和if n3>360errordlg('角度之和大于360，请重新输入','出错');end3.仿真程序。
MATLAB机械分析
图2 凸轮理论轮廓曲线
利用SolidWorks三维建模软件建立凸轮的三维模型，如图3所示。
图3凸轮三维模型
1.4
在凸轮设计时，经常需要设计一系列凸轮，这些凸轮之间有相位差。假设第二个凸轮与第一个凸轮之间的相位角为60°。利用MATLAB编制程序，将凸轮轮廓线上的每一个点坐标取出，乘以转换矩阵T，就得到凸轮二的轮廓线的各店的坐标。最后绘制凸轮二的理论轮廓线如图4。
图3关节角θ1和关节角θ2之间的角度约束关系
根据机械手的各连杆参数，建立连杆之间的D-H变换矩阵：
将表2中的数据代入可以得到各连杆坐标系之间的变换矩阵为
假设坐标系3中一点P，其坐标描述为 。该点在坐标系0中描述为 ，则有
取连杆3的关节点C为研究对象，点C在坐标系3中的坐标为 ，利用转换矩阵 ，可以找出点C在坐标系0（即固定坐标系）中的坐标 。则机械手工作头中心D在坐标系0中的坐标为 。
lunkuo(i+1,1)=h*cos(theta);
lunkuo(i+1,2)=h*sin(theta);
lunkuo(i+1,3)=0;
end
end
save 'F:\lunkuo.txt' lunkuo -ascii;
机械原理大作业—凸轮程序
附录：VB程序Option ExplicitDim s As DoubleDim v As DoubleDim a As DoubleDim h As DoubleDim f As DoubleDim f0 As DoubleDim f01 As DoubleDim fs As DoubleDim fs1 As DoubleDim w As DoubleDim fj1 As DoubleDim pi As DoubleDim pa As DoubleDim i As DoubleDim dsdf As Double Dim s0 As DoubleDim e As DoubleDim r0 As DoubleDim dydf As Double Dim dxdf As Double Dim x0 As DoubleDim y0 As DoubleDim x As DoubleDim y As DoubleDim x1 As DoubleDim y1 As DoubleDim x2 As DoubleDim y2 As DoubleDim dsdf2 As Double Dim fp As DoubleDim sp As DoubleDim dsdfp As Double Private Sub Form_Load() pi = 3.1415926pa = pi / 180h = 130f0 = 2 / 3 * pif01 = pi / 2fs = pi / 6fs1 = 2 / 3 * piw = 10End SubPrivate Sub Command1_Click()Picture1.Scale (-20, 200)-(380, -10)Picture1.Line (-20, 0)-(380, 0) 'XPicture1.Line (0, 200)-(0, -10) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture1.DrawStyle = 2Picture1.Line (i, 200)-(i, -10)Picture1.CurrentX = i - 10: Picture1.CurrentY = 0 Picture1.Print iNext iFor i = -200 To 500 Step 50 'Y轴坐标Picture1.DrawStyle = 2Picture1.Line (-20, i)-(380, i)Picture1.CurrentX = -20: Picture1.CurrentY = i + 7 Picture1.Print iNext iFor fj1 = 0 To 360 Step 0.1f = fj1 * paCall mPicture1.PSet (fj1, s)Next fj1End SubPrivate Sub Command2_Click()Picture2.Scale (-20, 1200)-(380, -2000)Picture2.Line (-20, 0)-(380, 0) 'XPicture2.Line (0, 1200)-(0, -2000) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture2.DrawStyle = 2Picture2.Line (i, 1200)-(i, -2000)Picture2.CurrentX = i - 10: Picture2.CurrentY = 0 Picture2.Print iNext iFor i = -2000 To 1200 Step 200 'Y轴坐标Picture2.Line (0, i)-(380, i)Picture2.CurrentX = -25: Picture2.CurrentY = i Picture2.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall mPicture2.PSet (fj1, v)Next fj1End SubPrivate Sub Command3_Click()Picture3.Scale (-40, 30000)-(380, -30000)Picture3.Line (-40, 0)-(380, 0) 'XPicture3.Line (0, 30000)-(0, -30000) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture3.DrawStyle = 2Picture3.Line (i, 30000)-(i, -30000)Picture3.CurrentX = i - 10: Picture3.CurrentY = 0 Picture3.Print iNext iFor i = -30000 To 30000 Step 5000 'Y轴坐标Picture3.Line (0, i)-(380, i)Picture3.CurrentX = -40: Picture3.CurrentY = i Picture3.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall mPicture3.PSet (fj1, a)Next fj1End SubPrivate Sub Command4_Click()Picture4.Scale (-200, 200)-(160, -160)Picture4.Line (-200, 0)-(160, 0) 'XPicture4.Line (0, 200)-(0, -160) 'YFor i = -200 To 160 Step 20 'X轴坐标Picture4.DrawStyle = 2Picture4.Line (i, 200)-(i, -160)Picture4.CurrentX = i - 10: Picture4.CurrentY = 0 Picture4.Print iNext iFor i = -160 To 200 Step 20 'Y轴坐标Picture4.Line (-200, i)-(160, i)Picture4.CurrentX = -15: Picture4.CurrentY = i + 5 Picture4.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall mPicture4.PSet (dsdf, s)Next fj1End SubPrivate Sub Command5_Click() '求点E的速度Picture5.Scale (-250, 150)-(150, -250)Picture5.Line (-250, 0)-(150, 0) 'XPicture5.Line (0, 150)-(0, -250) 'YFor i = -250 To 150 Step 50 'X轴坐标Picture5.DrawStyle = 2Picture5.Line (i, 150)-(i, -250)Picture5.CurrentX = i - 10: Picture5.CurrentY = 0 Picture5.Print iNext iFor i = -250 To 150 Step 50 'Y轴坐标Picture5.Line (-250, i)-(150, i)Picture5.CurrentX = -30: Picture5.CurrentY = i + 25 Picture5.Print iNextFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall mPicture5.PSet (-y, x)Picture5.PSet (-y0, x0)Picture5.PSet (x1, y1)Next fj1End SubPrivate Sub m()If f >= 0 And f <= 2 * pi / 3 Thens = 65 * (1 - Cos(3 / 2 * f))v = 97.5 * w * Sin(3 / 2 * f)a = 146.25 * w ^ 2 * Cos(1.5 * f)dsdf = 97.5 * Sin(3 / 2 * f)dsdf2 = 97.5 * 3 / 2 * Cos(3 / 2 * f)ElseEnd IfIf f >= 2 * pi / 3 And f <= 5 / 6 * pi Thens = 130v = 0a = 0dsdf = 0dsdf2 = 0ElseEnd IfIf f >= 5 / 6 * pi And f <= 13 / 12 * pi Thens = 130 - (1040 / (pi ^ 2)) * (f - 5 * pi / 6) ^ 2v = -2080 * w / (pi ^ 2) * (f - 5 * pi / 6)a = -2080 * w ^ 2 / (pi ^ 2)dsdf = -2080 / (pi ^ 2) * (f - 5 * pi / 6)dsdf2 = -2080 / (pi ^ 2)ElseEnd IfIf f >= 13 / 12 * pi And f <= 4 / 3 * pi Thens = 1040 / (pi ^ 2) * (pi * 4 / 3 - f) ^ 2v = -2080 * w / (pi ^ 2) * (pi * 4 / 3 - f)a = 2080 * w ^ 2 / (pi ^ 2)dsdf = -2080 / (pi ^ 2) * (pi * 4 / 3 - f)dsdf2 = 2080 / (pi ^ 2)ElseEnd IfIf f >= 4 / 3 * pi And f <= 2 * pi Thens = 0v = 0a = 0dsdf = 0dsdf2 = 0ElseEnd Iffp = 2 / 3 * Atn(3 / (2 * Tan(35 * pa)))sp = 65 * (1 - Cos(3 / 2 * fp))dsdfp = 97.5 * Sin(3 / 2 * fp)r0 = 100e = 0s0 = Sqr(r0 ^ 2 - e ^ 2)x = (s0 + s) * Cos(f) - e * Sin(f)y = (s0 + s) * Sin(f) + e * Cos(f)dxdf = -(s0 + s) * Sin(f) + dsdf * Cos(f) - e * Cos(f) dydf = (s0 + s) * Cos(f) + dsdf * Sin(f) - e * Sin(f) x0 = x - 20 * dydf / Sqr(dxdf ^ 2 + dydf ^ 2)y0 = y + 20 * dxdf / Sqr(dxdf ^ 2 + dydf ^ 2)x1 = r0 * Cos(f)y1 = r0 * Sin(f)End Sub。
matlab凸轮轮廓设计及仿真说明书.
偏置盘型凸轮创新课程设计课程名称：机械原理设计题目：偏置盘型凸轮设计院系：机电学院班级：09机41设计者：彭辉学号：09294040指导教师：王卫辰学校：江苏师范大学前言凸轮轮廓曲线的设计，一般可分为图解法和解析法．利用图解法能比较方便地绘制出各种平面凸轮的轮廓曲线．但这种方法仅适用于比较简单的结构，用它对复杂结构进行设计则比较困难，而且利用图解法进行结构设计，作图误差较大，对一些精度要求高的结构不能满足设计要求．解析法可以根据设计要求，通过推导机构中各部分之间的几何关系，建立相应的方程，精确地计算出轮廓线上各点的坐标，然后把凸轮的轮廓曲线精确地绘制出来．但是，当从动件运动规律比较复杂时，利用解析法获得凸轮的轮廓曲线的工作量比较大．而MATLAB软件提供了强大的矩阵处理和绘图功能，具有核心函数和工具箱．其编程代码接近数学推导公式，简洁直观，操作简易，人机交互性能好，且可以方便迅速地用三维图形、图像、声音、动画等表达计算结果、拓展思路口。
因此，基于MATLAB软件进行凸轮机构的解析法设计，可以解决设计工作量大的问题。
本此课程设计基于MATLAB软件进行凸轮轮廓曲线的解析法设计，并对的运动规律凸轮进行仿真，其具体方法为首先精确地计算出轮廓线上各点的坐标，然后运用MATLAB绘制比较精确的凸轮轮廓曲线和推杆的位移、速度及加速度曲线以及仿真。
目录前言 1第一章：工作意义 3 1.1本次课程设计意义3 1.2 已知条件4第二章：工作设计过程 5 2.1：设计思路 5 2.2：滚子从动件各个阶段相关方程 6 2.3：盘型凸轮理论与实际轮廓方程 7第三章：工作程序过程 7 3.1：滚子从动件各各阶段MATLAB程序编制 8 3.2：凸轮的理论实际运动仿真程序编制 12 第四章：运行结果 17 4.1：滚子运动的位移图 17 4.2：滚子运动的速度图 17 4.3：滚子运动的加速度图，局部加速度图 18 4.4：滚子运动的仿真图 19 4.5：滚子运动的理论与实际轮廓图 20第五章：设计总结 21 5.1：总结 21第六章：参考文献 226.1:参考文献 22第一章：工作意义1.1本次课程设计意义凸轮是一个具有曲线轮廓或凹槽的构件，一般为主动件，作等速回转运动或往复直线运动。
基于MATLAB软件的凸轮轮廓曲线设计_
基于ＭＡＴＬＡＢ软件的凸轮轮廓曲线设计摘要：以偏置移动从动件盘形凸轮为例，基于ＭＡＴＬＡＢ软件对凸轮轮廓曲线进展了解析法设计．绘制出轮廓曲线。
运行结果说明：在从动件运动规律确定的情况下，利用ＭＡＴＬＡＢ软件以很方便、快捷地得到凸轮的轮廓曲线。
关键词：凸轮机构；凸轮轮廓曲线；ＭＡＴＬＡＢ；解析法前言凸轮轮廓曲线的设计，一般可分为图解法和解析法．利用图解法能比拟方便地绘制出各种平面凸轮的轮廓曲线．但这种方法仅适用于比拟简单的构造，用它对复杂构造进展设计那么比拟困难，而且利用图解法进展构造设计，作图误差较大，对一些精度要求高的构造不能满足设计要求。
解析法可以根据设计要求，通过推导机构中各局部之间的几何关系，建立相应的方程，准确地计算出轮廓线上各点的坐标，然后把凸轮的轮廓曲线准确地绘制出来．但是，当从动件运动规律比拟复杂时，利用解析法获得凸轮的轮廓曲线的工作量比拟大．而ＭＡＴＬＡＢ软件提供了强大的矩阵处理和绘图功能，具有核心函数和工具箱．其编程代码接近数学推导公式，简洁直观，操作简易，人机交互性能好，且可以方便迅速地用三维图形、图像、声音、动画等表达计算结果、拓展思路［１］。
因此，基于ＭＡＴＬＡＢ软件进展凸轮机构的解析法设计，可以解决设计工作量大的问题。
本文基于ＭＡＴＬＡＢ软件进展凸轮轮廓曲线的解析法设计，利用?机械原理?课程的计算机辅助教学，及常用机构的计算机辅助设计．其具体方法为首先准确地计算出轮廓线上各点的坐标，然后运用ＭＡＴＬＡＢ绘制比拟准确的凸轮轮廓曲线。
1 设计的意义与条件1.1意义凸轮机构是由具有曲线轮廓或凹槽的构件，通过高副接触带动从动件实现预期运动规律的一种高副机构，它广泛地应用于各种机械，特别是自动机械、自动控制装置和装配生产线中，是工程实际中用于实现机械化和自动化的一种常用机构。
所以，在凸轮的加工中，准确确实定凸轮的轮廓，这对于保证凸轮所带动从动件的运动规律是尤为重要的。
1.2条件偏置移动从动件盘形凸轮设计条件〔图1〕:凸轮作逆时针方向转动,从动件偏置在凸轮轴心的右边从动件在推程作等加速/等减速运动,在回程作余弦加速度运动基圆半径rb = 40 mm，滚子半径rt = 10mm，推杆偏距e = 15 mm，推程升程h = 50 mm，推程运动角ft = 100度，远休止角fs = 60度回程运动角fh = 90度，推程许用压力角alp = 35度。
基于matlab的凸轮计算
基于matlab的凸轮计算>> r=10;>> r0=50;>> e=0;>> delt0=120*pi/180;>> h=50;>> i=100;>> s0=sqrt(r0^2-e^2);>> syms delt1 delt2 delt3;>> s1=h*((delt1/delt0)-sin(2*pi*delt1/delt0)/(2*pi)); >> v1=diff(s1);>> a1=diff(s1,2);>> s2=h*(2+sin(2*pi*delt2/delt0)/(2*pi)-(delt2/delt0)); >> v2=diff(s2);>> a2=diff(s2,2);>> delt1=linspace(0,delt0,i);>> delt2=linspace(delt0,240*pi/180,i);>> delt3=linspace(240*pi/180,2*pi,i);>> s3=0*delt3;>> delt=[delt1 delt2 delt3];>> s11=subs(s1,delt1);>> v11=subs(v1,delt1);>> a11=subs(a1,delt1);>> s22=subs(s2,delt2);>> v22=subs(v2,delt2);>> a22=subs(a2,delt2);>> s=[s11 s22 s3];>> v=[v11 v22 s3];>> a=[a11 a22 s3];>> plot(delt*180/pi,s,delt*180*pi,v,'-.',delt*180/pi,a,'--');>> title('凸轮位移,速度,加速度曲线');>> legend('位移曲线','速度曲线','加速度曲线');>> axis([0 360 -80 80]);>> grid on>> for j=1:3*ixx(j)=(s0+s(j))*sin(delt(j))-e*cos(delt(j));yy(j)=(s0+s(j))*cos(delt(j))+e*sin(delt(j));end>> for m=1:3*isyms deltxyx=(s0+s(m))*sin(deltxy)+e*cos(deltxy);y=(s0+s(m))*cos(deltxy)-e*sin(deltxy);sx=diff(x,deltxy)/(sqrt((diff(x,deltxy))^2+(diff(y,deltxy))^2));cx=diff(y,deltxy)/(sqrt((diff(x,deltxy))^2+(diff(y,deltxy))^2));deltxy=delt(m);ax(m)=subs(sx,deltxy);bx(m)=subs(cx,deltxy);xxx(m)=xx(m)-r*bx(m);yyy(m)=yy(m)-r*ax(m);endfigure>> plot(xxx,yyy)>> grid on>> title('凸轮实际轮廓曲线');>> cave=[xxx;yyy;0*xxx-5];>> fprintf('%10.6f',cave)0.000000 40.000000 -5.000000 1.269243 39.991388 -5.000000 2.538018 39.966909 -5.000000 3.806015 39.928587 -5.000000 5.073092 39.878420 -5.000000 6.339275 39.818364 -5.000000 7.604751 39.750315 -5.000000 8.869863 39.676091 -5.000000 10.135100 39.597420 -5.000000 11.401095 39.515925 -5.000000 12.668608 39.433108 -5.000000 13.938519 39.350338 -5.000000 15.211816 39.268838 -5.000000 16.489581 39.189678 -5.000000 17.772978 39.113760 -5.000000 19.063237 39.041812 -5.000000 20.361640 38.974382 -5.000000 21.669506 38.911829 -5.000000 22.988173 38.854322 -5.000000 24.318985 38.801835 -5.000000 25.663272 38.754145 -5.000000 27.022336 38.710835 -5.000000 28.397436 38.671291 -5.000000 29.789767 38.634708 -5.000000 31.200450 38.600093 -5.000000 32.630511 38.566270 -5.000000 34.080872 38.531891 -5.000000 35.552332 38.495438 -5.000000 37.045557 38.455237 -5.000000 38.561065 38.409469 -5.000000 40.099216 38.356180 -5.000000 41.660201 38.293294 -5.000000 43.244034 38.218627 -5.000000 44.850541 38.129905 -5.000000 46.479357 38.024773 -5.000000 48.129916 37.900818 -5.000000 49.801452 37.755579 -5.000000 51.492993 37.586571 -5.000000 53.203358 37.391297 -5.000000 54.931164 37.167264 -5.000000 56.674820 36.912008 -5.000000 58.432535 36.623105 -5.000000 60.202321 36.298187 -5.000000 61.981998 35.934966 -5.000000 63.769205 35.531242 -5.000000 65.561403 35.084926 -5.000000 67.355891 34.594048 -5.000000 69.149813 34.056777 -5.000000 70.940171 33.471432 -5.000000 72.723841 32.836494 -5.000000 74.497584 32.150618 -5.000000 76.258063 31.412642 -5.000000 78.001860 30.621598 -5.000000 79.725490 29.776719 -5.000000 81.425420 28.877443 -5.000000 83.098087 27.923421 -5.000000 84.739914 26.914517 -5.000000 86.347331 25.850814 -5.000000 87.916791 24.732613 -5.000000 89.444787 23.560431 -5.000000 90.927876 22.334999 -5.000000 92.362689 21.057259 -5.000000 93.745954 19.728360 -5.000000 95.074508 18.349647 -5.000000 96.34531916.922660 -5.000000 97.555495 15.449119 -5.000000 98.702302 13.930915 -5.000000 99.783179 12.370102 -5.000000100.795745 10.768880 -5.000000101.737815 9.129583 -5.000000102.607408 7.454669 -5.000000103.402758 5.746697 -5.000000104.122316 4.008319 -5.000000104.764763 2.242259 -5.000000105.329010 0.451300 -5.000000105.814201 -1.361736 -5.000000106.219718 -3.194005 -5.000000106.545177 -5.042653 -5.000000106.790431 -6.904840 -5.000000106.955562 -8.777752 -5.000000107.040882-10.658618 -5.000000107.046925-12.544727 -5.000000106.974436-14.433446 -5.000000106.824372-16.322227 -5.000000106.597882-18.208630 -5.000000106.296303-20.090331 -5.000000105.921147-21.965134 -5.000000105.474084-23.830985 -5.000000104.956934-25.685979 -5.000000104.371648-27.528371 -5.000000103.720293-29.356583 -5.000000103.005038-31.169209 -5.000000102.228137-32.965021 -5.000000101.391911-34.742971 -5.000000100.498728-36.502197 -5.000000 99.550994-38.242018 -5.000000 98.551128-39.961937 -5.000000 97.501545-41.661638 -5.000000 96.404646-43.340982 -5.000000 95.262794-45.000000 -5.000000 95.262794-45.000000 -5.000000 94.077722-46.638540 -5.000000 92.848837-48.255106 -5.000000 91.575076-49.847744 -5.000000 90.255505-51.414382 -5.000000 88.889338-52.952839 -5.000000 87.475945-54.460840 -5.000000 86.014863-55.936028 -5.000000 84.505807-57.375980 -5.000000 82.948676-58.778223 -5.000000 81.343563-60.140244 -5.000000 79.690755-61.459517 -5.000000 77.990742-62.733508 -5.000000 76.244218-63.959702 -5.000000 74.452080-65.135612 -5.000000 72.615433-66.258803 -5.00000070.735581-67.326904 -5.000000 68.814030-68.337626 -5.000000 66.852480-69.288782 -5.000000 64.852821-70.178297 -5.000000 62.817124-71.004227 -5.000000 60.747632-71.764774 -5.000000 58.646752-72.458299 -5.000000 56.517044-73.083332 -5.000000 54.361205-73.638590 -5.000000 52.182062-74.122982 -5.000000 49.982551-74.535624 -5.000000 47.765709-74.875841 -5.000000 45.534653-75.143179 -5.000000 43.292567-75.337411 -5.000000 41.042683-75.458537 -5.000000 38.788266-75.506791 -5.000000 36.532596-75.482641 -5.000000 34.278950-75.386788 -5.000000 32.030586-75.220167 -5.000000 29.790723-74.983939 -5.000000 27.562529-74.679493 -5.000000 25.349097-74.308435 -5.000000 23.153435-73.872583 -5.000000 20.978450-73.373958 -5.000000 18.826928-72.814773 -5.000000 16.701525-72.197424 -5.000000 14.604753-71.524474 -5.000000 12.538966-70.798644 -5.000000 10.506349-70.022795 -5.000000 8.508912-69.199916 -5.000000 6.548474-68.333103 -5.000000 4.626665-67.425549 -5.000000 2.744911-66.480521 -5.000000 0.904434-65.501344 -5.000000 -0.893750-64.491384 -5.000000 -2.648838-63.454030 -5.000000 -4.360234-62.392674 -5.000000 -6.027553-61.310693 -5.000000 -7.650616-60.211434 -5.000000 -9.229446-59.098191 -5.000000-10.764265-57.974193 -5.000000-12.255484-56.842582 -5.000000-13.703698-55.706402 -5.000000-15.109673-54.568580 -5.000000-16.474341-53.431914 -5.000000-17.798784-52.299057 -5.000000-19.084223-51.172507 -5.000000-20.332006-50.054594 -5.000000-21.543592-48.947468 -5.000000-22.720536-47.853097 -5.000000-23.864474-46.773252 -5.000000-24.977110-45.709505 -5.000000-26.060194-44.663222 -5.000000-27.115508-43.635563 -5.000000-28.144852-42.627476 -5.000000-29.150023-41.639701 -5.000000-30.132800-40.672768 -5.000000-31.094928-39.727000 -5.000000-32.038102-38.802521 -5.000000-32.963951-37.899254 -5.000000-33.874024-37.016935 -5.000000-34.769774-36.155117 -5.000000-35.652545-35.313181 -5.000000-36.523561-34.490347 -5.000000-37.383912-33.685683 -5.000000-38.234547-32.898120 -5.000000-39.076261-32.126466 -5.000000-39.909689-31.369417 -5.000000-40.735297-30.625575 -5.000000-41.553381-29.893465 -5.000000-42.364059-29.171544 -5.000000-43.167267-28.458227 -5.000000-43.962764-27.751896 -5.000000-44.750126-27.050921 -5.000000-45.528751-26.353673 -5.000000-46.297861-25.658546 -5.000000-47.056505-24.963968 -5.000000-47.803568-24.268420 -5.000000-48.537777-23.570450 -5.000000-49.257709-22.868689 -5.000000-49.961799-22.161864 -5.000000-50.648356-21.448813 -5.000000-51.315569-20.728494 -5.000000-51.961524-20.000000 -5.000000-51.961524-20.000000 -5.000000-52.584515-19.262731 -5.000000-53.183972-18.516841 -5.000000-53.759626-17.762665 -5.000000-54.311222-17.000538 -5.000000-54.838511-16.230804 -5.000000-55.341258-15.453805 -5.000000-55.819237-14.669890 -5.000000-56.272235-13.879410 -5.000000-56.700049-13.082719 -5.000000-57.102488-12.280172 -5.000000-57.479371-11.472129 -5.000000-57.830530-10.658953 -5.000000-58.155807 -9.841005 -5.000000-58.455057 -9.018654 -5.000000-58.728147 -8.192267 -5.000000-58.974953 -7.362213 -5.000000-59.195366 -6.528864 -5.000000-59.389287 -5.692594 -5.000000-59.556628 -4.853775 -5.000000-59.697316 -4.012785 -5.000000-59.811287 -3.169998 -5.000000-59.898489 -2.325793 -5.000000-59.958886 -1.480547 -5.000000-59.992448 -0.634639 -5.000000-59.999161 0.211554 -5.000000-59.979022 1.057652 -5.000000-59.932040 1.903277 -5.000000-59.858237 2.748049 -5.000000-59.757644 3.591592 -5.000000-59.630308 4.433528 -5.000000-59.476285 5.273479 -5.000000-59.2956436.111071 -5.000000-59.088465 6.945927 -5.000000-58.8548437.777675 -5.000000-58.594880 8.605942 -5.000000-58.3086949.430357 -5.000000-57.996413 10.250552 -5.000000-57.658176 11.066160 -5.000000-57.294134 11.876815 -5.000000-56.904452 12.682155 -5.000000-56.489302 13.481819 -5.000000-56.048872 14.275449 -5.000000-55.583357 15.062690 -5.000000-55.092966 15.843191 -5.000000-54.577920 16.616601 -5.000000-54.038447 17.382574 -5.000000-53.474791 18.140768 -5.000000-52.887202 18.890843 -5.000000-52.275944 19.632464 -5.000000-51.641291 20.365298 -5.000000-50.983526 21.089019 -5.000000-50.302944 21.803301 -5.000000-49.599849 22.507825 -5.000000-48.874557 23.202276 -5.000000-48.127392 23.886344 -5.000000-47.358687 24.559721 -5.000000-46.568788 25.222107 -5.000000-45.758047 25.873205 -5.000000-44.926828 26.512723 -5.000000-44.075503 27.140376 -5.000000-43.204451 27.755883 -5.000000-42.314065 28.358968 -5.000000-41.404741 28.949362 -5.000000-40.476887 29.526799 -5.000000-39.530917 30.091022 -5.000000-38.567257 30.641778 -5.000000-37.586335 31.178820 -5.000000-36.588593 31.701909 -5.000000-35.574476 32.210810 -5.000000-34.544438 32.705296 -5.000000-33.498939 33.185144 -5.000000-32.438449 33.650141 -5.000000-31.363441 34.100079 -5.000000-30.274397 34.534755 -5.000000-29.171804 34.953975 -5.000000-28.056156 35.357552 -5.000000-26.927951 35.745306 -5.000000-25.787695 36.117062 -5.000000-24.635898 36.472654 -5.000000-23.473075 36.811923 -5.000000-22.299747 37.134717 -5.000000-21.116440 37.440893 -5.000000-19.923681 37.730312 -5.000000-18.722007 38.002845 -5.000000-17.511953 38.258370 -5.000000-16.294062 38.496774 -5.000000-15.068879 38.717948 -5.000000-13.836952 38.921795 -5.000000-12.598833 39.108223 -5.000000-11.355075 39.277148 -5.000000-10.106235 39.428495 -5.000000 -8.852872 39.562197 -5.000000 -7.595547 39.678193 -5.000000 -6.334823 39.776431 -5.000000 -5.071264 39.856868 -5.000000 -3.805435 39.919467 -5.000000 -2.537903 39.964201 -5.000000 -1.269236 39.991049 -5.000000 -0.000000 40.000000 -5.000000>>。
机械原理课程设计Matlab编程
/*Matlab程序*/l1 = 59.1000;l2 = 263.9000;l3=120;l4=266.83;l5=180;l6=45;x2=170;y2=132.7289;w1=9.4248;N=42:10:402;ay=119:10:479a=2*l1*l3*sin(N/180*pi);b=2*l3*(l1*cos(N/180*pi)-l4);c=l2^2-l1^2-l3^2-l4^2+2*l1*l4*cos(N/180*pi);jiao3=2*atan((a- sqrt(a.^2+b.^2-c.^2))./ (b-c))/pi*180+77g=2*l1*l2*sin(N/180*pi);h=2*l2*(l1*cos(N/180*pi)-l4);m=l1^2+l2^2+l4^2-l3^2-2*l1*l4*cos(N/180*pi);jiao2=2*atan((g- sqrt(g.^2+h.^2-m.^2))./ (h-m))/pi*180+77lof=-sqrt(l6^2-x2^2-l5^2+2*l5*x2*cos((180+jiao3)/180*pi)+l5^2*sin(j iao3/180*pi).^2)+y2-l5* sin((180+jiao3)/180*pi)j12=N-(jiao2-77);j32=jiao3-jiao2;j13=(N-(jiao3-77));j23=(jiao2-jiao3);w3=(w1*l1*sin(j12/180*pi))./ (l3*sin(j32/180*pi))w2=(-1*w1*l1*sin(j13/180*pi))./(l2*sin(j23/180*pi))a3=(w1^2*l1*cos(j12/180*pi)+w2.^2*l2-(w3.^2).*(l3*cos(j32/180*pi)) )./ (l3*sin(j32/180*pi))a2=(-w1^2*l1*cos(j13/180*pi)-(w2.^2).*(l2*cos(j23/180*pi))+l3*w3.^ 2)./ (l3*sin(j23/180*pi))jiao4=acos((x2-l5*cos((180+jiao3)/180*pi))/l6)/pi*180w4=((-l5*sin((pi+jiao3)/180*pi)).*w3)./ (l6*sin(jiao4/180*pi))vof=((l5*sin((180+jiao3-jiao4)/180*pi)).*w3).* sin(jiao4/180*pi)aof=(l6*w4.^2+(l5*w3.^2).*(cos((180+jiao3-jiao4)/180*pi))+l5*a3.*si n((180+jiao3-jiao4)/180*pi))./sin(jiao4/180*pi)作图程序：/*F点的位移*/plot(N+77,lof,'-xk')xlabel('AB杆的角度'),ylabel('F点的位移/(mm)')title('F点的位移曲线图')text(100, 171.3339,'初始值= 171.3339')text(425, 171.3339,'最高值= 180')/*F点的速度*/plot(N+77,vof,'-xk')text(100,-285,'初始值= -347.0624')text(390,-430,'终点值= -347.0624')xlabel('AB杆的角度'),ylabel('F点的速度/(mm/s)')title('F点的速度曲线图')/*F点的加速度*/plot(N+77,aof,'-xk')xlabel('AB杆的角度'),ylabel('F点的加速度/(mm/s^2)')title('F点的加速度曲线图')text(100, 1.0e+004 *-0.9848,'初始值= 1.0e+004 *（-0.9048）') text(350, 1.0e+004 *-0.9048,'终点值= 1.0e+004 *（-0.9048）')/*结果*/ay =Columns 1 through 20119 129 139 149 159 169 179 189 199 209 219 229 239 249 259 269 279 289 299 309Columns 21 through 37319 329 339 349 359 369 379 389 399 409 419 429 439 449 459 469 479jiao3 =Columns 1 through 12152.7283 155.6055 159.2201 163.3851 167.9364 172.7308 177.6411 182.5508 187.3495 191.9295 196.1844 200.0109Columns 13 through 24203.3123 206.0053 208.0272 209.3409 209.9360 209.8245 209.0330 207.5944 205.5410 202.9003 199.6946 195.9435Columns 25 through 36191.6697 186.9095 181.7262 176.2288 170.5899 165.0593 159.9573 155.6378 152.4212 150.5209 150.0024 150.7929Column 37152.7283jiao2 =Columns 1 through 1297.9329 96.9493 96.1177 95.4427 94.9235 94.5564 94.3364 94.2587 94.3188 94.5127 94.8367 95.2869Columns 13 through 2495.8585 96.5454 97.3393 98.2294 99.2013 100.2373 101.3154 102.4095 103.4889 104.5187 105.4597 106.2696Columns 25 through 36106.9048 107.3227 107.4861 107.3682 106.9579 106.2651 105.3223 104.1834 102.9167 101.5959 100.2898 99.0551Column 3797.9329lof =Columns 1 through 12171.3339 162.4831 151.6215 139.2662 125.7535 111.3249 96.2283 80.8029 65.5172 50.9387 37.6489 26.1414Columns 13 through 2416.7408 9.5650 4.5415 1.4741 0.1412 0.3882 2.1781 5.5874 10.7654 17.8808 27.0707 38.3902Columns 25 through 3651.7607 66.9223 83.4132 100.6062 117.8015 134.3068 149.4275 162.3850 172.2940 178.3195 179.9923 177.4473Column 37171.3339w3 =Columns 1 through 122.0420 2.88413.54184.0289 4.3651 4.5686 4.6533 4.6279 4.4965 4.2598 3.9173 3.4690Columns 13 through 242.9187 2.2766 1.5612 0.7980 0.0174 -0.7506 -1.4803 -2.1528 -2.7570 -3.2875 -3.7434 -4.1259Columns 25 through 36-4.4354 -4.6697 -4.8213 -4.8744 -4.8021 -4.5650 -4.1167 -3.4196 -2.4746 -1.3428 -0.1370 1.0209Column 372.0420w2 =Columns 1 through 12-1.4343 -1.1068 -0.8180 -0.5658 -0.3428 -0.1403 0.0504 0.2372 0.4268 0.6248 0.8348 1.0577Columns 13 through 241.2907 1.5265 1.7527 1.95262.1078 2.20042.2158 2.1438 1.9799 1.7238 1.3800 0.9564Columns 25 through 360.4649 -0.0783 -0.6512 -1.2249 -1.7619 -2.2168 -2.5411 -2.6940 -2.6601 -2.4610 -2.1504 -1.7915Column 37-1.4343a3 =Columns 1 through 1252.5560 41.5031 31.2310 22.2213 14.4470 7.6608 1.5603 -4.1325 -9.6312 -15.0579 -20.4136 -25.5519Columns 13 through 24-30.1734 -33.8717 -36.2458 -37.0535 -36.3289 -34.3895 -31.7205 -28.8042 -25.9756 -23.3519 -20.8311 -18.1297Columns 25 through 36-14.8311 -10.4321 -4.3902 3.8075 14.5130 27.6704 42.4330 56.8069 67.8232 72.7326 70.6513 63.0454Column 3752.5560a2 =Columns 1 through 1242.6168 37.6998 32.7111 28.4933 25.3433 23.2851 22.2337 22.0665 22.6383 23.7616 25.1694 26.4765Columns 13 through 2427.1651 26.6283 24.2896 19.7685 13.0154 4.3334 -5.7233 -16.5098 -27.4206 -37.9407 -47.6216 -56.0112Columns 25 through 36-62.5725 -66.6194 -67.2958 -63.6367 -54.7645 -40.2872 -20.8999 1.0456 21.6515 36.8894 44.6857 45.7784Column 3742.6168jiao4 =Columns 1 through 1277.1498 82.2481 87.8236 93.1653 97.6943 100.9570 102.6405 102.6068 100.9154 97.8068 93.6522 88.8959Columns 13 through 2484.0113 79.4692 75.7071 73.0882 71.8537 72.0875 73.7150 76.5391 80.2895 84.6612 89.3256 93.9195Columns 25 through 3698.0213 101.1380 102.7329 102.3311 99.6946 94.9908 88.8525 82.3019 76.5687 72.8043 71.7241 73.3615Column 3777.1498w4 =Columns 1 through 12-3.4249 -4.2204 -4.2959 -3.7606 -2.7325 -1.3398 0.2606 1.8815 3.3352 4.4720 5.1962 5.4568Columns 13 through 245.2293 4.5113 3.3352 1.7918 0.0401 -1.7169 -3.2847 -4.5255 -5.3698 -5.7983 -5.8117 -5.4089Columns 25 through 36-4.5802 -3.3225 -1.6778 0.2193 2.1277 3.7480 4.7881 4.9961 4.2100 2.4944 0.2607 -1.8728Column 37-3.4249vof =Columns 1 through 12-347.0624 -492.8494 -603.7818 -681.3764 -732.8054 -766.8573 -789.4516 -800.4475 -793.1801 -757.7000 -686.9200 -582.3864Columns 13 through 24-455.6459 -323.7142 -201.3442 -95.0327 -1.9924 86.4566 179.8406 284.1862 399.4460 519.0511 631.6416 724.6701Columns 25 through 36788.9608 822.4766 830.9271 823.5243 805.0980 769.5568 700.9268 584.3674 420.0967 225.6063 22.9226 -171.8563Column 37-347.0624aof =1.0e+004 *Columns 1 through 12-0.9048 -0.6847 -0.5221 -0.4123 -0.3300 -0.2449 -0.1309 0.0224 0.2041 0.3874 0.5412 0.6412Columns 13 through 240.6754 0.6479 0.5799 0.5058 0.4597 0.4593 0.4981 0.5531 0.5991 0.6163 0.5913 0.5159Columns 25 through 360.3871 0.2120 0.0106 -0.1886 -0.3676 -0.5349 -0.7184 -0.9361 -1.1628 -1.3170 -1.3111 -1.1438Column 37-0.9048。