| 查看: 922 | 回复: 4 | ||
hailiang金虫 (正式写手)
|
[求助]
请教MATLAB中下面的提示是什么意思？
|
|
请教MATLAB中下面的提示是什么意思？ In an assignment A(I) = B, the number of elements in B and I must be the same. Error in hlduffing (line 5) dy(2)=f*cos(0.8*t)-0.2*y(2)-y(1)-y(1)^3; Error in odearguments (line 88) f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in ode45 (line 114) [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ... Error in hltest2013 (line 7) [t,y]=ode45('hlduffing',t0,[1;1]) |
回复此楼» 猜你喜欢
🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队
已经有0人回复
-
第一性原理计算方向2026级博士申请 PRB*1，四级484
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有85人回复
-
求助VISSIM破解版软件
已经有0人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
请教下大家6.66kpa 45度是什么意思啊！
已经有5人回复- 请教高斯公式，曲面积分和体积积分的转换。
已经有8人回复
- 请教个matlab中fft(X,n)，关于n取值的问题
已经有7人回复
- ＭＡＴＬＡＢ计算编程请教
已经有15人回复
- 请问一下matlab这个是什么问题
已经有10人回复
- matlab二重积分的错误提示，指的是什么
已经有53人回复
- 各位，请问下面：对于非线性项的处理是什么意思？
已经有5人回复
- 技术入门费 是什么意思啊？请教一下各位
已经有30人回复
- 以下提示是什么意思？请指教！
已经有3人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 请教大家一个反卷积的问题~~帮忙推导或者matlab编程计算~~
已经有5人回复
- 初学者，请教大家关于matlab解多元二次方程组
已经有4人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- 请教matlab反应动力学参数估计遇到的问题，谢谢
已经有15人回复
- 继续求助Matlab绘图！
已经有8人回复
- 大家帮忙看一下，关于模拟文献的
已经有6人回复
- VS2005下matlab调用
已经有8人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【分享】大家帮忙看一下这个Ni的峰归属于NiO还是Ni2O3
已经有7人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 计算强帖: 3
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2013-03-08 00:25:54
hailiang
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 3551.2
- 散金: 838
- 红花: 2
- 帖子: 866
- 在线: 792.5小时
- 虫号: 56375
- 注册: 2005-01-08
- 性别: GG
- 专业: 动力学与控制
|
麻烦帮我看看下面求解微分方程的程序问题出在哪？ 不对f循环的时候没有问题的，做循环就出问题了 function dy = hlduffing(t,y) global f dy=zeros(3,1); dy(1)=y(2); dy(2)=f*cos(0.8*t)-0.2*y(2)-y(1)-y(1)^3; end function hlduffing1_bifur_f_getmax % 最大值法求解分岔图 global f clear all t0=[0 100];%积分时间 %bifurcation for f=linspace(10,32,1000); % f=10:32:1000; %r的变化精度 [t,y]=ode45('hlduffing',t0,[1;1]); [Xmax]=getmax(y(:,1)); plot(f,Xmax,'k','markersize',1) hold on clear Xmax end |
3楼2013-03-08 10:29:33
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
hailiang: 金币+10, ★★★★★最佳答案 2013-04-07 20:37:40
csgt0: 金币+2, 谢谢 2013-04-09 10:44:23
hailiang: 金币+10, ★★★★★最佳答案 2013-04-07 20:37:40
csgt0: 金币+2, 谢谢 2013-04-09 10:44:23
|
这里有两个错误 ： 1. clear all 的位置 它应该在程序的最上方，而在你的程序里，它在 global f 之后，这样运行到 clear all 之后，全局变量 f 会被删除，你的子程序就无法获得 f 的值了 2. dy=zeros(3,1); 子程序里面，dy 的预设的错误，你的程序假设它是 （2, 1） : dy(1) / dy(2), 但是在程序一开始，却执行 dy=zeros(3,1); 这样就会出现 长度错误，正确的应该是 dy=zeros(2,1); 不过因为我没有你的 getmax 程序，无法运行，但是至少能解决你开始所提出的问题了 |
4楼2013-04-07 18:58:08
5楼2015-02-04 16:58:24