%%-----------------------------
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,依据图象下面的等式是对的。
全部