數(shù)值分析課程設(shè)計(jì)報(bào)告(編輯修改稿)
2025-02-14 22:31 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 d C(i,:)=C(i,:)/C(i,i)。 for j=i+1:n C(j,:)=C(j,:)C(i,:)*C(j,i)。 endendC(n,:)=C(n,:)/C(n,n)。%因?yàn)檠h(huán)結(jié)構(gòu)問(wèn)題，所以C[n,:]并沒(méi)有化一，所以在最后才把最后一行化一for i=n:1:2 for j=1:i1 C(j,:)=C(j,:)C(i,:)*C(j,i)。 endend%通過(guò)行變換，使得C的前n列為單位陣形式for i=n:1:1 t=flag(i)。 if t~=i p=C(:,t)。C(:,t)=C(:,i)。C(:,i)=p。 endend%通過(guò)行變化矩陣還原行的順序for i=1:n for j=1:n B(i,j)=C(i,n+j)。 endend%將C的后n列賦給B矩陣4． 程序運(yùn)行操作過(guò)程與輸出結(jié)果操作時(shí)間①給出所要求逆的矩陣A：②執(zhí)行inv_dm(A)③執(zhí)行inv(A)④檢驗(yàn)正確性顯然從數(shù)量級(jí)上我們可以容易判斷這是一個(gè)單位陣，負(fù)號(hào)的出現(xiàn)是由于計(jì)算機(jī)默認(rèn)誤差造成細(xì)小誤差。
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1