Using Heun's Method but Problem With Increment

5 ビュー (過去 30 日間)
Muhammad Hadyan Utoro
Muhammad Hadyan Utoro 2021 年 8 月 28 日
Hi everyone,
I'm really new to programming and try to use Heun's Method to this skydiving problem, but I kept getting this error, "Unable to perform assignment because the left and right sides have a different number of elements. Error in EulervsHeunStart (line 17)
vh(i) = vh;"
I know I have to use the increment to get the second slope (corrector). Here's my script:
while t<=tend
%save both velocities at time ti
vh(i) = vh;
ve(i) = ve;
tim(i)= t;
%good old Euler
a = g-c/m*ve.^2;
ve = ve+a*dt;
% x = x+v*dt;
%Heun: compute slope at ti
a = g-c/m*vh.^2;
%compute predicted velocity at ti+1
vh = vh+a*dt;
%compute second slope at v(ti+1)
a = g-c/m*vh(i).^2;
vh(i+1) = vh(i)+a*dt;
%compute average of both slopes
%use average slope for extrapolating v from ti to ti+1
%increments
t=t+dt;
i=i+1;
end
Thank you!
  2 件のコメント
darova
darova 2021 年 8 月 29 日
PLease atach the whole code
Muhammad Hadyan Utoro
Muhammad Hadyan Utoro 2021 年 8 月 29 日
Hi darova,
please find the attachment for the whole code, but they're pretty much the same as I put on the question and 've' and 'vh' variables are for velocity using euler and velocity using heun respectively. I'd like to compare the error.
Many thanks

サインインしてコメントする。

採用された回答

darova
darova 2021 年 8 月 29 日
I think you should use index (i) everywhere. vh variable has two elements after first iteration
You trying to assign two elements into one position in this line
vh(i) = vh;
  1 件のコメント
Muhammad Hadyan Utoro
Muhammad Hadyan Utoro 2021 年 8 月 30 日
I've used all index on every ve and vh, but it even said, "Unable to perform assignment because the left and right sides have a different number of elements.
Error in EulervsHeunStart (line 17)
vh(i) = vh;"
I don't understand how this works, but thanks again darova. I really appreciate your help

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeProgramming についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by