这里有两个关于未知数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都是零点几的数字。
这个问题困扰了我好多天,实在想不出来该怎么解了。
请问高手,我要解这两个方程该怎样改程式啊?谢谢