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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » MATLAB求解一个比较复杂的方程
152/2返回列表上一页12
查看: 3472  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主
引用回帖:
Originally posted by 信彼南山 at 2011-07-01 16:38:14:
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 ...

高手啊！
又学习了学习！

人才是万能的，方法都是想出来的，呵呵。
赞一下 回复此楼
明德厚学、求是创新
11楼2011-07-01 17:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyzhixue

铁虫 (初入文坛)
引用回帖:
Originally posted by 信彼南山 at 2011-07-01 16:38:26:
运行结果：
a =

-4.9837156546668308559867719464884-.34823945996364925341340692246585*i
-4.9837156546668308559867719464884+.34823945996364925341340692246585*i
  4.983715654666830855986771946 ...

哇塞，，你们都好强啊。。我这个菜鸟真是太佩服了。。向你们学习！！谢谢你们啦。。
赞一下 回复此楼
12楼2011-07-03 21:59:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyzhixue

铁虫 (初入文坛)
引用回帖:
Originally posted by xiegangmai at 2011-06-30 11:15:59:
呵呵，应该可以了，换用'levenberg-marquardt'算法，增加迭代次数。

用fsolve求解，初值选取有很大关系，当初值选为[ 100; 100; 0 ]时，
[ X, Fval, exitflag ] = fsolve( @cyzhixue, [ 100; 100; 0 ], optio ...

你好，谢谢你回答我的问题，我还想问一下就是这个初值的选择有什么依据吗？？我又换了个方程，选了好多都不行啊。。数学没学好 呵呵。。
赞一下 回复此楼
13楼2011-07-03 22:07:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Luckner

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by xiegangmai at 2011-06-30 11:03:09:
尝试了一下用符号求解，找不到解析解。


[code]syms a b c;
solve('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- ...

请问这样的问题怎么解决呀？谢谢
函数是：
function f=myhs1(x)
f=[(2*x(1)+3*x(2)+x(3))^2*(x(1)+x(2)-x(3))/((0.4125-x(1)-x(2))*(0.21*0.825-0.5*x(1))^0.5*(1.65+1.5*x(1)+2*x(2))^1.5)-3.389*exp(11);(2*x(1)+3*x(2)+x(3))^3*(x(1)+x(2)-x(3))/((0.4125-x(1)-x(2))*(0.4125-x(2)-x(3))*(1.65+1.5*x(1)+2*x(2))^2)-1.155;(2*x(1)+3*x(2)+x(3))*x(3)/((x(1)+x(2)-x(3))*(0.4125-x(2)-x(3)))-8.469];


运行如下：
>> x0=[0.1    0.1    0.1];
>> [r,val]=fsolve(@myhs1,x0,opt)

结果为：
Solver stopped prematurely.
r =

    0.3465    0.0658    0.2101


val =

  1.0e+003 *

    0.4843
    1.8293
   -0.0001

我试了好多组初值都不行，要么出现上面的情况，要么是“No solution found.”
请问这是怎么回事呀？还是哪儿出问题了？
不胜感激！
赞一下 回复此楼
物换星移，沧海桑田，人心却没有多少改变。
14楼2011-12-19 22:01:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lylpiao

木虫 (正式写手)
送红花一朵
想问一下楼主，会不会解这个方程：
cos(v-t)*(-k*sin(h)+c*i*cos(h))-(i*sin(h)*(n*t+n*pi-b+k*v)+i*a*cos(h))*sin(v-t)=0;其中a,b,c,h,i都是已知的，只有v和t未知，求v=f(t)的形式，要怎么求？
赞一下 回复此楼
把一切交给时间！
15楼2015-12-14 11:08:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cyzhixue 的主题更新
152/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭