| 查看: 1540 | 回复: 18 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
geegee铜虫 (小有名气)
|
[求助]
matlab模拟求助
|
|||
|
各常量的取值如下： a2=1.0*10^(-8); a5=4.0*10^(4); a7=1.0*10^(-7); a8=1.0*10^(-2); a10=1.0*10^(0); a11=2.0*10^(-1); A=5.00; B=0.36; 三个方程组如下： 0=0.5*a2*A+0.5*a3*(B-C)*A-0.5*a11*x0*y0+a7*a5^2*A^2*x0^2; 0=-a10*(B-C)*y0+a8*a5*A^2*x0-a11*x0*y0; 0=(-1/2)*a11*y0+2*a7*a5^2*A^2*x0-a10*(B-C)-a11*x0; 变量为C,a3. 要求当C取值为[0,0.36）时，做出a3与C的关系图。 |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有166人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab求助？？？
已经有6人回复
- MATLAB GUI 求助
已经有3人回复
- Matlab编程求助！
已经有3人回复
- 求助matlab写循环
已经有6人回复
- matlab计算面积求助！！难
已经有8人回复
- 求助 matlab交点问题
已经有19人回复
- 求助matlab问题
已经有5人回复
- 继续求助Matlab绘图！
已经有8人回复
- 求助啊......matlab
已经有10人回复
- matlab 积分求助
已经有8人回复
- 求助matlab问题
已经有5人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 求助大侠帮忙，有关matlab模拟数据
已经有29人回复
- 【求助】复杂精馏过程的MATLAB模拟
已经有3人回复
geegee
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1024.8
- 散金: 42
- 红花: 2
- 帖子: 217
- 在线: 121.3小时
- 虫号: 1264903
- 注册: 2011-04-14
- 性别: MM
- 专业: 能源化工
|
麻烦帮我重新编辑一下。 备注：（1）x0，y0，a3都要大于0！！！！ （2）最后输出x0，y0，a3，C的取值。用TXT文本。 各常量的取值如下： a2=1.0*10^(-8); a5=4.0*10^(4); a7=1.0*10^(-7); a8=1.0*10^(-2); a10=1.0*10^(0); a11=2.0*10^(-1); A=5.00; B=0.36; 三个方程组如下： 0=0.5*a2*A+0.5*a3*(B-C)*A-0.5*a11*x0*y0+a7*a5^2*A^2*x0^2; 0=-a10*(B-C)*y0+a8*a5*A^2*x0-a11*x0*y0; 0=(-1/2)*a11*y0+2*a7*a5^2*A^2*x0-a10*(B-C)-a11*x0; 变量为C,a3. 要求当C取值为[0,0.36）时，做出a3与C的关系图。 |
9楼2013-07-05 19:13:29
geegee
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1024.8
- 散金: 42
- 红花: 2
- 帖子: 217
- 在线: 121.3小时
- 虫号: 1264903
- 注册: 2011-04-14
- 性别: MM
- 专业: 能源化工
2楼2013-07-05 11:35:24
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68128.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-05 17:25:39
geegee: 金币+20, ★★★★★最佳答案 2013-07-05 17:37:12
感谢参与，应助指数 +1
csgt0: 金币+2, 谢谢 2013-07-05 17:25:39
geegee: 金币+20, ★★★★★最佳答案 2013-07-05 17:37:12
|
以下程序复制进一个m文件中，运行即可。我是按照C从0到0.36取20个点做的。 function solution_nonl_C clear all;clc format long global C %-------------根据计算结果，修改迭代初值------------------------------------ x0=[1 1 1]; C=0; n=21; %C的取值个数 for i=1:n xx=fsolve(@nonlin,x0); a3(i)=xx(1); x0(i)=xx(2); y0(i)=xx(3); C=C+0.36/(n-1); end [a3' x0' y0'] aa3=a3'; aaa3=[aa3(1:14);aa3(end)]; cspan=0:0.36/(n-1):0.36; cspan1=[cspan(1:14) cspan(end)]; figure,plot(cspan,aa3,'ko',cspan1,aaa3,'b'),legend('全部点','满足x0和y0大于0的点','Location','best'),... xlabel('C'),ylabel('a3'); function f= nonlin(x) global C a2=1.0*10^(-8); a5=4.0*10^(4); a7=1.0*10^(-7); a8=1.0*10^(-2); a10=1.0*10^(0); a11=2.0*10^(-1); A=5.00; B=0.36; a3=x(1); x0=x(2); y0=x(3); f(1)=0.5*a2*A+0.5*a3*(B-C)*A-0.5*a11*x0*y0+a7*a5^2*A^2*x0^2; f(2)=-a10*(B-C)*y0+a8*a5*A^2*x0-a11*x0*y0; f(3)=(-1/2)*a11*y0+2*a7*a5^2*A^2*x0-a10*(B-C)-a11*x0;  附图23.jpg |
3楼2013-07-05 14:30:54
geegee
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1024.8
- 散金: 42
- 红花: 2
- 帖子: 217
- 在线: 121.3小时
- 虫号: 1264903
- 注册: 2011-04-14
- 性别: MM
- 专业: 能源化工
4楼2013-07-05 17:35:28
10