【计算方法重要程序.docxVIP

第2章插值法一、Lagrange插值functionlagr_interpx0=[1 2 3 4 5];y0=[1 1.414 1.732 2 2.3];x=9;m=length(x0);s=0;for i=1:m t=1;for j=1:mif i~=j t=t*(x-x0(j))/(x0(i)-x0(j));endends=s+t*y0(i);ends第3章函数逼近与快速Fourier变换一、正交多项式拟合functionorth_iterx=(37:43);y=[3.40;3.00;2.10;1.53;1.80;1.90;2.90];n=length(x);S(:,1)=x.^0;for i=2:3 lambda(i)=((x.*S(:,i-1))*S(:,i-1))/((S(:,i-1))*S(:,i-1));if i2miu(i)=(S(:,i-1)*S(:,i-1))/(S(:,i-2)*S(:,i-2));elsemiu(i)=0;endif i3 S(:,i)=(x-lambda(i)*ones(n,1)).*S(:,i-1);else S(:,i)=(x-lambda(i)*ones(n,1)).*S(:,i-1)-miu(i)*S(:,i-2);endendfor i=1:3 a(i)=y*S(:,i)/(S(:,i)*S(:,i));enda=a二、矛盾方程组最小二乘拟合function fitx=(37:43);y=[3.4;3;2.1;1.53;1.8;1.9;2.9];X=[x.^0,x.^2,sin(x)];b=X\y第4章数值积分与数值微分一、复合梯形求积functionfuhetixinga=0;b=1;n=1e3;h=(b-a)/n;x=a+h:h:b-h;T=h/2*(f(a)+2*sum(f(x))+f(b))function y=f(x)y=4./(1+x.^2);二、复合Simpson求积Function fuheSimpsona=0;b=1;n=1e4;h=(b-a)/n;x1=a+h/2:h:b-h/2;x2=a+h:h:b-h;T=h/6*(f(a)+4*sum(f(x1))+2*sum(f(x2))+f(b))function y=f(x)y=4./(1+x.^2);三、Romberg求积function Romberga=0;b=1;eps=1e-15;R1=1e8;R0=1e10;i=1;T(1)=(b-a)/2*(f(a)+f(b));while abs(R1-R0)eps h=(b-a)/2^i; x=a+h:2*h:b-h;T(i+1)=T(i)/2+h*sum(f(x));S(i)=T(i+1)+(T(i+1)-T(i))/3;if i1C(i-1)=S(i)+(S(i)-S(i-1))/15;endif i2R(i-2)=C(i-1)+(C(i-1)-C(i-2))/63;endif i3 R0=R(i-3); R1=R(i-2);end i=i+1;ends=R0function y=f(x)y=4./(1+x.^2);第5章解线性方程组的直接方法一、列主元素Gauss消去法functionGaussA=[2,-1,3;4,2,5;1,2,0];b=[1;4;7];G=[A,b];n=length(A);[p,q]=max(abs(G(:,1)));t=G(q,:);G(q,:)=G(1,:);G(1,:)=t;for j=1:n-1for i=j+1:nm(i,j)=G(i,j)/G(j,j);G(i,:)=G(i,:)-m(i,j)*G(j,:);endendA=G(:,1:n);b=G(:,n+1);x(n)=b(n)/A(n,n);for i=n-1:-1:1 x(i)=(b(i)-A(i,i+1:n)*(x(i+1:n)))/A(i,i);endx=x二、列主元素Gauss – Jordan消去法function GJA=[2,-1,3;4,2,5;1,2,0];b=[1;4;7];G=[A,b];n=length(A);[p,q]=max(abs(G(:,1)));t=G(q,:);G(q,:)=G(1,:);G(1,:)=t;for j=1:n-1for i=j+1:nm(i,j)=G(i,j)/G(j,j);G(i,:)=G(i,:)-m(i,j)*G(j,:);endendfor j=n:-1:2for i=j-1:-1:1m(i,j)=G(i,j)/G(j,j); G(i,:)=G(i,:)/m(i,j)-G(j,:);endendfor i=1:nG(i,:)=G(i,:)/G(i,i);endx=G(:,n+1)三、三角分解法（LU分