无法使用fsolve之后的情节

2视图(30天)
jayash
jayash 2015年5月7日
回答: lei凯2015年5月8日
Ut = 0;
Ut < 20
cn = @ (q1、q2 n0) ((Ut * n0)。/ (1-cos(2π*。* q1) cos(2 *π。* q2)));
y = @ (q1、q2 n0) 0.25 . * ((1 + cn (q1、q2 n0)。^ (1/2) + (1 + cn (q1、q2 n0)。^ (1/2) 2);
一个= -0.5;
b = 0.5;
v = @ (n0) (quad2d (@ (q1、q2) y (q1、q2 n0), a, b, a, b));
简历= @ (n0) n0 + (0.5 v (n0)) * 1;
x0 = 0.1;
n = fzero(简历,x0);
情节(Ut, n,‘*’)
包含(“U / t”)
ylabel (“n0 / n”)
持有
Ut = Ut + 1;
结束
情节应该类似看起来像y = 1 / x但我只得到一点。有人请帮帮我。

答案(2)

沃尔特·罗伯森
沃尔特·罗伯森 2015年5月8日
你忘了告诉我们,你也得到一些输出,即
警告:到达最大数量的函数评估(2000)。测试结果失败全球错误。
>在quad2d在244
@ (n0) (quad2d (@ (q1、q2) y (q1、q2 n0), a, b, a, b))
@ (n0) n0 + (0.5 v (n0)) * 1
fzero在286
错误使用fzero(第309行)
函数开始猜价值必须是有限的和真实的。
当Ut = 1(第二次迭代)。集成是失败。你确定积分收敛吗?

lei凯
lei凯 2015年5月8日
有问题你的函数的简历,当Ut = 1, n = fzero(简历,x0);将错误

类别

标签

s manbetx 845

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!