| 查看: 1828 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
matlab建模中solver的问题
|
|||
|
%建立方程 function dy=fai(x,y) dy=2*(1/2*1.5*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*0.7*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*cos(y)))/(100+500) %运行matlab程序求解 [X,Y]=ode45('fai',[0 0.008],[0 1.5],y(0)=0) plot(x,y) 运行结果 ??? [X,Y]=ode45('fai',[0 0.008],[0 1.5],y(0)=0) | Error: The expression to the left of the equals sign is not a valid target for an assignment. 向各位大侠求助怎么改，谢谢！  |
回复此楼» 猜你喜欢
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有178人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有111人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有68人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关MS 的建模 优化及计算的一些问题
已经有6人回复
- 数学建模matlab算法大全
已经有720人回复
- Matlab中RBF建模，仿真是遇到的问题
已经有6人回复
- 【求助成功】关于ZnO理想晶胞建模的问题
已经有20人回复
- 【求助】建模问题
已经有12人回复
- 【讨论】数学建模的培训问题
已经有15人回复
- 【求助】ABINIT中的建模问题
已经有26人回复
- 【求助】MS建模导入lammps问题求救
已经有14人回复
- 【讨论】求助：数学建模问题
已经有8人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
【答案】应助回帖
驴肉火烧(金币+1): 十分感谢！！！ 2011-10-06 09:48:36
|
%建立方程 function dy=fai(t,y) global k1 k2 dy=2*(1/2*1.5*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*0.7*(0.00003*t^4-0.0002*t^3+0.008*t^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*cos(y)))/(100+500); %运行matlab程序求解 >> global k1 k2 >> k1=0; >> k2=0.008; >> [x,y]=ode45('fai',[0 15],0); >> plot(x,y) |
3楼2011-10-06 09:34:19
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2011-10-06 08:22:25
|
%运行matlab程序求解 [x,y]=ode45(@fai,[0 15],0,[],0,0.008); plot(x,y) %建立方程 function dy=fai(x,y,k1,k2) dy=2*(1/2*1.5*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*0.7*(0.00003*x^4-0.0002*x^3+0.008*x^2)^2+1/2*k2*0.008^2+(1.5+0.7)*9.8*0.008-100*0.008-1/2*k1*(0.1-0.1*cos(y)))/(100+500); 以上是change0618这个朋友通过站内消息帮忙改正后的正确答案，谢谢！ |
4楼2011-10-06 09:41:55
