应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 【求助】方程组求解
51/1返回列表
查看: 326  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lshuihui62

铜虫 (小有名气)

[交流] 【求助】方程组求解

非线性方程组（见图）：其中Dζ、Dγ、Dβ、Yζ、Yγ、Yβ为未知数，对应到程序里面依次为x(1)...x(6)，其余均为常数，我想问一下在MATLAB中，可不可以用解一般非线性方程组的方法解这个非线性方程组，下面是我的程序代码，但是没有解出来，我想请教一下我的问题在哪里？未知数的最小值是零，但是为确保分母不为零，初值就不能赋为零。
function F=myfun(x)
F=[((7.22-3.56)*sqrt(x(1))*exp(-(x(4))^2/(4*x(1))))/(erf((x(4))/(2*sqrt(x(1)))))-(sqrt(pi)*x(4)*(3.56-3.05)/2)-(((3.05-0.924)*sqrt(x(2))*exp(-(x(4))^2/(4*x(2))))/(erf((x(5))/(2*sqrt(x(2))))-erf((x(4))/(2*sqrt(x(2))))));
((0.263-0.016)*sqrt(x(3))*exp(-(x(5))^2/(4*x(3))))/(erf((x(5))/(2*sqrt(x(3)))-erf((x(6))/(2*sqrt(x(3))))))-(sqrt(pi)*x(5)*(0.924-0.263)/2)-(((3.05-0.924)*sqrt(x(2))*exp(-(x(5))^2/(4*x(2))))/(erf((x(4))/(2*sqrt(x(2))))-erf((x(5))/(2*sqrt(x(2))))))];

>> x0=[];
>> options=optimset('Display','iter');
>> [x,fval]=fsolve(@myfun,x0,options)
回复此楼

» 猜你喜欢
1楼 2009-06-01 15:54:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lshuihui62

铜虫 (小有名气)
我想问一下这个非线性方程组可不可以用一般的解方程命令求出其通解？
具体的用哪个命令比较好一些？
赞一下 回复此楼
4楼2009-06-07 10:55:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

q68

铁杆木虫 (著名写手)
这个方程看起来挺复杂的，等高手了！
赞一下 回复此楼
2楼2009-06-01 18:34:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭