利用matlab解方程组，并且画出方程组中变量随另一自变量变化的曲线

请问可以根据下面8个方程在matlab中得到FAB,FIJ,FEF,FCD随ob的变化曲线吗？这八个方程里面共8个未知参数，即oa,oc,o5,o6,FAB,FIJ,FEF,FCD。然后od=12°，ob=o°：90°，G=95000。

eq1=FAB*sin(oa-ob)+FIJ*sin(o5-ob)+FEF*(o6-ob)+FCD*sin(oc-ob);

eq2=FAB*cos(oa-ob)+FIJ*cos(o5-ob)+FEF*cos(o6-ob)+FCD*cos(oc-ob); eq3=G*2970.03*cos(ob+10/180*pi)+FIJ*400*sin(o5-ob)+FEF*1000*sin(o6+ob)+FCD*1408.86*sin(oc-ob); eq4=G*2576.41*cos(ob+11/180*pi)+FAB*400*sin(oa-ob)+FEF*600*sin(o6+ob)+FCD*1008.86*sin(oc-ob); eq5=G*1990.75*cos(ob+14/180*pi)+FAB*1000*sin(oa-ob)+FIJ*600*sin(o5-ob)+FCD*408.86*sin(oc-ob); eq6=G*1597.78*cos(ob+18/180*pi)+FAB*1408.86*sin(oa-ob)+FIJ*1008.86*sin(o5-ob)+FEF*408.86*sin(o6+ob);

eq7=931.82*sin(oa)+1408.86*sin(ob)-1785.94*sin(oc)-500*sin(od); eq8=931.82*cos(oa)+1408.86*cos(ob)-1785.94*cos(oc)-500*cos(od);

Matlab求解方程组有解析解法(符号解法)和数值解法之分，这里建议使用数值解法，循环给出ob，数值解法求出自变量，构造向量后绘图。