matlab计算抛物型偏微分方程的数值解(1)functionuxy=Ini2U(x,y)formatlong;uxy=exp(-10*x*x-10*y*y);functionu=peParabImp(c,dt,n,minx,maxx,lbu,rbu,M)formatlong;maxx=100;minx=1;n=10h=(maxx-minx)/(n-1);%u0(1)=lbu;u0(1)=0;%u0(n)=rbu;u0(n)=0;forj=2:n-1u0(j)=PrIniU(minx(j-1)*h);endu1=u0;M=100dt=0.01c=0.1lbu=...
functionuxy=Ini2U(x,y)formatlong;uxy=exp(-10*x*x-10*y*y);functionu=peParabImp(c,dt,n,minx,maxx,lbu,rbu,M)formatlong;maxx=100;minx=1;n=10h=(maxx-minx)/(n-1);%u0(1)=lbu;u0(1)=0;%u0(n)=rbu;u0(n)=0;forj=2:n-1u0(j)=PrIniU(minx(j-1)*h);endu1=u0;M=100dt=0.01c=0.1lbu=0rbu=0fork=1:MA=zeros(n-2,n-2);cb= -transpose(u0(2:(n-1)));cb(1)=cb(1)-dt*c*lbu/h/h;cb(n-2)=cb(n-2)-dt*c*rbu/h/h;A(1,1)=-2*dt*c/h/h-1;A(1,2)=dt*c/h/h;fori=2:n-3A(i,i-1)=dt*c/h/h;A(i,i)=-2*dt*c/h/h-1;A(i,i1)=dt*c/h/h ;endA(n-2,n-2)=-2*dt*c/h/h-1;A(n-2,n-3)=dt*c/h/h;u1(2:(n-1))=A\cb;u0=u1;endu=u1;formatshort;
本文档为【matlab计算抛物型偏微分方程的数值解(1)】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483