'-r');axis([-10^(-15),10^(-15),-1,2]);legend('第二种算法');运行结果:显然第二种算法结果不准确,是因为计算机中的舍入误差造成的,当时,,计算机进行舍入造成恒等于1,结果函数值恒为1。程序:秦九韶算法:QinJS.mfunctiony=QinJS(a,x)%y输出函数值%a多项式系数,由高次到零次%x给定点n=length(a);s=a(1);fori=2:ns=s*x+a(i);endy=s;计算p(x):test3.mclearall;clc;x=1.6:0.2:2.4;%x=2的邻域disp('x=2的邻域:');xa=[1-18144-6722016-40325376-46082304-512];p=zeros(1,5);fori=1:5p(i)=QinJS(a,x(i));enddisp('相应多项式p值:');pxk=1.95:0.01:20.5;nk=length(xk);pk=zeros(1,nk);k=1;fork=1:nkpk(k)=QinJS(a,xk(k));endplot(xk,pk,'-r');xlabel('x');ylabel('p(x)');运行结果:x=2的邻域:x=1.60001.80002.00002.20002.4000相应多项式p值:p=1.0e-003*-0.2621-0.000500.00050.2621p(x)在[1.95,20.5]上的图像程序:LU分解,.mfunction[L,U]=(A)%不带列主元的LU分解N=size(A);n=N(1);L=eye(n);U=zeros(n);fori=1:nU(1,i)=A(1,i);L(i,1)=A(i,1)/U(1,1);endfori=2:nforj=i:nz=0;fork=1:i-1z=z+L(i,k)*U(k,j);end
全文预览
