在數值計算中,某個函數f(x)的值可以僅在有限數量的點處獲得,其中x = x 1,x 2,...,x N 在這種情況下,當您想要的值不是您所擁有的值時,插值是有效的。 https://github.com/kbarbary/Dierckx.jl 按,]切換到Pkg模式並添加Dierckx。
=========程式===========
import Pkg Pkg.add("Dierckx") using Dierckx x = [0., 1., 2., 3., 4.] y = [-1., 0., 7., 26., 63.] # x.^3 - 1. spl = Spline1D(x, y) a=spl([1.5, 2.5]) # result = [2.375, 14.625] println(a) a=spl(1.5) # result = 2.375 println(a) a=derivative(spl, 1.5) # 在x = 1.5;值為5.75 println(a) a=integrate(spl, 0., 4.) # 從x = 0到x = 4的積分;值為60.0 println(a) a=roots(spl) # 找到0點。值為[1.0] println(a)
輸出畫面
Resolving package versions... Updating `~/.julia/Project.toml` [no changes] Updating `~/.julia/Manifest.toml` [no changes] [2.375, 14.625] 2.3750000000000004 6.7500000000000036 60.00000000000001 [1.0]
沒有留言:
張貼留言