全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
17323 2
2015-04-01
这里有两个关于未知数c1和c2的方程,其中一个方程只包含未知数c1,另一个方程包含c1和c2,我单独解第一个方程的话可以得到c1的解,但是第二个方程就解不出结果了,请高手指点一下。
方程一是:1.05157-1.222668+.0048063*(1-exp(-.479452*c1))+.0098652*(1-exp(-.9835616*c1))+.0147561*(1-exp(-1.479452*c1))+.0196046*(1-exp(-1.983562*c1))+.0242522*(1-exp(-2.479452*c1))+.0288602*(1-exp(-2.983562*c1))+0.4*.9673068*(1-exp(-2.983562*c1))=0
方程二是:1.12898-3.341505+.0076901*(1-exp(-.479452*c1))+.0157843*(1-exp(-.9835616*c1))+.0236098*(1-exp(-1.479452*c1))+.0313674*(1-exp(-1.983562*c1))+.0388035*(1-exp(-2.479452*c1))+.0461763*(1-exp(-2.983562*c1))+.0532934*(1-exp(-2.983562*c1-.4986301*c2))+.0603177*(1-exp(-2.983562*c1-1.00274*c2))+.0670641*(1-exp(-2.983562*c1-1.49863*c2))+.0737597*(1-exp(-2.983562*c1-2.00274*c2))+.0801891*(1-exp(-2.983562*c1-2.49863*c2))+.0865684*(1-exp(-2.983562*c1-3.00274*c2))+.0926922*(1-exp(-2.983562*c1-3.49863*c2))+.0987661*(1-exp(-2.983562*c1-4.00274*c2))+.1046266*(1-exp(-2.983562*c1-4.50137*c2))+.1104049*(1-exp(-2.983562*c1-5.005479*c2))+0.4*.8637214*(1-exp(-2.983562*c1-5.005479*c2))=0

下面是我的程式:>> c1=solve('1.05157-1.222668+.0048063*(1-exp(-.479452*c1))+.0098652*(1-exp(-.9835616*c1))+.0147561*(1-exp(-1.479452*c1))+.0196046*(1-exp(-1.983562*c1))+.0242522*(1-exp(-2.479452*c1))+.0288602*(1-exp(-2.983562*c1))+0.4*.9673068*(1-exp(-2.983562*c1))','c1')

c1 =

0.15438446386393466762440363354485
>> c2=solve('1.12898-3.341505+.0076901*(1-exp(-.479452*c1))+.0157843*(1-exp(-.9835616*c1))+.0236098*(1-exp(-1.479452*c1))+.0313674*(1-exp(-1.983562*c1))+.0388035*(1-exp(-2.479452*c1))+.0461763*(1-exp(-2.983562*c1))+.0532934*(1-exp(-2.983562*c1-.4986301*c2))+.0603177*(1-exp(-2.983562*c1-1.00274*c2))+.0670641*(1-exp(-2.983562*c1-1.49863*c2))+.0737597*(1-exp(-2.983562*c1-2.00274*c2))+.0801891*(1-exp(-2.983562*c1-2.49863*c2))+.0865684*(1-exp(-2.983562*c1-3.00274*c2))+.0926922*(1-exp(-2.983562*c1-3.49863*c2))+.0987661*(1-exp(-2.983562*c1-4.00274*c2))+.1046266*(1-exp(-2.983562*c1-4.50137*c2))+.1104049*(1-exp(-2.983562*c1-5.005479*c2))+0.4*.8637214*(1-exp(-2.983562*c1-5.005479*c2))','c2')
Warning: Explicit solution could not be found.
> In solve at 179

c2 =

[ empty sym ]

我试过了fsolve函数,但是计算出来的结果跟上面计算出来的c1差太多,理论上c1和c2都是零点几的数字。
这个问题困扰了我好多天,实在想不出来该怎么解了。
请问高手,我要解这两个方程该怎样改程式啊?谢谢
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2015-4-17 16:45:58
function y = myfunc( c )
y = 1.05157-1.222668+.0048063*(1-exp(-.479452*c(1)))+.0098652*(1-exp(-.9835616*c(1)))+...
    .0147561*(1-exp(-1.479452*c(1)))+.0196046*(1-exp(-1.983562*c(1)))+...
    .0242522*(1-exp(-2.479452*c(1)))+.0288602*(1-exp(-2.983562*c(1)))+...
    0.4*.9673068*(1-exp(-2.983562*c(1)));
end

命令窗口:
[x,fval] = fsolve(@myfunc,0.1)

结果:
x =
    0.1544
fval =
  -6.9046e-10
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2015-4-17 16:53:25
第二个函数, 当c1确定为0.15438446386393466762440363354485时, 该函数关于c2似乎没有零点, 当c2趋于无穷时, 该函数的似乎图像趋于-1. 你可以试着画一下第二个函数关于c2的图像.
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群