1. LU 分解：[L U］=lu(A)；2。 追赶法function x=zhuiganfa（A，b）［n,n]=size（A); for i=1：n if（i==1） l(i）=A(i,i）; y（i)=b（i）/l(i）； else l（i）=A（i，i）—A（i，i-1）＊u(i—1）； y（i）=(b(i)—y（i-1)＊A(i，i—1))/l（i); end if(i>n=101；c1=ones（1，n—1)；a1=diag（c1,—1)； 这个-1 说明行位置—1c2=12＊ones（1,n);a2=diag(c2）; c3=ones(1，n—1）;a3=diag（c3，1); a=a1+a2+a3;3. 拉格朗日插值function yh=lage(x,y,xh）n=length（x)；m=length（xh)；yh=zeros（1,m）；c1=ones（n—1，1）；c2=ones(1，m）;for i=1：nxp=x(［1：i—1 i+1：n]);yh=yh+y（i)＊prod((c1*xh—xp'＊c2）./（x(i）—xp'＊c2））；endend>〉 x=［11,12］;>〉 y=[2，4]；>〉 xh=［11。75]；〉〉 lage（x，y，xh）ans =3。50004 最小二乘法1.最小二乘的 xi 和 yi 为：xi19531964198219902000yi5。826.9510.0811。3412.66要拟合的函数为：y=a+bx—cxy 注意不是多项式2.编程函数为:function z = erchen（x，y）x1=ones（5,1）； A=［x1,x，-x。＊y]; 注意点乘z=A\y； 注意左除a=z（1）；b=z（2）；c=z(3）；end 输入：≻≻ x=[1953 1964 1982 1990 2000］’；≻≻ y=［5。82 6。95 10。08 11。34 12.66］’；≻≻erchen（x，y)ans = 2。9456 =a —0.0014 =b -0.0005 =c1.最小二乘的 xi 和 yi 为：xi00。250.50。751yi11。2841。64872。11702。7183要拟合的函数为：y=a+bx+cx2是多项式2.编程函数为：function z = erchen2（x，y）x1=ones(5,1); A=［x1，x，x.^2];z=A\y；a=z(1）；b=z(2)；c=z（3)；end输入：≻≻ x=［0 0。25 0。5 0.75 1.00］’；≻≻ y=［1。00 1。284 1.6487 2。1170 2。7183]’;≻≻ erchen2（x,y）ans = 1。0051 0。8642 0.8437最小二乘多项式拟合的简单函数方法：≻≻ x=[0 0。25 0。5 0.75 1.00］'；≻≻ y=［1。00 1.284 1.6487 2.1170 2.7183]’；≻≻ P=polyfit(x,y，2） 要拟合成 4...
