24小时热门版块排行榜    

查看: 2993  |  回复: 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-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[考研] 070300化学319求调剂 +4 锦鲤0909 2026-03-17 4/200 2026-03-17 18:21 by 重科小霸王
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 化学工程321分求调剂 +11 大米饭! 2026-03-15 14/700 2026-03-17 17:11 by ruiyingmiao
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 08工科 320总分 求调剂 +4 梨花珞晚风 2026-03-17 4/200 2026-03-17 13:38 by houyaoxu
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 一志愿中科院,化学方向,295求调剂 +4 一氧二氮 2026-03-11 4/200 2026-03-13 22:35 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
信息提示
请填处理意见