《matlab软件与基础数学实验》程序(部分)

《《matlab软件与基础数学实验》程序(部分)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、追击问题：constructionqualityacceptanceandassessmentRegulation(ProfessionalEdition)(DL/T5210.2-2009~DL/T5210.8-2009);1.9thequalitycheckoutandevaluationofelectricequipmentinstallationengineeringcode(DL/T5161.1-2002~5161.17-2002);1.10thenormsofconstructionsupervision,theelectricpowerconstru

2、ctionsupervisionregulations21constructionqualityacceptanceandassessmentRegulation(ProfessionalEdition)(DL/T5210.2-2009~DL/T5210.8-2009);1.9thequalitycheckoutandevaluationofelectricequipmentinstallationengineeringcode(DL/T5161.1-2002~5161.17-2002);1.10thenormsofconstructionsupervision

3、,theelectricpowerconstructionsupervisionregulations21constructionqualityacceptanceandassessmentRegulation(ProfessionalEdition)(DL/T5210.2-2009~DL/T5210.8-2009);1.9thequalitycheckoutandevaluationofelectricequipmentinstallationengineeringcode(DL/T5161.1-2002~5161.17-2002);1.10thenormso

4、fconstructionsupervision,theelectricpowerconstructionsupervisionregulations21一敌舰在某海域内以椭圆轨迹航行，其在时间t时刻的坐标为：x(t)=10+20costy(t)=20+5sint我方战舰恰位于原点处，我战舰向敌舰发射制导鱼雷，鱼雷的速率为20，其运行方向始终指向敌舰，试问敌舰航行在何处将被击中？2.若敌舰的运行轨迹变为x(t)=10+20costy(t)=20+20sint试问敌舰航行在何处将被击中？(无法击中)3.若敌舰的运行轨迹变为x(t)=10+20costy(t)=20

5、+20sint鱼雷速率提高至21，结果如何？%Matlab程序：clear;clch=0.01;%时间步长k=1;t(1)=0;x(1)=0;y(1)=0;%初始值r=10;whiler>=0.05%k<=250%m=(20+5*sin(t(k))-y(k))/(10+20*cos(t(k))-x(k)+1.e-10)+1.e-10;constructionqualityacceptanceandassessmentRegulation(ProfessionalEdition)(DL/T5210.2-2009~DL/T5210.8-2009);1.9thequa

6、litycheckoutandevaluationofelectricequipmentinstallationengineeringcode(DL/T5161.1-2002~5161.17-2002);1.10thenormsofconstructionsupervision,theelectricpowerconstructionsupervisionregulations21if10+20*cos(t(k))-x(k)>=0x(k+1)=x(k)+20*h/sqrt(1+m^2);elsex(k+1)=x(k)-20*h/sqrt(1+m^2);endif

7、20+5*sin(t(k))-y(k)>=0y(k+1)=y(k)+20*h/sqrt(1+1/m/m);elsey(k+1)=y(k)-20*h/sqrt(1+1/m/m);endr=(x(k)-10-20*cos(t(k)))^2+(y(k)-20-5*sin(t(k)))^2;r=sqrt(r);t(k+1)=h*k;k=k+1;plot(10+20*cos(t(k)),20+5*sin(t(k)),'r*')holdonaxis([-1032-330]);plot(x,y,'o')pause(0.02)endt=t(end),x=x(end),y=y(e

8、nd)construct