预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共30页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

数学建模,数学实验与Matlab的应用实例特拉法尔加海战TheBattleofTrafalgar1805年,英国海军上将纳尔逊勋爵指挥英国海军舰队赢得特拉法尔加海战的胜利。为纪念这场经典海战的200周年,来自35个国家的160余艘战舰和船只聚集在英格兰南部海域,重演当年海战。在2005年6月28日的庆典中,英女王伊丽莎白二世登上英国皇家海军船只,完成这次世界最大规模的海军舰队大检阅。纪念特拉法尔加海战英举行最大规模海军检阅 http://jczs.sina.com.cn2005年06月29日 特拉法尔加海战200周年的纪念系列活动,由英国皇家海军和英国政府共同举办。作为2005英国海洋年的先导活动,它将一直持续到明年10月。 拿破仑建立法兰西第一帝国后,为夺取英吉利海峡的制海权,命海军将领维尔纳夫与罗什福尔、布勒斯特的分舰队会师后远征英国本土。然而由于途中误会,维尔纳夫分舰队南下驶往西班牙南部的加的斯港。 英国海军名将纳尔逊得知后,乘“胜利”号旗舰从朴次茅斯出航,与英海军科林伍德分舰队会合。1805年10月20日,包括5艘西班牙主力舰在内的法西联合舰队(主力舰33艘、驱逐舰5艘,装备火炮300门,兵员3万人)从加的斯港出航,准备转移到地中海。次日晨,当法西舰队行进到特拉法尔加角附近时,被英国舰队(主力舰27艘,装备火炮2500门,兵员2万人)发现,于是双方展开大战。 英国海军最终取得了海战的胜利。法西舰队18艘舰船被击沉或俘获,7000人被俘,2600人死伤。而英军统帅纳尔逊被法国狙击手击中,伤重不治而殉国,英国人把他视为民族英雄,他的旗舰“胜利”号也被当作重要的历史文物而保存了下来。这次海战确立了英国在海上的霸权地位,而纳尔逊把帆船时代特别壮观的海战线式战术发展到最高峰,成为帆船时代海战的典范。167艘战舰抵英,经典海战重演.乐队在游行中英国演员装扮纳尔逊在座舰胜利号前重现当年海战的油画例特拉法尔加海战(TheBattleofTrafalgar) 1805年,纳尔逊海军上将率领的英国舰队与法国-西班牙联合舰队在特拉法尔加角进行了一场海战.当时法国-西班牙联合舰队有战舰F0=33艘,英国舰队有战舰B0=27艘.双方的战斗力与各自的战舰数成比例.影响战斗力的三个因素: 假设在每一个回合中,各方损失的舰只数是对方舰只数的10%,则有数学模型 Bk+1=Bk0.1Fk, Fk+1=Fk0.1Bk. 这正是前面已讨论过的正规战争模型.记B为,F为,,则有 用Matlab进行实验: X=[27,33]’; A=[1,-0.1;-0.1,1]; k=0; k=k+1,X=A*X 每运行一次X=A*X代表进行一次战斗,k记录了战斗次数.最后一行可用“”反复调用.如果双方都在战斗中投入全部兵力,X0=(27,33)T,则经过10个回合后, X10=A10X0=(2.6791,18.2416)T, 结果将是英国舰队仅剩3艘战舰,且至少有一艘遭到重创,而法国-西班牙舰队还剩18艘战舰.显然是英方战败.但是实际的结局是战斗以英方的胜利告终.纳尔逊决定采用各个击破(Divide-and-Conquer)的策略:先用13艘战舰去攻击对方的A组3艘战舰,留14艘作预备队.然后用战斗剩余的战舰与预备队合并去攻击对方的B组,最后用剩余的力量与对方剩余的力量决战.设每一个回合中,各方损失的舰只数是对方舰只数的5%.用数学模型作数学实验如下:从实验的结果来看,英方将获胜.事实上,纳尔逊率领的英国舰队确实获得了胜利,法西联合舰队在第三轮战斗中还剩13艘战舰时战败而去,但是纳尔逊本人也在战斗中阵亡.讨论题: 如果海战的双方分别是水面舰队和潜艇部队,那么应该用那一种战争模型?请自设数据和合理的参数进行数学实验并分析结果.设X为参战一方(水面舰队)的战舰数,Y为另一方(潜艇部队)的潜艇数,则可建立如下数学模型: Xk+1=XkYk, Yk+1=XkXkYk.例一只猎犬发现其正东方100米处有一只野兔,野兔以速度v向其正北方100米处的洞穴逃跑,猎犬向野兔追去,速度是2v.求:1.猎犬的运动轨迹;2.追上兔子的时间、地点. 模型1:消去s解得模型2:模型1的解: 求解二阶常微分方程用M文件ff26.m来定义此方程组,内容如下: functionz=ff26(x,y) z(1)=y(2); z(2)=sqrt(1+y(2).^2)./(2*(100-x)); z=[z(1)z(2)]';27模型2的解: 求解一阶微分方程组:解M文件ff26_1.m: functionz=ff26_1(t,y) v=1; z(1)=2*v*(100-y(1))/sqrt((v*t-y(2))^2+… (100-y(1))^2); z(2)=2*v*(v*t-y(2)