北科大_matlab第五次作业_多项式及其相关计算

《数学实验》报告

实验名称多项式及其相关计算学院计算机与通信工程学院专业班级

姓名

学号

2015年 5月

一、 【实验目的】

1、 学习多项式的拟合；

2、 学习二维多项式插值方法。

二、 【实验任务】

1.P130 T8

在钢线碳含量对于电阻的效应的研究中，得到以下数据，分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形。

2.P130 T12

用不同的方法对9

16z 2

2y x -=在（-3，3）上的二维插值效果进行比较。 三、 【实验程序】

1.P130 T8 clf

x=[0.10 0.30 0.40 0.55 0.70 0.80 0.95]; y=[15 18 19 21 22.6 23.8 26]; p1=polyfit(x,y,1) p3=polyfit(x,y,3); p5=polyfit(x,y,5);

disp('一阶拟合函数'),f1=poly2str(p1,'x') disp('三阶拟合函数'),f3=poly2str(p3,'x'); disp('五阶拟合函数'),f5=poly2str(p5,'x'); x1=0.1:0.01:0.95; y1=polyval(p1,x1); y3=polyval(p3,x1); y5=polyval(p5,x1);

plot(x,y,'rp',x1,y1,'g-',x1,y3,'k-',x1,y5,'m-');

legend('拟合点','一次拟合','三次拟合','五次拟合')

2.P130T12

clf

[x,y]=meshgrid(-3:0.5:3);

z=x.^2/16-y.^2/9;

[x1,y1]=meshgrid(-3:0.1:3);

z1=x1.^2/16-y1.^2/9;

figure(1)

subplot(1,2,1),mesh(x,y,z),title('数据点')

subplot(1,2,2),mesh(x1,y1,z1),title('函数图像')

[xi,yi]=meshgrid(-3:0.15:3);

zi1=interp2(x,y,z,xi,yi,'*nearest');

zi2=interp2(x,y,z,xi,yi,'*linear');

zi3=interp2(x,y,z,xi,yi,'*spline');

zi4=interp2(x,y,z,xi,yi,'*cubic');

figure(2)

subplot(2,2,1),mesh(xi,yi,zi1),title('最近点插值') subplot(2,2,2),mesh(xi,yi,zi2),title('线性插值') subplot(2,2,3),mesh(xi,yi,zi3),title('样条插值') subplot(2,2,4),mesh(xi,yi,zi4),title('立方插值')

四、【实验结果】

1.P130T8

2.P130T12

五、【实验总结】

1、学会了多项式的拟合的方法；

2、学会了二维多项式插值的方法。

哈工大机械原理大作业 凸轮机构设计 题

H a r b i n I n s t i t u t e o f T e c h n o l o g y 机械原理大作业二 课程名称： 机械原理 设计题目： 凸轮机构设计 一．设计题目 设计直动从动件盘形凸轮机构， 1.运动规律（等加速等减速运动） 推程 0450≤≤? 推程 009045≤≤? 2.运动规律（等加速等减速运动） 回程 00200160≤≤? 回程 00240200≤≤? 三．推杆位移、速度、加速度线图及凸轮s d ds -φ 线图 采用VB 编程，其源程序及图像如下： 1.位移： Private Sub Command1_Click() Timer1.Enabled = True '开启计时器 End Sub Private Sub Timer1_Timer() Static i As Single

Dim s As Single, q As Single 'i作为静态变量，控制流程；s代表位移；q代表角度 Picture1.CurrentX = 0 Picture1.CurrentY = 0 i = i + 0.1 If i <= 45 Then q = i s = 240 * (q / 90) ^ 2 Picture1.PSet Step(q, -s), vbRed ElseIf i >= 45 And i <= 90 Then q = i s = 120 - 240 * ((90 - q) ^ 2) / (90 ^ 2) Picture1.PSet Step(q, -s), vbGreen ElseIf i >= 90 And i <= 150 Then q = i s = 120 Picture1.PSet Step(q, -s), vbBlack ElseIf i >= 150 And i <= 190 Then q = i s = 120 - 240 * (q - 150) ^ 2 / 6400 Picture1.PSet Step(q, -s), vbBlue ElseIf i >= 190 And i <= 230 Then

机械原理大作业

机械原理大作业 This model paper was revised by the Standardization Office on December 10, 2020

机械原理大作业三 课程名称：机械原理 设计题目：齿轮传动设计 院系： 班级： 设计者： 学号： 指导教师： 设计时间： 1、设计题目 机构运动简图 机械传动系统原始参数

2、传动比的分配计算 电动机转速min /745r n =，输出转速m in /1201r n =，min /1702r n =， min /2303r n ，带传动的最大传动比5.2max =p i ,滑移齿轮传动的最大传动比4m ax =v i ，定轴齿轮传动的最大传动比4m ax =d i 。 根据传动系统的原始参数可知，传动系统的总传动比为： 传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。设带传动的传动比为5.2max =p i ，滑移齿轮的传动比为321v v v i i i 、、，定轴齿轮传动的传动比为f i ，则总传动比 令 4max 1==v v i i 则可得定轴齿轮传动部分的传动比为 滑移齿轮传动的传动比为 设定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为 3、齿轮齿数的确定 根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、7、8、9和10为角度变位齿轮，其齿数： 35,18,39,14,43,111098765======z z z z z z ；它们的齿顶高系数1=* a h ，径向间 隙系数25.0=*c ，分度圆压力角020=α，实际中心距mm a 51'=。

机械原理大作业2-齿轮机构分析

机械原理大作业2-齿轮机构分析

Harbin Institute of Technology 机械原理大作业三 题目：齿轮传动设计 院系：机电工程学院 班级： 姓名： 学号：

哈尔滨工业大学 1、设计题目 如图所示机械传动系统，运动由电动机1输入，经过机械传动系统变速后由圆锥齿轮16输出三种不同的转速，据下表中的原始数据，设计该传动系统。

2、传动比的分配计算 电动机转速n=745r/min，输出转速n1=23 r/min，n2=29 r/min，n3=35 r/min，带传动的最大传动比i pmax=2.8,滑移齿轮传动的最大传动比i vmax=4.5，定轴齿轮传动的最大传动比i dmax=4.5。 根据传动系统的原始参数可知，传动系统的总传动比为 i1=n/n1=745/35=21.286, i2=n/n2=745/29=25.690, i3=n/n3=745/23=32.391, 传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。 设带传动的传动比为i pmax=2.8，滑移齿轮的传动比为i v1, i v2 和i v3, 定轴齿轮传动的传动比为i f，则总传动比 i1= i pmax*i v1*i f, i2= i pmax*i v2*i f,

i3= i pmax*i v3*i f, 令i v3=i vmax=4.5，则可得定轴齿轮传动部分的传动比i f=i3/(i pmax*i vmax)= 32.391/(2.8*4.5)= 2.571, 滑移齿轮传动的传动比 i v1 =i1/(i pmax*i vmax) =21.286/(2.8*2.571)= 2.957 i v2 =i2/(i pmax*i vmax) =25.690/(2.8*2.571)= 3.569 定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为 id=3√i f= 3√2.571 =1.370 小于等于 i pmax = 4 3、设定齿轮齿数及基本参数 根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、7、8、9和10为角度变位齿轮，其齿数：z5 = 13，z6 = 38，z7 = 11，z8 =39，z9 = 9，z10 =40。它们的齿顶高系数h a* = 1，径向间隙系数c* = 0.25，分度圆压力角α = 20°，实际中心距a’= 51mm。 根据定轴齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮11、12、13和14为角度变位齿轮，其齿数：z11=z13=14，z12=z14=19。它们的齿顶高系数h a* =1，径向间隙系数c*=0.25，分度圆压力角α = 20°,实际中心距a’=51mm。 圆锤齿轮15和16选择为标准齿轮，其齿数：z15=17，z16=24。它们的齿顶高系数h a* =1，径向间隙系数c*=0.2，分度圆压力角α=20°。 4、滑移齿轮变速传动中每对齿轮的几何尺寸及重合度

机械原理大作业

机械原理大作业 二、题目（平面机构的力分析） 在图示的正弦机构中,已知l AB =100 mm,h1=120 mm,h2 =80 mm,W1 =10 rad/s(常数)，滑块2和构件3的重量分别为G2 =40 N和G3 =100 N,质心S2 和S3 的位置如图所示，加于构件3上的生产阻力Fr=400 N，构件1的重力和惯性力略去不计。试用解析法求机构在Φ1=60°、150°、220°位置时各运动副反力和需加于构件1上的平衡力偶M 。 b Array 二、受力分析图

三、算法 (1)运动分析 AB l l =1 滑块2 22112112/,/s m w l a s m w l v c c == 滑块3 21113113/cos ,sin s m l w v m l s ??== 212 113/sin s m w l a ?-= (2)确定惯性力 N w l g G a m F c 2 1122212)/(== N w l g G a m F 121133313sin )/(?-== (3)受力分析 i F F i F F x R D R x R C R 43434343,=-= j F j F F R R R 232323-==

j F i F j F i F F R x R y R x R R 2121121212--=+= j F F F y R x R R 414141+= 取移动副为首解副 ① 取构件3为分离体，并对C 点取矩 由0=∑y F 得 1323F F F r R -= 由0=∑x F 得 C R D R F F 4343= 由 ∑=0C M 得 2112343/cos h l F F R D R ?= ②取构件2为分离体 由0=∑x F 得 11212cos ?R x R F F = 由0 =∑y F 得 1123212sin ?F F F R y R -= ③取构件1为分离体，并对A 点取矩 由0=∑x F 得 x R x R F F 1241= 由0 =∑ y F 得 y R y R F F 1241= 由0=A M 得 1132cos ?l F M R b = 四、根据算法编写Matlab 程序如下： %--------------已知条件---------------------------------- G2=40; G3=100; g=9.8; fai=0; l1=0.1; w1=10; Fr=400; h2=0.8; %--------分布计算，也可将所有变量放在一个矩阵中求解------------------- for i=1:37 a2=l1*(w1^2); a3=-l1*(w1^2)*sin(fai); F12=(G2/g)*a2;

Matlab中多项式的运算

Matlab中多项式的运算 1:直接键入p=[ 1 2 3 4]系统就自动建立起多项式y=x3+2x2+3x+4 2:利用roots(p)命令就可以求解这个多项式的根,例如：>> p=[1 2 3 4 ]; >> roots(p) ans = -1.6506 -0.1747 + 1.5469i -0.1747 - 1.5469i 3：利用poly(a)命令可以由多项式的根求解多项式的系数，其中a=[ 2 3]的表示形式。例如：>> a=[2 3]; >> poly(a) ans = 1 -5 6 则这个多项式为x2+5x+6. 4：多项式的相加减：若干个多项式相加减时就是把它们的系数数组直接相加减，但是系数数组元素的个数必须相等，若不够，可以补0，例如：a=[1 2 3 4]; b=[2 2 1 2]; c=a+b c =[ 3 4 4 6]

5：利用conv(a,b)命令可以求解a,b两个多项式的乘积。例如：>> a=[1 2]; >> b=[1 -2]; >> c=conv(a,b) c = 1 0 -4.因为a=x+2,b=x-2,所以a,b的乘积为c=x2-4. 6:利用polyder(a)命令可以去、求多项式的微分（求导数），例如：>> a=[1 2 0 -5 6]; >> b=polyder(a) b =4 6 0 -5,其中a=x4+2x3-5x+6;b=4x3+6x2-5. 6;给出x的范围，利用polyval(a)命令可以求出x对应的多项式的值，例如：>> x=-1:0.1:2; >> a=[1 2 0 -5 6]; >> y=polyval(a,x); >> plot(x,y) >> grid 绘图结果如下：

哈工大机械原理大作业连杆

Harbin Institute of Technology 机械原理大作业一 课程名称： 机械原理 设计题目： 连杆机构运动分析 院 系： 机电工程学院 班 级： 设 计 者： 学 号： 指导教师： 设计时间： 1.运动分析题目 (11)在图所示的六杆机构中，已知： AB l =150mm, AC l =550mm, BD l =80mm, DE l =500mm,曲柄以等角速度1w =10rad/s 沿逆时针方向回转，求构件3的角速度、角加速度和构件5的位移、速度、加速度。 2.机构的结构分析 建立以点A 为原点的固定平面直角坐标系A-x, y,如下图： 机构结构分析 该机构由Ⅰ级杆组RR （原动件1）、Ⅱ级杆组RPR （杆2及滑块3）和Ⅱ级杆组RRP （杆4及滑块5）组成。 3.建立组成机构的各基本杆组的运动分析数学模型 原动件1（Ⅰ级杆组RR ） 由图所示，原动件杆1的转角a=0-360°，角速度1w =10rad/s ，角加速度1a =0,运动副A 的位置坐标A x =A y =0,速度

(A， A)，加速度 （A

， A ）， 原动件1的长度AB l =150mm 。 求出运动副B 的位置坐标（B x , B y ）、速度 （B

，B）和加速度 （B ， B）。

杆2、滑块3杆组（RPR Ⅱ级杆组） 已出运动副B 的位置（B x , B y ）、速度 （B ，B ） 和加速度

（B ， B ）， 已知运动副C 的位置坐标C x =0, C y =550mm,速度，加速度，杆长AC l =550mm 。 求出构件2的转角b,角速度2w 和角加速度2a . 构件二上点D 的运动

matlab实验五多项式和符号运算

实验五：Matlab多项式和符号运算 一、实验目的 1．掌握Matlab多项式的运算。 2．了解符号运算。 二、实验内容 1．将多项式()(2)(3)(7)(1) =-+-+化为x的降幂排列。 P x x x x x syms x; y=(x-2)*(x+3)*(x-7)*(x+1); expand(y) ans = x^4-5*x^3-19*x^2+29*x+42 2．求一元高次方程的根。 98765432 --++--++= 53015027313658204100576-28800 x x x x x x x x x syms x y; y=x^9-5*x^8-30*x^7+150*x^6-1365*x^4-820*x^3+410 0*x^2+576*x-2880; solve(y,x) ans = 6.81947687944124431946 1.42761488953013276419+.8192491831*i 2.865487219+2.49263348244446271927*i

-1.887673354+1.812452594*i -.9583509633 -5.922730991 -1.887673354-1.812452594*i 2.865487219-2.49263348244446271927*i 1.42761488953013276419-.8192491831*i 3．求一元高次方程的根，并画出左边多项式函数在[2,2] x∈-区间内的曲线。 42 -+= x x 210 a=[1 0 -2 0 1]; r=roots(a) syms x; x=-2:2; y=[1 0 -2 0 1]; plot(x,y) r = 1.0000 + 0.0000i 1.0000 - 0.0000i -1.0000 -1.0000

机械原理大作业

机械原理大作业三 课程名称: 机械原理 级: 者: 号: 指导教师: 设计时间: 1.2机械传动系统原始参数 设计题目： 系: 齿轮传动设计 1、设计题 目 1.1机构运动简图 - 11 7/7777777^77 3 UtH TH7T 8 'T "r 9 7TTTT 10 12 - 77777" 13 ///// u 2

电动机转速n 745r/min ,输出转速n01 12r/mi n , n02 17r /mi n , n°323r/min，带传动的最大传动比i pmax 2.5 ,滑移齿轮传动的最大传动比 i vmax 4，定轴齿轮传动的最大传动比i d max 4。 根据传动系统的原始参数可知，传动系统的总传动比为： 传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实 现。设带传动的传动比为i pmax 2.5，滑移齿轮的传动比为9、心、「3，定轴齿轮传动的传动比为i f，则总传动比 i vi i vmax 则可得定轴齿轮传动部分的传动比为 滑移齿轮传动的传动比为 设定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为 3、齿轮齿数的确定 根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、 7、8 9和10为角度变位齿轮，其齿数: Z5 11,Z6 43,Z7 14,Z8 39,Z9 18,乙。35 ；它们的齿顶高系数0 1，径向间隙

系数c 0.25，分度圆压力角200，实际中心距a' 51mm。 根据定轴齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮11、12、13和14为角度变位齿轮，其齿数：Z11 z13 13，乙 2 z14 24。它们的齿顶高系数d 1，径向间隙系数c 0.25，分度圆压力角200,实际中心距 a' 46mm。圆锥齿轮15和16选择为标准齿轮令13,乙 6 24，齿顶高系数 h a 1，径向间隙系数c 0.20，分度圆压力角为200（等于啮合角’）。 4、滑移齿轮变速传动中每对齿轮几何尺寸及重合度的计算 4.1滑移齿轮5和齿轮6

哈工大机械原理大作业

连杆的运动的分析 一．连杆运动分析题目 图1-13 连杆机构简图 二．机构的结构分析及基本杆组划分 1.。结构分析与自由度计算 机构各构件都在同一平面内活动，活动构件数n=5, PL=7,分布在A、B、C、E、F。没有高副，则机构的自由度为 F=3n-2PL-PH=3*5-2*7-0=1 2．基本杆组划分 图1-13中1为原动件，先移除，之后按拆杆组法进行拆分，即可得到由杆3和滑块2组成的RPR II级杆组，杆4和滑块5组成的RRP II级杆组。机构分解图如下：

图二 图一 图三 三．各基本杆组的运动分析数学模型 图一为一级杆组， ? c o s l A B x B =， ? sin lAB y B = 图二为RPR II 杆组， C B C B j j B E j B E y y B x x A A B S l C E y x S l C E x x -=-==-+=-+=0000 )/a r c t a n (s i n )(c o s )(?? ? 由此可求得E 点坐标，进而求得F 点坐标。 图三为RRP II 级杆组， B i i E F i E F y H H A l E F A l E F y y l E F x x --==+=+=111)/a r c s i n (s i n c o s ??? 对其求一阶导数为速度，求二阶导数为加速度。

lAB=108; lCE=620; lEF=300; H1=350; H=635; syms t; fai=(255*pi/30)*t; xB=lAB*cos(fai); yB=lAB*sin(fai); xC=0; yC=-350; A0=xB-xC; B0=yB-yC; S=sqrt(A0.^2+B0.^2); zj=atan(B0/A0); xE=xB+(lCE-S)*cos(zj); yE=yB+(lCE-S)*sin(zj); a=0:0.0001:20/255; Xe=subs(xE,t,a); Ye=subs(yE,t,a); A1=H-H1-yB; zi=asin(A1/lEF); xF=xE+lEF*cos(zi); vF=diff(xF,t); aF=diff(xF,t,2); m=0:0.001:120/255; xF=subs(xF,t,m); vF=subs(vF,t,m); aF=subs(aF,t,m); plot(m,xF) title('位移随时间变化图像') xlabel('t(s)'),ylabel(' x') lAB=108; lCE=620; lEF=300; H1=350; H=635; syms t; fai=(255*pi/30)*t; xB=lAB*cos(fai); yB=lAB*sin(fai); xC=0;

哈工大机械原理大作业二凸轮机构设计(29)

设计说明书 1 设计题目 如图所示直动从动件盘形凸轮机构，其原始参数见下表，据此设计该凸轮机构。 2、推杆升程、回程运动方程及位移、速度、加速度线图 2.1凸轮运动理论分析 推程运动方程： 01cos 2h s π?????=-?? ?Φ???? 1 00sin 2h v πωπ??? = ?ΦΦ?? 22 12 00cos 2h a πωπ???= ?ΦΦ?? 回程运动方程： ()0' 1s s h ?-Φ+Φ?? =- ??Φ ? ? 1'0 h v ω=- Φ 0a = 2.2求位移、速度、加速度线图MATLAB 程序 pi= 3.1415926; c=pi/180; h=140; f0=120; fs=45; f01=90; fs1=105; %升程 f=0:1:360; for n=0:f0

s(n+1)=h/2*(1-cos(pi/f0*f(n+1))); v(n+1)=pi*h/(2*f0*c)*sin(pi/f0*f(n+1)); a(n+1)=pi^2*h/(2*f0^2*c^2)*cos(pi/f0*f(n+1)); end %远休程 for n=f0:f0+fs s(n+1)=140; v(n+1)=0; a(n+1)=0; end %回程 for n=f0+fs:f0+fs+f01 s(n+1)=h*(1-(f(n+1)-(f0+fs))/f01); v(n+1)=-h/(f01*c); a(n+1)=0; end %近休程 for n=f0+fs+f01:360; s(n+1)=0; v(n+1)=0; a(n+1)=0; end figure(1);plot(f,s,'k');xlabel('\phi/\circ');ylabel('s/mm');grid on;title('推杆位移线图') figure(2);plot(f,v,'k');xlabel('\phi/\circ');ylabel('v/(mm/s)');grid on;title('推杆速度线图') figure(3);plot(f,a,'k');xlabel('\phi/\circ');ylabel('a/(mm/s2');grid on;title('推杆加速度线图') 2.3位移、速度、加速度线图

matlab多项式运算和方程组的求解

二、多项式 （1）多项式的表达式和创建 MATLAB中使用一维向量来表示多项式，将多项式的系数按照降幂次序存放在向量中。 例如：多项式2X4+3X3+5X2+1可以用向量[2 3 5 0 1]来表示。 例2-1，输入多项式3x4-10x3+15x+1000 在命令窗口输入： p=[3 -10 0 15 1000] 输出结果如下： （2）多项式求根 1、多项式的根 找出多项式的根，即使多项式为零的值，MATLAB提供了特定的函数roots求解多项式的根。 例2-2，求解多项式3x4-10x3+15x+1000的根。 在命令窗口输入：

输出的结果如下： 2、由根创建多项式 在MATLAB中，无论是一个多项式，还是它的根，都是以向量形式存储的，按照惯例，多项式是行向量，根是列向量。因此当我们给出一个多项式时，MATLAB 也可以构造出相应的多项式，这个过程需要使用函数poly。 例2-3 输入及结果 （3）多项式四则运算 1，多项式的加法 MATLAB并未提供一个特别的函数，如果两个多项式向量大小相同，那么多项

式相加时就和标准的数组加法相同。 例2-4 在命令窗口输入： a=[1 3 5 7 9];b=[1 2 4 6 8]; c=a+b 输出结果： C（x）=2x4+5x3+9x2+13x+17 2、多项式的乘法运算 在MATLAB中，函数conv支持多项式乘法（运算法则为执行两个数组的卷积）。例2-5 在命令窗口输入： a=[1 3 5 7 9]; b=[1 2 4 6 8]; c=conv(a,b) 输出的结果如下：

C(x)=x8+5x7+15x6+35x5+69x4+100x3+118x2+110x+72 PS：conv指令只能进行两个多项式的乘法，两个以上的多项式的乘法需要重复使用conv。 3、多项式的除法运算 在MATLAB中，由函数deconv完成的。 例2-6 在命令窗口输入： c=[1 5 15 35 69 100 118 110 72];b=[1 2 4 6 8]; [a,r]=deconv (c,b) 输出的结果： （4）多项式微分

哈工大机械原理大作业凸轮机构设计题

Har bi n I nst i t ute of Technol ogy 械原理大作业二课程名称：机械原理 设计题目：凸轮机构设计 凸轮推杆运动规律 1.运动规律（等加速等减速运动） 推程 0 450 推程 450900 2.运动规律（等加速等减速运动） 回程16002000 回程20002400 ds s 三．推杆位移、速度、加速度线图及凸轮d线图 采用VB编程，其源程序及图像如下： 1.位移： Private Sub Command1_Click（） Timer1.Enabled = True ' 开启计时器 End Sub Private Sub Timer1_Timer() Static i As Single

表角度 Picture1.CurrentX = 0 Picture1.CurrentY = 0 1 = i + 0.1 If i <= 45 Then q = i s = 240 * (q / 90) ^ 2 Picture1.PSet Step(q, -s), vbRed ElseIf i >= 45 And i <= 90 Then q = i s = 120 - 240 * ((90 - q) ^ 2) / (90 ^ 2) Picture1.PSet Step(q, -s), vbGreen ElseIf i >= 90 And i <= 150 Then q = i s = 120 Picture1.PSet Step(q, -s), vbBlack ElseIf i >= 150 And i <= 190 Then q = i s = 120 - 240 * (q - 150) ^ 2 / 6400 Picture1.PSet Step(q, -s), vbBlue Dim s As Single, q As Single 'i 作为静态变量，控制流程； s 代表位移； q 代

机械原理大作业

机械原理大作业 课程名称：机械原理 设计题目：连杆机构运动分析 院系：机械工程院 班级： xxxx 学号： xxxxx 设计者： xx 设计时间：2016年6月

一、题目 1-12：所示的六连杆机构中，各构件尺寸分别为：lAB =200mm，lBC=500mm，lCD=800mm，xF=400mm，xD=350mm，yD=350mm，w1=100rad/s，求构件5上的F点的位移、速度和加速度。 二、数学模型 1.建立直角坐标系 以F点为直角坐标系的原点建立直角坐标系X-Y,如下图所示。

2.机构结构分析 该机构由I级杆组RR（原动件AB）、II级杆组RRR（杆2、3）、II级杆组PRP （杆5、滑块4）组成。 3.各基本杆组运动分析 1.I级杆组RR（原动件AB） 已知原动件AB的转角

φ＝0－2Π 原动件AB的角速度 w=10rad/s 原动件AB的角加速度 α=0 运动副A的位置 xA=-400,yA=0 运动副A的速度 vA=0,vA=0 运动副A的加速度 aA=0,aA=0 可得： xB=xA+lAB*cos(φ) yB=yA+lAB*sin(φ) 速度和加速度分析： vxB=vxA-wl*AB*sin(Φ) vyB=vyA+w*lAB*sin(φ) axB=axA-w2*lAB*cos(φ)-e*lAB*sin(φ) ayB=ayA-w2*lAB*sin(φ)+e*lAB*cos(φ)

2.II级杆组RRR（杆2、3） 杆2的角位置、角速度、角加速度 lBC=500mm，lCD=800mm，xD=350mm，yD=350mm， ψ2=arctan﹛[Bo+﹙Ao2+Bo2-Co2﹚?]／﹙Ao+Bo﹚﹜ ψ3=arctan[﹙yC-yD)／(xC-xD)] Ao=2*LBC(xD-xB) Bo=2*LBC(yD-yB) lBD2=(xD-xB)2+(yD-yB)2 Co=lBC2+lBD2-lCD2 xC=xB+lBC*cos(ψ2) yC=xB+lBC*sin(ψ2) 求导可得C点的角速度和角加速度。

哈工大机械原理大作业——凸轮——2号

哈工大机械原理大作业——凸轮——2号

————————————————————————————————作者：————————————————————————————————日期：

Harbin Institute of Technology 机械原理大作业 课程名称：机械原理 设计题目：凸轮机构设计

一、设计题目 (1)凸轮机构运动简图： (2)凸轮机构的原始参数 序号升程升程运 动角 升程运 动规律 升程许 用压力 角 回程运 动角 回程运 动规律 回程许 用压力 角 远休止 角 近休 止角 14 90°120°余弦加 速度 35°90°3-4-5 多项式 65°80°70° (1) 推杆升程、回程运动方程如下： A.推杆升程方程： 设为1rad s ω= 升程位移为： ()() 1cos451cos1.5 2 h s π ψψψ ?? ?? =-=- ?? ? Φ ?? ?? 2 3 ψπ ≤≤升程速度为： ()() 1 1 00 sin67.5sin1.5 2 h v πωπ ψψωψ ?? == ? ΦΦ ?? 2 3 ψπ ≤≤升程加速度为： ()() 22 2 1 1 00 cos101.25cos1.5 2 h a πωπ ψψωψ ?? == ? ΦΦ ?? 2 3 ψπ ≤≤ B.推杆回程方程：

回程位移为： ()()345 111110156s h T T T ψ??=--+?? 1029 918 ψπ≤≤ 回程速度为： ()()2211110 3012h v T T T ωψ=- -+'Φ 1029 918ψπ≤≤ 回程加速度为： ()()22 11112 60132h a T T T ωψ=--+'Φ 1029918ψπ≤≤ 其中：() 010 s T ψ-Φ+Φ= 'Φ 1029 918 ψπ≤≤ (2) 利用Matlab 绘制推杆位移、速度、加速度线图 A. 推杆位移线图 clc clear x1=linspace(0,2*pi/3,300); x2=linspace(2*pi/3,10*pi/9,300); x3=linspace(10*pi/9,29*pi/18,300); x4=linspace(29*pi/18,2*pi,300); T1=(x3-10*pi/9)/(pi/2); s1=45*(1-cos(1.5*x1)) s2=90; s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5)); s4=0; plot(x1,s1,'r',x2,s2,'r',x3,s3,'r',x4,s4,'r') xlabel('角度ψ/rad'); ylabel('位移s/mm') title('推杆位移线图') grid axis([0,7,-10,100]) 得到推杆位移线图：

哈工大机械原理大作业24题

班级 1013102 学号 6 机械原理大作业说明书 题目 1、连杆机构运动分析 2、凸轮机构设计 3、齿轮传动设计 学生姓名

1连杆机构运动分析1.设计题目：

一、先建立如下坐标系： 二、划分杆组如下，进行结构分析： 该机构由I级杆组RR（如图1）、II级杆组RPR（如图2、3）和II级杆组RRP（如图4）组成。 （1）（2） （3）（4）

三、运动分析数学模型： （1）同一构件上点的运动分析： 如右图所示的原动件1,已知杆1的角速度=10/rad s ω，杆长1l =170mm,A y =0，A x =110mm 。可求得下图中B 点的位置B x 、B y ，速度xB v 、yB v ，加速度xB a 、yB a 。 θcos 1l xB =，θsin 1l yB = θωυsin 1l xB -=，θωυcos 1l yB =， 222B 2==-cos =-B xB i d x a l x dt ω?ω 222 2 ==-sin =-B yB i B d y a l y dt ω?ω。 （2）RPRII 级杆组的运动分析： a. 如右图所示是由2个回转副和1个移 动副组成的II 级组。已知两个外运动副C 、B 的位置（B x 、B y 、c x =110mm 、C y =0）、速度（xB υ，yB υ， xC υ=0， yC υ=0）和加速度 （0,0,,==yC xC yB xB a a a a ）。可确定下图中D 点的位置、速度和加速度。确定构件3的角位移1?、角速度1ω、角加速度1α。 1sin 31..??l x dt dx C B -= 1sin 131cos 13.....2????l l x dt x d C B --= 1cos 31..??l y dt dy C B += 1cos 131sin 13.....2????l l y dt y d C B +-= 根据关系：1111d 122..11. α??ω??====dt d dt ， 故可得出： D x =)1cos( 4β?++l x C

机械原理大作业二-凸轮机构设计

机械原理大作业二 课程名称：机械原理 设计题目：凸轮机构设计 院系：机械设计制造及其自动化 班级：1208104 完成者：郑鹏伟 学号：1120810416 指导教师：林琳刘福利 设计时间：2014年6月4日 哈尔滨工业大学

一、 设计题目： 凸轮的机构运动简图如下图所示 : 序 号 升程 （mm ） 升程运动角（°） 升程运 动规律 升程 许用 压力角 （°） 回程运 动角 （°） 回程运动规律 回程许用压力角（°） 远休止角（°） 近休止角（°） 14 90 120 余弦 加速 度 35 90 等减等加速 65 75 75 二、 凸轮推杆升程、回程运动方程及推杆位移，速度加速度线图： （1）凸轮推杆升程运动方程： 根据题意知： 001207590 75s s Φ=Φ='Φ='Φ= （1）从动件升程运动方程（设为1rad s ω=）

1 2 22 12 s [1cos()]2sin()2cos( ) 2h h h a π?πωπ υ?πωπ ?=-Φ= ΦΦ= ΦΦ （2）从动件远休止运动方程 在远休止s Φ段，即213 3 12π?π ≤≤时，90,0,0s h mm a υ====。 （3）从动件回程运动方程 升程段采用等减等加运动规律，运动方程为： ①当回程 0002s s ?'ΦΦ+Φ≤≤Φ+Φ+ 134()123π?π≤≤时： 20s 20 1 02 2 12 2[-+]4[()]4s h s h h h a ?ωυ?ω=- ΦΦ'Φ=- -Φ+Φ'Φ=-'Φ（） ②当回程 0002s s ?'Φ'Φ+Φ+ ≤≤Φ+Φ+Φ419()312π?π≤≤时： 2 002 01 0020 2 12 2[)]4[)]4s s h s h h a ?ωυ?ω'= Φ+Φ+Φ-'Φ'=- Φ+Φ+Φ-'Φ='Φ（（ （4）从动件近休止运动方程 在近休止s 'Φ段，即19 212π?π ≤≤时，s 0,0,0a υ===。 （2）推杆位移 %t 表示转角 s 表示位移 t=0:0.01:2/3*pi; %升程阶段

MATLAB数据分析与多项式计算_习题答案

第6章 MATLAB数据分析与多项式计算 习题6 一、选择题 1．设A=[1,2,3,4,5;3,4,5,6,7]，则min(max(A))的值是（）。B A．1 B．3 C．5 D．7 2．已知a为3×3矩阵，则运行mean(a)命令是（）。B A．计算a每行的平均值B．计算a每列的平均值 C．a增加一行平均值D．a增加一列平均值 3．在MATLAB命令行窗口输入下列命令： >> x=[1,2,3,4]; >> y=polyval(x,1); 则y的值为（）。D A．5 B．8 C．24 D．10 4．设P是多项式系数向量，A为方阵，则函数polyval(P,A)与函数polyvalm(P,A)的值（）。D A．一个是标量，一个是方阵B．都是标量 C．值相等D．值不相等 5．在MATLAB命令行窗口输入下列命令： >> A=[1,0,-2]; >> x=roots(A); 则x(1)的值为（）。C A．1 B．-2 C．D． 6．关于数据插值与曲线拟合，下列说法不正确的是（）。A A．3次样条方法的插值结果肯定比线性插值方法精度高。 B．插值函数是必须满足原始数据点坐标，而拟合函数则是整体最接近原始数据点，而不一定要必须经过原始数据点。 C．曲线拟合常常采用最小二乘原理，即要求拟合函数与原始数据的均方误差达到极小。 D．插值和拟合都是通过已知数据集来求取未知点的函数值。 二、填空题 1．设A=[1,2,3;10 20 30;4 5 6]，则sum(A)= ，median(A)= 。 [15 27 39]，[4 5 6[ 2．向量[2,0,-1]所代表的多项式是。2x2-1

机械原理大作业一

连杆机构的运动分析 一.题目 如图所示是曲柄摇杆机构，各构件长度分别为a，b，c，d，试研究各构件长度的变化对机构急回特性的影响规律。 二.机构分析 四连杆机构可分为如下两个基本杆组 Ⅰ级杆组 RRRⅡ级杆组 AB为曲柄，做周转运动；CD为摇杆，做摆动运动； BC为连杆；AB，CD均为连架杆，AB为主动件。

三.建立数学模型 θ为极位夹角，φ为最大摆角 必须满足条件为：1.a≤b，a≤c，a≤d（a为最短杆）； 2.L min+L max≤其他两杆之和。 下面分析杆长和极位夹角的关系： 在△AC2B中， =； 在△AC1B中， =。 θ=- K=

最后分以下四种情况讨论： 1.机架长度d变化 令a=5，b=30，c=29 d由6开始变化至54，步长为1 输出杆长a，b，c，d和K。 2.连杆长度b变化 令a=5，b=29，d=30 b由6开始变化至54，步长为1 输出杆长a，b，c，d和K。 3.摇杆长度c变化 令a=5，b=29，d=30 c由6开始变化至54，步长为1 输出杆长a，b，c，d和K。 4.曲柄长度a变化 令b=29，c=28，d=30 a由5开始变化至27，步长为1 输出杆长a，b，c，d和K。

四．MATLAB计算编程a=5;b=30;c=29; d=6:1:54; m=(d.^2-216)./(50.*d); n=(384+d.^2)./(70.*d); p=acos(m); q=acos(n); w=p-q; o=(w.*180)/3.14; K=(180+o)./(180-o); fprintf('%.6f\n',K); plot(d,K,'b') xlabel('机架长度d变化时 '); ylabel('极位夹角/度'); tilte('极位夹角变化图'); ———————————————————————————————————— ——— a=5;d=30;c=29; b=6:1:54; m=((b-5).^2+59)./(60.*(b- 5)); n=(59+(b+5).^2)./(60.*(b+ 5)); p=acos(m); q=acos(n); w=p-q; o=(w.*180)/3.14; K=(180+o)./(180-o); fprintf('%.6f\n',K); plot(b,K,'b') xlabel('连杆长度b变化时'); ylabel('极位夹角/度'); tilte('极位夹角变化图');

哈工大机械原理大作业三上传版

哈尔滨工业大学 机械原理大作业三 齿轮传动系统设计说明书 题目：（3） 课程名称：机械原理 学院：外国语学院 姓名：XX 班号：XXX 学号：XXX

一：设计题目 二：传动比的分配计算 根据传动系统的原始参数可知，传动系统的总传动比为： 667.9615 14501 3===n n i 048.692114502 2===n n i 769.552614503 1===n n i 传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。设带传动的传动比为5.2max =p i ，滑移齿轮的传动比为321,v v v i i i 和，定轴齿轮传动的传动比为f i ，则总传动比： f v p i i i i 1max 1 = f v p i i i i 2max 2 = f v p i i i i 3max 3= 令： 4max 3 ==v v i i 则可得定轴齿轮传动部分的传动比： 677.94 5.2677 .96max max 3=×== v p f i i i i 滑移齿轮传动的传动比： 305.2677 .95.2769.55max 11 =×==f p v i i i i

854.2677 .95.2048 .69max 22 =×= = f p v i i i i 定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为 4≤131.2677.9max 33====d f d i i i 三：齿轮齿数的确定 根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、7和8为角度变位齿轮，其齿数分别为：22,51,19,54它们的齿顶高系数为1径向间隙系数 25.0*=c ，齿轮9与10齿顶高系数为0.8，C=0.3，采用短齿。分度圆压力角α=20°， 实际中心距75' =a 。 根据定轴齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮11、12、13和14为角度变位齿轮，其齿数：42,20,42,2014 131211 ====z z z z 。它们的齿顶高系数* a h =1， 径向间隙系数* c =0.25，分度圆压力角α=20°,实际中心距93' =a mm 。圆锤齿轮15和16 选择为标准齿轮42,202515 ==z z ，齿顶高系数 * a h =1，径向间隙系数* c =0.2，分度圆压力 角α=20°（等于啮合角' α）。 四：滑移齿轮变速传动中每对齿轮几何尺寸及重合度的计算 表1 滑移齿轮5、6参数 序号 项目 代号 计算公式及计算结果 1 齿数 齿轮5 Z 5 22 齿轮6 Z 6 51 2 模数 m 2 3 压力角 α 200 4 齿顶高系数 *a h 1