Question

Matlab 列的方程无法求出正确的解

程序如下clc,clearx1 = [-1.2352 -1.2081 -1.1813 -1.1546 -1.1281 -1.1018 -1.0756 -1.0496 -1.0237 -0.998 -0.9724 -0.947 -0.9217 -0.8965 -0.8714 -0.8464 -0.8215 -0.7967 -0.7719 -0.7473 -0.7227];y1 = [0.173 0.189 0.2048 0.2203 0.2356 0.2505 0.2653 0.2798 0.294 0.308 0.3218 0.3354 0.3488 0.3619 0.3748 0.3876 0.4001 0.4124 0.4246 0.4366 0.4484];x2 = [1.1637 1.2212 1.2791 1.3373 1.396 1.4552 1.5148 1.575 1.6357 1.697 1.7589 1.8215 1.8848 1.9488 2.0136 2.0792 2.1457 2.2131 2.2815 2.3508 2.4213];y2 = [3.336 3.3299 3.3242 3.3188 3.3137 3.3091 3.3048 3.3007 3.2971 3.2937 3.2907 3.2881 3.2859 3.284 3.2824 3.2813 3.2805 3.2801 3.2801 3.2804 3.2812];l1 = (x1.^2+y1.^2).^0.5;l2 = (x2.^2+y2.^2).^0.5;T1=12.683:0.05:13.683;T2=13.15:0.05:14.15;syms alpha beta H;eq1=cosd(alpha)*sind(T1(1)+(beta-120)/15-12)/sind(atand(y1(1)/x1(1)))-l1(1)/(H^2+l1(1)^2)^0.5;eq2=cosd(alpha)*sind(T1(2)+(beta-120)/15-12)/sind(atand(y1(2)/x1(2)))-l1(2)/(H^2+l1(2)^2)^0.5;eq3=cosd(alpha)*sind(T1(3)+(beta-120)/15-12)/sind(atand(y1(3)/x1(3)))-l1(3)/(H^2+l1(3)^2)^0.5;[alpha,beta,H]=vpasolve([eq1==0,eq2==0,eq3==0],[alpha,beta,H])求出来的解很奇怪，而且肯定不符合要求附图里有方程的表示，其中alpha就是φ，beta和β，H和H是对应的，这三个是未知数；通过三组T1、y1、x1、l1来列出三个等式求三个未知数至于T2,y2,x2,l2则是另一组数据，处理方法同上

Answer 1

由于给出的方程比较复杂，用matlab的solve函数是很难得到x(a,c)的表达式。如a、c为已知值，则可以（用vpasolve,fsolve等）得到x的数值解。

追问

我用的已经是vpasolve了

Answer 2

.2081 -1.1813 -1.1546 -1.1281 -1.1018 -1.0