matlab怎么显示插值得到的函数自己用MATLAB编了牛顿插值的函数 放在三个m文件里面 但是不知道怎么输出最后得到的那个多项式function u = chashang(x,y)N = length(x);if N == 2u = (y(2)-y(1))/(x(2)-x(1));elsex0 =
来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/19 17:45:21
matlab怎么显示插值得到的函数自己用MATLAB编了牛顿插值的函数 放在三个m文件里面 但是不知道怎么输出最后得到的那个多项式function u = chashang(x,y)N = length(x);if N == 2u = (y(2)-y(1))/(x(2)-x(1));elsex0 =
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');
matlab怎么显示插值得到的函数自己用MATLAB编了牛顿插值的函数 放在三个m文件里面 但是不知道怎么输出最后得到的那个多项式function u = chashang(x,y)N = length(x);if N == 2u = (y(2)-y(1))/(x(2)-x(1));elsex0 =
p=poly2str(x0,'x');
是这样吗?望采纳