數(shù)值分析上機(jī)(四)資料

1、Author by prs,版權(quán)所有Ps：題目均來自數(shù)值分析第五版作者：李慶揚(yáng)，王能超，易大義 編 出 版 社：清華大學(xué)出版社誤差分析 誤差分析問題：求下列方程的實(shí)根（1） 2 x 3 2 0 x x e ? ? ? ?（2） 3 2 2 +10x 20 0 x x ? ? ?要求：（1）設(shè)計(jì)一種不動點(diǎn)迭代法，要使迭代序列收斂，然后再用斯特芬森加速迭代，計(jì)算到 為止。 （2）用牛頓迭代，同樣計(jì)算到 ，輸出迭 81 | x x | 10

2、k k?? ? ? 81 | x x | 10 k k?? ? ?代初值及各次迭代值和迭代次數(shù) k，比較方法的優(yōu)劣。代碼部分：/**函數(shù)**/function y = fun(x)y=x^3+2*x^2+10*x-20;endfunction y = fun1( x)y=x^2-3*x+2-exp(x);% y=2*log(x)+log(3);endfunction [y,k]= niudun(x0)%NUIDUN Summary of

3、 this function goes here% Detailed explanation goes herex(1)=x0;k=1;des=1;while des>1.0e-8x(k+1)=x(k)-fun1(x(k))/dfun1(x(k));des=abs(x(k+1)-x(k));k=k+1;endy=x(k);k=k;endfunction [y,k]= sitefensen(x0,f)%SITEFENSEN S

4、ummary of this function goes here% Detailed explanation goes here%x0 為初值，n 為迭代次數(shù)，f 為迭代函數(shù)x(1)=x0;Author by prs,版權(quán)所有disp('不動點(diǎn)迭代解->')fprintf('%f\n',x)disp('迭代次數(shù)->')fprintf('%d\n',k)di

5、sp('誤差->')fprintf('%f\n',abs(0-fun1(x)))%斯蒂芬森加速f='(x^2+2-exp(x))/3';f=inline(f);[yy,kk]=sitefensen(0.5,f);disp('斯蒂芬森加速解->')fprintf('%f\n',yy)disp('迭代次數(shù)->')fprintf(

6、'%d\n',kk)disp('誤差->')fprintf('%f\n',abs(0-fun1(yy)))%牛頓迭代法[yy1,kk1]=niudun(0.5);disp('牛頓迭代法解->')fprintf('%f\n',yy1)disp('迭代次數(shù)->')fprintf('%d\n',kk1)disp(&