10弧度是一个不寻常的角度,但这可能是对你的问题不是那么重要。
找到一个变量与其他变量的函数方程系统使用象征性的工具箱
12个视图(30天)
显示旧的评论
我有3个方程4变量之间的连接,假设:
信谊a b c d;
r = b + d;
e =π/ 2;
eq1 = c - * sin (d) + b * cos (d) = = 0;
eq2 = r - c * tan (e) - * cos (d) - b * sin (d) = = 0;
eq3 = b = = e * r + * tan (e);
现在我想找到c (d),意味着,我想找到什么是c作为另一个变量的函数,在这种情况下- d。
然后执行其他操作c (d) differntiate等。
什么是最好的/最简单的方法使用符号工具箱?
谢谢。
接受的答案
沃尔特·罗伯森
2020年4月4日
信谊a b c d;
r = b + d;
e =符号(π)/ 10;%不π/ 2,没有办法!
eq1 = c - * sin (d) + b * cos (d) = = 0;
eq2 = r - c * tan (e) - * cos (d) - b * sin (d) = = 0;
eq3 = b = = e * r + * tan (e);
索尔=解决([eq1、eq2 eq3], [a, b, c],“returnconditions”,真正的);
c = symfun(简化(sol.c), d);
这将是有效的
除了
在四组d值的值可以通过分析推导出sol.conditions(在每种情况下添加一个* 2π的整数倍数被另一个禁止值)
2的评论
沃尔特·罗伯森
2020年4月4日
双(rad2deg (theta_max (3))
ans =
53.0352986075461 - 6.4856523980211 e - 71 - i
复杂的部分是由于舍入误差。
不过,我建议你重新编写代码来减少重复数字常量,你切换到符号模式更大的准确性。使用用方程解决()和浮点常量总是一个类别错误:为准确的解决方案,解决()和浮点常量只根据定义近似(除了那些整数)。金宝搏官方网站
Z = @ (v)信谊(v);
%未知变量
信谊Pa N Tθ;
%给变量
factor507 = Z (5.07);
factor93 = Z (0.93);
L = factor507 /罪(θ);
δ=函数(Z (10));
c = Z (10);
γ= Z (20);
函数φ= (Z (20));
Q = c * factor507 / tan(θ);
W = ((factor507 / tan(θ))* factor93 + factor507 * (factor507 / tan(θ))* Z(0.5)) *伽马;
%方程
eq1 = Pa-N * sin(θ)+ T * cos(θ)= = 0;
eq2 = Q + W-Pa * tan(δ)- n * cos(θ)- t * sin(θ)= = 0;
eq3 = T = = c * L + N * tan(φ);
%的解决方案
索尔=解决([eq1、eq2 eq3], [N T Pa],“returnconditions”,真正的);
Pa_theta = symfun(简化(sol.Pa),θ);
eq4 = diff (Pa_thetaθ)= = 0;
theta_max =解决(eq4θ);
rad2deg (theta_max)
答案(1)
埃米尔哈姆萨
2020年4月4日
试试这个例子
信谊a b c d;
r = b + d;
e =π/ 2;
eq1 = c - * sin (d) + b * cos (d) = = 0;
eq2 = r - c * tan (e) - * cos (d) - b * sin (d) = = 0;
eq3 = b = = e * r + * tan (e);
索尔=解决([eq1 eq2 eq3], [a b c]);
C = sol.c;% c函数的d
%现在发现它的导数
dC_dd = diff (C, d);% dc / dd
%你也可以将其集成
IC = int (C, d);
