| 查看: 789 | 回复: 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； |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有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
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.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
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.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