24小时热门版块排行榜    

查看: 1969  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

foreverolia

新虫 (初入文坛)

[求助] matlab求解非线性方程组报错,请各位大神指点! 已有1人参与

clear
>> syms n beta;
>> f1=sym(\'(f01-4*(2-n)*n*f01^(n-1)*sigma1^2*beta*yita1)/(1+2*n*(n-1)*f01^(n-2)*sigma1^2*beta*yita1)=fc1\')

f1 =

(f01 + 4*f01^(n - 1)*n*sigma1^2*yita1*(n - 2)*beta)/(2*f01^(n - 2)*n*sigma1^2*yita1*(n - 1)*beta + 1) = fc1

>> f2=sym(\'(f02-4*(2-n)*n*f02^(n-1)*sigma2^2*beta*yita2)/(1+2*n*(n-1)*f02^(n-2)*sigma2^2*beta*yita2)=fc2\')

f2 =

(f02 + 4*f02^(n - 1)*n*sigma2^2*yita2*(n - 2)*beta)/(2*f02^(n - 2)*n*sigma2^2*yita2*(n - 1)*beta + 1) = fc2

>> fo1=4.92;
>> fc1=3.76;
>> fo2=2.45;
>> fc2=2.31;
>> sigma1=1.50;
>> sigma2=0.635;
>> [n,beta]=solve(f1,f2)
Warning: Explicit solution could not be found.
> In solve at 81

n =

[ empty sym ]


beta =

     []
f1和f2是所求的方程组,算法总是报错,不知道该如何求解。。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

foreverolia

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by NicoleLas at 2014-09-09 10:18:27
看来你应该求下数值解,解析解MATLAB说找不到,或者你换用Maple或者Mathematica求下解析解。

我用plot绘图显示曲线和y=0有两个交点,理论上应该有两个解,可是为什么只能输出一个解。。。高手能不能帮忙写一下代码。。。
4楼2014-09-09 13:05:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

foreverolia

新虫 (初入文坛)

小女有一组数据要计算,但是大部分都显示帖子中的错误,不知道该如何纠正。。。求各位高手指点
2楼2014-09-09 09:19:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

看来你应该求下数值解,解析解MATLAB说找不到,或者你换用Maple或者Mathematica求下解析解。
Be water, my friend.
3楼2014-09-09 10:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NicoleLas

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
为图方便,给你用1stOpt编程算了下。
代码:
CODE:
Title "notsure";
Parameters n[0,1000,0],betta,yita1,yita2;
Function (4.92 + 4*4.92^(n - 1)*n*1.50^2*yita1*(n - 2)*betta)/(2*4.92^(n - 2)*n*1.50^2*yita1*(n - 1)*betta + 1) = 3.76;
         (2.45 + 4*2.45^(n - 1)*n*0.635^2*yita2*(n - 2)*betta)/(2*2.45^(n - 2)*n*0.635^2*yita2*(n - 1)*betta + 1) = 2.31;

采用L-M+UGO结果:
可能结果1:
n: 1
betta: 0.0879576864290957
yita1: 1.46535105823627
yita2: 0.986840117382106

可能结果2:
n: 1
betta: 0.235202634137847
yita1: 0.547990839309009
yita2: 0.369044224009309

上述结果都不一定正确,因为我不知道你的变量的取值范围,而且是否是连续或者整数变量,我假设n为非负整数变量,其余均为实数域连续变量了,你根据你的问题修改下吧,也可以采取其它算法试一下。
Be water, my friend.
5楼2014-09-09 18:18:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:40 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 04:32 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:29 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:23 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:56 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 01:24 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 00:59 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:12 by Z7qfCsdIiD1T
[考博] 申请环境,微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
信息提示
请填处理意见