fsolve解方程的结果和初值一样，好像没带进去一样… - 程序语言 - MATLAB/Mathematica

版块导航: 正在加载中...

客户端APP下载

复试调剂辅导

登录注册

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

24小时热门版块排行榜石溪大学接受考研调剂申请>

>论坛更新日志 (7108)
>虫友互识 (564)
>文献求助 (424)
>硕博家园 (411)
>导师招生 (373)
>考博 (279)
>招聘信息布告栏 (261)
>考研 (217)
>休闲灌水 (216)
>论文投稿 (160)
>博后之家 (95)
>教师之家 (80)
>基金申请 (66)
>海外博后 (62)
>有机交流 (60)
>电化学 (54)

【调剂】北京石油化工学院2024年16个专业接受调剂

返回列表

jia544

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 309.1
帖子: 110
在线: 3.3小时
虫号: 26939297

已领完

fsolve解方程的结果和初值一样，好像没带进去一样…: 领取红包 (小木虫手机app专属红包)

扫一扫，下载小木虫客户端

function f=myfun3(b)
x(1)=0.386126923;x(2)=0.638873077;x(3)=8.17E-09;x(4)=0.00792481;x(5)=0.243680566;x(6)=0.067260778;x(7)=0.118;
f(1)=(3.55e-10)-x(2).*x(3).*exp(2.*b(2).*x(5)+2.*b(1).*x(1)).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5))./(x(1).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5)+2.*b(3).*x(5)+2.*b(1).*x(2)));
f(2)=(4.60e-07)-x(3).*x(5).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5)).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5))./(0.004941889.*1);
f(3)=(4.66e-11)-x(3).*x(6).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5)).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5))./(x(5).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5)));
f(4)=(1.0577e-14)-x(3).*x(4).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5)).*exp(-0.5547.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5./(1+1.2.*((x(1)+x(5)+4.*x(6)+x(3)+x(4))./2).^0.5));
f(5)=x(2)+x(1)-1.025;
f(6)=x(7).*1.025-x(5)-x(6)-0.004941889;
f(7)=x(1)+x(3)-x(4)-x(5)-2.*x(6);
end

OPT=optimset;
OPT.LargeScale='off';
fun=@myfun3;
b0=[0.1893,0.141,0.2037];
[b,Y,c,d]=fsolve(fun,b0,OPT);

无论b输入什么样的初值，解出来的结果都和初值一样是咋回事呢，郁闷，求大家帮我看看是哪里的问题，应该咋解决呢

发自小木虫Android客户端

回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2022-06-18 10:26:52

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

cupb0

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1282.2
帖子: 267
在线: 24小时
虫号: 24932156

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

先手工化简一下，里面常数的量级差别太大，f(5)看起来是不满足的

发自小木虫Android客户端

赞一下

回复此楼

31楼2022-06-20 16:05:53

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

hzlhm

至尊木虫 (著名写手)

应助: 372 (硕士)
金币: 16128.8
帖子: 2560
在线: 571.9小时
虫号: 1517335

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

题主不知怎么考虑的，只要三个变量，却列了七个方程，而其中六个方程是废弃方程，没有联系。

赞一下

回复此楼

32楼2022-06-20 16:46:37

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

hzlhm

至尊木虫 (著名写手)

应助: 372 (硕士)
金币: 16128.8
帖子: 2560
在线: 571.9小时
虫号: 1517335

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

把已有的x(1).x(2)....去掉，看成未注变量。这样的话，得到如下结果

123.jpg

赞一下

回复此楼

33楼2022-06-20 16:59:20

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

jia544

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 309.1
帖子: 110
在线: 3.3小时
虫号: 26939297

引用回帖:

32楼: Originally posted by hzlhm at 2022-06-20 16:46:37
题主不知怎么考虑的，只要三个变量，却列了七个方程，而其中六个方程是废弃方程，没有联系。

对，是我的疏忽，把方程写错了，已经解决啦，谢谢～

发自小木虫Android客户端

赞一下

回复此楼

35楼2022-06-22 21:09:47

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

简单回复

tzynew2楼

2022-06-18 10:30 回复

jia544(金币+1): 谢谢参与

[论文投稿] 无under review，直接DIP，咋么办。 +3	lizhengke06 2024-04-19	4/200	2024-04-19 19:24 by lizhengke06
[考研] 找一个允许我去清华联合培养的研究生导师 +7	lekinna 2024-04-19	14/700	2024-04-19 17:58 by ld_cauchy
[考研] 考研英语带学 +6	choice666 2024-04-18	6/300	2024-04-19 17:46 by charles-c
[考博] 山东大学化学与化工学院刘磊课题组博士研究生招生启事 +3	sdorganic 2024-04-17	5/250	2024-04-19 17:04 by 368ghnf
[基金申请] 申请省自然科学基金，研究区能否是省外区域 100+3	喜欢兔兔的我 2024-04-15	12/600	2024-04-19 16:40 by 淀粉搬运工
[基金申请] 基金和生小孩 +33	Ausy 2024-04-15	35/1750	2024-04-19 16:21 by feng6531
[教师之家] 美国教授开除中国研究生：我就不该录取你！ +7	苏东坡二世 2024-04-15	7/350	2024-04-19 16:20 by feng6531
[考博] 25申博记录贴 +3	我属驴核动力驴 2024-04-18	4/200	2024-04-19 10:53 by 安塔瓦拉多
[基金申请] 国自然青年基金只有一篇能第一标注，其他都是第二标注能结题吗？ 10+4	lice_1987 2024-04-18	6/300	2024-04-19 08:29 by kudofaye
[有机交流] 兄弟们帮我看看这两个结构怎么合成 +5	xl2088131 2024-04-17	5/250	2024-04-19 08:09 by 西瓜霜华
[考研] 274求调剂 +7	四个大字，你滴� 2024-04-14	10/500	2024-04-17 22:38 by 华北刘兵
[论文投稿] 可以打电话问编辑部是否可以先发录用通知吗 +7	双倍好运锦鲤 2024-04-14	10/500	2024-04-17 13:38 by cjzhu
[考研] 求调剂 +4	桃岸雪 2024-04-15	5/250	2024-04-15 18:49 by mthwyj
[考研] 347求调剂 +4	安谷先生 2024-04-15	4/200	2024-04-15 18:27 by huao连
[考研] 材料与化工267求调剂 +5	工科材料267 2024-04-15	8/400	2024-04-15 14:11 by wyx颖颖
[材料综合] 请问哪里可以测试低压吸氢的PCT曲线 +3	yunshengcd 2024-04-14	5/250	2024-04-15 13:59 by zqdsb
[考研] 专硕调剂招生 +3	电致发光 2024-04-15	4/200	2024-04-15 07:34 by ashorewmj
[考研] 290求调剂 +3	杨yhr 2024-04-14	5/250	2024-04-14 21:50 by coco1981
[考研] 284生物与医药求调剂（考数学没过A区国家线） +5	胖男孩 2024-04-12	5/250	2024-04-13 18:38 by yezixiaop
[考研] 290，一志愿报考重庆大学生物与医药专硕，求调剂 +4	GPX4 2024-04-13	5/250	2024-04-13 17:19 by lincunhui

24小时热门版块排行榜 石溪大学接受考研调剂申请>

jia544

» 猜你喜欢

» 抢金币啦！回帖就可以得到:查看全部散金贴

cupb0

hzlhm

hzlhm

jia544

24小时热门版块排行榜石溪大学接受考研调剂申请>

» 抢金币啦！回帖就可以得到:
查看全部散金贴