| 查看: 962 | 回复: 5 | ||
chem.Jlt新虫 (著名写手)
|
[求助]
我现在有如下①②两个方程，③是限制条件，求各位大佬给看一下怎么写这个matlab语言 已有1人参与
|
|
①(3a+b)^3*（a-b)*1.2^2=0.267*(1-a)*(6-a-b)(7+2a)^2； ②(3a+b)*b=1.368(a-b)(6-a-b)； ③0<a<1;a>b； |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有14人回复
-
这样的filecode谁见过
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
fileCode有新解读？
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
chem.Jlt
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 8985.2
- 散金: 898
- 红花: 2
- 帖子: 2403
- 在线: 80.4小时
- 虫号: 8714408
- 注册: 2018-05-10
- 性别: GG
- 专业: 能源化工
2楼2018-11-30 08:10:37
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
前两个方程求出的结果是唯一的，a=0.9394,b= 0.6019 .你的第三个式子的意义何在？ function p=example_fsolve p0=[1;1]; %%% 给定初始值 options=optimset('MaxFunEvals',20000,'MaxIter',2000); [p,fval,exitflag]=fsolve(@f,p0,options); function F=f(p) %%% 待求方程 F=[(3*p(1)+p(2))^2*(p(1)-p(2))*1.2^2-0.267*(1-p(1))*(6-p(1)-p(2))*(7+2*p(1))^2; (3*p(1)+p(2))*p(2)-1.368*(p(1)-p(2))*(6-p(1)-p(2))]; end end %% ****** P0=[]; P=[]; exitflag=0; while p(1)<=0||p(1)>=1||p(1)<=p(2) [p,p0,fval,exitflag]=example_fsolve; if exitflag==1; P0=[P0;p0]; P=[P;p]; end end |
3楼2018-11-30 09:45:02
chem.Jlt
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 8985.2
- 散金: 898
- 红花: 2
- 帖子: 2403
- 在线: 80.4小时
- 虫号: 8714408
- 注册: 2018-05-10
- 性别: GG
- 专业: 能源化工
4楼2018-11-30 20:51:18
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
chem.Jlt5楼2018-11-30 21:00:08
chem.Jlt
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 8985.2
- 散金: 898
- 红花: 2
- 帖子: 2403
- 在线: 80.4小时
- 虫号: 8714408
- 注册: 2018-05-10
- 性别: GG
- 专业: 能源化工
6楼2018-11-30 23:28:59
20