小木虫

当前位置：首页 > 计算模拟 >想用matlab解一个方程组谁能帮帮我呢。

想用matlab解一个方程组谁能帮帮我呢。

数学小白，现遇到一个方程组不会解，希望能得到高人指点，必有重谢。

IMG_0607.JPG 返回小木虫查看更多

今日热帖

yaohonghao

消去A，算出B，C
sjt19941213
过氧酰硝酸酯

引用回帖:
2楼: Originally posted by yaohonghao at 2017-04-17 20:24:50
消去A，算出B，C

不要开玩笑，哪有那么简单
dingd

1stOpt试试：
CODE:
Constant y=[4.05,3.87,3.77],x=[1.0613,1.20116,1.3392];
Function For(i=1:3,x,y)(A/(1+B*exp(C/(x*10^(-21))))=y);

a: 3.43657094484967
b: -0.0113334833952189
c: 2.75151179795343E-21
过氧酰硝酸酯

引用回帖:
5楼: Originally posted by dingd at 2017-04-17 22:06:46
1stOpt试试：

Constant y=,x=;
Function For(i=1:3,x,y)(A/(1+B*exp(C/(x*10^(-21))))=y);

a: 3.43657094484967
b: -0.0113334833952189
c: 2.75151179795343E-21...

你好你能加下我QQ么我好详细咨询一下您
593989702
smqh5207

[ 164.87993147 41.31307048 83.37959768]
from math import exp
from scipy.optimize import fsolve
def f(x):
x0 = float(x[0])
x1 = float(x[1])
x2 = float(x[2])
return [
      4.05-x0/(1+x1*exp(x2/1.0631*1e-21)),
      3.87-x0/(1+x1*exp(x2/1.20116*1e-21)),
      3.77-x0/(1+x1*exp(x2/1.3392*1e-21))
      ]
result = fsolve(f, [1,1,1])
print result，
smqh5207

三个误差：[0.15333332902489438, -0.026666670975105333, -0.12666667097510542]。可以考虑其他算法，使得误差最小，比如marqut，wegstein加速，牛拉算法等，多了去了