| 查看: 1743 | 回复: 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的关系图。 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有208人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
16楼2013-07-05 20:35:04
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
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.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
30