作业帮 > 综合 > 作业

matlab怎么显示插值得到的函数

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/25 16:34:59
matlab怎么显示插值得到的函数
自己用MATLAB编了牛顿插值的函数 放在三个m文件里面 但是不知道怎么输出最后得到的那个多项式
function u = chashang(x,y)
N = length(x);
if N == 2
u = (y(2)-y(1))/(x(2)-x(1));
else
x0 = [x(1:N-2) x(N)];
y0 = [y(1:N-2) y(N)];
Temp1 = chashang(x0,y0);
Temp2 = chashang(x(1:N-1),y(1:N-1));
u = (Temp1-Temp2)/(x(N)-x(N-1));
end
function p = newton(x0,y0)
N = length(x0);
if N == 1
p = y0
else
p =@(x) y0(1);
temp = @(x) 1;
for i=1:N-1
x00 = x0(1:i+1);
y00 = y0(1:i+1);
temp = @(x) temp(x)*(x-x0(i));
p =@(x) p(x) + chashang(x00,y00)*temp(x);
end
end
x0 = [0.2 0.4 0.6 0.8 1.0]
y0 = [0.98 0.92 0.81 0.64 0.38];
fh = newton(x0,y0)
figure;
plot(x0,y0,'*r')
hold on;
fplot(fh,[0,1.5],'g');
p=poly2str(x0,'x');
是这样吗?望采纳
再问: 什么意思,不太懂。
我的问题是要把我自己插值得到的函数的具体表达式显示出来 但是不知道应该怎么操作
再答: 可以啊,就是一步就可以了。