Question

如何用MATLAB编写下面函数

其中Fa=22; Z=17; Db=6.85; k=2.222x10^-4;q=0.04;a0=15
求解a在（0~90）的数值解？

Answer 1

Fa=22;
Z=17;
Db=6.85;
k=2.222*10^-4;
q=0.04;
a0=15;
f=sym('(cosa0/cosx-1)*sinx^(2/3)=k/q*(Fa/ZDb)^(2/3)');
x=solve(f)

 结果：

x =
 
exp(-pi*m*3*i)*((cosx*k*(Fa/ZDb)^(2/3))/(cosa0*q - cosx*q))^(3/2)

追问

我想要得到它的数值解，不知道怎么可以算出来？

追答

f=sym('(cos(15)/cos(x)-1)*sin(x)^(2/3)-2.222*10^-4/0.04*(22/(17*6.85))^(2/3)');
x=solve(f);

结果

x =
 
98.09592137875194398171547628755 - 0.0024808294545748299441784137043216*i