how i can write periodic boundary condition

2 次查看（过去 30 天）

显示更早的评论

lulu 2020-1-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/502970-how-i-can-write-periodic-boundary-condition

评论： Mahesh Taparia 2020-2-3

clear

clc

xmin=0;

xmax=1;

N=10;

dt=0.001;

t=0;

tmax=0.5;

a=0.1;

L=2;

dx=(xmax-xmin)/N;

x=xmin-dx:dx:xmax+dx;

u0=exp(-(x-L/2).^2);

u=u0;

unp1=u0;

v0=exp(-(x-L/2).^2);

v=v0;

vnp1=v0;

nstep=tmax/dt;

for n=1:nstep

u(1)=u(3);

u(3+N)=u(1+N);

v(3)=v(1);

v(N+1)=v(N+3);

for i=2:N+2

unp1(i)= -a*u(i-1);

vnp1(i)=a*v(i-1);

end

t=t+dt;

u=unp1;

v=vnp1;

end

plot(u,'b');

hold on

plot(v);

hold off

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

Mahesh Taparia 2020-2-3

Can you give more details what you want to do?

请先登录，再进行评论。

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

回答（0 个）

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

类别

MATLAB Mathematics Graph and Network Algorithms Construction Directed Graphs

在 Help Center 和 File Exchange 中查找有关 Directed Graphs 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by