数值计算常见编程汇总20.docxVIP

秦九韶法 利用秦九韶算法求多项式p(x)=x^5-3x^4+4x^2-x+1在x=3时的值。 x=input(请输入x：); n=input(请输入多项式最高幂数n:); a=input(请按照升幂顺序写出系数矩阵：); v0=a(n+1); for k=1:n v1=v0*x+a(n+1-k); v0=v1; end v=v0 请输入x：3 请输入多项式最高幂数n:5 请按照升幂顺序写出系数矩阵：[1 -1 4 0 -3 1] v =34 二分法 用二分法求方程x^3-x-1=0在[1,2]内的近似根，要求误差不超过10-3 a=1; b=2; fa=a^3-a-1; fb=b^3-b-1; c=(a+b)/2; fc=c^3-c-1; if fa*fb0 ,break ,end while abs(fc)0.001 c=(a+b)/2; fc=c^3-c-1; if fa*fc0; b=c; fb=fc; else a=c; fa=fc; end end format long fx=fc, x=c fx = -4.659488331526518e-05 x = 1.324707031250000 拉格朗日差值法： 已知点（1 0）、（2 7）、（3 16）求x=4的值。 x=input(以行向量形式输入x的取值); y=input(以行向量形式输入y的取值); n=input(多项式的次数); x1=input(计算在x处的函数值,x1=); for k=1:(n+1) l(k)=1; end for k=1:(n+1) for i=1:k-1 l(k)=l(k)*(x1-x(i))/(x(k)-x(i)); end for i=(k+1):(n+1) l(k)=l(k)*(x1-x(i))/(x(k)-x(i)); end end p=0; for k=1:n+1 p=p+l(k)*y(k); end disp(p); 以行向量形式输入x的取值[1 2 3] 以行向量形式输入y的取值[0 7 16] 多项式的次数2 计算在x处的函数值,x1=4 27 分段插值法： 已知点（0 0.9）、（4 3.4）、（6 5.8）（8 8.9）（10 12）求x=7.5的值。 x = [0 4 6 8 10]; y = [0.9 3.4 5.8 8.9 12]; n = length(x); x0=input(输入待求x0的值：); for i=1:n-1 if x0 x(1) || x0 x(n) disp(不符合条件，无法进行分段插值); break; end if x(i) = x0 x0 x(i+1) f = y(i) + ( y(i+1) - y(i)) / ( x(i+1)-x(i) ) * ( x0 - x(i) ); end end F 输入待求x0的值：7.5 f =8.1250 牛顿差值 已知点（2 3），（3 7），（5 9），（8 12）求x=7的值 x0=2;x1=3;x2=5;x3=8; y0=3;y1=7;y2=9;y3=12; x=7; f11=(y1-y0)/(x1-x0); f12=(y2-y0)/(x2-x0); f13=(y3-y0)/(x3-x0); f22=(f12-f11)/(x2-x1); f23=(f13-f11)/(x3-x1); f33=(f23-f22)/(x3-x2); N=y0+f11*(x-x0)+f22*(x-x1)*(x-x0)+f33*(x-x2)*(x-x1)*(x-x0); p=N p = 9.6667 分段三次插值： 已知：p(0)=0,p(1)=3,p(2)=12,p(3)=33,p，(0)=2，p，(1)=5，p，(2)=14，p，(3)=29，用分段三次插值法求解p(1.2) clc; x=input(已知x的取值); y=input(已知y的取值); yy=input(已知y的一阶导数的取值); n=input(‘请输入x值的个数：’); x0=input(x0的取值);y0=0; for i=1:n-1 if(x0x(1)||x0x(n)) disp(输入的值不在所给范围内); break; end if(x(i)=x0x0x(i+1)) h=x(i+1)-x(i); t=(x0-x(i))/h; y