数学建模实验报告

发布 2023-05-17 21:14:28 阅读 6730

成员。专业:12级信息与计算科学。

学号。指导老师。

试用摆角的角加速度的三种表达式, 即式(1.11)、(1.23)和(1.24),取步长为,,,的值如前,计算当变化时角加速度的值,并列表加以比较。

实验程序:function m1_2(t)

r=100;l=300;w=240/60*2*pi;

a0=-r*w^2*sin(t)*(l^2-r^2)./l^2-r^2*sin(t).^2).^3/2)

a1=-w^2*r*sin(t)/l

a2=-w^2*(r*sin(t)/l+r^3*(sin(t).^3-sin(2*t).*cos(t))/2*l^3))

> m1_2([0:pi/12:pi])

运行结果如下:

a0 =columns 1 through 9

columns 10 through 13

a1 =columns 1 through 9

columns 10 through 13

a2 =columns 1 through 9

columns 10 through 13

从结果中可以看出误差的大小,取决于近似表达式的精度,在利用泰勒公式求近似模型时,如果展开的精度越高,则误差就越小,在数据表中也可以看出, 取得精度比高,所以结果与真实值相差的更小。

实验任务二:

应用数学软件matlab对问题(4.12)~(4.14)进行数值计算,先运用euler法,与表4.

2以及表4.3的数据比较,并以更小的步长计算结果。再用改进euler法计算(步长与euler法步长相同)。

euler方法:

实验程序:function m_1(n)

h=120;

h=h/n;

lamda=90/450;

x(1)=0;p(1)=0;

y=0:h:h;

for i=0:n-1

x(i+2)=x(i+1)+h*p(i+1);

p(i+2)=p(i+1)+h*(lamda*sqrt(1+p(i+1)^2)/(h-y(i+1)))

endx;p]’

l=x(n+1)

t=x(n+1)/90

>m_1(4)ans =

l =t =书中表(4.1)

实验程序:function m_2(n)

k=1;for n=n

h=120;

h=h/n;

lamda=90/450;

x0=0;p0=0;

for i=0:n-1

x1=x0+h*p0;

p1=p0+h*(lamda*sqrt(1+p0^2)/(h-i*h));

x0=x1;

p0=p1;

endl(k)=x1;

t(k)=x1/90;

k=k+1;

endn;l;t]'

>m_2([4,8,12,24,48,96,120,240])ans =

书中表(4.2)

实验程序:function m_3(t)

h=120;ve=90;vw=450;

x(1)=0;y(1)=0;t(1)=0;

for i=1:10e6

m=(ve*t(i)-x(i))/h-y(i));

x(i+1)=x(i)+vw*t/sqrt(1+1/m.^2);

y(i+1)=y(i)+vw*t/sqrt(1+m.^2);

t(i+1)=t+t(i);

if y(i+1)>=h

break;

endend

t;x;y]'

l=x(i+1)

t=x(i+1)/ve

>m_3(0.1)

in m_3 at 6ans =

l =t =书中表(4.3)

>m_3(0.05)

warning: divide by zero.

in m_3 at 6ans =

l =t =书中表(4.4)

改进的euler方法:

用改进的euler方法取步长为0.1和0.05时(与euler 法步长相同)计算。

function m_5(t)

h=120;ve=90;vw=450;

x(1)=0;y(1)=0;t(1)=0;

for i=1:10e6

m=(ve*t(i)-x(i))/h-y(i));

x1(i+1)=x(i)+vw*t/sqrt(1+1/m.^2);

y1(i+1)=y(i)+vw*t/sqrt(1+m.^2);

t(i+1)=i*t;

x(i+1)=0.5*(x1(i+1)+x(i)+vw*t/sqrt(1+((h-y1(i+1))/ve*t(i+1)-x1(i+1)))2));

y(i+1)=0.5*(y1(i+1)+y(i)+vw*t/sqrt(1+((ve*t(i+1)-x1(i+1))/h-y1(i+1)))2));

if y(i+1)>=h

break;

endend

t;x;y]'

l=x(i+1)

t=x(i+1)/ve

> m_5(0.1)

warning: divide by zero.

in m4_5 at 6ans =

l =t =

数学建模 实验报告

解 根据题意我们可以建立目标函数和约束条件。目标函数 约束条件 均大于0且为整数。所以周一开始工作到从周日开始工作的雇员数分别为 2 人 2 人 4 人 3 人 3 人 0 人 8 人时,总聘用费用最少为 11000 元。直接将区间范围改为 5,5 时,结果如下 最小值点取在x 1.5715处,然而...

数学建模实验报告

桂林电子科技大学2017 2018学年第1学期。数学建模。一 实验目的。学会使用matlab软件解决线性规划问题的最优值求解问题。学会使用matlab软件解决非线性规划问题的最优值求解问题。学会将实际问题归结为线性规划问题或者非线性规划问题用matlab软件建立巧当的数学模型来求解。二 实验内容。用...

数学建模实验报告

湖南科技学院实验报告。实验目的。1 学会操作matlab。2 应用软件建立层次分析模型,选择最佳方案。实验内容 算法 程序 步骤和方法 输入 输出 实验结果 实验结果。分析 问题 某企业由于生产效益较好,年底取得一笔利润,领导决定拿出一部分资金分别用于 1 为企业员工发年终奖金 2 扩建集体福利设施...