برنامه متلب حل معادله موج یک بعدی درجه یک به روش اختلاف بالادست

برنامه فوق برای حل معادله موج یک بعدی درجه یک به روش اختلاف بالادست می باشد. این برنامه شکل موج را به صورت انیمیشن ترسیم می کند و دارای شرط مرزی دریکله در ابتدا و انتها می باشد. شرط اولیه را تابع sin(x) در نظر گرفتم. این برنامه فاقد گزارش کار می باشد.

clear all

close all

clc

dx=(pi/12);% interval

dt=0.1; % time interval

c=1; %velocity of wave

v=(c*dt/dx);% noo number

.

.

.

for i=1:Jx

u1(i)=sin(x(i));

end

.

.

.