和happy教授结果略有出入啊, 不知道有没有问题以及更简便的方法:@)
k=linspace(0,20);
f = @(w,k)0.01*w/100.*sum(1./(w-0.5*sin(k/2).^2));
fplot(@(w)f(w,k),)
[ 本帖最后由 ChaChing 于 2010-4-28 00:16 编辑 ] 结果的不同主要是计算的误差造成的
你可以改变步长看看，呵呵，差别可能更大
[ 本帖最后由 ChaChing 于 2010-4-28 00:04 编辑 ] 多谢各位的讨论，这个问题是我刚开始用Matlab时（没有当版主之前）的一个问题，问题本是二维(kx,ky)的，我稍微简化了一下拿来讨论，现在看来，是我不小心才没有画出奇异性。------能引起大家的讨论也算是没有白问了，相当感谢各位的热心。
[ 本帖最后由 xjzuo 于 2008-6-4 16:10 编辑 ]
页:
1
[2]