数值分析方程求根实验matlab,matlab数值分析(非线性方程求根)上机实验报告.doc

数值分析试验报告

非线性方程求根

二分法解方程

题目：用二分法求方程的所有根

方法：二分法

程序

function x=erfenfa(a,b)

if (a*a*a-3*a-1)*(b*b*b-3*b-1)>0

disp('区间错误，请重新调试！');

break

else

x=(a+b)/2.0;

k=0;

while (k<=100&abs(x*x*x-3*x-1)>0.0001&b-a>0.0001)

if (x*x*x-3*x-1)*(a*a*a-3*a-1)<0

b=x;

else

a=x;

end

x=(a+b)/2.0;

k=k+1;

k,x

if k==100

display('计算方法错误，请重新调试！');

end

end

end

disp(['k=',num2str(k)])

结果

>> format compact

>> erfenfa(1,2)

k =

1

x =

1.75000000000000

k =

2

x =

1.87500000000000

k =

3

x =

1.93750000000000

k =

4

x =

1.90625000000000

k =

5

x =

1.89062500000000

k =

6

x =

1.88281250000000

k =

7

x =

1.87890625000000

k =

8

x =

1.88085937500000

k =

9

x =

1.87988281250000

k =

10

x =

1.87939453125000

k=10

ans =

1.87939453125000

>> erfenfa(-1,1)

k =

1

x =

-0.50000000000000

k =

2

x =

-0.25000000000000

k =

3

x =

-0.37500000000000

k =

4

x =

-0.31250000000000

k =

5

x =

-0.34375000000000

k =

6

x =

-0.35937500000000

k =

7

x =

-0.35156250000000

k =

8

x =

-0.34765625000000

k =

9

x =

-0.34570312500000

k =

10

x =

-0.34667968750000

k =

11

x =

-0.34716796875000

k =

12

x =

-0.34741210937500

k =

13

x =

-0.34729003906250

k=13

ans =

-0.34729003906250

>> erfenfa(-2,-1)

k =

1

x =

-1.75000000000000

k =

2

x =

-1.62500000000000

k =

3

x =

-1.56250000000000

k =

4

x =

-1.53125000000000

k =

5

x =

-1.54687500000000

k =

6

x =

-1.53906250000000

k =

7

x =

-1.53515625000000

k =

8

x =

-1.53320312500000

k =

9

x =

-1.53222656250000

k =

10

x =

-1.53173828125000

k =

11

x =

-1.53198242187500

k =

12

x =

-1.53210449218750

k=12

ans =

-1.53210449218750

5，拓展：在程序编写完成调试过程中发现如果给出错误区间程序仍能给出一个错误答案，于是找同学讨教，讨论得出了正确的方法

不动点解方程

题目： 解方程的根

方法：不动点迭代法

程序

function x=budongdian(x0)

p=x0;

for k=1:10

x0=p;

p=20/(x0*x0+4*x0+10);

if a