找回密码
 注册
Simdroid-非首页
查看: 53|回复: 0

方程求解

[复制链接]
发表于 2011-9-5 15:36:54 | 显示全部楼层 |阅读模式 来自 江苏南京
方程求解如下:
clc;
clear;
format long;
XY=0.164984207E+07;
%XY=0.205921093E+07;
%XY=0.257842948E+07;

x1=1E20;
y1=(20.384E+06)*(x1.^(-0.2))+(735.392E+06)*(x1.^(-0.98))-XY;
x2=1;
y2=(20.384E+06)*(x2.^(-0.2))+(735.392E+06)*(x2.^(-0.98))-XY;
x0=(x1+x2)/2;
y0=(20.384E+06)*(x0.^(-0.2))+(735.392E+06)*(x0.^(-0.98))-XY;
%****************************************************
while(abs(y0)>1E-10)
    if(y0>0)
        x2=x0;
        y2=(20.384E+06)*(x2.^(-0.2))+(735.392E+06)*(x2.^(-0.98))-XY;
    else
        x1=x0;
        y1=(20.384E+06)*(x1.^(-0.2))+(735.392E+06)*(x1.^(-0.98))-XY;
    end
   
    x0=(x1+x2)/2;
    y0=(20.384E+06)*(x0.^(-0.2))+(735.392E+06)*(x0.^(-0.98))-XY;
end
%****************************************************
x0

在方程求解中,对应XY的值
XY=0.164984207E+07;
XY=0.205921093E+07;
XY=0.257842948E+07;
第一个和第三个都能求解出x0,而第二个却不能求解。Matlab中始终显示“busy”,请帮忙指出错误所在!谢谢!

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 15:03 , Processed in 0.025986 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表