《MATLAB语言》课程论文[u1]
MATLAB在力学分析中的应用[u2]
姓名：
学号：
专业：
班级：
指导老师：
学院：
完成日期：[u3]
MATLAB在力学分析中的应用[u4]
（ 姓名 12008243764 20##级2班）[u5]
[摘要][u6] 大学物理力学中涉及许多复杂的数值计算问题，例如非线性问题，对其手工求解较为复杂，而MATLAB语言正是处理非线性问题的很好工具，既能进行数值求解，又能绘制有关曲线，非常方便实用。另外，利用其可减少工作量，节约时间，加深理解，同样可以培养应用能力 。[u7]
[关键词][u8] 力学 重力场 阻尼振动 MATLAB语言 图形绘制[u9]
一、问题的提出[u10]
MATLAB 语言是当今国际上科学界 (尤其是自动控制领域) 最具影响力、也是最有活力的软件。它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、便捷的与其他程序和语言接口的功能。MATLAB 语言在各国高校与研究单位起着重大的作用.它是一种集数值计算、符号运算、可视化建模、仿真和图形处理等多种功能……[u11]
二、重力场中小球落点问题
在物理课程的学习中我们可以明确的得到解决落体运动的方程：
[u12] （1）[u13]
例：一弹性球，初始高度 h=10m,向上初速度 v0=15m/s, 与地相碰的速度衰减系数 k=0.8,计算任意时刻球的速度和位置。
分析：用传统计算方法解决时我们需要列出传统方程，
我们明显可以感觉到，这样的计算不仅繁琐费时，而且没有图示很难给以直观的感受，现在我们用MATLAB语言来对此例题做以下解析：
MATLAB程序如下：
clear all %有衰减弹性小球运动程序
v0=15; h=10; %初速度、高度
g=-9.8; k=0.8; % 重力加速度 衰减系数
[u14] 运行结果如图1所示。
[u15]
图1 小球落地速度及位置曲线[u16]
通过以上程序对小球落地速度、位置以及运动过程的坐标描述，我们就会发现其在此类问题中直观的表述，那么现在我们来解决另外一个问题。
三、解决阻尼振动与受迫震动图像问题
1、阻尼振动方程[u17]
红线—简谐振动，蓝线
的阻尼振动，绿线
的阻尼振动，阻尼振动周期比自由振动要长，当
时，振幅按指数迅速缩减。
四、结论
从以上利用MATLAB语言对3种基本力学模型的分析我们不难的出以下结论：
五、课程体会
经过一学期紧张而有序的课程学习，在忙碌之余也得到了颇多的收获。我深深体会到MATLAB语言相对于同类程序语言更方便更简洁易懂，……
[参考文献][u18]
[1] 刘卫国.MATLAB程序设计与应用（第二版）[M].北京：高等教育出版社，2006.
[2] 马文蔚.物理学(上册)（第四版）[M],北京：高等教育出版社，1999.
说明：
1.页面为A4，页边距上下均为2.5厘米，左右均为2.2厘米。
2.行距为单倍行距。
3.页码不显示首页。
4.注意参考文献的格式。
[u1]楷体，4号，加粗，左对齐
[u2]宋体，1号，加粗，居中
[u3]楷体，4号，加粗，居中
[u4]宋体，3号，加粗，居中
[u5]楷体，5号，居中
[u6]宋体，小4号，加粗
[u7]仿宋，5号
[u8]宋体，小4号，加粗
[u9]仿宋，5号，3~5个
[u10]一级标题，宋体，4号，加粗
[u11]正文，宋体，小4号
[u12]公式用数学公式编辑，居中，编号
[u13]公式编号，新罗马，小4号
[u14]程序，从MATLAB文件窗口直接粘贴，每个语句都必须加解释
[u15]图片，居中
[u16]图要有编号，名称，宋体，5号，加粗
[u17]二级标题，宋体，小4号，加粗
[u18]宋体，小4号，加粗
第二篇：MATLAB与数值分析实验报告模板
电子科技大学 电子工程 学院
标 准 实 验 报 告
（实验）课程名称 MATLAB与数值分析
学生姓名： 代鹏
学 号： 2013020912022
指导教师：程建
一、实验名称
二、实验目的
三、实验内容
四、实验数据及结果分析：
五、总结及心得体会：
clear all
clc
x=[],y=[];
a=input('a=')
b=input('b=')
m=input('m=')
x(1)=input('x(1)=')
y(1)=input('y(1)=')
for n=1:m
x(n+1)=a*x(n)-b*(y(n)-x(n)^2);
y(n+1)=b*x(n)+a*(y(n)-x(n)^2);
end
plot(x,y)
xlabel('x'),ylabel('y')
clear all
clc
a=input('a=');
m=length(a);
for b=1:m
n=m+1-b;
for j=1:(n-1)
if a(j)>a(j+1)
t=a(j);
a(j)=a(j+1);
a(j+1)=t;
end
end
end
a
a=[6 3 7 9 1 5]
a =
1 3 5 6 7 9
clear all
clc
r=input('r=');
o1=6;
o2=6;
m=[0:0.01:2*pi];
x=[],y=[];
n=length(m);
for k=1;3
o1=o1+3/2*r*(k-1);
for i=1:n
x(i)=o1+r*cos(m(i));
y(i)=o2+r*sin(m(i));
end
plot(x,y