作业帮 > 数学 > 作业

matlab 解方程 遇到的困难求帮助

来源:学生作业帮 编辑:拍题作业网作业帮 分类:数学作业 时间:2024/04/30 13:25:53
matlab 解方程 遇到的困难求帮助
已知 g=9.8;U=5;d=100;theta=0:2*pi; 求k0
g=9.8;U=5;d=100;theta=0:2*pi;
fk=@(k) g*tanh(k)*d/U^2./(cos(theta)).^2;
fsolve(fk,pi*ones(size(theta))
解全为零
再问: 不可能解全为零啊!后面还需要用到这个数据呀
再答: 我用contour画出的图,看起来就是全为零的。。。你查查你的公式吧
再问: 你那里是tanh(k*d)才对,但是我运行你的程序怎么出错了呀? fsolve(fk,pi*ones(size(theta)) | Error: Expression or statement is incorrect--possibly unbalanced (, {, or [.
再答: 就算是tanh(k*d)也不可能啊 1/cos(x)^2 这个在分母,不可能为零, g/U^2 这个是常数且不为零 tanh(x) 单调 所以你的方程的解只有在tanh(x)=0时有解,而tanh(x)=0的解就是x=0 fsolve(fk,pi*ones(size(theta)))%少贴了半个括号
再问: 没有理解 这方程左边不是有一个K吗,是求解非齐次方程哇,很感谢你详细的解答。
再答: g=9.8;U=5;d=100;theta=0:2*pi; fk=@(k) g*tanh(k*d)/U^2./(cos(theta)).^2-k; fsolve(fk,ones(size(theta)) 这样就是你要的解了,但是你的份子上有1/0的问题,在theta=pi/2与1.5pi附近的解会很不准