| 查看: 1807 | 回复: 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. 向各位大侠求助怎么改，谢谢！  |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有195人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关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
28