随笔分类 - 数值分析
摘要:A=[1,-1,1,-4;5,-4,3,12;2,1,1,11;2,-1,7,-1] L=eye(length(A)) %开始消元过程 for k=1:(length(A)) a=A(k,k) for i=k+1:(length(A)) c=-A(i,k) L(i,k)=-c./a for j=1: (length(A)) A(i,j)=A(i,j)+c....
阅读全文
摘要:A=[1,-1,1,-4;5,-4,3,12;2,1,1,11;2,-1,7,-1] Adet=1 %开始消元过程 for k=1:(length(A)) a=A(k,k) Adet = Adet.*a for i=1:(length(A)) A(k,i)=A(k,i)/a end for i=k+1:(length(A)) c=-A(i,k...
阅读全文
摘要:e=sym('e') E=[1,0,0,0;0,1,0,0;0,0,1,0;0,0,0,1] E=e.*E A=[1,-1,1,-4;5,-4,3,12;2,1,1,11;2,-1,7,-1] A=A-E Adet=1 %开始消元过程 for k=1:(length(A)) a=A(k,k) Adet = Adet.*a for i=1:(length(A)) A(...
阅读全文
摘要:A=[1,-1,1,-4;5,-4,3,12;2,1,1,11;2,-1,7,-1] E=[1,0,0,0;0,1,0,0;0,0,1,0;0,0,0,1] b=[2,4,3,0]' %开始消元过程 for k=1:(length(A)) a=A(k,k) for i=1:(length(A)) A(k,i)=A(k,i)/a E(k,i)=E(k, i)/a...
阅读全文
摘要:a=1.0 b=1.5 f=sym('x^4-x-2') e=eps(realmin).*10 while (b-a)>0.0001 c=(a+b)/2 if sign(subs(f,'x',c))==sign(subs(f,'x',a)) a=c else b=c end end a = 1 b = 1.5000 f = x^4 - x - ...
阅读全文

浙公网安备 33010602011771号