【运动学】基于matlab飞机往返运动（相对运动速度.zip

%飞机往返的时间 clear %清除变量 u=0:0.05:0.8; %空气速率与飞机速率之比的向量 theta=0:30:90; %角度向量 th=theta*pi/180; %化为弧度 [U,TH]=meshgrid(u,th); %速率和角度矩阵 T1=1/2./(sqrt(1-(U.*sin(TH)).^2)+U.*cos(TH)); T2=1/2./(sqrt(1-(U.*sin(TH)).^2)-U.*cos(TH)); T=sqrt(1-(U.*sin(TH)).^2)./(1-U.^2); %时间矩阵 figure %创建图形窗口 %plot(u,T,'LineWidth',2) %画往返时间曲线族 plot(u,T(1,:),'o-',u,T(2,:),'d-',u,T(3,:),'s-',u,T(4,:),'v-','LineWidth',2)%画往返时间曲线族 hold on %保持图像 % plot(u,T1,'--','LineWidth',2) %画A到B的时间曲线族 % plot(u,T2,':','LineWidth',2) %画B到A的时间曲线族 plot(u,T1(1,:),'o--',u,T1(2,:),'d--',u,T1(3,:),'s--',u,T1(4,:),'v--','LineWidth',2)%画A到B的时间曲线族 plot(u,T2(1,:),'o:',u,T2(2,:),'d:',u,T2(3,:),'s:',u,T2(4,:),'v:','LineWidth',2)%画B到A的时间曲线族 grid on %加网格 fs=16; %字体大小 title('飞机往返的时间与风速的关系','FontSize',fs)%标题 xlabel('\itu/v','FontSize',fs) %标记横坐标 ylabel('\itt/t\rm_0','FontSize',fs) %标记纵坐标 n=length(theta); %角度个数 legend([repmat('\it\theta\rm=',n,1),num2str(theta'),repmat('\circ',n,1)],2)%图例 text(0,1.2,'\itt\rm_0=2\itl/v','FontSize',fs)%标记时间单位 text(0,2,'虚线是A到B的时间,点虚线是B到A的时间','FontSize',fs)%说明单向时间