计算方法实验一，方程求根

分别用Matlab和C写了一下，初学Matlab，如有不足还请指正。

实验内容：
1）在区间[0,1]内用二分法求方程 $e^x+10*x-2$ 的近似根，要求误差不超过 $0.5\times10^{-3}$ 。
2）取初值 $x_0=0$ ，用迭代公式 $x_{k+1}=\frac{2-e^{x_k}}{10},(k=0,1,2,...)$ 求方程 $e^x+10*x-2$ 的近似根，要求误差不超过 $0.5\times10^{-3}$ 。
3）取初值 $x_0=0$ ，用牛顿迭代法求方程 $e^x+10*x-2$ 的近似根，要求误差不超过 $0.5\times10^{-3}$ 。
牛顿迭代公式： $x_{k+1}=x_k-\frac{f(x_k)}{f'(x_k)}$ 。

1. C

#include <stdio.h>
#include <Windows.h>
#include <math.h>

#define WUCHA (0.5*pow(10, -4))

double f_x(double x);
double erfenfa();
double diedai();
double newtown();

int cnt = 0;

int main()
{
	// 二分法
	printf("二分法：%lf\n", erfenfa());
	printf("计算次数%d\n", cnt);

	// 迭代法
	printf("迭代法：%lf\n", diedai());
	printf("计算次数%d\n", cnt);

	// 牛顿法
	printf("牛顿法：%lf\n", newtown());
	printf("计算次数%d\n", cnt);

	system("pause");
	return 0;
}


// f(x)
double f_x(double x)
{
	return (exp(x) + 10 * x - 2);
}

// f'(x)
double fd_x(double x)
{
	return (exp(x) + 10);
}

// 牛顿法
double newtown()
{
	cnt = 0;

	double x0 = 0.0;
	double ret = f_x(x0);
	double xk = x0 - f_x(x0) / fd_x(x0);

	for (int k = 0; ; k++)
	{
		cnt++;
		if (fabs(ret - 0) <= WUCHA)
		{
			break;
		}

		xk = xk - f_x(xk) / fd_x(xk);
		ret = f_x(xk);
	}

	return xk;
}

// 迭代法
double diedai()
{
	cnt = 0;
	double x0 = 0.0;
	double ret = f_x(x0);
	double xk = x0;

	//int k = 0;
	for (int k = 0; ; k++)
	{
		cnt++;
		if (fabs(ret - 0) <= WUCHA)
		{
			break;
		}

		xk = ((2 - exp(xk)) / 10);
		ret = f_x(xk);
	}

	return xk;
}

// 二分法
double erfenfa()
{
	cnt = 0;
	double left = 0.0;
	double right = 0.1;
	double x = (left + right) / 2;
	double ret = f_x(x);

	while (fabs(ret - 0) >= WUCHA)
	{
		cnt++;
		if (ret > 0)
		{
			right = x;
		}
		if (ret < 0)
		{
			left = x;
		}
		x = (left + right) / 2;
		ret = f_x(x);
	}
	return x;
}

分界线

2. Matlab

2.1 二分法

function x = erfen(f,a,b,e)
% f是函数名 a,b是区间端点，e为精度
fa = feval(f,a); % feval是求函数值
fb = feval(f,b);
if fa*fb>0
    error('此函数在[%d,%d]区间无解',a,b);
end

k = 0;
x = (a+b)/2;
while(b-a)>(2*e)
    fx=feval(f,x);
    if fa*fx<0
        b=x;
        fb=fx;
    else
        a=x;
        fb=fx;
    end
    k=k+1;
    x=(a+b)/2;
    
    fprintf("%d:x=%f\n", k, x);
    
end

2.2 不动点迭代法

function x = diedai(f,diedaif,x0,e)

N = 10; %最大迭代次数

ret = feval(diedaif,x0);
xk = x0;
k = 0;

while abs(x0-ret)>e & k<N
    
    xk = feval(diedaif,xk);
    ret = feval(diedaif,xk);
    
    k = k+1;
    fprintf("xk=%f\n", xk);
    %if
    %    break;
    %end
    
end

if k==N
    fprintf('超过最大迭代次数限制\n');

end

2.3 牛顿法

function x = newton(f,fd,x0,e)
% f->f(x)
% fd->f'(x)
% x0=0
% e=0.0005

N = 100; %最大迭代次数100
x = x0;
x0 = x+2*e;
k = 0;
while abs(x0-x)>e & k<N
    k = k+1;
    x0=x;
    x=x0-feval(f,x0)/feval(fd,x0);
    fprintf("x=%f\n", x);
end

if k==N
    warning('超过最大迭代次数限制');

end