计算方法matlab上机报告(流程图、源代码等).docxVIP

上机报告 共轭梯度法 （1）计算过程如下： 第一步：取初始向量, 计算 第步：计算 （2） syms n; %定义一个变量n A=input(请输入矩阵A) b=input(请输入矩阵b) n=size(A) X=zeros(n); D=zeros(n); R=zeros(n); x=X(:,1); %将矩阵X的第一列赋值给x作为初始向量 r=b-A*x; %将x代入求得初始非零残向量 R(:,1)=r; d=r; %求得初始方向向量 D(:,1)=d; for i=1:n; %利用循环进行迭代求得各向量 AF(i)=R(:,i)*R(:,i)/(D(:,i)*A*D(:,i)); X(:,i+1)= X(:,i)+AF(i)*D(:,i); R(:,i+1)=b-A*X(:,i+1); BT(i)=norm(R(:,i+1))^2/norm(R(:,i))^2; D(:,i+1)=R(:,i+1)+BT(i)*D(:,i) ; end x=X(:,i+1) %输出计算结果  龙贝格积分 s=input(请输入被积函数表达式：f(x)=,s) %键盘输入被积函数表达式 f=inline(s) a=input(请输入积分下限a=)%输入积分下限a b=input(请输入积分上限b=)%输入积分上限b f0=f(a);%求下限值f(a) f1=f(b);%求上限值f(b) j=zeros(8,4);%定义一个矩阵用来存放T,S,C,R值 j(1,1)=(b-a)*0.5*(f0+f1);%计算出T1 for i=2:8;%对i循环赋值 t=2^(i-2); fj=zeros(t,1);%构造一个矩阵为接下来存放f(a+(2i-1)*(b-a)/2^(k+1))的值做准备 for m=1:t; fj(m,1)=f(a+(2*m-1)*(b-a)/2^(i-1));%将计算得到的f(a+(2i-1)*(b-a)/2^(k+1))的值赋值给对应的矩阵元素 ff=sum(fj);%对矩阵所有元素求和 end j(i,1)=0.5*j(i-1,1)+(1/2^(i-1))*ff;%得到所有的T值 j(i,2)=j(i,1)+(1/3)*(j(i,1)-j(i-1,1));%得到所有的S值 j(i,3)=j(i,2)+(1/3)*(j(i,2)-j(i-1,2));%得到所有的C值 j(i,4)=j(i,3)+(1/3)*(j(i,3)-j(i-1,3));%得到所有的R值 s=j(i,4)-j(i-1,4); end j %输出计算结果表 If=vpa(j(i,4),7) %得到精确的积分结果  三样条插值 X=zeros(1,6) ; %定义一些下面将用到的矩阵 y=zeros(1,6); f=zeros(6); l=zeros(1,6); for i=1:1:6 X(i)=((i-1)*2)/5-1; %将区间等分取点 y(i)=1/(1+25*X(i)*X(i)); %得到对应点的函数值 end for j=2:6 f(1,1)= y(1); f(j, 1)=y(j); for k=2:j f(j,k)=(f(j,k-1)-f(j-1,k-1))/(X(j)-X(j-k+1)); %利用循环求差商 end end f %得到差商表，其中对角线上的为需要的差商值 syms x; l=[x,x,x,x,x,x]; g=l-X; n=g; for t=2:6 n(1,t)= n(1,t-1)* n(1,t); end N=zeros(1,1); N=f(1,1); for r=1:5 N= N +n(1,r)*f(r+1,r+1); %利用循环求5次牛顿插值多项式 end N %得到5次牛顿插值多项式 P=zeros(1,101); X1=[-1:0.02:1]; %取间距为0.02的等分点作为作图的横坐标 for u=1:101 x=X1(u); subs(N) ; p(1,u)= subs(N(1,1)) ; %用上面所求出的5次牛顿插值多项式算得所取等分点的函数值作为纵坐标 end plot(X1,p, -c) %绘制5次牛顿插值曲线 hold on XX=zeros(1,11) ; %定义一些下面将用到的矩阵 yy=zeros(1,11); F=zeros(11); L=zeros(1,11); for I=1:1:11