- 由于我不想用matlab的bode()函数,所以就想自己写一下简单的bode()函数;
- 碰到的第一个问题是:十倍频程的横坐标怎么形成?我查了一下发现有一个函数是logspace(),它的说明是这样的。
y = logspace(a,b,n)
generates n points between decades 10^a and 10^b. 我试了一下
n=logspace(0,4,5);
结果n=[1 10 100 1000 10000]
,因此横坐标可以做到十倍频程了,指定n可以把十倍频程的密度变大。代码是
%
% 1
% H(s) = -------
% 1+s
%
f=logspace(-2,2,1e5);
w=2*pi*f;
Theta=atan(-w).