版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

驴肉火烧

金虫 (小有名气)

应助: 16 (小学生)
金币: 1345.5
帖子: 224
在线: 109.2小时
虫号: 1012862
注册: 2010-05-06
专业: 传动机械学

[求助] 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人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有186人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2011-10-05 21:17:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

计算强帖: 1
应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

【答案】应助回帖

%建立方程
function dy=fai(t,y,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程序求解
[x,y]=ode45('fai',[0 15],0,[],0,0.008)
plot(x,y)

赞一下

回复此楼

2楼2011-10-06 08:22:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

驴肉火烧

金虫 (小有名气)

应助: 16 (小学生)
金币: 1345.5
帖子: 224
在线: 109.2小时
虫号: 1012862
注册: 2010-05-06
专业: 传动机械学

引用回帖:

2楼: Originally posted by change0618 at 2011-10-06 08:22:25:
%建立方程
function dy=fai(t,y,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*c ...

%运行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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主驴肉火烧的主题更新

返回列表