| 查看: 966 | 回复: 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； |
回复此楼» 猜你喜欢
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有8人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
关于Filecode分析方法
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
独孤神宇
版主 (知名作家)
- 应助: 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
- 专业: 能源化工
2楼2018-11-30 08:10:37
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
