MATLAB论文 1245346 李星辰Word格式.docx
《MATLAB论文 1245346 李星辰Word格式.docx》由会员分享,可在线阅读,更多相关《MATLAB论文 1245346 李星辰Word格式.docx(17页珍藏版)》请在冰豆网上搜索。
按系统受力情况分,有自由振动、衰减振动和受迫振动。
按弹性力和阻尼力性质分,有线性振动和非线性振动。
振动又可分为确定性振动和随机振动,后者无确定性规律,如车辆行进中的颠簸。
振动是自然界和工程界常见的现象。
振动的消极方面是:
影响仪器设备功能,降低机械设备的工作精度,加剧构件磨损,甚至引起结构疲劳破坏;
振动的积极方面是:
有许多需利用振动的设备和工艺(如振动传输、振动研磨、振动沉桩等)。
振动分析的基本任务是讨论系统的激励(即输入,指系统的外来扰动,又称干扰)、响应(即输出,指系统受激励后的反应)和系统动态特性(或物理参数)三者之间的关系。
20世纪60年代以后,计算机和振动测试技术的重大进展,为综合利用分析、实验和计算方法解决振动问题开拓了广阔的前景。
在物体振动这样一个看似简单但又包含着很多复杂计算的运动中,在人为的计算时是很难精确的实现,而通过MATLAB可以处理诸多科学中的许多问题,利用它来研究物理学中的机械振动,不仅特别方便还非常有效。
例如,通过它来对单个波的振动情况以及同一方向上两列频率相同的波的合成情况都可以进行研究。
再者,在求解复杂问题时也会涉及求解方程组、阻尼振动以及受破振动等问题。
上述讨论中遇到的问题,在MATLAB中都可以快速而准确解决,为我们解决问题上都会有很大的帮助,和带来更多方便。
而且MATLAB可以多次使用,多次改变数值,对一些有规律性的问题的验证也是一大帮助。
下面就来简单的介绍一下MATLAB在物体振动方面的具体应用。
二、实例应用
例一:
关于物体震动的计算的应用
质量为0.10kg的物体,以振幅
作简谐运动,其最大加速度为4.0m/s.求:
(1)震动的周期;
(2)通过平衡位置时的动能;
(3)总能量;
(4)物体在何处其动能和势能相等;
(1)
通过平衡位置时的速度最大,所以得:
(2)
(3)
当
时,可得x的位置
即:
(4)
程序如下:
m=0.01;
%m为物体的质量
amax=4.0;
%amax为最大加速度
A=1.0*10^-2;
%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);
%求动能和势能相等时的位移
T,
Ekmax,
E,
x
T=
0.3142
Ekmax=
2.0000e-004
E=
x=
0.0071
例二:
振子的计算
以质量为0.01kg的物体作简谐运动,其振幅为0.08m,周期为4s,起始时刻物体在x=0.04m处,向OX轴负方向运动,试求:
画出此时刻的0到4
的振动图形。
解题思路:
(5)
A=0.08m(6)
T=4s(7)
(8)
T=0时,x=0.04;
得0.04=0.08cosy
得到
(9)
程序如下:
t=0:
pi/200:
4*pi;
y=0.08*cos(pi/2*t+pi/3);
plot(t,y)
运行结果如图1所示
图1振子计算运行结果
例三:
关于阻尼振动方面计算的应用。
有一单摆在空气(室温为20℃)中来回摆动,其摆线长s=1.0m,摆锤是一半径r=
的铅球,求:
(1)摆动周期;
(2)振幅减小10%所需的时间;
(已知铅球密度为
,20℃时空气的粘滞阻力
)
解题步骤如下:
粘滞阻力为:
(10)
得:
C=6
rn(11)
阻尼系数
(12)
(13)
有阻尼的情况下,单摆的振幅:
(14)
(15)
得:
(16)
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=-6rnv=-Cr
m=4/3*pi*r^3*p;
%求小球的质量
k=C/2/m;
%求阻尼系数
%求单摆周期
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
6.3470
t1=
1.7429e-004
例四:
相互垂直的简谐振动的合成李萨如图
如果两个简谐振动分别在x轴和y轴上进行,他们的振动方程分别为
合成后,可得质点的轨迹为椭圆方程
若两分振动有不同的频率,且两频率之比为有理数时,则合成后的质点运动具有稳定、封闭的轨迹。
称其为李萨如图形。
t=1:
0.001:
10;
%设定时间的范围从1到10
a1=input('
振幅1='
w1=input('
频率1='
phi1=input('
初相位1='
%要求用户输入x轴上简谐振动1的状态参量
a2=input('
振幅2='
w2=input('
频率2='
phi2=input('
初相位2='
%要求用户输入y轴上简谐振动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=20
频率1=>
>
4
初相位1=pi
振幅2=20
频率2=2
初相位2=20
图二李萨如图形
拍现象:
振幅1=20频率1=4初相位1=pi振幅2=20频率2=2初相位2=20
例五:
关于平面简谐波和简谐振动
一余弦波在弦上传播,其波函数为
式中
和
的单位为
,t的单位为
1、试求其振幅、波长、频率、周期和波速。
2、分别画出对应
俩时刻弦上的波形图。
解:
由已知的波函数求波动的特征量,我们一般采用比较系数法,由
(17)
上式说明此简谐波向
正方向传播,将它与波函数的标准形式相比较得:
程序如下
x=-0.25:
0.25;
%设定x的取值范围
fori=0:
2%用循环语句令i分别取0,1,2
t=0.0025*i;
y=0.002*cos(pi*5*x-200*pi*t);
%此为时刻t的波函数
ifi==0;
%用选择语句分别用不同的颜色线型画不同时刻的波形图
plot(x,y,'
k-'
)%用黑色实线画t=0时刻的波形图
holdon%保存图形的命令,否则后一幅图会覆盖前一幅
gridon%绘制网格
elseifi==1;
r--'
)%用红色虚线画t=0.0025时刻的波形图
elseifi==2;
b-.'
)%用蓝色虚点线画t=0.005时刻的波形图
end%结束if语句
end
end%结束for语句
图三平面简谐波和简谐振动
例六:
弹簧振子在策动力、弹性力和阻尼力的作用下做受迫振动,其振动方程为
令
,
,则上式可写成
同样先将受迫振动的微分方程分解为一阶微分方程:
(18)
则原方程化为
(19)
f=input('
单位质量物体所受的最大策动力fo='
%输入函数所需的参数
w=input('
策动力的频率w='
wo=input('
振动物体的固有频率wo='
beita=input('
阻尼因子beita='
[t,x]=ode23('
shoupozd'
[0:
10],[6,8],[],f,w,wo,beita);
%调用ode23,[0:
10]为t的积分区间,
%[6,8]为初始条件(这里只是任取的),即x=6,x'
=8,[]说明以后是输入参数,f,w,wo,beita是参数名,其顺序必须与所调用的%ode文件的参数顺序相同
plot(t,x(:
1))%取x
(1),绘制t-x
(1)位移曲线
xlabel('
t'
ylabel('
title('
受迫振动'
单位质量物体所受的最大策动力fo=10
策动力的频率w=5
振动物体的固有频率wo=15
阻尼因子beita=0.1
图四受迫运动
三、结论
振动是物体的一种很普通的运动形式,所谓机械振动是物体在一定位置附近所作的周期性往复运动。
例如,心脏的跳动、钟摆的摆动、活塞的往复运动、固体原子的振动等等。
这种运动都是在某一数值附近作往复的周期性运动,而在我们所学过的知识中,我们仅仅能解决的只是一些非常理想的振动——无阻尼振动,在处理这些问题时,遇到的问题都是很容易解决的,不需要很复杂的变换,以及涉及一些特殊角度问题,而在余弦函数中特殊角度很容易解决,而一般的角度却是无法精确计算的,这对一些求精确值上从在的一系列问题,在MATLAB中都可以解决,对于总的结果的影响也会相对小很多,为计算以及解决其他问题上也是一个很大的前进。
同时,以上给出的例子,运用了MATLAB关于计算和绘图方面的功能,MATLAB还有许多功能有待利用,从这些例子中可以看出只要具备一般的高等数学知识,便可以轻松地利用MATLAB解决这些问题,从其他计算机高级语言解决这些问题问题却要繁琐得多,用图像和图形来表征时间和空间变化的物理量或物理现象,有助于把抽象思维化为形象思维,开拓学生的视野,激发学习兴趣,启迪思维,从而获得一种解决问题的方法,把MATLAB运用于理工科特别是普通物理的教学中,在理工科教学引入计算机辅助教学,改革传统的理工科教学模式,有很大的应用前景,这中间还有很多问题有待于探索,我认为用MATLAB处理物理方面的应用是很有必要的,可以提高我们的学习效率,增强我们对理论的理解,提高了我们的分析,解决问题的能力,为今后的学习工作打下基础。
四课程体会
在第一章中,讲的是一些matlab的入门技术,除了一些基本操作与介绍之外,还初步认识了简单指令的编制,认识了一些matlab的特殊符号,例如运算用到的加减乘除。
(+—*/)对我来说还是比较容易上手的。
但在第二章开始,开始有点难度了,在第二章我接触到两种数据类型,一个是double,还有就是char,另外还介绍了赋值语句,数学计算,常用函数,输入输出语句和数据文件。
其中让我最容易混淆的就是运算的优先等级,当所需要运算的公式较长时,常常因没有弄好运算优先级而把程序编错,这一章也接触了更多的特殊符号,在多次看书之后,还是顺利把老师所布置的作业完成了。
在第三章中,我接触到基本的matlab选择结构,还有控制这个结构的关系运算符和逻辑运算符。
也就是if结构。
这个结构对运算起到很大作用,跟elseif配合使用的话,就可以将很图方法。
还有,我们学习如何控制画图的附加功能,例如线的宽度和符号的颜色。
这些属性可由指定的“propertyname”和值Value决定,“propertyname”和值Value将出现在plot命令的数据后。
在这个学期的学习中,主要还是围绕1到3章学习,学了一些初步的矩阵运算以及画图方法。
在书上的练习中也得到了实践,各种矩阵的运算也可以运用到现在所学的电路计算中。
活学活用,的确为我们提供了不少的方便。
但现今所学到的matlab技术还只是很初步,但是我已经了解到了matlab的实用性,所以还是会在今后继续学习这门课程的。
学习了MATLAB这门课程,我了解该软件的基本功能,也知道了该软件在我们生活中的重要地位。
随着社会的不断发展,科技的不断进步,计算机的普及,它也被应用在越来越多的方面。
MATLAB的基本单位是矩阵,它的指令表达式与数学,工程中常用的形式十分相似,,故用MATLAb来解决问题要比C语言等完成相同的事情简捷得多,MATLAB的最突出的特点就是简洁。
MATLAB相对于其他的一些编程软件有许多的优点:
一语言简洁,使用方便灵活,库函数极其丰富。
二运算符丰富。
三MATLAB既具有结构化的控制语句(如for循环while循环),又有面对对象编程的特性。
四语言限制不严格,程序设计自由度大。
五程序的可移植性很好,基本上不做修改就可以在各种型号的计算机和操作系统上运行。
六MATLAB的图形功能强大。
用matlab创建矩阵时,方法有两种:
第一可以直接依次输入矩阵各行各列的元素,但矩阵元素必须用[]括住,矩阵元素必须用逗号或空格分隔,在[]内矩阵的行与行之间必须用分号分隔。
第二用MATLAB函数创建矩阵。
MATLAB可以进行矩阵的加减、乘除的元素,求可逆矩阵、转置矩阵、求矩阵的特征值,求线性方程组等等、
MATLAB商务功能是非常强大的,MATLAB不仅有强大的运算功能,它还有强大的绘图功能,我对它的了解也仅仅就是一点点,或许说还没有入门,比如说它含有丰富的内建函数,例如数学函数中的三角函数,复函数,多项式函数,数据分析函数的求平均值,最大最小值排序等,以及逻辑/选择函数如if-else等,还有用来模拟随机发生事件的随机函数。
这些我都不了解。
虽学习MATLAB的时间虽然很短,但却让我了解到了它的强大和它功能。
我想就算时间足够,老师也不能把所有的都讲解我们,因为一个软件的功能需要我们自己不断去探索,老师也不可能知道所有。
老师只是个指路的明灯,最终的学习还是要靠自己。
而且在摸索过程中,我们能够发现和体会学习的快乐!
MATLAB是一种专业的计算机程序,主要用于工程科学的矩阵数学运算。
但在之后的几年里,它又渐渐发展成为一种非常灵活的计算体系,可以解决更多技术上的问题。
在解决工程技术问题方面,matlab比其他任何计算机语言都简单高效,对于我们学电气信息工程专业的学生,这门课程给我们日后本专业的工作提供了很大的方便。
[参考文献]
[1]陈奋策,Matlab在物理实验中的应用[M],厦门,厦门大学出版社,2009。
[2]刘卫国,Matlab程序设计与应用(第二版)[M],北京,高等教育出版社,2003。
[3]杨述武,普通物理实验[M],北京,高等教育出版社,1992。
[4]马文蔚,物理学(第四版),东南大学等七所工科院校编,高等教育出版社。