MATLAB问题 已知速度v=1，和速度方向变化函数n1（如下），求该运动的轨迹

MATLAB问题 已知速度v=1，和速度方向变化函数n1（如下），求该运动的轨迹

clc
clear
close all
R_L=[0.2 0.4 0.8 1.2];
n=10;%rpm
omega=2*pi*n/60;
t=0:0.01:24;
for j=1:length(R_L)
        s=0;
    for i=1:length(t),
        n1(j,i)=60/2/pi*R_L(j)*omega*cos(omega*t(i))/(1+R_L(j)^2*sin(omega*t(i))^2);
        s0=n1(j,i)*0.01;
        if i==1;
            S(j,i)=s0;
        else
            S(j,i)=S(j,i-1)+s0;
        end
    end
end
subplot(221)
[Ax,h1,h2]=plotyy(t,n1(1,:),t,S(1,:));
set(get(Ax(2),'Ylabel'),'string','routing S')
xlabel('Time [s]');ylabel('Speed n1 [rpm]')
title('角速度分析图 R/L=0.2,n=10rpm ')
subplot(222)
[Ax,h1,h2]=plotyy(t,n1(2,:),t,S(2,:));
set(get(Ax(2),'Ylabel'),'string','routing S')
xlabel('Time [s]');ylabel('Speed n1 [rpm]')
title('角速度分析图 R/L=0.4,n=10rpm ')
subplot(223)
[Ax,h1,h2]=plotyy(t,n1(3,:),t,S(3,:));
set(get(Ax(2),'Ylabel'),'string','routing S')
xlabel('Time [s]');ylabel('Speed n1 [rpm]')
title('角速度分析图 R/L=0.8,n=10rpm ')
subplot(224)
[Ax,h1,h2]=plotyy(t,n1(4,:),t,S(4,:));
set(get(Ax(2),'Ylabel'),'string','routing S')
xlabel('Time [s]');ylabel('Speed n1 [rpm]')
title('角速度分析图 R/L=1 2,n=10rpm ')

追问求的不是改变R_L的参数，是轨迹，角速度是n1，速度为1，的轨迹追答不好意思哈，程序有误，已更正，详见附件

如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！

点此我要举报以上问答信息