修改matlab程序
%%-----------------------------
syms t w dt real %定义符号变量，real表示符号变量定义在实数范围内。
dt=pi/180,w=1*5*pi/9, %%--角度要变换成弧度，w有三个值，改一下，计算一次。
for k=1:40 %由于本方程不容易得到解析解，所以，用循环变量，k从1增加到50，每步默认增加1。
f=k, %把变化的循环变量赋值给f舵机频响。
c=1/(4*t...全部
%%-----------------------------
syms t w dt real %定义符号变量，real表示符号变量定义在实数范围内。
dt=pi/180,w=1*5*pi/9, %%--角度要变换成弧度，w有三个值，改一下，计算一次。
for k=1:40 %由于本方程不容易得到解析解，所以，用循环变量，k从1增加到50，每步默认增加1。
f=k, %把变化的循环变量赋值给f舵机频响。
c=1/(4*t*f), %------------------------------------------
s=w*t*((c-1)+exp((-1*c)))-dt/2^0。5, %%%把关系式列出。
ss=subs(s), % 把s中的已知的变量（包括具体数值和符号表达式），代入。
ff=solve(ss,'t') , %解方程，求ss=0的解，把舵机时间常数t看作变量解方程。
t1=double(ff(1)),% w=100时，循环到36，即f=36，方程无解，循环被终止，这时，只得到35个值，直接plot（）画图即可。
%t2=double(ff(2)), %由于上式解方程中，有两值，所以，分开赋值给t1、t2并把它们变换为双精度数值，以防止以分数形式给出。
fx(k)=f,tx1(k)=t1,%tx2(k)=t2, %把变量赋值给能循环记录的数组。
end %循环计算完毕，结束。
fx=double(fx),
hold on, %固定，以后的图象都在一张图上。
plot(fx,tx1), %画出ff一个解的图象。
▲▲▲▲▲▲▲▲▲▲▲▲
程序已经修改完毕，花点时间自己看懂了，以后就可以随你修改了。（4。4。9）中的有错误，上面的等式中分母少了一个4，依据图象下面的等式是对的。收起
