syms t f w dt real %定义符号变量,real表示符号变量定义在实数范围内。
dt=pi/180,w=1*5*pi/9, %%--角度要变换成弧度,w有三个值,改一下,计算一次。
for k=1:50 %由于本方程不容易得到解析解,所以,用循环变量,k从1增加到50,每步默认增加1。
t=k, %把变化的循环变量赋值给舵机时间常数t。
c=1/(4*t*f),
s=w*t*((c-1)+exp((-1*c)))-dt/2^0。5, %%%把关系式列出。
ss=subs(s), % 把s中的已知的变量(包括具体数值和符号表达式),代入。
ff=solve(ss,'f') , %解方程,求ss=0的解,把f舵机频响看作变量解方程。
f1=double(ff(1)),f2=double(ff(2)), %由于上式解方程中,有两值,所以,分开赋值给f1、f2并把它们变换为双精度数值,以防止以分数形式给出。
tx(k)=t,fx1(k)=f1,fx2(k)=f2, %把变量赋值给能循环记录的数组。
end %循环计算完毕,结束。
plot(tx,fx1), %画出ff一个解的图象。
hold on, %固定,以后的图象都在一张图上。
plot(tx,fx2) %化出第二个解的图象。
%-----------------------------------------
%%图象画好后,可以直接编辑修改,比如,删除其中不要的曲线,或者改变线型、文字等等。
。
全部