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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab建模中solver的问题
51/1返回列表
查看: 1626  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

驴肉火烧

金虫 (小有名气)

[求助] 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.

向各位大侠求助怎么改，谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-10-05 21:17:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

驴肉火烧

金虫 (小有名气)
引用回帖:
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的回帖
查看全部 4 个回答

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

%建立方程
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): 十分感谢！！！ 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的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 329求调剂 +4 钮恩雪 2026-03-25 4/200 2026-03-26 20:53 by sanrepian
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +6 晨昏线与星海 2026-03-20 6/300 2026-03-26 19:20 by macy2011
[考研] 279求调剂 +6 红衣隐官 2026-03-21 6/300 2026-03-26 18:32 by 不吃魚的貓
[考研] 考研调剂 +8 小蜡新笔 2026-03-26 8/400 2026-03-26 16:18 by dick_runner
[考研] 275求调剂 +9 Micky11223 2026-03-25 11/550 2026-03-26 15:54 by 不吃魚的貓
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +6 哈哈157349 2026-03-21 6/300 2026-03-26 10:39 by 醉在风里
[考研] 299求调剂 +4 15188958825 2026-03-25 4/200 2026-03-25 22:56 by 418490947
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 材料求调剂 +4 .m.. 2026-03-25 4/200 2026-03-25 21:30 by peike
[考研] 材料学硕，求调剂 6+4 糖葫芦888ll 2026-03-22 9/450 2026-03-25 11:19 by greychen00
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 材料/农业专业，07/08开头均可，过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪！？！
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
信息提示
关闭
请填处理意见
关闭