数值分析第五章答案Word下载.docx
- 文档编号:14989359
- 上传时间:2022-10-26
- 格式:DOCX
- 页数:13
- 大小:18.02KB
数值分析第五章答案Word下载.docx
《数值分析第五章答案Word下载.docx》由会员分享,可在线阅读,更多相关《数值分析第五章答案Word下载.docx(13页珍藏版)》请在冰豆网上搜索。
p4=vpa(collect((sum(d))),5)%p4即为4次牛顿插值多项式,并保留小数点后5位数pp=csape(x1,y1,variational);
%调用三次样条函数
q=pp.coefs;
disp(三次样条函数);
fori=1:
4
s=q(i,:
)*[(x-x1(i))^3;
(x-x1(i))^2;
(x-x1(i));
1];
s=vpa(collect(s),5)
x2=0.2:
0.08:
1.08;
dot=[121112];
figure
ezplot(p4,[0.2,1.08]);
holdon
y2=fnval(pp,x2);
x=x2(dot);
y3=eval(p4);
y4=fnval(pp,x2(dot));
plot(x2,y2,r,x2(dot),y3,b*,x2(dot),y4,co);
title(4次牛顿插值及三次样条);
结果如下:
4次牛顿插值多项式
p4=-0.52083*x^4+0.83333*x^3-1.1042*x^2+0.19167*x+0.98三次样条函数
x∈[0.2,0.4]时,s=-1.3393*x^3+0.80357*x^2-0.40714*x+1.04x∈[0.4,0.6]时,s=0.44643*x^3-1.3393*x^2+0.45*x+0.92571x∈[0.6,0.8]时,s=-1.6964*x^3+2.5179*x^2-1.8643*x+1.3886x∈[0.8,1.0]时,s=2.5893*x^3-7.7679*x^2+6.3643*x-0.80571输出图如下
2-3
(1)
clc;
x1=[01491625364964];
y1=[012345678];
%插值点
a=ones(n,2);
a(:
2)=-x1;
c=1;
n
c=conv(c,a(i,:
));
q=zeros(n,n);
r=zeros(n,n+1);
[q(i,:
),r(i,:
)]=deconv(c,a(i,:
%wn+1/(x-xk)
end
dw=zeros(1,n);
dw(i)=y1(i)/polyval(q(i,:
),x1(i));
%系数
p=dw*q;
symsxl8;
l8(i)=p(n-i+1)*x^(i-1);
disp(8次拉格朗日插值);
l8=vpa(collect((sum(l8))),5)
xi=0:
64;
yi=polyval(p,xi);
plot(xi,yi,x1,y1,r*);
title(8次拉格朗日插值);
8次拉格朗日插值
l8=-3.2806e-10*x^8+6.7127e-8*x^7-5.4292e-6*x^6+0.00022297*x^5-0.0049807*x^4+0.060429*x^3-0.38141*x^2+1.3257*x
输出图如下:
第五章
4-1(3)
y=@(x)sqrt(x).*log(x);
a=0;
b=1;
tol=1e-4;
p=quad(y,a,b,tol);
fprintf(采用自适应辛普森积分结果为:
%d\n,p);
采用自适应辛普森积分结果为:
-4.439756e-01
第九章
9-1
(a)程序:
a=1;
b=2;
%定义域
h=0.05;
%步长
n=(b-a)/h;
y0=1;
%初值
f=@(x,y)1/x^2-y/x;
%微分函数
xn=linspace(a,b,n+1);
%将定义域分为n等份yn=zeros(1,n);
%结果矩阵
yn
(1)=y0;
%赋初值
%以下根据改进欧拉公式求解
xn=xn(i);
xnn=xn(i+1);
yn=yn(i);
yp=yn+h*f(xn,yn);
yc=yn+h*f(xnn,yp);
yn=(yp+yc)/2;
yn(i+1)=yn;
xn=yn;
%以下根据经典四阶r-k法公式求解
k1=f(xn,yn);
k2=f(xn+h/2,yn+h/2*k1);
k3=f(xn+h/2,yn+h/2*k2);
k4=f(xn+h,yn+h*k3);
yn=yn+h/6*(k1+2*k2+2*k3+k4);
yn(i+1)=yn;
disp(改进欧拉法四阶经典r-k法);
disp([xnyn])
改进欧拉法四阶经典r-k法11
0.998870.99885
0.995770.9978
0.991140.99694
0.985320.99634
0.978570.99603
0.971110.99606
0.963110.99645
0.95470.99723
0.945980.99841
0.937051
0.927981.002
0.918831.0044
0.909641.0073
0.900451.0106
0.891291.0143
0.882181.0184
0.873151.0229
0.864211.0278
0.855381.0331
0.846651.0388
(b)程序:
h=[0.10.0250.01];
y0=1/3;
f=@(x,y)-50*y+50*x^2+2*x;
%微分函数xi=linspace(a,b,11);
y=1/3*exp(-50*xi)+xi.^2;
%准确解ym=zeros(1,11);
forj=1:
3
【篇二:
数值分析(第五版)计算实习题第五章作业】
题:
lu分解法:
建立m文件
functionh1=zhijielu(a,b)%h1各阶主子式的行列式值
[nn]=size(a);
ra=rank(a);
ifra~=n
disp(请注意:
因为a的n阶行列式h1等于零,所以a不能进行lu分解。
a的秩ra如下:
)ra,h1=det(a);
return
ifra==n
forp=1:
h(p)=det(a(1:
p,1:
p));
h1=h(1:
n);
ifh(1,i)==0
因为a的r阶主子式等于零,所以a不能进行lu分解。
a的秩ra和各阶顺序主子式h1依次如下:
)
h1;
ra
ifh(1,i)~=0
因为a的r阶主子式都不等于零,所以a能进行lu分解。
u(1,j)=a(1,j);
fork=2:
forj=2:
l(1,1)=1;
l(i,i)=1;
ifij
l(2,1)=a(2,1)/u(1,1);
l(i,1)=a(i,1)/u(1,1);
l(i,k)=(a(i,k)-l(i,1:
k-1)*u(1:
k-1,k))/u(k,k);
else
u(k,j)=a(k,j)-l(k,1:
k-1,j);
ra,u,l,x=inv(u)*inv(l)*b
输入:
a=[10-701;
-32.09999962;
5-15-1;
2102];
b=[8;
5.900001;
5;
h1=zhijielu(a,b)
输出:
请注意:
ra=
4
u=
10.0000-7.000001.0000
02.10006.00002.3000
00-2.1429-4.2381
0-0.0000012.7333
l=
1.0000000
-0.30001.000000
0.50001.19051.0000-0.0000
0.20001.14293.20001.0000
x=
-0.2749
-1.3298
1.2969
1.4398
h1=
10.0000-0.0000-150.0001-762.0001
列主元高斯消去法:
function[ra,rb,n,x]=liezhu(a,b)
b=[ab];
n=length(b);
rb=rank(b);
zhicha=rb-ra;
ifzhicha0
因为ra~=rb,所以方程组无解)
warningoffmatlab:
return_outside_of_loop
ifra==rb
因为ra=rb,所以方程组有唯一解)
x=zeros(n,1);
c=zeros(1,n+1);
n-1
[y,j]=max(abs(b(p:
n,p)));
c=b(p,:
b(p,:
)=b(j+p-1,:
b(j+p-1,:
)=c;
fork=p+1:
m=b(k,p)/b(p,p);
b(k,p:
n+1)=b(k,p:
n+1)-m*b(p,p:
n+1);
b=b(1:
n,n+1);
a=b(1:
n,1:
x(n)=b(n)/a(n,n);
forq=n-1:
1
x(q)=(b(q)-sum(a(q,q+1:
n)*x(q+1:
n)))/a(q,q);
因
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数值 分析 第五 答案