| 查看: 1553 | 回复: 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的关系图。 |
回复此楼» 猜你喜欢
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有198人回复
-
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有28人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有22人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
geegee
管理员
- 应助: 0 (幼儿园)
- 金币: 1024.8
- 散金: 42
- 红花: 2
- 帖子: 217
- 在线: 121.3小时
- 虫号: 1264903
- 注册: 2011-04-14
- 性别: MM
- 专业: 能源化工
15楼2013-07-05 20:30:15
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