基于matlab的马格努斯力的球类三维空气动力学仿真

function air_resistance global A; %matlab自动将a作为数组处理 A=[]; m=input('m=');%数组的个数 c=colormap(jet(m)); %曲线颜色 n=input('n=');%数组的长度 for i=1:m for j=1:n fprintf('a%.0f=',j); x=input('');%分别输入各个数的值 A=[A,x]; end [t,x]=ode45(@f,[0:0.01:15],[A(1),A(2),A(3),A(4)]); plot(x(:,1),x(:,3),'color',c(i,:), 'LineWidth',4); xlim([0 160]) % sets the x limits ylim([0 160]) % sets the y limits xlabel({'X方向位移'},'FontSize',12,'FontWeight','bold'); ylabel({'Y方向位移'},'FontSize',12,'FontWeight','bold'); grid on; hold on; A=[]; end function dx=f(t,x) dx = zeros(4,1); dx(1)=x(2); dx(3)=x(4);