| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 2631 人关注过本帖
标题:求助:fsolve求解非线性方程组的问题
只看楼主 加入收藏
zjjwlx
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2006-9-3
收藏
 问题点数:0 回复次数:4 
求助:fsolve求解非线性方程组的问题

需要求解的方程组如下:
function y=fc(x)
pr=15e-6
ps=17e-6
ec=100e3
ef=200
tf=150e-7
ti=2e-7
q=1.602e-19
k=1.380e-23
t=300
eo=8.854e-14
esi=11.9
nd=1e16
ni=1.5e10
ei=3.9
dat=ec/log((1+pr/ps)/(1-pr/ps))
bat=q/(k*t)
ld=sqrt(esi*eo/(q*nd*bat))

y(1)=-sqrt(2)*esi*eo/(bat*ld)*sqrt((ni^2)/(nd^2)*(exp(-bat*x(1))+bat*x(1)-1)+exp(bat*x(1))-bat*x(1)-1)+ps*tanh((x(2)-ec)/(2*dat))
y(2)=x(1)+sqrt(2)*esi*eo/(bat*ld)*sqrt((ni^2)/(nd^2)*(exp(-bat*x(1))+bat*x(1)-1)+exp(bat*x(1))-bat*x(1)-1)*ti/(ei*eo)+x(2)*tf-5

y=[y(1) y(2)]


其中x(1)的范围为1左右,x(2)的范围在3e5左右。

我解的时候是这样的:
>>x0=[1 3e5];
>>fsolve('fc',x0)
运行以后,最后显示的是:
Optimization terminated: first-order optimality is less than options.TolFun.

ans =

1.0e+005 *

0.0000 1.0003

结果明显是不对的。方程是没有问题的。大家能够帮我解决一下这个问题吗??我很着急啊~~

搜索更多相关主题的帖子: 方程 fsolve 非线性 求解 
2006-09-07 10:02
sxs77
Rank: 1
等 级:新手上路
帖 子:5
专家分:0
注 册:2006-9-7
收藏
得分:0 

兄弟,我和你遇到同样的问题了,这个问题我也琢磨了两天了,解决不了,还望高手赐教!

2006-09-07 15:55
zjjwlx
Rank: 1
等 级:新手上路
帖 子:3
专家分:0
注 册:2006-9-3
收藏
得分:0 

盼望中。。。

2006-09-07 18:52
冰山一角
Rank: 1
等 级:新手上路
帖 子:385
专家分:0
注 册:2006-9-5
收藏
得分:0 
?????

做程序员太乏味?来这里www..cn试试吧,你肯定能找到乐趣!
2006-09-08 11:30
zt_chem
Rank: 1
等 级:新手上路
帖 子:2
专家分:0
注 册:2011-6-29
收藏
得分:0 
我也遇到类似的问题,期待求解
2011-06-29 16:26
快速回复:求助:fsolve求解非线性方程组的问题
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.029458 second(s), 9 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved