How to code the linear Splines function?

2018 11 12
1 个回答
更新时间：2018 11 12
6 次查看（30 天）

0 个投票

Hi all, I am getting stuck with the code of linear Splines function. The error is: "Unable to perform assignment because the indices on the left side are not compatible with the size of the right side" I appreciate for any help. Thanks!
if true
% code
function out=LSplines(x,pointx,pointy)
n=length(pointx);
for i=2:n
t1(i)= (pointx(i)-x);
b1(i)= (pointx(i)-pointx(i-1));
t2(i)= (x-pointx(i-1));
b2(i)= (pointx(i)-pointx(i-1));
end
P==0;
for i=2:n
P(i)= P(i)+pointy(i-1)*t1(i)/b1(i)+pointy(i)*t2(i)/b2(i)
end
end
end

请先登录，再回答此问题。

回答（1 个）

madhan ravi
madhan ravi 2018-11-12
编辑：madhan ravi 2018-11-12

0 个投票

function out=LSplines(x,pointx,pointy)
n=length(pointx);
for i=2:n
t1(i)= (pointx(i)-x);
b1(i)= (pointx(i)-pointx(i-1));
t2(i)= (x-pointx(i-1));
b2(i)= (pointx(i)-pointx(i-1));
end
P=cell(1,n);
P{1}=0;
for i=2:n
P{i}= P{i-1}+pointy(i-1)*t1(i)/b1(i)+pointy(i)*t2(i)/b2(i);
end
out=P;
end

8 个评论
显示 6更早的评论 隐藏 6更早的评论

Ni Su
Ni Su 2018-11-12
Thank you! Madhan, Still having error, "Error in LSplines (line 5)
t1(i)= (pointx(i)-x);" This part cause to be incompatible of right and left side.
Ni Su
Ni Su 2018-11-12
编辑：madhan ravi 2018-11-12
Data is: Write a Mat lab function which will take in n data points and an x value and produce the y value corresponding to the linear spline approximation. Plot the linear spline through the points:(-2,1),(-1.5,3),(-1,0),(-0.5,2),(0,3),(0.5,2),(1,-4),(1.5,7),(2,-1). Thank you!
madhan ravi's reply: so what's x , pointx and pointy in your function?
Ni Su
Ni Su 2018-11-12
编辑：madhan ravi 2018-11-12
pointx: is a vector of x value through the given points:[-2 -1.5 -1 -0.5 0 0.5 1 1.5 2],and same for pointy:[1 -3 0 2 3 2 -4 7 -1],
madhan ravi's reply : what's x?
Ni Su
Ni Su 2018-11-12
x is the values that we use to plot it out (help the line smoother): for example:x=[-3:0.01:3], and y is the corresponding values of x with the Splines function.
Ni Su
Ni Su 2018-11-12
this code is working, some issues happen when I try to plot it out.
function P=LSplines(x,pointx,pointy)
n=length(pointx);
P=cell(1,n);
P{1}=0;
for i=2:n
P{i}= P{i-1}+pointy(i-1)*(pointx(i)-x)/(pointx(i)-pointx(i-1))+pointy(i)*(x-pointx(i-1))/(pointx(i)-pointx(i-1));
end
end
Ni Su
Ni Su 2018-11-12
编辑：Ni Su 2018-11-12
The main code to plot out and the error is:"Invalid data argument"
if true
%Maincode
pointx = input('Enter the vector of x points: \n');
pointy = input('Enter the vector of y points: \n');
x = [min(pointx)-1:0.1:max(pointx)+1];
y = LSplines(x,pointx,pointy);
plot(pointx,pointy,'r*');
hold on
plot(x,y,'b')
end

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Geometric Transformation and Image Registration 的更多信息

标签

提问:

Ni Su
Ni Su
2018-11-12

编辑：

Ni Su
Ni Su
2018-11-12

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by