24小时热门版块排行榜    

查看: 3374  |  回复: 12

newtonnapo

木虫 (正式写手)

毛毛虫

[求助] matlab中dydt的编写要求

我编的一段程序的一部分,yeild是一个行向量,他老是说??? Error: File: F:\ivpodefun.m Line: 20 Column: 14Missing MATLAB operator.

Line: 20 Column: 14就是下面加粗的地方
for i=3:1:39
    for j=1:1:38
    if i>j+1
        sum=sum+1;
    dyeild(i)dt=dyeild(i)dt-k(sum)*yeild(i);
    dyeild(j)dt=dyeild(j)dt+k(sum)*yeild(i);
    dyeild(39+i-j)dt=dyeild(39+i-j)dt+k(sum)*yeild(i);
    else
    end
    end
end
回复此楼

» 猜你喜欢

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

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

newtonnapo

木虫 (正式写手)

毛毛虫

我是个新手,看了半天了也没弄明白,希望大虾不吝指导啊
攒人品。
2楼2011-06-08 15:58:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝

【答案】应助回帖


臭水沟(金币+1): 谢谢参与~~ 2011-06-08 22:25:41
引用回帖:
Originally posted by newtonnapo at 2011-06-08 15:57:11:
我编的一段程序的一部分,yeild是一个行向量,他老是说??? Error: File: F:\ivpodefun.m Line: 20 Column: 14Missing MATLAB operator.

Line: 20 Column: 14就是下面加粗的地方
for i=3:1:39
    f ...

dyeild与yeid是两个不同的数组,请问你的dyeild是否定义?
[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
3楼2011-06-08 21:13:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

newtonnapo

木虫 (正式写手)

毛毛虫

臭水沟: 回复别人帖子时,点击“引用回复”,这样对方才能看到你的回复! 2011-06-09 09:14:27
啊?dyeild是对yeid的求导,yeild是一个很长的东西
攒人品。
4楼2011-06-09 08:49:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

newtonnapo

木虫 (正式写手)

毛毛虫

引用回帖:
Originally posted by lds588 at 2011-06-08 21:13:59:
dyeild与yeid是两个不同的数组,请问你的dyeild是否定义?

整个程序是这样的,有什么不对的吗

function dyeilddt=ivpodefun(t,yeild)
format long
sum=0;
feed=cell(1,2);
feed(1,1)={[1 1 0 0;2 1 0 0;3 1 0 0;4 1 0 0;5 1 0 0;6 1 0 0;7 1 0 0;8 1 0 0;9 1 0 0;10 1 0 0;
    11 1 0 0;12 1 0 0;13 1 0 0;14 1 0 0;15 1 0 0;16 1 0 0;17 1 0 0;18 1 0 0;19 1 0 0;20 1 0 0;
    21 1 0 0;22 1 0 0;23 1 0 0;24 1 0 0;25 1 0 0;26 1 0 0;27 1 0 0;28 1 0 0;29 1 0 0;30 1 0 0;
    31 1 0 0;32 1 0 0;33 1 0 0;34 1 0 0;35 1 0 0;36 1 0 0;37 1 0 0;38 1 0 0;39 1 0 0]};
feed(1,2)={[2 0 0 0;3 0 0 0;4 0 0 0;5 0 0 0;6 0 0 0;7 0 0 0;8 0 0 0;9 0 0 0;10 0 0 0;
    11 0 0 0;12 0 0 0;13 0 0 0;14 0 0 0;15 0 0 0;16 0 0 0;17 0 0 0;18 0 0 0;19 0 0 0;20 0 0 0;
    21 0 0 0;22 0 0 0;23 0 0 0;24 0 0 0;25 0 0 0;26 0 0 0;27 0 0 0;28 0 0 0;29 0 0 0;30 0 0 0;
    31 0 0 0;32 0 0 0;33 0 0 0;34 0 0 0;35 0 0 0;36 0 0 0;37 0 0 0;38 0 0 0;39 0 0 0]};
yeild=xlsread('F:\crack\1','Sheet1','F6:F82');
k=xlsread('F:\crack\1','Sheet1','A2:A667');
dyeilddt=zeros(77,1);
for i=3:1:39
    for j=1:1:38
    if i>j+1
        sum=sum+1;
    dyeild(i)dt=dyeild(i)dt-k(sum)*yeild(i);
    dyeild(j)dt=dyeild(j)dt+k(sum)*yeild(i);
    dyeild(39+i-j)dt=dyeild(39+i-j)dt+k(sum)*yeild(i);
    else
    end
    end
end
攒人品。
5楼2011-06-09 13:40:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝

引用回帖:
Originally posted by newtonnapo at 2011-06-09 13:40:16:
整个程序是这样的,有什么不对的吗

function dyeilddt=ivpodefun(t,yeild)
format long
sum=0;
feed=cell(1,2);
feed(1,1)={[1 1 0 0;2 1 0 0;3 1 0 0;4 1 0 0;5 1 0 0;6 1 0 0;7 1 0 0;8 1 0 0;9 1 0  ...

dydt?请问您是想对时间求导数?
[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
6楼2011-06-09 13:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

newtonnapo

木虫 (正式写手)

毛毛虫

引用回帖:
Originally posted by lds588 at 2011-06-09 13:57:06:
dydt?请问您是想对时间求导数?

en,是的
攒人品。
7楼2011-06-09 18:25:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

【答案】应助回帖

★ ★
xiegangmai(金币+2): 谢谢应助! 2011-06-09 23:37:29
dyeild(i)dt---------------dyeilddt(i)
dyeild(j)dt---------------dyeilddt(j)
dyeild(39+i-j)dt---------------dyeilddt(39+i-j)
8楼2011-06-09 22:31:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lds588

铁杆木虫 (著名写手)

愚钝

★ ★
xiegangmai(金币+2): 谢谢参与!这些确实很基础。 2011-06-09 23:37:54
引用回帖:
Originally posted by newtonnapo at 2011-06-09 18:25:05:
en,是的

matlab 没有你想的那么智能,求导要自己写,只写dydt它是不会工作的,建议用diff函数求导,dyeild(i)dt在matlab中只能看做一个变量名,而且是错的变量名,语法有些不太合适,建议你能看点基础的书~~~
可能我这样说有些冒昧,请你原谅
[url=http://weibo.com/u/2577306830?s=6uyXnP][img]http://service.t.sina.com.cn/widget/qmd/2577306830/fa2fdb00/1.png[/img][/url]
9楼2011-06-09 23:11:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

newtonnapo

木虫 (正式写手)

毛毛虫

引用回帖:
Originally posted by lds588 at 2011-06-09 23:11:04:
matlab 没有你想的那么智能,求导要自己写,只写dydt它是不会工作的,建议用diff函数求导,dyeild(i)dt在matlab中只能看做一个变量名,而且是错的变量名,语法有些不太合适,建议你能看点基础的书~~~
可能我这 ...

我就是在看书的啊,边看边写的。。。。
攒人品。
10楼2011-06-10 09:11:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 newtonnapo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 9/450 2026-09-06 17:48 by myrtle
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见