| 查看: 3340 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
MATLAB求解一个比较复杂的方程
|
|||
|
a/(1.0+exp(-b*(0.003-c)))-a/(1.0+exp(b*c))=1.460; a/(1.0+exp(-b*(0.012-c)))-a/(1.0+exp(b*c))=4.388; a/(1.0+exp(-b*(0.030-c)))-a/(1.0+exp(b*c))=5.455; 有谁知道上面的三个方程怎么解，初学MATLAB，还不会求解这种方程，谁能指导下我，谢谢~ 实际我是想用一个函数去拟合我的三个点。最后转化成上面的求解。 |
回复此楼» 本帖已获得的红花（最新10朵）
lylpiao» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab求解非线性方程组
已经有16人回复
- matlab用多元线性最小二乘法求解多元线性方程的问题
已经有11人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 求助一个差分方程的MATLAB实现
已经有3人回复
- 【资源】用matlab求解超越方程
已经有13人回复
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
★ ★ ★
xiegangmai(金币+3): 辛苦了！用符号算法可以得到全部解！ 2011-07-01 16:53:22
xiegangmai(金币+3): 辛苦了！用符号算法可以得到全部解！ 2011-07-01 16:53:22
|
eq1='a/(1+x*y)-a/(1+x)=1.460'; eq2='a/(1+x*y^4)-a/(1+x)=4.388'; eq3='a/(1+x*y^10)-a/(1+x)=5.455'; [a,x,y]=solve(eq1,eq2,eq3) b=-log(y)/0.003 c=log(x)./b s1=a./(1.0+exp(-b.*(0.003-c)))-a./(1.0+exp(b.*c)) s2=a./(1.0+exp(-b.*(0.012-c)))-a./(1.0+exp(b.*c)) s3=a./(1.0+exp(-b.*(0.030-c)))-a./(1.0+exp(b.*c)) |
9楼2011-07-01 16:38:14
snowwinter10
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 4180.3
- 散金: 2012
- 红花: 1
- 帖子: 607
- 在线: 155.6小时
- 虫号: 930945
- 注册: 2009-12-20
- 性别: GG
- 专业: 教学论
★ ★
xiegangmai(金币+2): 谢谢参与！鼓励讨论交流！ 2011-06-30 11:04:10
xiegangmai(金币+2): 谢谢参与！鼓励讨论交流！ 2011-06-30 11:04:10
|
试着写了一下，但是求解的结果很奇怪，不知道是数值参数还是方法的问题。希望有所参考价值。本人对matlab也是新手，我用的是非线性方程组求解函数solve，源码和结果如下。 clc clear all s1='a/(1.0+exp(-b*(0.003-c)))-a/(1.0+exp(b*c))=1.460'; s2='a/(1.0+exp(-b*(0.012-c)))-a/(1.0+exp(b*c))=4.388'; s3='a/(1.0+exp(-b*(0.030-c)))-a/(1.0+exp(b*c))=5.455'; [a b c]=solve(s1,s2,s3) 结果 a = 0. -0. b = (NaN)+(NaN)*i (NaN)+(NaN)*i c = (NaN)+(NaN)*i (NaN)+(NaN)*i |
2楼2011-06-30 10:53:22
【答案】应助回帖
★ ★ ★ ★ ★
臭水沟(金币+5): 谢谢应助~~ 2011-07-01 19:55:02
臭水沟(金币+5): 谢谢应助~~ 2011-07-01 19:55:02
|
尝试了一下用符号求解，找不到解析解。 得到： Warning: Explicit solution could not be found. > In solve at 81 ans = [ empty sym ] 改用数值解： （1）先写方程函数： （2）在命令窗口输入： 得到： Solver stopped prematurely. fsolve stopped because it exceeded the function evaluation limit, options.MaxFunEvals = 300 (the default value). X = 6.1816 5.8800 0.0390 Fval = -1.4331 -4.2800 -5.1839 exitflag = 0 求解失败。 可能方程组非线性太强了。 [ Last edited by xiegangmai on 2011-6-30 at 11:11 ] |
3楼2011-06-30 11:03:09
【答案】应助回帖
★ ★
cyzhixue(金币+10): 真的是太感谢你了。。崇拜中。。真的和我的点拟合的一模一样啊。 2011-06-30 17:21:47
臭水沟(金币+2): 谢谢应助~~ 2011-07-01 19:55:13
cyzhixue(金币+10): 真的是太感谢你了。。崇拜中。。真的和我的点拟合的一模一样啊。 2011-06-30 17:21:47
臭水沟(金币+2): 谢谢应助~~ 2011-07-01 19:55:13
|
呵呵，应该可以了，换用'levenberg-marquardt'算法，增加迭代次数。 用fsolve求解，初值选取有很大关系，当初值选为[ 100; 100; 0 ]时， [ X, Fval, exitflag ] = fsolve( @cyzhixue, [ 100; 100; 0 ], options ) 所得结果： X = 10.8784 182.8811 0.0001 Fval = 1.0e-009 * -0.1962 -0.2380 0.0061 exitflag = 1 应该能满足精度要求了。 |
4楼2011-06-30 11:15:59
