文件名称:matlab最简单的代码-juliafem:二维Poisson方程的简单FEM代码
文件大小:178KB
文件格式:ZIP
更新时间:2024-06-11 05:29:36
系统开源
matlab最简单的代码该代码是对“大约50行Matlab的注释:简短的有限元实现”的Matlab代码的重写。 用茱莉亚语言。 对于matlab版本,请参见 去做 添加四边形的绘图解决方案。 Matlab代码使用trisurf,可以绘制三角形和四边形。 我在PyPlot中找不到等效的功能。 我通过PyPlot的tricontour函数在Julia代码中使用等高线图。 如果您在“样本”中运行示例,它将通过对它们进行三角剖分来绘制四边形。 如何在“正方形”中运行示例 您需要安装PyPlot软件包以绘制结果。 import Pkg Pkg.add("PyPlot") 首先,生成网格 bash> cd square 通过在Matlab中运行square.m生成网格(julia版本未完成) matlab> square(30) matlab> quit 运行fem代码 bash> julia run.jl 你应该得到这样的解决方案 构造稀疏矩阵 初始化空的nxn稀疏矩阵的最简单方法是 julia> A = spzeros(n,n) 这利用了64个字节的整数和浮点数。 对于较小的问题,使用32个字
【文件预览】:
juliafem-master
----sample()
--------run.jl(68B)
--------elements4.dat(102B)
--------coordinates.dat(195B)
--------dirichlet.dat(64B)
--------elements3.dat(44B)
--------Data.jl(265B)
--------neumann.dat(32B)
----src()
--------fem.jl(3KB)
----square()
--------run.jl(68B)
--------square.jl(1KB)
--------square.m(1KB)
--------Data.jl(300B)
----README.md(2KB)
----square.png(175KB)