When I plot f towards time amplitude is increasing I don't know why. It should be constant. what's wrong?

1 次查看（过去 30 天）

Nikoloz Nikuradze 2020-9-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/591637-when-i-plot-f-towards-time-amplitude-is-increasing-i-don-t-know-why-it-should-be-constant-what-s-w

评论： Nikoloz Nikuradze 2020-9-10

f(1)=10;    %% Ratation angle (Rad)
w(1)=0;    %% Angular velocity (Rad/sec)
E(1)=0;    %% Angular acceleration
t(1)=0;    %% Time vector (sec)
dt=0.001;
N=10000;
%% The main code
for i= 1:N
    t(i+1)=t(i)+dt;
    
    E(i)=-f(i); 
    w(i+1)=w(i)+E(i)*dt;
    f(i+1)=f(i)+w(i)*dt;
end
plot(t,f);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Matt J 2020-9-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/591637-when-i-plot-f-towards-time-amplitude-is-increasing-i-don-t-know-why-it-should-be-constant-what-s-w#answer_492544

编辑：Matt J 2020-9-10

在 MATLAB Online 中打开

I don't see why f should be constant, but nor is it monotonically increasing. Once I fix your initialization code, I see a sine-like wave:

clear
N=10000;
[f,w,E,t]=deal(zeros(1,N));
f(1)=10;
dt=0.001;
%% The main code
for i= 1:N
    t(i+1)=t(i)+dt;
    
    E(i)=-f(i); 
    w(i+1)=w(i)+E(i)*dt;
    f(i+1)=f(i)+w(i)*dt;
end
plot(t,f);