24小时热门版块排行榜    

查看: 2974  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 9/450 2026-03-16 17:35 by 沐霖12138
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 本人考085602 化学工程 专硕 +12 不知道叫什么! 2026-03-15 14/700 2026-03-16 16:45 by 我的船我的海
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 306求调剂 +4 唐薏薏 2026-03-09 4/200 2026-03-14 01:19 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 321求调剂 +3 CUcat 2026-03-10 3/150 2026-03-14 00:25 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
请填处理意见