我想使用函数的输出作为构建多项式的函数的输入:
这是我的代码:
function c = interpolation(x, y)
n = length(x);
V = ones(n);
for j = 2:n
V(:,j) = x.*V(:,j-1);
end
c = V \ y;
disp(V)
for i = 0:n-1
fprintf('c%d= %.3f\n', i, c(i+1));
end
polynome(c);
function p = polynome(x)
n = length(x);
for l= 0:n-1
polynome = polynome * x^l;
end
end
单独的第一个功能,有效。这意味着我的代码可以工作,如果我从第13行开始注释到结束,并且我得到c值,其数量取决于开头输入的x向量的长度。
我想使用那个c值来构建这种形式的多项式:p(x) = c0 + c^1*x1 + c2^x2 + .... + c(n-1)^x(n-1)
并绘制该多项式,其中通过2个向量在开始时给出的点xi,yi作为函数插值的输入。
有人可以帮助我吗?
答案 0 :(得分:0)
制作单独的函数多项式,例如
function y=polynome(x,c)
y=sum(c.*x.^(0:length(c)-1));
或只是使用
y=sum(c.*x.^(0:length(c)-1));
计算系数c的多项式。
如果您有多个x值,例如
x=[0:.1:3]';
y=repmat(x,1,numel(c)).^repmat((0:numel(c)-1),numel(x),1)*c';
应该给你多项式的值