Question

matlab solve解方程结果不对

matlab solve解方程结果不对r=1.9;
d=0.175;
ac=11.5;
bc=6.5;
l=11.66;
cc=pi;

aa=solve('(bc*(1-cos(aa))）^2+(ac+r*cos(cc)-bc*sin(aa))^2+(r*sin(cc))^2-l^2','aa');
x=subs(aa);
x=double(aa);
%把得出来的两个解带进去算不对，这怎么解决，求大神帮忙。

Answer 1

对于你的问题，最好这样来处理。

r=1.9;d=0.175;ac=11.5;bc=6.5;l=11.66;cc=pi;

syms aa real

aa=solve((bc*(1-cos(aa)))^2+(ac+r*cos(cc)-bc*sin(aa))^2+(r*sin(cc))^2-l^2);

vpa(aa)

ans =

 -0.32170999022918155334949100099991

   2.2729436161889969111963218680334

将得到的两个解回带验算，结果正确，等于0

追问

如果我的cc是一个从0增加到pi的数组，但是我想求的结果要在±pi╱2范围内

他这个绿色红色的分别是两个解的图像，下面蓝色的是cc的正弦函数，我想要的是在0附近波动的解，这怎么办啊。

追答

可以用循环语句来实现。
clc,clear all,close all
cc1=-pi/2:pi/10:pi/2;
n=length(cc1);
i=1;
while i<=n
r=1.9;d=0.175;ac=11.5;bc=6.5;l=11.66;
syms aa
cc=cc1(i);
aa=solve((bc*(1-cos(aa)))^2+(ac+r*cos(cc)-bc*sin(aa))^2+(r*sin(cc))^2-l^2);
aa(i)=aa;
i=i+1;
end
real(double(aa))

追问

你跑过吗？我跑出来是这样

追答

是这样的

追问

中间这么多0不对吧

我想要就是前面红绿线组成的在0附近波动的线

就是标出来的线，可以吗

谢谢啦，我用循环加if函数解决了问题