实例介绍
【实例简介】matlab偏微分方程的求解代码 含流体力学计算模型 一维波动差分 一维热传导有限元等
matlab偏微分方程的求解代码
并包含一些运行样例,真题运行结果与分析-2013微分方程数值解编程试题
有限元方法的程序
有限差分方法的程序
二维空间的偏微分方程可以用PDE函数组求解,该组是用有限元方法实现的,也可以不用编程,使用MATLAB的交互程序PDETOOL求解.
三维空间的偏微分方程求解可以使用ANSYS,包括一些流体力学计算模型.
【实例截图】
【核心代码】
function ex2p2
clear all;
M=100;N=256;
a=0;b=1;T0=0;T1=2*pi;
h=(b-a)/M;
k=(T1-T0)/N;
r=k/h^2;
x=a:h:b;
A=-r/2*sparse(2:M-1,1:M-2,ones(1,M-2),M-1,M-1) (1 r)*sparse(1:M-1,1:M-1,ones(1,M-1),M-1,M-1)-r/2*sparse(1:M-2,2:M-1,ones(1,M-2),M-1,M-1);
B=r/2*sparse(2:M-1,1:M-2,ones(1,M-2),M-1,M-1) (1-r)*sparse(1:M-1,1:M-1,ones(1,M-1),M-1,M-1) r/2*sparse(1:M-2,2:M-1,ones(1,M-2),M-1,M-1);
K=inv(A)*B;
u=zeros(M 1,N 1);
u(:,1)=x'.*(1-x)';
u0=u(2:M,1);
for i=1:N
u0=K*u0;
u(2:M,i 1)=u0;
end
mesh(u);
.
├── matlab偏微分方程的求解代码 含流体力学计算模型 一维波动差分 一维热传导有限元等.zip
└── 例子
├── 2013pdeold
│ ├── 12.jpg
│ ├── HyperLF.m
│ ├── HyperLW1.m
│ ├── HyperUP.m
│ ├── Test12.m
│ ├── Test1ParabGuDianXian.m
│ ├── Test2ParabGuDianYin.m
│ ├── Test3EllpsDirichlet.m
│ ├── Test4EllpsNewmann.m
│ ├── Test5.asv
│ ├── Test56.asv
│ ├── Test56.jpg
│ ├── Test56.m
│ ├── Test5HyperCN.m
│ ├── Test6HyperLW.m
│ ├── Test7Hyper2RXian.m
│ ├── Test8Hyper2RYin.m
│ ├── hyper2rxian50x100-1.jpg
│ ├── hyper2rxian50x100.jpg
│ ├── hyper2rxian50x100t=0p5line.jpg
│ ├── hyper2rxian50x100x=0p5line.jpg
│ ├── hyper2ryin50x100.jpg
│ ├── hyper2ryin50x100t=0p5.jpg
│ ├── hyper2ryin50x100t=0p5line.jpg
│ ├── hyper2ryin50x100x=0p5line.jpg
│ ├── hyper3D.jpg
│ ├── hyperbolic-x=0p5line.jpg
│ ├── hyperbolic-x=0p5line1.jpg
│ ├── hyperboliccn-t=0p5line.jpg
│ ├── hyperboliccn-x=0p5line.jpg
│ ├── hyperboliccn.jpg
│ ├── hyperboliccn100x100-t=0p5line.jpg
│ ├── hyperboliccn100x100-x=0p5line.jpg
│ ├── hyperboliccn100x100.jpg
│ ├── hyperboliccn100x100t=0p5.jpg
│ ├── hyperboliccn100x100x=0p5.jpg
│ ├── hyperboliccn100x200-t=0p5line.jpg
│ ├── hyperboliccn100x200-x=0p5line.jpg
│ ├── hyperboliccn100x200.jpg
│ ├── hyperboliccn300x300.jpg
│ ├── hyperboliccn300x300t=0p25line.jpg
│ ├── hyperboliccn300x300t=0p5line.jpg
│ ├── hyperboliccn300x300x=0p25line.jpg
│ ├── hyperboliccn300x300x=0p5line.jpg
│ ├── hyperboliclw-t=0.5line.jpg
│ ├── hyperboliclw-x=0.5line.jpg
│ ├── hyperboliclw.jpg
│ ├── hyperboliclw50x1000-t=0p5line.jpg
│ ├── hyperboliclw50x1000-x=0p5line.jpg
│ ├── hypercn100x100.jpg
│ ├── hypercs1x2pi.jpg
│ ├── t1.m
│ ├── test3.m
│ └── test4.m
├── 2013微分方程数值解编程试题
│ ├── 2013微分方程数值解编程试题.doc
│ ├── Test12.m
│ ├── Test1ParabGuDianXian.m
│ ├── Test2ParabGuDianYin.m
│ ├── Test3EllpsDirichlet.m
│ ├── Test4EllpsNewmann.m
│ ├── Test56.m
│ ├── Test5HyperCN.m
│ ├── Test6HyperLW.m
│ ├── Test7Hyper2RXian.m
│ ├── Test8Hyper2RYin.m
│ ├── ~$13微分方程数值解编程试题.doc
│ └── 新建文件夹
├── EllpsDirichlet-5point.m
├── EllpsNewmann.m
├── Hyper2RXian.m
├── Hyper2RYin.m
├── HyperCN.m
├── HyperLF.m
├── M-Files
│ ├── BeamAssemble.m
│ ├── BeamElementForces.m
│ ├── BeamElementMomentDiagram.m
│ ├── BeamElementShearDiagram.m
│ ├── BeamElementStiffness.m
│ ├── BilinearQuadAssemble.m
│ ├── BilinearQuadElementArea.m
│ ├── BilinearQuadElementPStresses.m
│ ├── BilinearQuadElementStiffness.m
│ ├── BilinearQuadElementStiffness2.m
│ ├── BilinearQuadElementStresses.m
│ ├── GridAssemble.m
│ ├── GridElementForces.m
│ ├── GridElementLength.m
│ ├── GridElementStiffness.m
│ ├── LinearBarAssemble.m
│ ├── LinearBarElementForces.m
│ ├── LinearBarElementStiffness.m
│ ├── LinearBarElementStresses.m
│ ├── LinearTriangleAssemble.m
│ ├── LinearTriangleElementArea.m
│ ├── LinearTriangleElementPStresses.m
│ ├── LinearTriangleElementStiffness.m
│ ├── LinearTriangleElementStresses.m
│ ├── PlaneFrameAssemble.m
│ ├── PlaneFrameElementAxialDiagram.m
│ ├── PlaneFrameElementForces.m
│ ├── PlaneFrameElementLength.m
│ ├── PlaneFrameElementMomentDiagram.m
│ ├── PlaneFrameElementShearDiagram.m
│ ├── PlaneFrameElementStiffness.m
│ ├── PlaneFrameInclinedSupport.m
│ ├── PlaneTrussAssemble.m
│ ├── PlaneTrussElementForce.m
│ ├── PlaneTrussElementLength.m
│ ├── PlaneTrussElementStiffness.m
│ ├── PlaneTrussElementStress.m
│ ├── PlaneTrussInclinedSupport.m
│ ├── QuadTriangleAssemble.m
│ ├── QuadTriangleElementArea.m
│ ├── QuadTriangleElementPStresses.m
│ ├── QuadTriangleElementStiffness.m
│ ├── QuadTriangleElementStresses.m
│ ├── QuadraticBarAssemble.m
│ ├── QuadraticBarElementForces.m
│ ├── QuadraticBarElementStiffness.m
│ ├── QuadraticBarElementStresses.m
│ ├── QuadraticQuadAssemble.m
│ ├── QuadraticQuadElementArea.m
│ ├── QuadraticQuadElementPStresses.m
│ ├── QuadraticQuadElementStiffness.m
│ ├── QuadraticQuadElementStresses.m
│ ├── SpaceFrameAssemble.m
│ ├── SpaceFrameElementAxialDiagram.m
│ ├── SpaceFrameElementForces.m
│ ├── SpaceFrameElementLength.m
│ ├── SpaceFrameElementMomentYDiagram.m
│ ├── SpaceFrameElementMomentZDiagram.m
│ ├── SpaceFrameElementShearYDiagram.m
│ ├── SpaceFrameElementShearZDiagram.m
│ ├── SpaceFrameElementStiffness.m
│ ├── SpaceFrameElementTorsionDiagram.m
│ ├── SpaceTrussAssemble.m
│ ├── SpaceTrussElementForce.m
│ ├── SpaceTrussElementLength.m
│ ├── SpaceTrussElementStiffness.m
│ ├── SpaceTrussElementStress.m
│ ├── SpringAssemble.m
│ ├── SpringElementForces.m
│ ├── SpringElementStiffness.m
│ ├── TetrahedronAssemble.m
│ ├── TetrahedronElementPStresses.m
│ ├── TetrahedronElementStiffness.m
│ ├── TetrahedronElementStresses.m
│ └── TetrahedronElementVolume.m
├── OneDHeatEquFEM0.m
├── ParabGuDianXian.m
├── ParabGuDianYin.m
├── allpde
│ ├── 1DHeatEquFEM.m
│ ├── Adi2Dpara.m
│ ├── Adischem.m
│ ├── EllpsDirichlet.m
│ ├── HyperbolicFDM.m
│ ├── HyperbolicFDM1.m
│ ├── ParabGuDianXian.m
│ ├── ParabGuDianYin.m
│ ├── Parabgudianxian.jpg
│ ├── ch3ex3-1.jpg
│ ├── ch3ex3-2.jpg
│ ├── ch3ex3-3.jpg
│ ├── ch3ex3-pdetool.m
│ ├── ellpsdiricjlet.jpg
│ ├── ex2p1-20X256.jpg
│ ├── ex2p1-30X256.jpg
│ ├── ex2p1-40X256.jpg
│ ├── ex2p1-8X256.jpg
│ ├── ex2p2.m
│ ├── laplacezero.m
│ ├── parabGudianxian1.jpg
│ ├── pdepeex1-1.jpg
│ ├── pdepeex1-2.jpg
│ ├── pdepeex1.jpg
│ ├── pdepeex2-1.jpg
│ ├── pdepeex2-2.jpg
│ ├── pdepeexe1.m
│ ├── pdepeexe2.m
│ ├── t=0.5parab.jpg
│ └── x=0.5parab.jpg
├── ddex1.m
├── exe1.m
├── exe2.jpg
├── exe2.m
├── exe21-1.jpg
├── exe21.jpg
├── exe21.m
├── exe3.m
├── exe31.jpg
├── exe32.jpg
├── heattransfer.m
├── mmbvp.m
├── myddex.m
├── parabolic.m
├── pdeex1.jpg
├── pdex1.m
├── simpleellipse.m
├── whlupbvp4c.m
├── 均匀网格2DSEEP
│ ├── InitMesh.m
│ ├── LinearTriELementRightEE.m
│ ├── LinearTriELementStiffnessAE.m
│ ├── LinearTriELementStiffnessBE.m
│ ├── LinearTriRightAssembleE.m
│ ├── LinearTriStiffnessAssembleAB.m
│ ├── PDEPlot.m
│ ├── QUAD0.08x687x2009925
│ │ ├── QUAD0.08Mesh.jpg
│ │ └── QUADLog.txt
│ ├── QUAD0.8x11x2009925
│ │ ├── QUAD0.8Mesh.jpg
│ │ └── QUADLog.txt
│ ├── QUAD10.08x565x2009925
│ │ ├── QUAD10.08Mesh.jpg
│ │ └── QUAD1Log.txt
│ ├── QUAD10.08x567x2009925
│ │ ├── QUAD10.08Mesh.jpg
│ │ └── QUAD1Log.txt
│ ├── QUAD10.08x694x20091123
│ │ ├── QUAD10.08Mesh.jpg
│ │ ├── QUAD1Log.txt
│ │ └── QUAD1hcont694-0.jpg
│ ├── QUAD10.08x694x2009925
│ │ ├── QUAD10.08Mesh.jpg
│ │ ├── QUAD1Log.txt
│ │ ├── QUAD1hcont694-0.jpg
│ │ ├── QUAD1hcont694-1.jpg
│ │ ├── QUAD1hcont694-2.jpg
│ │ ├── QUAD1hcont694-3.jpg
│ │ └── QUAD1hcont694-4.jpg
│ ├── QUAD10.08x742x2009925
│ │ ├── QUAD10.08Mesh.jpg
│ │ └── QUAD1Log.txt
│ ├── QUAD10.3x59x20091123
│ │ ├── QUAD10.3Mesh.jpg
│ │ ├── QUAD1Log.txt
│ │ ├── QUAD1hcont59-0.jpg
│ │ ├── QUAD1hcont59-1.jpg
│ │ └── QUAD1hcont59-2.jpg
│ ├── Qfun.m
│ ├── RECT0.08x571x2009925
│ │ ├── RECT0.08Mesh.jpg
│ │ └── RECTLog.txt
│ ├── RKSolution.m
│ ├── RegionBoundaryNum.m
│ ├── SEEPRKUNLSQR.m
│ ├── StageSolutionOne.m
│ ├── TRAPEZE0.08x830x2009925
│ │ ├── TRAPEZE0.08Mesh.jpg
│ │ └── TRAPEZELog.txt
│ ├── fun1.m
│ ├── fun2.m
│ ├── fun3.m
│ ├── fun4.m
│ ├── funinit.m
│ └── hs_err_pid184.log
├── 一维波动差分
│ ├── HyperbolicFDM.m
│ └── HyperbolicFDM1.m
└── 一维热传导有限元
└── 1DHeatEquFEM.m
19 directories, 249 files
小贴士
感谢您为本站写下的评论,您的评论对其它用户来说具有重要的参考价值,所以请认真填写。
- 类似“顶”、“沙发”之类没有营养的文字,对勤劳贡献的楼主来说是令人沮丧的反馈信息。
- 相信您也不想看到一排文字/表情墙,所以请不要反馈意义不大的重复字符,也请尽量不要纯表情的回复。
- 提问之前请再仔细看一遍楼主的说明,或许是您遗漏了。
- 请勿到处挖坑绊人、招贴广告。既占空间让人厌烦,又没人会搭理,于人于己都无利。
关于好例子网
本站旨在为广大IT学习爱好者提供一个非营利性互相学习交流分享平台。本站所有资源都可以被免费获取学习研究。本站资源来自网友分享,对搜索内容的合法性不具有预见性、识别性、控制性,仅供学习研究,请务必在下载后24小时内给予删除,不得用于其他任何用途,否则后果自负。基于互联网的特殊性,平台无法对用户传输的作品、信息、内容的权属或合法性、安全性、合规性、真实性、科学性、完整权、有效性等进行实质审查;无论平台是否已进行审查,用户均应自行承担因其传输的作品、信息、内容而可能或已经产生的侵权或权属纠纷等法律责任。本站所有资源不代表本站的观点或立场,基于网友分享,根据中国法律《信息网络传播权保护条例》第二十二与二十三条之规定,若资源存在侵权或相关问题请联系本站客服人员,点此联系我们。关于更多版权及免责申明参见 版权及免责申明
网友评论
我要评论