如何用matlab求解如图方程组
1个回答
展开全部
根据题主给出的关系式，经简化得到是二元一次方程，即
(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-10-21 广告
上海华然企业咨询有限公司专注于AI与数据合规咨询服务。我们的核心团队来自头部互联网企业、红圈律所和专业安全服务机构。凭借深刻的AI产品理解、上百个AI产品的合规咨询和算法备案经验，为客户提供专业的算法备案、AI安全评估、数据出境等合规服务，...
点击进入详情页
本回答由上海华然企业咨询提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
