偏微分方程数值解某实验报告材料.docxVIP

偏微分方程数值解 上 机 实 验 报 告 （一）实验一 上机题目： 用线性元求解下列边值问题的数值解： -y+π24y=π22sinπ2x，0x1 y(0)=0，y1=0 实验程序： function S=bz x=fzero(@zfun,1); [t y]=ode45(@odefun,[0 1],[0 x]); S.t=t; S.y=y; plot(t,y) xlabel(x:′从0一直到1) ylabel(y) title(线性元求解边值问题的数值解) function dy=odefun(x,y) dy=[0 0]; dy(1)=y(2); dy(2)=(pi^2)/4*y(1)-((pi^2)/2)*sin(x*pi/2); function z=zfun(x); [t y]=ode45(@odefun,[0 1],[0 x]); z=y(end)-0; 实验结果： 1．以步长h=0.05进行逐步运算，运行上面matlab程序结果如下： 2．在0x1区间上，随着x的不断变化，x，y之间关系如下图所示： （二）实验二 上机题目： 求解Helmholtz方程的边值问题： ，于 ，于 其中k=1,5,10,15,20 实验程序： （采用有限元方法，这里对[0,1]*[0,1]采用n