Ngày tải lên :
23/01/2014, 06:20
... (1)
vớiđiềukiệnbên:
403
CHƯƠNG 9: PHƯƠNG TRÌNH VI PHÂN ĐẠO HÀM RIÊNG
§1.KHÁINIỆMCHUNG
Phương trình vi phân đạo hàm riêng( PDE)làmộtlớpcác phương trình
vi phân cósốbiếnđộclậplớnhơn1.Trongchươngnàytasẽkhảosátcác
phương trình vi
phân đạo hàm ... 22
qq
q
qq
⎡
⎤
=
⎢
⎥
⎣
⎦
1
2
g
g
g
2.SửdụngPDETOOL: PDETOOLgiải phương trình vi phân đạo hàm riêng
bằng cách dùng phương phápFEM.Đểgiải phương trình tatheocácbước
sau:
) Nhập lệnh pdetool ...
432
Mx=16;
dx=(xf‐x0)/Mx;
xi=x0+[0:Mx]*dx;
My=16;
dy=(yf‐y0)/My;
yi=y0+[0:My]*dy;
fori=1:length(xi)
forj=1:length(yi)
fors=1:Ns
ifinpolygon(xi(i),yi(j),N(S(s,:),1),N(S(s,:),2))>0
Z(i,j)=U(s,:)*[1xi(i)
yi(j)]ʹ;%Pt.(4.5b)
break;
end
end
end
end
figure(2);
clf;
mesh(xi,yi,Z)
%desosanh
bx0=inline(ʹ0ʹ);
bxf=inline(ʹ0ʹ);
by0=inline(ʹ0ʹ);
byf=inline(ʹ0ʹ);
D=[x0xfy0yf];
[U,x,y]=poisson(f,g,
bx0,bxf,by0,byf,D,Mx,My,1e‐6,50);
figure(3)
clf;
mesh(x,y,U)
§6.GUICỦAMATLABĐỂGIẢIPDE
1.Các phương trình cóthểgiảiđượcbằngPDETOOL:CôngcụPDETOOL
củaMATLABcóthểdùngđểgiảicác loại phương trình sau:
a. Phương trình elliptic:Tasẽgiải phương trình elliptic
...