上海大学课程实验-数值代数教程.docx

课程实验报告COURSE PAPER课程名称：数值代数与计算方法课程号： 授课教师：学号：姓名：所属：计算机科学与工程打印时间：2015评语：题目一：算法：对于问题一：对于问题二：直接编写递归函数程序，算出三个差分方程的10个近似值程序:1）%main.mclc;clear all;a=1;b=-2;c=-3;[x1,x2]=roots(a,b,c)%roots.mfunction [x1,x2]=roots(a,b,c)d=sqrt(b*b-4*a*c);if d0 x1=(-2*c)/(b+d); x2=(-b-d)/(2*a);elseif d0 x1=(-b+d)/(2*a); x2=(-2*c)/(b-d);end2）%solu.m function [X,R,P,Q]=solu(X0,R0,P0,P1,Q0,Q1)X(1)=X0;R(1)=R0;P(1)=P0;P(2)=P1;Q(1)=Q0;Q(2)=Q1;for i=1:9 X(i+1)=X(i)/2; X(i)=X(i+1);endfor i=1:9 R(i+1)=R(i)/2; R(i)=R(i+1);endfor i=3:10 P(i)=3/2*P(i-1)-1/2*P(i-2); P(i-1)=P(i); P(i-2)=P(i-1);endfor i=3:10 Q(i)=5/2*Q(i-1)-Q(i-2); Q(i-1)=Q(i); Q(i-2)=Q(i-1);end%Xn.mclc;clear all;X0=1;R0=0.994;P0=1;P1=0.497;Q0=1;Q1=0.497;[X,R,P,Q]=solu(X0,R0,P0,P1,Q0,Q1);for i=1:10 x(i)=i;endplot(x,X-R,r*);hold on结果：1）x1 = 3x2 =-12）题目二：算法：1、此题我们以第一个式子为例，求出g(Xn+1)= (2/(2+3*Xn-Xn.^3))^(1/2)2、利用二分法：利用牛顿迭代法，并根据题目提示的立方根算法：此题我们以第三个式子为例：，设A=71/3 ,求出x的3次方便是近似值。程序：1、%fixpt.mfunction [k,p,err,P]=fixpt(p0,tol,max1)P(1)=p0;for k=2:max1 P(k)=g(P(k-1)); err=abs(P(k)-P(k-1)); relerr=err/(abs(P(k))+eps); p=P(k); if (errtol)||(relerrtol), break; endendif k==max1 disp(maximum number of iterations exceeded)endP=P;%g.mfunction y=g(x)y=sqrt(2/(2+3*x-x.^3));%iteration.mclc;clear all;p0=0.5;tol=1.0e-9;max1=20;[k,p,err,P]=fixpt(p0,tol,max1);2、二分法：%bisect.mfunction [c,err,yc,k,an,cn,bn,ycn]=bisect(a,b,delta)ya=f(a);yb=f(b);an(1)=a;bn(1)=b;if ya*yb0 disp(No Root in this district!);endmax1=1+round((log(b-a)-log(delta))/log(2));for k=1:max1 c=(a+b)/2; cn(k)=c; yc=f(c); ycn(k)=yc; if yc==0 a=c; an(k+1)=a; b=c; bn(k+1)=b; elseif yb*yc0 b=c; bn(k+1)=b; yb=yc; else a=c; an(k+1)=a; ya=yc; end if b-adelta break; endendc=(a+b)/2;cn(max1)=c;err=abs(b-a);yc=f(c);ycn(max1)=yc;for i=1:9 k(i)=i-1;end%f.mfunction h=f(x)h=x*sin(x)-1;%main.mclc;clear all;a=0;b=2;de