版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
YN結(jié)束本次作業(yè)所使用的編譯環(huán)境為:Visualc++6.0doubleA[10][10],rr[10][10],qq[10][10],rq[10][1//對矩陣A進(jìn)行初始化,賦值{{A[i-1][j-1]=sin(0.5*i+0.2A[i-1][i-1]=cos(i+1.2*i)*}//定義函數(shù)對矩陣A進(jìn)行擬上三角化{intsgn(doublea);doubled,c,h,t,tr,u[10],p[10],q[10],w[10];{for(i=r+2;i<10;i++)t=t+(A[i]{{}}for(i=0;i<10;i++)sum=sum+p[i{{A[i][j]=A[i][j]-w[i]*u[j]-u[i]}}}}}}//sgn子程序intsgn(doublea){}//計(jì)算A的特性值的函數(shù)(書本P6311步)inttezhengzhi(){voidqmk(doublemk[10][10],intm,doublems,doublemt);//申明求mk的函數(shù)voidQR(doublemk[10][10],intm);//申明對mk進(jìn)行qr分解的函數(shù)detdk=A[m-1][m-1]*A[m][m]-A[m-1][fb=-1.0*(A[m-1][m-1]+A[m{}{s1_r=(-1.0*fb+sqrt(s2_r=(-1.0*fb-sqrt(}}{}{if(-1e-12<A[m-1][m-2]&&{{}elseif(m==-1)}{{mt=A[m-1][m-1]*A[m][m]-}}{doubletmp_sum;{{for(k=0;k<=m;k++)tmp_sum+=A[i][k]*}}for(i=0;i<=m;i++)mk[i][}//對mk做QR分解并得到A(k+1)的函數(shù){doubleur[10],vr[10],pr[10],qr[10],wr[10],tmp_dr,tmp_vr,tmp_pr,tmp_qr,tmp_tr,dr,ctmp_ir=tmp_ir+(mk[i]{234for(i=r;i<=m;i++)tmp_dr=tmp_dr+mk[i][r]*mfor(i=r+1;i<=m;i++)ur[i]=mktmp_vr=tmp_vr+mk[j][itmp_pr=tmp_pr+A[j][itmp_qr=tmp_qr+A[i][j}}for(i=0;i<=m;i++)tmp_tr=tmp_tr+pr[i]for(i=0;i<=m;i++)wr[i]=qr[i]-tr{{A[i][j]=A[i][]-wr[i]*ur[j]-ur[i]}}}{}{}voidtranslation(doub{{{}}}//列主元的高斯消元法求解特性向量{for(k=0;k<9;k++)/*選主元并消元*/{mk=k;//記錄最大組員的行號for(i=k;i<10;i++)//*選擇第K列主元素{{}if(mk!=k)/*將第K列主元素?fù)Q行到對角線上*/{{}}{A[i][j]=A[i][j]-mik*A[k}}/*回代求解*/x[9]=b[9]/A[10-1][1{{}}}{{{}{}{}}while(fabs(beita2-beita1)>=fabs(beita2)*exp(-1for(m=0;m<10;m++)xy[m]}//主程序{{{}}printf("矩陣A各個(gè)特性值的實(shí)部及虛部為:\n\n");{}printf("經(jīng)帶雙步位移的QR分解后的矩陣為:\n");{{if(-1.0e-12<A[i][j]&&A[i][j}}{{{}{}}}1.通過上三角化后的矩陣經(jīng)過上三角化后的矩陣為255889387184e-8016.83259944534e-11.51886895646939998e+008-6.388578124485e-001-4.15487338868395-0819.783997%21285e-0822.557763574.08836·B.88e·8-1.21695343e+8-4.11683513396+E?1+e?-41.3W?583821e+e?1.883691775ke-1-4.246335358369e-8817.988955239384-0821.68819928869e-.08368e-B.8e:·88.88838e-81.5681262522-·88.12584937515e-14.421756832232-1-1.588616128136718-01-2.7365958588920-801-7.32446835-081-6.7099254814498-0812.5446199298821893682e-0811.19587191495e-0811.9298618e-014.9882594688e-011.1376916837762.0089008-88.8808908-80.00930+08B.0808+088.0008008-08B.08900-088.0808908-087.031672187e-B?11.582188688475e-BP13.86259461422.88888+8B.88B8908+88.88888+08B.808+88.088-8B.8388e+88.088e+B?84.843887682783e-B?13.992777995-8.e-001-9.e-002-1.e+000-7.e-0011.e-001-2.e+0002.e+0001.8e+0003.e-0012.e-0012.e+0001.8e-0011.e+000-6.e0.e+0001.e+000-1.5e+000-1.e+000-6.e-001-2.e+000.e+0000.e+000-1.e+000-0.e+0000.e+0000.e+0001.e+0008.e-0014.e-001-3.e-0024.e-0010.e+0000.e+0000.e+0000.e+0000.e+000-7.e-001-2.e-002-90.e+0000.e+0000.e+0000.e+0000.e+0000.e+000-7.e-001第1個(gè)特征值為:(3.383839617436e+880,a.098888888888e+888i)第2個(gè)特征值為:(-2.323496218212e+880,8.938405177208e-881i)第3個(gè)特征值為:(-2.323496218212e+880,-8.9304851?7200e-881i)第4個(gè)特征值為:(1.577548557113e+088,a.000808808880e+080i)第5個(gè)特征值為:(-1.484839822259e+080,a.888808808888e+888i)第6個(gè)特征值為:(-9.805309562982e-881,1.139489127438e-801i)第?個(gè)特征值為:(-9.805309562902e-Q81,-1.139489127430e-Q01i)第8個(gè)特征值為:(9.355889878188e-O?1,0.000eene?8800e+000i)第9個(gè)特征值為:(6.360627875746e-Q81,a.80080880888Ce+080i)第10個(gè)特征值為:(5.650488993581e-882,a.0008Be008088e+000i)5218517e-802-1.882578788827e+088-4.886268122088BP?e·-2.184744213e+8-2.3615279847Ke+8-1.4561256863c-02-4.7641314:21.81639532183e·8-2.318975622He5516812783e-081-6.537876947738e-B?13.227152128265524e-BE]-4.154391785514-082-1.861198124719-801.04875636749-01-22887e-BR1-5.962368585233e-081-4.9115922718343047e-822.98478806855e-081-3.57W39888185223188e-811.877411851833e-8011.36825651900ke+Bee6.368511809728e-012.737834413363e-0810+02.45769935855-05.6516496536請證值3.3309617406e-對應(yīng)的特征向量為:(2.3838426916e4.363067W26e-+82.96232952853.41鉗鋤4924-3.1922812991e+3.4685172特征值1.57548559113e+對應(yīng)的特征向量為:(5.7552638375e-12.5274889591+088.788449279494-1.77693466823e-8884418782+001.95284971187+4E82.8319869547a+0081.94129特征值-1.錫403982259++00對應(yīng)的特征向量為:(-2.4848382252+085.091259892-01-1.482753484281+8-1.48S8532881e833998777+0閥-1.484361479621+088-1.銅058559854+08-1.48482423030+000-1.4398331697%+085.64354827408-H?2)特征值9.35897818e-謝對團(tuán)的特征向量為:(1.935588251875+1.9附42405381-839.26681862717-謝9.227378151e-813592N?4+01+8981.06133463+-l8.416198235e-17.528673539e-4.86929持癥值5.58489935Mle:對應(yīng)的特證向?yàn)椋?1.8564124521特性值9.e-001對應(yīng)的特性向量為:(1.e+0001.e-0039.e-0019.e-001特性值6.e-001對應(yīng)的特性向量為:(1.e+000-2.e-0016.e-0015.e-0016.e-0018.e-0017.e
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 校外課程設(shè)計(jì) 英語
- 2022年廣東韶關(guān)新豐縣一年級上學(xué)期六單元語文試卷及答案
- 2021-2022學(xué)年江蘇省南京市棲霞區(qū)小學(xué)三年級上冊語文期末試題及答案
- 2020年浙江省杭州市椒江區(qū)六年級下冊期末語文試卷及答案
- 2024年標(biāo)準(zhǔn)墊資服務(wù)合作合同樣本版B版
- 2024年度知識產(chǎn)權(quán)抵押貸款合同大全3篇
- 2025年高考?xì)v史復(fù)習(xí)熱搜題速遞之秦漢時(shí)期(2024年7月)
- 2025年小升初復(fù)習(xí)之小題狂練300題(多句選詞):名詞(10題)
- 2022-2023學(xué)年遼寧沈陽鐵西區(qū)五年級上冊語文期末試卷及答案
- 游泳腳步課程設(shè)計(jì)
- 系統(tǒng)運(yùn)行維護(hù)方案
- 轉(zhuǎn)子找靜平衡方法
- 2025年九省聯(lián)考新高考 政治試卷(含答案解析)
- 終極戰(zhàn)略規(guī)劃指南:深度剖析Cross SWOT分析、市場洞察與內(nèi)部能力優(yōu)化的綜合行動(dòng)方案
- 中國偏頭痛診治指南(第一版)2023解讀
- 湖北省武漢市黃陂區(qū)2024年數(shù)學(xué)六年級第一學(xué)期期末學(xué)業(yè)質(zhì)量監(jiān)測模擬試題含解析
- 關(guān)于開展2024年度保密自查自評專項(xiàng)檢查工作的實(shí)施方案
- 商場反恐防暴應(yīng)急預(yù)案演練方案
- 2024年天津市西青經(jīng)濟(jì)開發(fā)集團(tuán)限公司公開招聘工作人員高頻500題難、易錯(cuò)點(diǎn)模擬試題附帶答案詳解
- 數(shù)據(jù)庫設(shè)計(jì)規(guī)范標(biāo)準(zhǔn)
- 2023年全國職業(yè)院校技能大賽賽項(xiàng)-ZZ019 智能財(cái)稅基本技能賽題 - 模塊三-答案
評論
0/150
提交評論