使用quad的错误(第75行)integrand函数必须返回与输入向量长度相同的输出向量。
8次浏览(最近30天)
显示旧的注释
我已经根据答案编辑了代码。
谁能解释一下错误背后的原因,并帮助我纠正以下代码:
U = 1;
E = @(q) 2*(1-cos(2* *q));
hq = @(q,n0) ((e(q)).^2+2*U*n0*(e(q))).^0.5;
y = @(q,n0) ((e(q))+(U*n0))/hq(q,n0))-1;
A = -0.5;
B = 0.5;
V = @(n0) quad(@(q) y(q,n0),a,b);
CV =@(n0) n0+(0.5*v(n0))-1;
而U < 20
N0 = 0.1;
Options = optimset(“显示”,“通路”);
N = fsolve(cv(n0),n0,0.1,options);
情节(U, n)
持有在
U = U + 1;
结束
错误:
错误使用quad(第75行)
的被积函数必须返回与输入向量长度相同的输出向量。
错误@ (n0)四(@ (q) y (q, n0), a, b)
错误@ (n0) n0 + (0.5 * v (n0)) 1
错误在simulv1中(第12行)
N = fsolve(cv(n0),n0,0.1,options);