24小时热门版块排行榜    

查看: 1970  |  回复: 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的回帖

NicoleLas

木虫 (正式写手)

引用回帖:
6楼: Originally posted by foreverolia at 2014-09-10 20:18:39
高手,想跟您请教一下,plot函数画图的时候,为什么x=1:1:10和x=1:0.1:10得到的函数图像走势不同?应该是一样的呀。。。程序代码如下:clear all;
>> n=1:1:10;
>> y=(4.92-3.72)/(2.45-2.34)-((1.5 ...

大体趋势是一样的。
   不同处在于步长设置,第二个步长更细,采样点更多,反映的情况更加真实,你可以尝试以0.01为步长,更加接近真实情况。
Be water, my friend.
7楼2014-09-10 21:53:23
已阅   回复此楼   关注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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:59 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文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:10 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一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 22:39 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见