Can anyone pls tell me why I am getting the same plot instead of putting different values of H and B? The code is attached with this question for your kind perusal

Question

Gaurav Sharma 2018-10-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/422747-can-anyone-pls-tell-me-why-i-am-getting-the-same-plot-instead-of-putting-different-values-of-h-and-b

编辑： Walter Roberson 2019-2-2

PROB1.m

I am not getting the different plots for different valued of the H and B. Every row of the subplot has the same plot as of the first one of that row. can you Help me with this? the values I am putting are a=1,c=1,m=1,omega=5, T1=10, x0=1,v0=0.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

KSSV 2018-10-8

It has two functions RK and IEM which are not provided to check.

请先登录，再进行评论。

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

Answer 1

Gaurav Sharma 2018-10-8

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/422747-can-anyone-pls-tell-me-why-i-am-getting-the-same-plot-instead-of-putting-different-values-of-h-and-b#answer_340264

编辑：Walter Roberson 2019-2-2

在 MATLAB Online 中打开

RK:

function [OUT3,OUT4] = RK(A,M,B,C,OMEGA,T1,H,X0,V0) %RK Summary of this function goes here % Detailed explanation goes here
Y1=zeros(1,T1/H+1); Y2=zeros(1,T1/H+1); w=zeros (1,T1/H+1); Y1(1)=X0; Y2(1)=V0; w(1)=0;
    for t=1:T1/H
      K11=Y2(t);
      K12=(-A/M)* abs(Y1(t))^B* (Y1(t)/abs (Y1(t)))+ C/M* cos (OMEGA*w(t));
      K21= Y2(t)+ (H/2)* K12;
      K22= (-A/M)* abs (Y1(t)+ (H/2)* K11)^B*((Y1(t)+ (H/2)* K11)/ abs (Y1(t)+ (H/2)* K11))+ (C/M)* cos(OMEGA *w(t));
      K31= Y2(t)+ (H/2)* K22;
      K32= (-A/M)* abs (Y1(t)+ (H/2)* K21)^B*((Y1(t)+ (H/2)* K21)/ abs (Y1(t)+ (H/2)* K21))+ (C/M)* cos(OMEGA *w(t));
      K41= Y2(t)+ (H)* K32;
      K42=(-A/M)* abs (Y1(t)+ (H/2)* K31)^B*((Y1(t)+ (H/2)* K31)/ abs (Y1(t)+ (H/2)* K31))+ (C/M)* cos(OMEGA *w(t));
      Y1(t+1)= Y1(t)+H*((K11/6)+(K21/3)+(K31/3)+(K41/6));
      Y2(t+1)= Y2(t)+H*((K12/6)+(K22/3)+(K32/3)+(K42/6));
      w(t+1)=w(t)+H;
    end
OUT3= Y1; OUT4=Y2; plot(w,Y1,w,Y2);
end