南京郵電大學(xué)數(shù)值計算實踐報告

1、數(shù)值計算實踐I、方程求根、方程求根1、實驗?zāi)康膶嶒災(zāi)康氖煜ず驼莆誑ewton法割線法拋物線法的方法思路，并能夠在matlab上編程實現(xiàn)2、問題描述問題描述(1).給定一個三次方程分別用Newton法割線法拋物線法求解.方程的構(gòu)造方法:(a)根:方程的根為學(xué)號的后三位乘以倒數(shù)第二位加1再除以1000.假設(shè)你的學(xué)號為B06060141則根為141(41)1000=0.564(b)方程:以你的學(xué)號的后三位數(shù)分別作為方程的三次項二次項一次項的系

2、數(shù)根據(jù)所給的根以及三個系數(shù)確定常數(shù)項.例如:你的學(xué)號是B06060141則你的方程是x34x2xa0=0的形式.方程的根為0.564因此有0.564340.56420.564a0=0于是a0=2.015790144你的方程為x34x2x2.015790144=0.(2)假設(shè)方程是sinx4x2xa0=0的形式（三個系數(shù)分別是學(xué)號中的數(shù)字），重新解決類似的問題(3)構(gòu)造一個五次方程完成上面的工作.四次方程的構(gòu)造:將三次多項式再乘以(xp)

3、2得到對應(yīng)的五次多項式(p為已經(jīng)確定的方程的根顯然得到的五次方程有重根).(4)將（2）中的方程同樣乘以(xp)得到一個新的方程來求解注:(1)Newton法取0.5為初值割線法以01為初值拋物線法以00.51為初值(2)計算精度盡量地取高.終止準(zhǔn)則：根據(jù)來終止????||1nnpp(3)可供研究的問題：while(abs(b)eps1x(i))i=i1x(i)=x(i1)n_f(anx(i1))n_df(anx(i1))b=x(i)x

4、(i1)if(inn)errreturnendendy=x(i)建立n_f.m的源程序來求待求根的實數(shù)代數(shù)方程的函數(shù)，源程序代碼為：functiony=n_f(anx)%待求根的實數(shù)代數(shù)方程的函數(shù)y=0.0fi=1:(n1)y=ya(i)x^(n1i)end建立n_df.m的源程序來方程的一階導(dǎo)數(shù)的函數(shù)，源程序代碼為：functiony=n_df(anx)%方程的一階導(dǎo)數(shù)的函數(shù)y=0.0fi=1:ny=ya(i)(n1i)x^(ni)e