CP080_计算物理常微分方程解.ppt

下载文档

12
0
约 46页
2017-04-21 发布于上海
举报
版权申诉
保障服务

CP080_计算物理常微分方程解.ppt

1、本文档共46页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

CP080_计算物理常微分方程解

常微分方程的数值解法;绪论;常微分方程描写的物理现象;常微分方程;常微分方程数值解基本思想;8.1 Euler方法;8.1 Euler方法;8.1 Euler方法-梯形公式;8.2 改进Euler方法;8.2 改进Euler方法;1;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;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(欧拉法的曲线,改进欧拉法的曲线, 理论曲线) function y=f(x) global E R C; y=-1*x./(R*C);;K;Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd.;function main %欧拉法：二阶常微分方程 global E R C L; E=10;R=100;C=0.01;L=10; Q0=0;I0=0;h=0.01;t=[0:h:10]; Q(1)=Q0; I(1)=I0; for i=2:length(t) Q(i)=Q(i-1)+h*I(i-1);%Q I(i)=I(i-1)+h*f(Q(i-1),I(i-1));%I end plot(t,Q,r,t,I,m);%欧拉法的曲线 function y=f(Q,I) global E R C L; y=(E-Q/C-I*R)/L;;（用改进的Euler法解）：;function main %改进欧拉法：二阶常微分方程 global E R C L; E=10;R=100;C=0.01;L=10; Q0=0;I0=0;h=0.01;t=[0:h:10];Q(1)=Q0;I(1)=I0; for i=2:length(t) %1、预报 Q(i)=Q(i-1)+h*I(i-1); I(i)=I(i-1)+h*f(Q(i-1),I(i-1)); %2、计算Q k1=h*I(i-1); k2=h*I(i); Q(i)=Q(i-1)+1/2*(k1+k2); %3、计算I k1=h*f(Q(i-1),I(i-1)); k2=h*f(Q(i),I(i)); I(i)=I(i-1)+1/2*(k1+k2); end plot(t,Q,r,t,I,m);%改进欧拉法的曲线 function y=f(Q,I) global E R C L; y=(E-Q/C-I*R)/L;;8.3 龙格-库塔（R-K）方法;例：用R-K方法解例题8.1.1;function main %龙格库塔：一阶常微分方程 global E R C; E=10;R=10;C=0.01;h=0.01;t=[0:h:1];Q(1)=E/C; for i=2:length(t) k1=h*f(Q(i-1)); k2=h*f(Q(i-1)+k1/2); k3=h*f(Q(i-1)+k2/2); k4=h*f(Q(i-1)+k3); Q(i)=Q(i-1)+1/6*(k1+2*k2+2*k3+k4); end plot(t,Q,‘b’);hold on%龙格库塔方法的曲线 plot(t,Q(1)*exp(-1*t./(R*C)),.);%理论曲线 function y=f(x) global E R C; y=-1*x./(R*C);;例：用R-K方法解例题8.1.2;function main %龙格库塔：二阶常微分方程 global E R C L; E=10;R=100;C=0.01;L=10; Q0=0;I0=0;h=0.01;t=[0:h