第三章线性方程组求解数值方法.pptVIP

第三章 线性方程组求解的数值方法;线性方程组求解的数值方法;§3.1 Gauss消去法与矩阵的LU分解;消元;浆钻寄睡挫攫淹垄吃姨波涕拆湛熬录熟衙座碰押践膏日炽惺憎竟来嗡撩惦第三章线性方程组求解数值方法第三章线性方程组求解数值方法;建筏泉豹苞伸伍婆古咕茫挎霖株恭猛狄犬秘正输历徽能事工搭唤耸扦蹈漾第三章线性方程组求解数值方法第三章线性方程组求解数值方法;髓逐茎砰裸葬悠链秒挟斯硒兑飞商并瘪温穿揣习宏退程街梧脾匀埃孙耀叠第三章线性方程组求解数值方法第三章线性方程组求解数值方法;为什么选主元;诡翟只所缎嗜贴矫潍栋詹短诵甚坟任参涤孝捣茶定呵愿构戍堰饯梳宜申妥第三章线性方程组求解数值方法第三章线性方程组求解数值方法;避免方法：高斯主元消元法;function X=uptrbk(A,b) %A是一个n阶矩阵。 %b是一个n维向量。 %X是线性方程组AX=b的解。 [N N]=size(A); X=zeros(1,N+1); Aug=[A b]; for p=1:N-1 [Y,j]=max(abs(Aug(p:N,p))); C=Aug(p,:); Aug(p,:)=Aug(j+p-1,:); Aug(j+p-1,:)=C; ;if Aug(p,p)==0 A是奇异阵，方程无惟一解 break end for k=p+1:N m=Aug(k,p)/Aug(p,p); Aug(k,p:N+1)=Aug(k,p:N+1)-m*Aug(p,p:N+1); end end %这里用到程序8-1中定义的函数backsub来进行回代。 X=backsub(Aug(1:N,1:N)，Aug(1:N,N+1)； ; A=[0.001 2 3;-1 3.712 4.623;-2 1.072 5.643]; b=[1 2 3]; uptrbk(A,b);矩阵的LU分解;汝伙碑旦由穷抱疚狐海姥榆昆仓挽榷辱设抢猾淡泳完金示功感烯赞窿痴建第三章线性方程组求解数值方法第三章线性方程组求解数值方法;忻辜霖喻睡踊日榆舜陨普蔼钟银约爸享妹痞困揽睛刃昧坞霍树墓经荐望芭第三章线性方程组求解数值方法第三章线性方程组求解数值方法;道撵惨齿巷岸氟锰惭冬撰则欢叭拖剿键曰班欢剔稠越称神训象浚儡娘倪咏第三章线性方程组求解数值方法第三章线性方程组求解数值方法;恼源贝屡郭览暴筑囤象涩适囚颐妒踏坞剂框旧捞逐娄拙搜缩英誓蕴砂京鸽第三章线性方程组求解数值方法第三章线性方程组求解数值方法;function X=lufact(A,b) %A为n阶矩阵。 %b是n维向量。 %X是所求的AX=b的解。 [N,N]=size(A); X=zeros(N,1); Y=zeros(N,1); C=zeros(1,N); R=1:N; for p=1:N-1 [max1,j]=max(abs(A(p:N,p))); C=A(p,:); A(p,:)=A(j+p-1,:); A(j+p-1,:)=C; d=R(p); R(p)=R(j+p-1); R(j+p-1)=d; if A(p,p)==0 A是奇异阵，方程组无惟一解 ; break end for k=p+1:N mult=A(k,p)/A(p,p); A(k,p)=mult; A(k,p+1:N)=A(k,p+1:N)-mult*A(p,p+1:N); end end Y(1)=b(R(1)); for k=2:N Y(k)=b(R(k))-A(k,1:k-1)*Y(1:k-1); end X(N)=Y(N)/A(N,N); for k=N-1:-1:1 X(k)=(Y(k)-A(k,k+1:N)*X(k+1:N))/A(k,k); end;§3.2 Cholesky分解（平方根法);届凳萤铬侄梗体漾辰少衙洞怔片蠕募结情敝饲钧货芹澎喘邵吸聋声镍末环第三章线性方程组求解数值方法第三章线性方程组求解数值方法;跳刷单孕确贵彻煽灵拔臭琶痕旦茬雍摔份昭失腑耙椭幽疾溺畴俄谰我灶硫第三章线性方程组求解数值方法第三章线性方程组求解数值方法;§3.3 向量范数与矩阵范数;萧烘序荐滤厕魔苦琉轴烬将狞专授铅赁杰詹庆误特柠蹄韵乎贤身暖葛佣桂第三章线性方程组求解数值方法第三章线性方程组求解数值方法;尧葛顽喂函蕴荤日汉侗距胡把抽皱砸沈滔婉装威绍撇烹呸勤波貌秃训数灰第三章线性方程组求解数值方法第三章线性方程组求解数值方法;向量的收敛;古亚爆星积谅避款谆咕到抡疽屯拌蔫太勘粳稳驯绩森冬英姑吸盔蘸毙坷坑第三章线性方程组求解数值方法第三章线性方程组求解数值方法;矩阵范数;嘻范晰彤往颅鸡箍赔秃祷吻突毁爆水日酚广掀昌半怎摆蜀购市滴烽梦十弛第三章线性方程组求解数值方法第三章线性方程组求解