请MATLAB高手帮我解这道隐函数方程

[Lixiney]

急！请MATLAB高手帮我解这道隐函数方程
已知角1相应度数，求角3相应度数，我解了很久都解不出来，好急！
麻烦高手解围！


其中L1==70cm，L2=130cm，L3=90cm，L4=150cm, W=70rad/s

[hrt]

可利用fsolve函数求解；
先编写函数M文件，如下：
function y=f(x,a) %x是数组，存储theta1、theta2，a是theta1的值
y(1)=x(1)-a;
l1=0.7;l2=1.3;l3=0.9;l4=1.5;
y(2)=l1*cos(x(1))+l2*sqrt(1-((l3*sin(x(2))-l1*sin(x(1)))/l2)^2)-l3*cos(x(2))-l4;

再在matlab命令行中输入：
>> x0=[0,0.2];
>> [x,fv,ef]=fsolve(@f,x0,[],0)
即可得到：
x =

0 1.4033


fv =

1.0e-007 *

0 0.5511

即：theta1=0时，theta2=1.4033（rad）
其它情况依此即可，要注意初值x0的选取，否则可能不收敛。
~ ~

[Lixiney]

非常的谢谢您！很好，很开心！

[laosam280]

我把周期一个周期内计算了101个点，并绘制了图像。应该比较直观。txt里有这些点的结果和误差。