matlab作业 第13-14题

《matlab作业 第13-14题》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、13、计算这一年已经逝去的天数（包括当天）。在平年中，它的取值范围为1到365，在闰年中，它的取值范围1到366。编写一个MATLAB程序，输入年，月，日，输出为对应的该天是当年的第几天。闰年是这样规定的1）.能被400整除的年为闰年2）能被100整除但不能被400整除的年不为闰年3）.能被4整除但不能被100整除年为闰年4）.其余的年份不为闰年代码：clearall;closeall;year=input('请输入年：');month=input('请输入月：');day=input('请输入日：');date=[312831303130313130313031];sum=0;

2、forn=1:month-1sum=sum+date(n);endsum=sum+day;ifmonth>2ifmod(year,400)==0

3、mod(year,4)==0&mod(year,100)~=0sum=sum+1;endendfprintf('该天是当年的第%d几天。',sum);运行结果：请输入年：2000请输入月：3请输入日：6该天是当年的第66几天。请输入年：2001请输入月：3请输入日：6该天是当年的第65几天。请输入年：2008请输入月：3请输入日：6该天是当年的第66几天。14、画出小球的轨迹如果我们假设处于真空中，且忽略地球的曲率。在地球任意一点向空中

4、抛出一小球将会产生类似于图（a）所示的曲线。球在时刻t的高度将会遵守公式。yyvvyxθvxxorigin(a)impact(b)(a)如果物体上抛，将会产生一个抛特线轨迹。(b)速度可以分解为水平速度和竖直速度，水平速度和合速度之间的夹角为θ。12y(t)=y+vt+gt0yo2其中y0是初始高度，vy0代表初速度，g代表重力加速度。水平位移的计算公式为x(t)=x0+vx0tx0代表初始位移，vx0代表这个物体的水平初速度，竖直初速度与合初速度之间的关系为vx0=v0cosθvy0=v0sinθ假设一个小球的初始位置为(x0，y0)为(0，0)，初速度为20m/s，水平速度和

5、合速度之间的夹角为θ度，编写一个程序，1）画出这个小球的轨迹，并计算小球再次落地与初始位置之间的距离。这程序应当能画出多个抛物线，θ的取值从5到85度，每隔10度取一次，计算水平位移..2）θ的取值从0到90度，每隔1度取一次。最终应当确定那一个θ值使得水平位移最大。还有打印不同的抛物线时要用不同的颜色。地球的重力加速度为9.81m/s2。（1）代码：clearall;closeall;g=9.81;a=5:10:85;t=0:0.01:4.1;Vx=20.*cos(pi*a/180);Vy=20.*sin(pi*a/180);holdon;axis([045025]);fori

6、=1:1:9x=Vx(i).*t;y=Vy(i).*t-0.5.*g.*(t.^2);ifmod(i,5)==0plot(x,y,'r');elseifmod(i,5)==1plot(x,y,'g');elseifmod(i,5)==2plot(x,y,'b');elseifmod(i,5)==3plot(x,y,'y');elseifmod(i,5)==4plot(x,y,'k');endd(i)=2*Vy(i)*Vx(i)/g;fprintf('第%d条抛物线的水平位移为：%.3f米。',i,d(i));end运行结果：25第1条抛物线的水平位移为：7.080米。20第2

7、条抛物线的水平位移为：20.387米。第3条抛物线的水平位移为：31.235米。第4条抛物线的水平位移为：38.316米。15第5条抛物线的水平位移为：40.775米。第6条抛物线的水平位移为：38.316米。10Y/m第7条抛物线的水平位移为：31.235米。第8条抛物线的水平位移为：20.387米。5第9条抛物线的水平位移为：7.080米。0051015202530354045X/m（2）代码：clearall;closeall;g=9.81;a=0:1:90;t=0:0.01:4.1;Vx=20.*cos(pi*a/180);Vy=20.*sin(pi*a/180);hol

8、don;axis([045025]);fori=1:1:91x=Vx(i).*t;y=Vy(i).*t-0.5.*g.*(t.^2);ifmod(i,5)==0plot(x,y,'r');elseifmod(i,5)==1plot(x,y,'g');elseifmod(i,5)==2plot(x,y,'b');elseifmod(i,5)==3plot(x,y,'y');elseifmod(i,5)==4plot(x,y,'k');endd(i)=2*Vy(i)*Vx(i)/g;