matlab求解五式方程组 matlab求解四元一次方程组

matlab求解五式方程组matlab求解四元一次方程组

用MATLAB解五元一次方程

>> x=[775 800 825 850 875 ]';

y=[47.04 47.27 47.39 47.49 47.64 ]';

A=[ones(5,1) x x.^2 x.^3 x.^4];

xx=A\y

Warning: Matrix is close to singular or badly scaled.

Results may be inaccurate. RCOND = 1.439940e-018.

xx =

1.0e+003 *

-1.4332

0.0066

-0.0000

0.0000

-0.0000

>> rref([A y])

ans =

1.0e+003 *

0.0010 0 0 0 0 -1.4332

0 0.0010 0 0 0 0.0066

0 0 0.0010 0 0 -0.0000

0 0 0 0.0010 0 0.0000

0 0 0 0 0.0010 -0.0000

改用数值解：

1、编写函数：

function f=fun2(in)

x=in(1);y=in(2);z=in(3);t=in(4);

f(1)=((500-x)^2+(3300-y)^2+z^2)/(21.15-t)^2-((3200-x)^2+(3100-y)^2+z^2)/(17.95-t)^2;

f(2)=((300-x)^2+(200-y)^2+z^2)/(19.4833-t)^2-((3400-x)^2+(100-y)^2+z^2)/(16.8167-t)^2;

f(3)=((800-x)^2+(1600-y)^2+z^2)/(14.85-t)^2-((2500-x)^2+(1900-y)^2+z^2)/(10.2333-t)^2;

f(4)=((1400-x)^2+(2200-y)^2+z^2)/(13.2833-t)^2-((2300-x)^2+(2800-y)^2+z^2)/(14.7833-t)^2;

f(5)=((1700-x)^2+(700-y)^2+z^2)/(11.7667-t)^2-((2900-x)^2+(900-y)^2+z^2)/(11.7667-t)^2;

2、调用求解：

fsolve(@(x)fun2(x),x0);% x0可以根据实际确定出大致范围即可

是有根号解，其中一个。

其余的传到我的空间吧。

>> syms a % 定义 a 为符号变量 f=4*cos(a)-2*sin(a)-6*exp(-2*a); % 方程关系式 solve(f) % 解方程 ans = .86889441616212189601300889832946+0.*i