文件名称:数值分析_第五单元解析
文件大小: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)