利用MATLAB计算函数值求解

clc;clear;
% syms x ;
% y = 1/sqrt((1-x^2/500)^2+0.0003978*x^2);
% y1=matlabFunction(y)
% x = [6.28 12.57 18.85 20.11 20.86 21.36 21.87 25.13 34.36 43.98];
% y2=y1(x)
% l1=20*log10(y2)

% syms w  ww;
% phi1=(-atan(0.01968*w/(1-w^2/500)))*180/pi
% phi2=-(180-(atan(0.01968*ww/(ww^2/500-1)))*180/pi)
% 
% phi1m=matlabFunction(phi1)
% phi2m=matlabFunction(phi2)
% 
% w = [6.28 12.57 18.85 20.11 20.86 21.36 21.87 ]
% ww=[25.13 34.36 43.98]
% 
% phi1v=phi1m(w)
% phi2v=phi2m(ww)

syms p t;
a1=p/t
a2=matlabFunction(a1)
p=[2160 2690 4035 4347 4487 4616 4649 4056 1507 703]
t=2000;
A=a2(p,t)
L=20*log10(A)

syms w  ww;
phi1=(-atan(0.01968*w/(1-w^2/21.36^2)))*180/pi
phi2=-(180-(atan(0.01968*ww/(ww^2/21.36^2-1)))*180/pi)

phi1m=matlabFunction(phi1)
phi2m=matlabFunction(phi2)

w = [6.28 12.57 18.85 20.11 20.86 21.36 ]
ww=[21.87 25.13 34.36 43.98]

phi1v=phi1m(w)
phi2v=phi2m(ww)

输出结果：

a1 =
 
p/t
 

a2 =

  包含以下值的 function_handle:

    @(p,t)p./t

p =

        2160        2690        4035        4347        4487        4616        4649        4056        1507         703

A =

    1.0800    1.3450    2.0175    2.1735    2.2435    2.3080    2.3245    2.0280    0.7535    0.3515

L =

    0.6685    2.5744    6.0963    6.7432    7.0185    7.2647    7.3266    6.1414   -2.4583   -9.0815

 
phi1 =
 
(180*atan((123*w)/(6250*((8796093022208*w^2)/4013213922945191 - 1))))/pi
 
 
phi2 =
 
(180*atan((123*ww)/(6250*((8796093022208*ww^2)/4013213922945191 - 1))))/pi - 180
 

phi1m =

  包含以下值的 function_handle:

    @(w)(atan((w.*1.968e-2)./(w.^2.*2.191782743480761e-3-1.0)).*1.8e+2)./pi

phi2m =

  包含以下值的 function_handle:

    @(ww)(atan((ww.*1.968e-2)./(ww.^2.*2.191782743480761e-3-1.0)).*1.8e+2)./pi-1.8e+2

w =

    6.2800   12.5700   18.8500   20.1100   20.8600   21.3600

ww =

   21.8700   25.1300   34.3600   43.9800

phi1v =

   -7.7045  -20.7283  -59.1921  -73.9822  -83.5696  -90.0000

phi2v =

  -96.4060 -127.8382 -156.9300 -165.0409