| 查看: 1383 | 回复: 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]) |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有115人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教下大家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 (博后)
- 金币: 12915.6
- 红花: 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
