资源描述:
《矩阵与数值分析上机实习》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、矩阵与数值分析上机实习1.设,其精确值为.(1)编制按从大到小的顺序,计算的通用程序(2)编制按从小到大的顺序,计算的通用程序(3)按两种顺序分别计算并指出有效位数(编制程序时用单精度)(4)通过本上机题,你明白了什么从小到大,代码:%1---SN=%N=input('pleaseinputanumber(N>=2)')if(N<2)disp('wrongnumber')elseS=0;forj=2:1:NS=S+1/(j^2-1);enddisp('S:')disp(S)end结果pleaseinputanumber(N>=2)10^2N=100S:7.
2、4005e-001>>clearpleaseinputanumber(N>=2)10^4N=10000S:7.4990e-001>>clearpleaseinputanumber(N>=2)10^6N=1000000S:7.5000e-001>>从大到小代码:%1---SN=%eps('single')N=input('pleaseinputanumber(N>=2)')if(N<2)disp('wrongnumber')elseS=0;forj=N:-1:2S=S+1/(j^2-1);enddisp('S:')disp(S)end结果pleaseinpu
3、tanumber(N>=2)10^2N=100S:7.4005e-001>>clearans=1.1921e-007pleaseinputanumber(N>=2)10^4N=10000S:7.4990e-001>>clearans=1.1921e-007pleaseinputanumber(N>=2)10^6N=1000000S:7.5000e-001(4)计算的顺序影响结果。1.秦九韶算法。已知n次多项式,用秦九韶算法编写通用的程序计算函数在点的值,并计算在点的值.(提示:编写程序时,输入系数向量和点,输出结果,多项式的次数可以通过向量的长度来判断)代
4、码:A=input('请输入系数,由高次幂开始');n=input('请输入计算变量的值');len=length(A);val=zeros(len);val(1)=A(1);%%printf('len=%c',len)fori=2:1:len%disp(val(i-1))%disp(n)val(i)=val(i-1)*n+val(i);end%printf('¼ÆËã½á¹û£º%f',val(len))val(len)结果:请输入系数,由高次幂开始[73-511]请输入计算变量的值23ans=851691.分别用Gauss消元法和列主元消去法编程求解方
5、程组Ax=b,其中,.高斯消去法代码:A=[31-13000-10000-1335-90-1100000-931-100000000-1079-30000-9000-3057-70-500000-747-300000000-3041000000-50027-2000-9000-229];B=[-15,27,23,0,-20,12,-7,7,10];B=B.';%disp(A)%disp(B)C=[A,B];n=size(A,1);ra=rank(A);rc=rank(C);x=zeros(1,n);ifra~=rcdisp('nosolution')els
6、efori=1:1:(n-1)forj=i+1:1:ntemp=(-C(j,i))/C(i,i);fork=i:1:(n+1)C(j,k)=C(i,k)*temp+C(j,k);endendendendfori=n:-1:1ifi==ndisp(C(i,(n+1)))disp(C(i,(n+1))-sum(C(i,(i+1):n).*x((i+1):n)))disp(C(i,i))endx(i)=(C(i,(n+1))-sum(C(i,(i+1):n).*x((i+1):n)))/C(i,i);enddisp('solution')disp(x)结果:-2
7、.2051e-0029.4726e-0011.0410e+0007.4481e-002-2.6617e-0012.0020e-001-2.4242e-0022.3844e-0013.8439e-001列主元消去法代码:%t4---max%A=[31-13000-10000-1335-90-1100000-931-100000000-1079-30000-9000-3057-70-500000-747-300000000-3041000000-50027-2000-9000-229];B=[-15,27,23,0,-20,12,-7,7,10];B=B.';
8、%disp(A)%disp(B)C=[A,B];n=size(A,