如何用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-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大（共振点）；正弦振动在任一瞬间只包含一种频率的振动，而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率，所以样品上的共振点会同时激发...
点击进入详情页
本回答由富港检测技术（东莞）有限公司_提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
