... tr1xyTaxâydựng hàm wave2D()đểthựchiệnthuậttoántrên:function[u,x,y,t]=wave2D(a,D,T,it0,i1t0,bxyt,Mx,My,N)%giaia(u_xx+u_yy)=u_ttvoiD(1)<=x<=D(2),D(3)<=y<=D(4),0<=t%<=T%dieukiendau:u(x,y,0)=it0(x,y),u_t(x,y,0)=i1t0(x,y)%dieukienbien:u(x,y,t)=bxyt(x,y,t)voi(x,y)trenbien%Mx/My‐sokhoangchiatrentrucx/y%N‐sokhoangchaitheotdx=(D(2)‐D(1))/Mx;x=D(1)+[0:Mx]*dx;dy=(D(4)‐D(3))/My;y=D(3)+[0:My]ʹ*dy;dt=T/N;t=[0:N]*dt;%khoiganu=zeros(My+1,Mx+1);ut=zeros(My+1,Mx+1);forj=2:Mxfori=2:Myu(i,j)=it0(x(j),y(i));ut(i,j)=i1t0(x(j),y(i));endendadt2=a*dt*dt;rx=adt2/(dx*dx);ry=adt2/(dy*dy); 403CHƯƠNG 9: PHƯƠNG TRÌNH VI PHÂN ĐẠO HÀM RIÊNG §1.KHÁINIỆMCHUNG Phươngtrìnhviphân đạo hàm riêng( PDE)làmộtlớpcácphươngtrìnhviphâncósốbiếnđộclậplớnhơn1.Trongchươngnàytasẽkhảosátcácphươngtrìnhviphân đạo hàm ... t)vàđiềukiệnđầuu(x,y,0)=i0(x,y)Tathay đạo hàm bậc1theotởvếphảibằngsaiphân3điểmtạiđiểmgiữa(tk+1+tk)/2nhưphươngphápCrank‐Nicholson.Tacũngthaythếmộttrongcác đạo hàm bậchaiuxxvàuyybằngxấpxỉ3điểmtạithờiđiểmtkvà đạo hàm kiatạitk+1vàcó:++−+−⎛⎞−+ ... hết ta tạo các hàm cơ sở bằng cách dùng hàm fembasisftn() và vẽ một hàm trong số đó bằng cách dùng lệnh MATLABmesh()nhưhìnha.Thứhai,khôngtạora hàm cơsở,tadùnglệnhMATLABtrimesh()đểvẽcác hàm hìnhdạngchocácnútn=2,3,4và5...