| 查看: 1061 | 回复: 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]) |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有221人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教下大家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
10