24小时热门版块排行榜    

查看: 1620  |  回复: 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.

向各位大侠求助怎么改,谢谢!
回复此楼

» 猜你喜欢

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

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

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的回帖

驴肉火烧

金虫 (小有名气)

引用回帖:
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的回帖
相关版块跳转 我要订阅楼主 驴肉火烧 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 277材料科学与工程080500求调剂 +7 自由煎饼果子 2026-03-16 7/350 2026-03-22 22:40 by ACS Nano——
[考研] 263求调剂 +4 yqdszhdap- 2026-03-22 4/200 2026-03-22 21:20 by 1144970272
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 0856材料专硕353求调剂 +4 NIFFFfff 2026-03-20 4/200 2026-03-22 09:49 by 2026paper
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
信息提示
请填处理意见