【文件属性】:
文件名称:数值分析_第五单元解析
文件大小:360KB
文件格式:PDF
更新时间:2022-07-09 22:04:30
国科大 数值分析
clear all;
clc
N = 3;
format long;
x = -1:0.1:1;
for i=1:N
for j=1:N
H(i,j)=2/(2*i+2*j-3);
end
end
b = [1;1/2;1/3];
p = nagauss(H,b);%调用nagauss求解多项式系数 p
y = abs(x);
y1 = p(1) + p(2)*x.^2 + p(3)*x.^4;
plot(x,y,'o',x,y1,'-');
axis([-1 1 -0.5 1.5]);
title('函数最佳平方趋近多项式')
xlabel('x');
ylabel('y');
legend('原始数据','拟合曲线')
p
syms t;
yy = t^2;
cc = int(yy,t,-1,1)-p'*b
ccc = sqrt(cc)