| 查看: 1501 | 回复: 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人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
物理学I论文润色/翻译怎么收费?
已经有262人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有1人回复
-
投稿chemical physical letters不送审？
已经有2人回复
-
南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生
已经有21人回复
-
杂志要求提供official language-editing certificate
已经有4人回复
-
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有15人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
geegee: 金币+15, ★★★★★最佳答案 2013-07-05 22:29:07
geegee: 金币+15, ★★★★★最佳答案 2013-07-05 22:29:07
|
刚没传上附件，现在好了 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : solution_nonl_C1.m
2013-07-05 20:12:39, 843 bytes
14楼2013-07-05 20:13:00
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
- 金币: 68123.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