matlab系统频率特性实验.ppt

CT信号的频域分析及CT系统的频率响应;1.傅立叶变换;MATLAB提供了能直接求解傅立叶变换及逆变换的函数fourier()、ifourier(),调用格式如下：

傅立叶变换：

(1)F=fourier(f):表示求符号函数f的FT,默认返回关于ω的函数；

(2)F=fourier(f,v):表示求符号函数f的FT,返回关于v的函数；;傅立叶逆变换：

(1)f=ifourier(F):函数F的IFT,默认返回关于t的函数

(2)f=ifourier(F,u):符号函数F的IFT,返回关于u的函数

说明：在调用函数fourier()及ifourier()之前，要用sym命令对所用到的变量进行说明。;注意：;例1：;subplot(2,1,1)

ezplot(exp(-2*t)*(1/2+1/2*sign(t)))

subplot(2,1,2)

ezplot(abs(F));例2.单位阶跃信号;例2.单位阶跃信号;例：用matlab画出u(t+3)-u(t)的波形;例2：单位阶跃信号;例3：直流信号;例4：单位冲激信号;例5：符号函数;例6：用MATLAB计算的FT;例6：矩形脉冲信号（时宽为2）;h(t);MATLAB提供了专门对连续系统的频率响应H(jω)进行分析的函数freqs(),该函数可求出H(jω)的数值解，并能绘出系统的幅频特性和相频特性曲线。

调用格式：

(1).h=freqs(b,a,ω1:p:ω2)，其中

b=[bm,bm-1,…b2，b1]，a=[an,an-1,…a2，a1]

ω1:p:ω2：ω1为频率起始值；ω2为频率终止值

p为频率的取样间隔。

该调用格式将系统频率响应的样值返回给h。

;(2).[h,w]=freqs(b,a)，

该调用格式将计算默认频率范围内200个频率点上系统频率响应的样值返回给h，并将这200个频率点记录在w中。

(3).[h,w]=freqs(b,a,n)，

该调用格式将计算默认频率范围内n个频率点上系统频率响应的样值返回给h，并将这n个频率点记录在w中。;(4).freqs(b,a)

该调用格式不返回系统频率响应的样值，而是以对数坐标的方式绘出系统的幅频特性和相频特性曲线。

例1：a=[10.41];b=[100];

freqs(b,a)

gridoff;3.低通滤波器;该电路的频率响应H(jw)为：;b=[01];

a=[11];

[h,w]=freqs(b,a,100);

h1=abs(h);

h2=angle(h);

subplot(2,1,1);

plot(w,h1);

grid

注：angle(h)表示求复数h的相角，单位为弧度。;xlabel(w)

ylabel(幅度)

title(幅频特性);

subplot(2,1,2);

plot(w,h2*180/pi);

grid

xlabel(w)

ylabel(相位)

title(相频特性);;RC=0.1;RC=10;习题;

;3.用MATLAB计算Sa(t)的FT并画出波形;习题4：下图是一个简单的高通滤波器，画出其幅频特性和相频特性曲线。(分别令RC=1，RC=0.001及RC=0.00001比较其变化)