شبیه ساز معادله موج و قانون هوک با روش اویلر

در نرم افزار متلب (Matlab)، میخواهیم برنامه ای نوشته شود تا معادله موج مربوط به یک فنر را که از قانون هوک پیروی می کند از روش اویلر حل کند.

clear;clc;

k=1;m=1;

w=sqrt(k/m);

h=0.1;

t=0:h:20;

n=length(t);

x=zeros(1,n);

r=zeros(1,n);

x(1)=1;r(1)=0;

for q=1:n-1

x(q+1)=x(q)+h*r(q);

r(q+1)=r(q)-h*w^2*x(q);

end

subplot(211);plot(t,x);

subplot(212);plot(r,t);

در خط اول تمام حافظه متلب (Matlab) با دستور clear و محیط خط فرمان با دستور clc پاک می شود.

از خط دوم تا چهارم، پارامترهای فیزیکی و گام تعیین می شوند.

در خط پنجم محدوده t مشخص می گردد.

در خط ششم طول بردار t با دستور length به دست آمد.

در خط هفتم و هشتم با دستور zeros بردار مکان مقدار اولیه می گیرند.

در خط نهم شرایط اولیه نوشته شده است.

از خط دهم تا سیزدهم، الگوریتم مساله با حلقه for پیاده سازی شده است.

در دو خط آخر، پاسخ مساله با دستور plot و دستور subplot رسم شده است.

معادله موج و قانون هوک روش اویلر دستور clear دستور clc دستور length دستور zeros حلقه for دستور plot دستور subplot متلب

یکشنبه 2 خرداد 1395 ساعت 08:57 ق.ظ

Peducate.ir