数值分析(第五版)计算实习题第三章.pdfVIP

数值分析计算实习题第三章

第二次作业：

题一：

x=-1:0.2:1;y=1./(1+25.*x.^2);

f1=polyfit(x,y,3)

f=poly2sym(f1)

y1=polyval(f1,x)

x2=linspace(-1,1,10)

y2=interp1(x,y,x2)

plot(x,y,r*-,x,y1,b-)

holdon

plot(x2,y2,k)

legend(数据点,3次拟合曲线,3次多项式插值)

xlabel(X),ylabel(Y)

输出：f1=

0.0000-0.57520.00000.4841

f=

y1=

-0.09110.11600.27710.39210.46110.48410.27710.1160-0.0911

(3305*x^2)/5746++0.46110.3921

x2=

-1.0000-0.7778-0.5556-0.3333-0.11110.11110.33330.55560.77781.0000

y2=

0.03850.06340.12220.30000.72220.72220.30000.12220.06340.0385

题二：

X=[0.00.10.20.30.50.81.0];

Y=[1.00.410.500.610.912.022.46];

p1=polyfit(X,Y,3)

p2=polyfit(X,Y,4)

Y1=polyval(p1,X)

Y2=polyval(p2,X)

plot(X,Y,r*,X,Y1,b-.,X,Y2,g--)

p3=polyfit(X,Y,2)

Y3=polyval(p3,X)

f1=poly2sym(p1)

f2=poly2sym(p2)

f3=poly2sym(p3)

plot(X,Y,r*,X,Y1,b-.,X,Y2,g--,X,Y3,m--)

legend(数据点,3次多项式拟合,4次多项式拟合,2次多项式拟合)xlabel(X轴

),ylabel(Y轴)

输出：

p1=

-6.622112.8147-4.65910.9266

p2=

2.8853-12.334816.2747-5.29870.9427

Y1=

0.92660.58220.45440.50340.97302.01032.4602

Y2=

0.94270.56350.43990.50821.00051.98602.4692

p3=

3.1316-1.24000.7356

Y3=

0.73560.64290.61280.64540.89841.74772.6271

f1=

f2=

f3=

(18733*x^2)/5982-(74179*x)/59820+73337/99700

题三：

建立三角插值函数的m文件

function[A,B,Y1,Rm]=sanjiaobijin(X,Y,X1,m)%AB分别是m阶三角多项式Tm

（x）的系数aj，bj（j=1,2,...,m）的系数矩阵，Y1是Tm（x）在X1处的值，XY

数据点,Rm为均方误差

n=length(X)-1;max1=fix((n-1)/2);

ifmmax1

m=max1;

end

A=zeros(1,m+1);B=zeros(1,m+1);

Ym=(Y(1)+Y(n+1))/2;Y(1)=Ym;

Y(n+1)=Ym;A(1)=2*sum(Y)/n;

fori=1:m

B(i+1)=sin(i*X)*Y;

A(i+1)=cos(i*X)*Y;

end

A=2*A/n;B=2*B/n;

A(1)=A(1)/2;Y1=A(1);

fork=1:m

Y1=Y1+A(k+1)*cos(k*X1)+B(k+1)*sin(k*X1);

Tm=A(1)+A(k+1).*cos(k*X)+B(k+1).*sin(k*X);k=k+1;

end

Y,Tm,Rm=(sum(Y-Tm).^2)/n

输出：X=-pi:2*pi/33:pi;

Y=X.^2.*cos(X);

[A,B,Y1,Rm]=sanjiaobijin(X,Y,X1,16)

输出：

A=

1至12列

-0.13974.4002-2.83261.2355-0.91280.6982-0.67730.6635-