| 查看: 3333 | 回复: 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 |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有109人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
明天放榜？
已经有4人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab-常微分方程参数估计
已经有12人回复
- 【分享】MATLAB在智能网络管理信息处理中的应用
已经有16人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- matlab中如何将双坐标系图中的另外一个坐标系设置为当前坐标系呢？
已经有10人回复
- 【教程】MATLAB的S-Function编写指导
已经有103人回复
- 用matlab统计图片中某种颜色成分占总图片中所有颜色成分的比率
已经有4人回复
- 如何在matlab中绘制三维立体图
已经有13人回复
- 【求助】用MATLAB画三维图程序编写的问题
已经有6人回复
- 【求助】matlab 中，几个m文件调用
已经有14人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【原创】分享一个matlab编写的模糊KNN的例子，经过验证效果不错。
已经有12人回复
- matlab的fig中插入特殊字符或中文或公式的办法？要求eps格式
已经有5人回复
newtonnapo
木虫 (正式写手)
毛毛虫
- 应助: 4 (幼儿园)
- 金币: 2058
- 帖子: 334
- 在线: 93.7小时
- 虫号: 967058
- 注册: 2010-03-10
- 性别: GG
- 专业: 能源化工
4楼2011-06-09 08:49:02
newtonnapo
木虫 (正式写手)
毛毛虫
- 应助: 4 (幼儿园)
- 金币: 2058
- 帖子: 334
- 在线: 93.7小时
- 虫号: 967058
- 注册: 2010-03-10
- 性别: GG
- 专业: 能源化工
2楼2011-06-08 15:58:21
lds588
铁杆木虫 (著名写手)
愚钝
- 应助: 15 (小学生)
- 金币: 10949.6
- 散金: 2018
- 红花: 5
- 帖子: 1444
- 在线: 206.3小时
- 虫号: 890720
- 注册: 2009-11-01
- 专业: 凝聚态物性I:结构、力学和
3楼2011-06-08 21:13:59
newtonnapo
木虫 (正式写手)
毛毛虫
- 应助: 4 (幼儿园)
- 金币: 2058
- 帖子: 334
- 在线: 93.7小时
- 虫号: 967058
- 注册: 2010-03-10
- 性别: GG
- 专业: 能源化工
|
整个程序是这样的，有什么不对的吗 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
