如何用matlab求解这个方程,试过了牛顿迭代和二分法,没成功,求高手指点一个程序! 5

72*cos(phic)-sin(x-phic)*(rp+a/sin(x))+0.0744708325271972-63.7356217792314=0其中,rp=(pi... 72*cos(phic)-sin(x-phic)*(rp+a/sin(x))+0.0744708325271972-63.7356217792314=0
其中,rp=(pi*9/2-(tan(20/180*pi)+tan(30/180*pi))*(9(1+0.25)))/(sec(20/180*pi)+sec(30/180*pi)-tan(20/180*pi)-tan(30/180*pi));
a=9*(1+0.25)-rp;
bc=pi*9/2;
phic=(2*(a*cot(x)+bc))/(9*16);
展开
 我来答
我行我素850915
2014-06-03 · TA获得超过8241个赞
知道大有可为答主
回答量:5788
采纳率:77%
帮助的人:1733万
展开全部
包含三角函数的方程是多值解,一般把所有的已知量全代进去,化为变量最少的方程,再求解。
上面的方程,经简化为:72*cos(0.1177*cot(x)+0.1963)-sin(x-(0.1177*cot(x)+0.1963))*(2.7766+8.4734/sin(x))-63.6612=0
则可求解:
>> syms x
f=72*cos(0.1177*cot(x)+0.1963)-sin(x-(0.1177*cot(x)+0.1963))*(2.7766+8.4734/sin(x))-63.6612;
x=solve(f)
x =
-225.6517942299019031056165710969
追问
补充:x的范围(0,pi/2),solve函数是matlab自带求解器,可否给出二分法或者牛顿迭代法的程序呢?
问题延伸了哈:72*cos(0.1177*cot(x)+0.1963)-sin(x-(0.1177*cot(x)+0.1963))*(2.7766+8.4734/sin(x))-C=0,C在(59.4351801148340,63.7356217792314)这个范围内变化,求指教!
追答
for c=linspace(59.4351801148340,63.735621779231,10);
syms x
f=72*cos(0.1177*cot(x)+0.1963)-sin(x-(0.1177*cot(x)+0.1963))*(2.7766+8.4734/sin(x))-c;
x=solve(f)
end
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式