Hello People
I have been working on heat equation simulation using Julia Software. I was stuck with creating a surface plot (3d plotting)
I have values of x,t and T (is a matrix) how to make 3d plot
Please help me
1 Like
Do you need the following examples?
using Plots
f(t, x) = exp(-x^2/(2t))/√(2π*t)
u(t, x) = f(t+0.1, x+2) + 2f(t+0.2, x-2)
t = range(0, 3; length=101)
x = range(-5, 5; length=101)
T = u.(t', x)
surface(t, x, T; colorbar=false,
xlabel="t", ylabel="x", zlabel="T",
camera=(30, 70), color=reverse(cgrad(:RdYlBu_11)))
t = range(0, 3; length=51)
x = range(-5, 5; length=51)
T = u.(t', x)
wireframe(t, x, T; colorbar=false,
xlabel="t", ylabel="x", zlabel="T",
camera=(30, 70))
1 Like
Hello Thanks for sharing the examples. But when I am trying to run this code I am finding this errors please help me to clear them
using LinearAlgebra, SparseArrays, Plots, DifferentialEquations
using Surrogates
using Plots
gr()
L=20
N=100
k=50
x=range(0, stop=100, length=100);
dx=x[2]-x[1];
dt=0.5*dx^2/(2*k);
t=0:dt:10
# print("the vlues of t",t)
T=zeros(length(t),length(x))
s=size(T)
print("size of matrix",s)
T[1,2]=T[1,1]=100 #complete 1st row
T[end,end]=150 #last row of a matrix
for i in 2:1960
for j=1:99
T[i,j+1]=T[i,j]+k*dt/(dx^2)*(T[i+1,j]-2*T[i,j] + T[i-1,j])
end
end
surface(x,t,T)
This is my code
I am trying to simulate the function. But i am finding the errors
I’m using Windows 10, and the output of the following code is
VERSION = v"1.6.2"
Status `D:\.julia\environments\v1.6\Project.toml`
[91a5bcdd] Plots v1.20.1
Status `D:\.julia\environments\v1.6\Manifest.toml`
[28b8d3ca] GR v0.58.1
Code:
@show VERSION
using Pkg
Pkg.status("Plots")
Pkg.status("GR"; mode=PKGMODE_MANIFEST)
If you have older versions, I suggest you upgrade and try again.
VERSION = v"1.6.0"
Status `C:\Users\manis\.julia\environments\v1.6\Project.toml`
[91a5bcdd] Plots v1.20.0 `https://github.com/JuliaPlots/Plots.jl.git#master`
Status `C:\Users\manis\.julia\environments\v1.6\Manifest.toml`
[28b8d3ca] GR v0.57.4
I too have the same version
Corrected version of your code:
#using LinearAlgebra, SparseArrays, Plots, DifferentialEquations
#using Surrogates
using Plots
#gr()
L=20
N=100
k=50
x=range(0, stop=100, length=100);
dx=x[2]-x[1];
dt=0.5*dx^2/(2*k);
t=0:dt:10
# print("the vlues of t",t)
#T=zeros(length(t),length(x))
T=zeros(length(x),length(t)) # t <--> x
s=size(T)
print("size of matrix",s)
#T[1,2]=T[1,1]=100 #complete 1st row
#T[end,end]=150 #last row of a matrix
@. T[:, 1] = exp(-(x-50)^2/(2*5^2)) # initial values
#for i in 2:1960
for j in 1:1960 # i -> j
# for j=1:99
for i=2:99 # j -> i, 1:99 -> 2:99
T[i,j+1]=T[i,j]+k*dt/(dx^2)*(T[i+1,j]-2*T[i,j] + T[i-1,j])
end
end
#surface(x,t,T)
surface(t,x,T) # x <--> t
2 Likes
Can you please explain me this line. I havent understood this one as I am new to the Julia
To discover what a symbol means you can type a ? in the REPL. This will get you into the help mode of the RPEL. Once you are in there, you can type the symbol you want to undersantd and click enter.
doing ? @. enter gives me:
@. expr
Convert every function call or operator in expr into a "dot call" (e.g.
convert f(x) to f.(x)), and convert every assignment in expr to a "dot
assignment" (e.g. convert += to .+=).
If you want to avoid adding dots for selected function calls in expr, splice
those function calls in with $. For example, @. sqrt(abs($sort(x))) is
equivalent to sqrt.(abs.(sort(x))) (no dot for sort).
(@. is equivalent to a call to @__dot__.)
So your expression is equivalent to T[:, 1] .= exp.(-(x .- 50) .^ 2 ./ ( 2 .* 5 .^2)) # initial values
This basically assigns to each element of T[:,1] the value of the function evaluated at each point in the vector x.
The dot notation is broadcasting. You can understand it if you go through this section of the maual. Mathematical Operations and Elementary Functions · The Julia Language
1 Like
Under the condition T=zeros(length(x),length(t)), the code @. T[:, 1] = exp(-(x-50)^2/(2*5^2)) is equivalent to the following for loop:
for i in eachindex(x)
T[i, 1] = exp(-(x[i]-50)^2/(2*5^2))
end
Code for verification:
using Plots
k = 50
x = range(0, 100; length=100);
dx = step(x)
dt = 0.5*dx^2/(2k)
t = 0:dt:10
T = zeros(length(x), length(t))
U = zeros(length(x), length(t))
@. T[:, 1] = exp(-(x-50)^2/(2*5^2))
for i in eachindex(x)
U[i, 1] = exp(-(x[i]-50)^2/(2*5^2))
end
@show T ≈ U
plot(x, T[:, 1]; label="T[:,1]", lw=2)
plot!(x, U[:, 1]; label="U[:,1]", lw=2, ls=:dash)
Output:
T ≈ U = true
You could verify the results of the following code with REPL.
julia> T[:, 1]
julia> plot(T[:, 1])
julia> ?
help?> @.
It will be easier to use Julia if you always ask REPL❗️
Jupyter notebook: https://github.com/genkuroki/public/blob/main/0015/Surrogates/3D%20plot%20example%20-%20heat%20eq.ipynb
Edit: Add Jupyter notebook.
1 Like