| 查看: 1630 | 回复: 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. 向各位大侠求助怎么改，谢谢！  |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有67人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关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
- 专业: 化学反应工程
2楼2011-10-06 08:22:25
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
|
%运行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
