我假衡键坦亮银设你的x和x(1)--x(n)是不同的变量,那么我用a来表示你的x,用x数组代表你咐桐的x(1)--x(n)
我随便给这些数赋值
a=10
x=1:3:30
k=5
temp = true(size(x))
temp(k) = false
y = prod(a - x) / prod(x(k) - x(temp))
a-x是生成数组[a-x(1), a-x(2), ..., a-x(n)],prod()是吧这个数组相乘
temp是一个逻辑数组,和x相同大小,并且只有第k个数是false。x(temp)就是生成[x(1), x(2),..., x(k-1), x(k+1), ..., x(n)]
其他的应该容易理解
补充:拉格朗日插值多项式可也这么写:
%x0,y0是已知各点坐标,共n点
x0 = [1, 2, 3]
y0 = [1, 8, 27]
n = length(x0)
%x是要求的点
x = 4
y = 0
for k = 1:n
temp = true(size(x0))
temp(k) = false
y = y + y0(k) * prod(x - x0(temp)) / prod(x0(k) - x0(temp))
end
Matlab-2018a软件+教程百度网盘免费资源在线学习消御键
链接: https://pan.baidu.com/s/1GbXCK-5bfGG2Abvv34fR6g
?pwd=yc8c 提取码: yc8cMatlab-2018a软件+教程 内拿巧容截图展示 64位拆搜电脑建议下这个 32位电脑只能下这个 matlab 2014a(支持32位、64位).zip
2014A安装教程.mp4 Matlab-2018a破解版安装包.zip
3.png 2.png 1.png
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)