北京科技大学matlab大作业
《数学实验》报告
实验名称MATLAB在研究物体振动方面的应用学院
专业班级
姓名
学号
2015年 1月
一、【实验目的】
物体振动这样一个看似简单但又包含着很多复杂计算的运动中，在人为的计算时是很难精确的实现，而通过MATLAB可以处理诸多科学中的许多问题，利用它来研究物理学中的机械振动，不仅特别方便还非常有效。
二、【实验任务】
本列举振动的一些实例，用matlab语言编制计算机程序进行仿真以达到研究简谐振动以及振动的合成，振动的计算以及受迫振动。
三、【实验程序】
（一）简谐振动介绍
最简单和最基本的振动是简谐振动．任何复杂的振动，都可以看成为许多简谐振动的合成．
1．特点
质点作简谐振动的条件是：在任何时候所受到的力与质点离开平衡位置的位移成正比，其指向与位移相反，始终指向平衡位置．所受的力与位移的关系表示为
（1）
式中为正的常数．对于弹簧振子，就是弹簧劲度系数
2．运动的微分方程及其解
根据牛顿第二定律，作简谐振动的质点的微分方程写成
即
（2）
式中。如下面的（3）和（4）所示，是简谐振动的圆频率。
微分方程（2）的解是
（3）
或
（4）
式（7.3）也可以表为复数形式
（5）
但要约定取其实数部分．
利用三角公式，很容易导出A ，和B，C之间的关系
即（6）
3．速度和加速度
作简谐振动的质点，它的速度和加速度很容易得到．只要将（7.3）对时间分别求导一次和求导两次即可，
（7）
（8）
式（1）、（2）、（3）、（4）、（5）都是判别一个系统是否作简道振动的依椐．
4．圆频率、周期和频率之间的关系
，，（9）
，，三者不是独立的，只要知道其中一个，就可以由（7.9）求出其余两个。它们是由振动系统的固有性质决定，常称为固有圆频率，固有周期和固有频率．
5．振幅和初周相
（3）中和是两个积分常数，可由初始条件决定．将初始条件：“，，”代入（3）和（7），得
（10）
解得
（11）
求解质点作简谐振动的具体运动情况，也就是要确定（7.3）中的，，三个值．其中和由
初始条件决定，因此一般来说，首先必须确定初始值和，而根据（7.10）或（7.11）求出和值．至于（或或），它是由系统固有性质决定的，与初始情况无关．例如对于弹簧振子，，完全由弹簧劲度系数和物体质量所决定．弹簧的大（即所谓硬的弹簧），振动的圆频率也就
大。而物体的质量m大，就小．
6．简谐振动系统的能量
作简谐振动的质点动能为
（12）
振动系统弹性势能为
（13）
因此系统总机械能为
（14）
系统的动能和势能各随时间作周期性变化，在振动过程中动能和势能互相转换，而总机械能保持不变．这是简谐振动的一个特性．总机械能E与振动的振幅平方A 2，振动的圆频率平方成正比．
动能和势能在一个周期内对时间的平均值分别是
（15）
注意和在一周期内对时间的平均值均等于1／2．这样，
（16）
7.振动的合成
一个质点同时参与两个振动方向相同、频率相同的简谐振动，
合振动仍为简谐振动
（17）
利用振幅矢量图示法容易求得
（18）
二个振动方向相同、频率略有差异的简谐振动，其合振动不为简谐振动，产生“拍”现象．拍频为
（
，
为两分振动频率）． （19）
二个振动方向互相垂直的简谐振动合成：
（1）若二振动频率相同，合振动轨迹一般为一椭圆.
（2）若二振动频率成整数比，合振动轨迹为规则的稳定的闭合曲线，称利萨如图．但若不成整数比，轨迹为不闭合的复杂曲线.
（二）实际运用
例1.关于物体振动的计算的应用
质量为1kg 的物体，以振动1x10^-4m 作简谐运动，其最大加速度为4.0m/s.求：（1）震动的周期；（2）通过平衡位置时的动能；（3）总能量；（4）物体在何处其动能和势能相等；
2
max a A ω=
max
a A ω=
2T πω=
通过平衡位置时的速度最大，所以得：
max 222
max 11
22k E mv m A ω==? max
k E E = 当
k p
E E =时，可得x 的位置
即： 2
22p
E x m ω=
程序如下：
m=1; %m 为物体的质量
amax=4.0; %amax 为最大加速度 A=1.0*10^-4; %A 为振幅 W=sqrt(amax/A); %求角速度 T=2*pi/W; %求周期
Ekmax=1/2*m*W*W*A*A; %求最大动能 E=Ekmax; %求总能量 Ep=1/2*E; %求势能
x=sqrt(2*Ep/m/W/W); %求动能和势能相等时的位移
例2 振动图
以质量为0.01kg 的物体作简谐运动，其振幅为0.08m ，周期为4s ，起始时刻物体在x=0.04m 处，向OX 轴负方向运动，试求：画出此时刻的0到8π的振动图形。
解: cos()x A t y ω=+ A=0.08m T=4s
2T π
ω=
T=0时，x=0.04；
得 0.04=0.08cosy
(0.08)cos[]
23x t ππ
=+ 程序： t=0:pi/200:8*pi;
x=0.08*cos(pi/2*t+pi/3); plot(t,x,'rp')
legend('x=0.08*cos(pi/2*t+pi/3') 结果如下图：
例3：阻尼振动
有一单摆在空气（室温为20℃）中来回摆动，其摆线长s=1.0m ，摆锤是一半径r=3
5.010m -?的铅球，求：
（1）摆动周期；（2）振幅减小10%所需的时间；（已知铅球密度为33
2.6510p kg m -=?，20℃时空气的
粘滞阻力5
1.7810n Pa s -=?）
解：
0g
l ω=
6r r F rnv C π=-=-
得:C=6πrn
阻尼系数
2c
k m =
02T π
ω=
有阻尼的情况下，单摆的振幅：
11
k A Ae t = （
10.9
A A =
得：1
1
ln(
)
0.9
t
k
=
程序如下：
g=input('输入g的值'); %g为重力加速度其值为
l=input('输入l的值'); %l为摆线长
p=input('输入p的值'); %p为铅球密度
n=input('输入n的值'); %20℃时空气的粘滞阻力
r=input('输入r的值'); %r为小球的半径
W=sqrt(g/l); %求角速度
C=6*pi*r*n; %由粘滞阻力Fr=-6 rnv=-Cr
m=4/3*pi*r^3*p; %求小球的质量
k=C/2/m; %求阻尼系数
T=2*pi/W; %求单摆周期
t1=log(1/0.9)/k; %振幅减小10%所需的时间
T
t1
输入g的值0.98
输入l的值1.0
输入p的值2.65*10^-3
输入n的值1.78*10^-5
输入r的值5.0*10^-3
例4：相互垂直的简谐振动的合成
程序：t=1:0.001:15
a1=input('振幅1=');w1=input('频率1=');phi1=input('初相位1='); a2=input('振幅2=');w2=input('频率2=');phi2=input('初相位2='); x=a1*cos(w1*t+phi1);
y=a2*cos(w2*t+phi2);
subplot(2,2,1),plot(t,x),title('x轴上谐振1')
subplot(2,2,4),plot(y,t),title('y轴上谐振2')
subplot(2,2,3),comet(x,y),ylabel('y'),xlabel('x'),title('李萨如图形')
振幅1=10
频率1=2
初相位1=pi
振幅2=10
频率2=4
相位2=10
例5：关于平面简谐波和简谐振动
一余弦波在弦上传播，其波函数为
式中和的单位为，t的单位为
1、试求其振幅、波长、频率、周期和波速。
2、分别画出对应和俩时刻弦上的波形图。
解：由已知的波函数求波动的特征量，我们一般采用比较系数法，由
上式说明此简谐波向正方向传播，将它与波函数的标准形式相比较得：,,
程序：x=-0.25:0.001:0.25; %设定x的取值范围
for i=0:2 %用循环语句令i分别取0,1,2
t=0.0025*i;
y=0.002*cos(pi*5*x-200*pi*t); %此为时刻t的波函数
if i==0; %用选择语句分别用不同的颜色线型画不同时刻的波形图plot(x,y,'k-') %用黑色实线画t=0时刻的波形图
hold on %保存图形的命令，否则后一幅图会覆盖前一幅
grid on %绘制网格
else if i==1;
plot(x,y,'r--') %用红色虚线画t=0.0025时刻的波形图
else if i==2;
plot(x,y,'b-.') %用蓝色虚点线画t=0.005时刻的波形图
end %结束if语句
end
end
end %结束for语句
四．【实验结果】
机械原理大作业
机械原理大作业 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、滑移齿轮变速传动中每对齿轮的几何尺寸及重合度
matlab期末大作业
电气学科大类 Modern Control Systems Analysis and Design Using Matlab and Simulink Title: Automobile Velocity Control Name: 巫宇智 Student ID: U200811997 Class：电气0811
电气0811 巫宇智 Catalogue Preface (3) The Design Introduction (4) Relative Knowledge (5) Design and Analyze (6) Compare and Conclusion (19) After design (20) Appendix (22) Reference (22)
Automobile Velocity Control 1.Preface: With the high pace of human civilization development, the car has been a common tools for people. However, some problems also arise in such tendency. Among many problems, the velocity control seems to a significant challenge. In a automated highway system, using the velocity control system to maintain the speed of the car can effectively reduce the potential danger of driving a car and also will bring much convenience to drivers. This article aims at the discussion about velocity control system and the compensator to ameliorate the preference of the plant, thus meets the complicated demands from people. The discussion is based on the simulation of MATLAB. Key word: PI controller, root locus
哈工大机械原理大作业凸轮 - 黄建青
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 机械原理大作业二 课程名称：机械原理 设计题目：凸轮机构设计 院系：能源学院 班级： 1302402 设计者：黄建青 学号： 1130240222 指导教师：焦映厚陈照波 设计时间： 2015年06月23日
凸轮机构设计说明书 1. 设计题目 设计直动从动件盘形凸轮机构，机构运动简图如图1，机构的原始参数如表1所示。 图1 机构运动简图 表1 凸轮机构原始参数
计算流程框图： 2. 凸轮推杆升程，回程运动方程及推杆位移、速度、加速度线图 2.1 确定凸轮机构推杆升程、回程运动方程 设定角速度为ω=1 rad/s (1) 升程：0°<φ<50° 由公式可得 )]cos(1[20 ?π Φh s -=
)sin( 20 1 ?π ωπΦΦh v = )cos(20 2 2 12?π ωπΦΦh a = (2) 远休止：50°<φ<150° 由公式可得 s = 45 v = 0 a = 0 (3) 回程：150°<φ<240° 由公式得： ()()22 0000200000002200000 0,2(1)(1)1,12(1)(1),2(1)s s s s s s s s s Φhn s h ΦΦΦΦΦΦn Φn ΦΦn h n s h ΦΦΦΦΦΦn Φn n ΦΦΦn hn s ΦΦΦΦΦn Φn ??????'?=---+<≤++?'-? ???''-? =----++ <≤++???'-??? ?'---?'=-++<≤++'-?? 201 00000010002001 000 00n (),(1)(1)n ,(1)(1)n (1),(1)s s s s s s s s Φh v ΦΦΦΦΦΦn Φn ΦΦn h v ΦΦΦΦn Φn n ΦΦΦn h v ΦΦΦΦΦn ΦΦn ω??ω??ω??'=- --+<≤++?'-? ?''-? =- ++<≤++?'-? ?'---'?=--++<≤++''-??
春MATLAB仿真期末大作业
MATLAB仿真 期末大作业 姓名：班级：学号：指导教师：
2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。
仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。
哈工大机械原理大作业凸轮机构第四题
Harbin Institute of Technology 机械原理大作业二 课程名称：机械原理 设计题目：凸轮机构设计 姓名：李清蔚 学号：1140810304 班级：1408103 指导教师：林琳
一.设计题目 设计直动从动件盘形凸轮机构，其原始参数见表 1 表一：凸轮机构原始参数 升程（mm ) 升程 运动 角（o） 升程 运动 规律 升程 许用 压力 角（o） 回程 运动 角（o） 回程 运动 规律 回程 许用 压力 角（o） 远休 止角 （o） 近休 止角 （o） 40 90 等加 等减 速30 50 4-5-6- 7多 项式 60 100 120
二.凸轮推杆运动规律 （1）推程运动规律（等加速等减速运动） 推程F0=90° ①位移方程如下： ②速度方程如下： ③加速度方程如下： （2）回程运动规律（4-5-6-7多项式） 回程，F0=90°，F s=100°，F0’=50°其中回程过程的位移方程，速度方程，加速度方程如下：
三.运动线图及凸轮线图 本题目采用Matlab编程，写出凸轮每一段的运动方程，运用Matlab模拟将凸轮的运动曲线以及凸轮形状表现出来。代码见报告的结尾。 1、程序流程框图 开始 输入凸轮推程回 程的运动方程 输入凸轮基圆偏 距等基本参数 输出ds,dv,da图像 输出压力角、曲率半径图像 输出凸轮的构件形状 结束
2、运动规律ds图像如下： 速度规律dv图像如下： 加速度da规律如下图：
3.凸轮的基圆半径和偏距 以ds/dfψ-s图为基础，可分别作出三条限制线（推程许用压力角的切界限D t d t,回程许用压力角的限制线D t'd t'，起始点压力角许用线B0d'')，以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点，即可满足压力角的限制条件。 得图如下：得最小基圆对应的坐标位置O点坐标大约为（13，-50）经计算取偏距e=13mm，r0=51.67mm.
机械原理大作业
机械原理大作业 二、题目（平面机构的力分析） 在图示的正弦机构中,已知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;
期末大作业报告
期末大作业报告 课程名称：数字图像处理 设计题目：车牌识别 学院：信息工程与自动化学院 专业：计算机科学与技术 年级：xxxxx 学生姓名：xxxxxxx（学号xxxxxxxxxxxxx） 指导教师：xxxx 日期：20XX.6.10 教务处制 车牌识别 摘要:数字图像处理技术是20世纪60年代发展起来的一门新兴学科，随着图像处理理论和方法的进一步完善，使得数字图像处理技术在各个领域得到了广泛应用，并显示出广阔的应用前景。MATLAB既是一种直观、高效的计算机语言,同时又是一个科学计算平台。它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。根据它提供的500多个数学和工程函数，工程技术人员和科学工作者可以在它的集成环境中交互或编程以完成各自的计算。MATLAB中集成了功能强大的图像处理工具箱。由于MA TLAB语言的语法特征与C语言极为相似，而且更加简单，更加符合科技人员对数学表达式的书写格式，而且这种语言可移植性好、可扩展性强，再加上其中有丰富的图像处理函数，所以MA TLAB在图像处理的应用中具有很大的优势。车牌识别技术是智能交通系统的重要组成部分，在近年来得到了很大的发展。本文从预处理、边缘检测、车牌定位、字符分割、字符识别五个方面，具体介绍了车牌自动识别的原理。并用MATLAB软件编程来实现每一个部分,最后识别出汽车牌照。 关键词:车牌识别、数字图像处理、MATLAB
一、设计原理 车辆牌照识别系统的基本工作原理为：将摄像头拍摄到的包含车辆牌照的图像通过视频卡输入到计算机中进行预处理，再由检索模块对牌照进行搜索、检测、定位，并分割出包含牌照字符的矩形区域，然后对牌照字符进行二值化并将其分割为单个字符，然后输入JPEG或BMP 格式的数字，输出则为车牌号码的数字。牌照自动识别是一项利用车辆的动态视频或静态图像进行牌照号码、牌照颜色自动识别的模式识别技术。其硬件基础一般包括触发设备、摄像设备、照明设备、图像采集设备、识别车牌号码的处理机等，其软件核心包括车牌定位算法、车牌字符分割算法和光学字符识别算法等。某些牌照识别系统还具有通过视频图像判断车辆驶入视野的功能称之为视频车辆检测。一个完整的牌照识别系统应包括车辆检测、图像采集、牌照识别等几部分。当车辆检测部分检测到车辆到达时触发图像采集单元，采集当前的视频图像。牌照识别单元对图像进行处理，定位出牌照位置，再将牌照中的字符分割出来进行识别，然后组成牌照号码输出。 二、设计步骤 1. 提出总体设计方案: （1）车牌图像预处理方法 因为车牌图像都是在室外拍摄的，所以不可避免地会受到光照、气候等因素的影响，而且拍摄者的手部抖动与车辆的移动会造成图像的模糊。要去除这些干扰就得先对车牌图像进行预处理。由于当前数码相机的像素较高，原始图像的数据一般比较大，输入的彩色图像包含大量颜色信息，会占用较多的存储空间，且处理时也会降低系统的执行速度。因此对图像进行识别等处理时，常将彩色图像转换为灰度图像，以加快处理速度。对图像进行灰度化处理后常用的方法是图像二值化、去除背景图像、增强处理、边缘检测、滤波等处理等。
《科学计算与MATLAB》期末大作业
杭州电子科技大学信息工程学院《科学计算与MATLAB》期末大作业
给出程序、图、作业分析，程序需加注释。 1. 试编写名为fun.m 的MATLAB 函数，用以计算下述的值： ?? ? ??-<->=t t n t t t n t f 的)4/sin()(si 对所有)4/sin(其他情况)sin(的)4/sin()(si 对所有)4/sin()(ππππ 绘制t 关于函数f(t)的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。 function y=fun()%定义函数 % t=-6*pi:pi/10:6*pi; %定义变量范围 y = (sin(pi/4)).*(sin(t)>sin(pi/4))+(sin(-pi/4)).*(sin(t)
2.解以下线性方程组 ??? ??=+=++=--3 530 42231 321321x x x x x x x x A=[2 -1 -1;1 1 4;3 0 5];%输入矩阵 B=[2;0;3]; %输入矩阵 X = A\B %计算结果 3.已知矩阵? ? ??? ???? ???=44434241 3433323124232221 14131211A 求： (1)A(2:3,2:3) (2)A(:,1:2) (3)A(2:3,[1,3]) (4)[A,[ones(2,2)；eye(2)]]
A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];%输入矩阵A(2:3,2:3) %输出矩阵 A(:,1:2) %输出矩阵 A(2:3,[1,3]) %输出矩阵 [A,[ones(2,2);eye(2)]] %输出矩阵
机械原理大作业
机械原理大作业三 课程名称: 机械原理 级: 者: 号: 指导教师: 设计时间: 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
MATLAB期末大作业模板
MATLAB应用技术 期末大作业 专业： 姓名： 学号： 分数
一、在一个图形窗口中以子图形式同时绘制正弦、余弦、正切、余切曲线。请写下完整代码，展示图形结果。（请标注题图和坐标轴，用不同颜色和不同线型分别绘制以上曲线）。（15分） 二、某公司员工的工资计算方法如下。 （1）工作时数超过120小时者，超过部分加发15%。 （2）工作时数低于60小时者，扣发700元。 （3）其余按每小时84元发。 根据员工的工时数，计算应发工资。请写下完整的程序代码，并任意输入一工时数（使用input 函数），将结果展示（使用disp 函数）利用该代码进行计算工资，请写下计算结果。（15分） 三、编写一个函数文件，使其能够产生如下的分段函数： ?? ? ??≥<<≤-=66225.0,25.05.15.0)(x x x x x x f 请编写完整的函数文件（保存函数文件名为hanshu.m ），并编写脚本文件代码，任意输入x 值（使用input 函数），在脚本文件中调用函数文件求)(x f ，展示结果（使用disp 函数），请写下计算结果。（15分） 四、将5个学生的6门功课的成绩存入矩阵P 中，进行如下处理： （1）分别求每门课的最高分、最低分及相应学生的序号。 （2）分别求每门课的平均分和标准差。 （3）5门课总分的最高分、最低分及相应学生序号。 （4）将5门课总分按从大到小顺序存入score 中，相应学生序号存入num 。 请将各小题的运行代码完整写下来，并写下运行结果。（20分） 五、请利用所学的MATLAB 知识，自主设计一个图形用户界面，请完整记录它的设计过程，需提供文字、代码和图片，以充分说明设计的图形用户界面可实现
哈工大机械原理大作业
连杆的运动的分析 一．连杆运动分析题目 图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;
matlab综合大作业(附详细答案)
m a t l a b综合大作业(附详细 答案) 标准化文件发布号：（9312-EUATWW-MWUB-WUNN-INNUL-DQQTY-
《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）： 1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5) 实验结果： A = 2)将矩阵A按列拉长得到矩阵B； 实验程序：B=A(:) 实验结果： B =
3)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的 矩阵C； 实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)] 实验结果： C = 4)寻找矩阵A中大于0的元素；] 实验程序：G=A(find(A>0)) 实验结果： G = 5)求矩阵A的转置矩阵D； 实验程序：D=A' 实验结果： D = 6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E； 实验程序：E=flipud(fliplr(A)) 实验结果： E =
7)删除矩阵A的第2列和第4列得到矩阵F； 实验程序：F=A; F(:,[2,4])=[] 实验结果： F = 8)求矩阵A的特征值和特征向量； 实验程序：[Av,Ad]=eig(A) 实验结果： 特征向量Av = + - + - - + + - 特征值Ad = 0 0 0 0 0 + 0 0 0 0 0 - 0 0 0 0 0 0 0 0 0 0 9)求矩阵A的每一列的和值； 实验程序：lieSUM=sum(A) 实验结果： lieSUM = 10)求矩阵A的每一列的平均值； 实验程序：average=mean(A) 实验结果： average = 2．符号计算（10分，每小题5分）： 1)求方程组20,0 uy vz w y z w ++=++=关于,y z的解； 实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z'); y= S. y, z=S. z
MATLAB期末作业(MATLAB应用)
Matlab在汽车振动分析中的应用 XXX （昆明理工大学交通工程学院，昆明650500） 摘要：在明确汽车振动的产生原因及其危害的基础上，对汽车振动进行了理论分析、MATLAB编程计算以及试验研究。结果表明：将MATLAB强大的数据处理和可视化技术应用于汽车振动分析与控制中，既可以验证理论分析的结果，又可以预测汽车响应，具有很高的实用价值。 关键字：MATLAB 汽车振动激励阻尼自由度 1 汽车振动的产生原因及其危害 机械振动是一种特殊形式的运动，激励、质量、弹性和阻尼是振动系统的四大要素。如果把汽车作为一个系统来研究，汽车本身就是一个具有质量、弹簧和阻尼的振动系统。由于汽车内部各部分的固有频率不同，汽车在行驶中常因路面不平、车速和运动方向的变化，车轮、发动机和传动系统的不平衡，以及齿轮的冲击等各种外部和内部的激振作用而极易产生整车和局部的强烈振动。汽车的这种振动使汽车的动力性得不到充分的发挥，使经济性变坏。同时，还要影响汽车的通过性、操纵稳定性和平顺性，使乘员产生不舒服和疲乏的感觉，甚至损坏汽车的零部件和运载的货物，缩短汽车的使用寿命[1]。因此，研究汽车振动的目的主要有两方面：一是降低振动对汽车零部件的损伤、对汽车使用性能的危害；二是试图利用振动为汽车设计服务，利用振动机理设计制造振动机械以减轻劳动强度，提高工作效率是不乏先例的。 2汽车振动的理论分析 2.1建立振动的力学模型 当一个实际振动系统较复杂时，建立的模型越复杂，越接近实际情况，模拟越逼真，但往往使分析困难；相反地，建立模型时，分析越容易，但得到的结果可能不精确，因此，在建立振动系统力学模型中，总是在求得简化表达和逼真模拟二者之间的折衷[2]。振动分析的关键就是：根据研究的内容和要求，把所研究的对象以及外界对它的作用简化为一个既简单又能在动态特性方面与原来的研究对象等效的力学模型。 汽车是由多个系统组成的复杂的振动系统，每个系统都存在振动问题。主要包括发动机、传动系统、制动系统、转向系统、悬架系统、车身和车架系统存在的振动问题。研究汽车这样一个复杂的振动系统，要根据所分析的问题进行简化，具体简化方案有以下几种 [3]： 1)当汽车对称于其纵轴线时，汽车车身只有垂直振动和俯仰振动对平顺性影响最大。此时，将汽车简化成如图1a)所示的四个自由度的平面模型，因轮胎阻尼较小，在此予以忽略。在这个模型中，车身质量m2，主要考虑垂直和俯仰两个自由度前、后车轴质量m1f，m1r。有两个垂直自由度。 2)当汽车前、后轴悬架质量分配达到一定值时，前、后悬架系统的垂直振动几乎是独立的。于是可以将汽车进一步简化为如图1b)所示的车身和车轮两个自由度振动系统。 M2：为簧载(车身)质量，m1为非簧载(车轮)质量。分析平顺性时，只考虑两个质量的垂直自由度。
机械原理大作业
机械原理大作业 课程名称：机械原理 设计题目：连杆机构运动分析 院系：机械工程院 班级： 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点的角速度和角加速度。
哈工大机械原理大作业凸轮机构设计题
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 代
MATLAB期末大作业模版
《MATLAB》期末大作业 学院土木工程与建筑学院 专业 班级 姓名 指导教师李琳 2018 年 5 月16 日
明 作业内容题目2:问题描述：在[0 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先划定x的范围与间距，再列出y的表达式，利用plot函数绘制二维曲线。 （2）软件说明及源代码 >> x = 0:pi/10:2*pi; >>y = cos(5*x).*sin(x); >>plot(x,y) （3）实验结果 题目4:问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c (2)求f=0的解 （1）问题分析 这是一个符号函数显示以及符号函数的求解问题，第一问先定义常量与变量，在写出f表达式，利用pretty函数显示f。第二问利用solve函数求解f=0时的解。 （2）软件说明及源代码
第一问 >> syms a b c x; >> f=a*x^2+b*x+c; >> pretty(f) 第二问 >>syms a b c x; >>f=a*x^2+b*x+c; >> solve(f) （3）实验结果 1、 2、 题目5:问题描述：求积分 （1）问题分析 这是一个利用符号函数求积分的问题，先定义变量x，再列出I1表达式，利用int函数求在范围0到Pi/2上的积分。 （2）软件说明及源代码 >> syms x; >> I1=(1-2*sin(2*x))^0.5; >> int(I1,0,0.5*pi) （3）实验结果 题目6:问题描述：分别随机产生一个6×6的整数矩阵（元素可在［－20,20］之间），求该随机阵的秩，特征值和特征向量。 （1）问题分析 这是一个矩阵运算问题，先利用rand函数产生一个6*6的元素在-20到20
哈工大-机械原理大作业3-齿轮-23题完整
1、设计题目 1.1机构运动简图 1.2机械传动系统原始参数 2、传动比的分配计算 电动机转速n=970r/min，输出转速n1=41 r/min，n2=37 r/min，n3=33 r/min，带传动的最大传动比i pmax=2.5,滑移齿轮传动的最大传动比i vmax=4，定轴齿轮传动的最大传动比i dmax=4。 根据传动系统的原始参数可知，传动系统的总传动比为 i1=n n1=970/41=23.659i2=n n2 =970/37=26.216i3=n n3 =970/33=29.394 传动系统的总传动比由带传动、滑移齿轮传动和定轴齿轮传动三部分实现。设带传动的传动比为i pmax=2.5，滑移齿轮的传动比为i v1、i v2和i v3，定轴齿轮传动的传动比为i f，则总传动比 i1=i pmax i v1i f i2=i pmax i v2i f
i 3=i pmax i v3i f 令i v3=i vmax =4 则可得定轴齿轮传动部分的传动比为i f = i 3 i pmax ×i vmax = 29.3942.5×4 =2.939 滑移齿轮传动的传动比i v1 = i 1 i pmax ×i f = 23.659 2.5×2.939 =3.220 i v2=i 2i pmax ×i f =26.216 2.5×2.939 =3.568 定轴齿轮传动由3对齿轮传动组成，则每对齿轮的传动比为 i d = i f 3= 2.9393 =1.432≤i dmax =4 3、齿轮齿数的确定 根据滑移齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮5、6、7、8、9和10为角度变位齿轮，其齿数：z 5=12，z 6=38，z 7=11，z 8=39，z 9=10，z 10=40；它们的齿顶高系数h a ?=1，径向间隙系数c ?=0.25，分度圆压力角α=20°，实际中心距a ＇=52mm 。 根据定轴齿轮变速传动系统中对齿轮齿数的要求，可大致选择齿轮11、12、13和14为角度变位 齿轮，其齿数：z 11=z 13=12，z 12=z 14=17。它们的齿顶高系数h a ? =1，径向间隙系数c ?=0.25， 分度圆压力角α=20°,实际中心距a ＇=45mm 。圆锥齿轮15和16选择为标准齿轮z 15=17，z 16=25， 齿顶高系数h a ?=1，径向间隙系数c ?=0.2，分度圆压力角α=20°（等于啮合角α＇） 。
MATLAB期末大作业
学号：姓名： 《Matlab/Simulink在数学计算与仿真中的应用》大作业1.假设地球和火星绕太阳运转的半径分别为r和2r，利用comet指令动画显示从地球 到火星的转移轨迹（r可以任意取值，要求实时显示探测器、太阳、地球和火星的位置）。 解函数function comet(varargin) [ax,args,nargs] = axescheck(varargin{:}); error(nargchk(1,3,nargs,'struct')); % Parse the rest of the inputs if nargs < 2, x = args{1}; y = x; x = 1:length(y); end if nargs == 2, [x,y] = deal(args{:}); end if nargs < 3, p = 0.10; end if nargs == 3, [x,y,p] = deal(args{:}); end if ~isscalar(p) || ~isreal(p) || p < 0 || p >= 1 error('MATLAB:comet:InvalidP', ... 'The input ''p'' must be a real scalar between 0 and 1.'); End 指令 %particle_motion t = 0:5:16013; r1=6.7e6;%随便给定参数 %--------------------------- r2=2*r1; g=9.8; R=6.378e6; m=g*R^2; %内轨道 v_inner=sqrt(m/r1); w_inner=v_inner/r1; x_inter=r1*cos(w_inner*t); y_inter=r1*sin(w_inner*t); %外轨道 v_outer=sqrt(m/r2); w_outer=v_outer/r2;