如何用matlab求解如图方程组
展开全部
根据题主给出的关系式，经简化得到是二元一次方程，即
(0.5*(sin(θ)-cos(θ))*h-sin(θ)*h)^2+(sqrt(3*（(sin(θ)-cos(θ))*h）^2+(cos(θ))^2*h^2-2*sin(θ)*（(sin(θ)-cos(θ))*h）*h)+1.5*（(sin(θ)-cos(θ))*h）)^2=4*（(sin(θ)-cos(θ))*h）^2+h^2
是关于h（θ）的函数。所以该只能得到数值解。
也就是说，当θ=0～π时，利用for循环语句和vpasolve函数求得对应θ值的h值。
实现代码：
x0=0:pi/10:pi;
for i=1:length(x0)
syms h
%x=θ
x=x0(i);
r=(sin(x)-cos(x))*h;
eq1=(0.5*r-sin(x)*h)^2+(sqrt(3*r^2+(cos(x))^2*h^2-2*sin(x)*r*h)+1.5*r)^2;
eq2=4*r^2+h^2;
h=vpasolve(eq1==eq2);
h0(i)=h;
end
[x0' h0'] %显示θ和h值
(0.5*(sin(θ)-cos(θ))*h-sin(θ)*h)^2+(sqrt(3*（(sin(θ)-cos(θ))*h）^2+(cos(θ))^2*h^2-2*sin(θ)*（(sin(θ)-cos(θ))*h）*h)+1.5*（(sin(θ)-cos(θ))*h）)^2=4*（(sin(θ)-cos(θ))*h）^2+h^2
是关于h（θ）的函数。所以该只能得到数值解。
也就是说，当θ=0～π时，利用for循环语句和vpasolve函数求得对应θ值的h值。
实现代码：
x0=0:pi/10:pi;
for i=1:length(x0)
syms h
%x=θ
x=x0(i);
r=(sin(x)-cos(x))*h;
eq1=(0.5*r-sin(x)*h)^2+(sqrt(3*r^2+(cos(x))^2*h^2-2*sin(x)*r*h)+1.5*r)^2;
eq2=4*r^2+h^2;
h=vpasolve(eq1==eq2);
h0(i)=h;
end
[x0' h0'] %显示θ和h值
更多追问追答
追答
请题主核查相关表达式是否正确。感觉有点问题。
追问
代码运行后的结果h都为0，我改了一下方程，改成这样的方程该怎么修改代码。
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
东莞大凡
2024-11-19 广告
板格标定棋盘是我们东莞市大凡光学科技有限公司在精密光学测量领域的重要工具。它采用高精度设计，确保每一个格板都达到严格的校准标准。通过使用板格标定棋盘，我们能够有效地对光学测量系统进行校准，从而提升测量的准确性和可靠性。这一工具在光学仪器的研...
点击进入详情页
本回答由东莞大凡提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
