CP080-计算物理常微分方程解.pptVIP

常微分方程的数值解法 绪论 在工程和科学计算中，所建立的各种常微分方程的初值或边值问题，除很少几类的特殊方程能给出解析解，绝大多数的方程是很难甚至不可能给出解析解的，其主要原因在于积分工具的局限性。因此，人们转向用数值方法去解常微分方程，并获得相当大的成功，讨论和研究常微分方程的数值解法是有重要意义的。 常微分方程描写的物理现象 镭的衰变规律 单摆的运动 RLC振荡电路 物理场计算 常微分方程 常微分方程数值解基本思想 8.1 Euler方法 8.1 Euler方法 function main global E R C; E=10;R=10;C=0.01; Q0=E/C; h=0.01; t=[0:h:1]; QE(1)=Q0; QEG(1)=Q0; for i=2:length(t) QE(i)=QE(i-1)+h*f(QE(i-1));%欧拉法 k1=h*f(QEG(i-1));%改进欧拉法 k2=h*f(QEG(i-1)+k1); QEG(i)=QEG(i-1)+1/2*(k1+k2); end plot(t,QE,r);%欧拉法的曲线 hold on plot(t,QEG,b);%改进欧拉法的曲线 plot(t,Q0*exp(-1*t./(R*C)),‘.’);%理论曲线 legend(欧拉法的曲线,改进欧拉法的曲线