信号与系统实验四___MATLAB

傅里叶变换的 MATLAB 求解

symbolic Math Toolbox 中函数 fourier() 及 ifourier() 可直接求取函数的Fourier Transformation
正变换

F = fourier( f );%函数fFourier变换返回默认为ω的函数

F = fourier( f，v );%此时返回v的函数

反变换

f = ifourier( F );% F 的 fourier 逆变换，默认变量 w，默认返回为x的函数

f = ifourier( F, u );%默认返回为u的函数

例4-1Fourier Transformation

syms t
Fw = fourier( exp( -2 * abs( t ) ) )

例4-2Inverse Fourier Transform

syms t w
ft = ifourier( 1 / ( 1 + w^2 ) )

连续时间信号的频谱图

例4-3求取频谱函数

sym t
ft = 4cos(2pi6t)*(heaviside(t+1/4)-heaviside(t-1/4));
Fw = simplify( fourier( ft ) )
subplot( 1, 2, 1 )
ezplot( ft,[ -0.5, 0.5 ] )
grid on;
subplot( 1, 2, 2 )
ezplot( abs( Fw ), [ -24 * pi, 24 * pi ] );
grid on;

例4-4数值计算法求取Fourier变换

R = 0.02;t = -2 : R : 2;
f = Heaviside( t + 1 ) - Heaviside( t - 1 );
W1 = 2 * pi * 5;
N = 500;
k = 0: N;
W = k * W1 / N;
F = f * exp( -j * t’ * W) * R;
F= real( F );
W = [ -fliplr( W ), W( 2 : 501 ) ];
F = [ fliplr( F ), F( 2 : 501 ) ];
subplot( 2, 1, 1 ); plot( t, f );
xlabel( ‘t’ );ylabel( ‘f(t)’ );
title( ‘f(t)=u(t+1)-u(t-1)’ );
subplot( 2, 1, 2 );plot( W, F );
xlabel( ‘w’ );ylabel( ‘F(w)’ );
title( ‘f(t)的付氏变换 F(w)’ );

利用 MATLAB 分析 LTI 系统的频率特性

** H = freqs(b,a,w) ** 其中，a 和 b 分别是 H(jw)的分母和分子多项式的系数向量，w 为形如 w1:p:w2 的向量，定义系统频率响应的频率范围，w1 为频率起始值，w2 为频率终止值,p 为频率取样间隔。
H 返回 w 所定义的频率点上，系统频率响应的样值。

a=[1 2 1]; 						%分母多项式
b=[0 1];						%分子多项式
h=freqs(b,a,0:0.5:2*pi)			%0到2pi,抽取0.5

例4-5低通滤波器的频率响应

w = 0 : 0.025 : 5;
b = [ 1 ]; a = [ 1, 2, 2, 1 ];
H = freqs( b, a, w );
subplot( 2, 1, 1 );
plot( w, abs( H ) );
grid on;
xlabel( ‘\omega(rad/s)’ );
ylabel( ‘|H(j\omega)|’ );
title( ‘H(jw)幅频特性’ );
subplot(2,1,2);
plot(w,angle (H));grid;
xlabel(’\omega(rad/s)’);
ylabel(’\phi(\omega)’);
title(‘H(jw)的相频特性’);

取模abs() 取相角angle()

例4-6电路中的系统响应

RC=0.04;
t=linspace(-2,2,1024);
w1=5;w2=100;
H1=jw1/(jw1+1/RC);
H2=jw2/(jw2+1/RC);
f=cos(5t)+cos(100t);
y=abs(H1)cos(w1t+angle(H1))+ abs(H2)cos(w2t+angle(H2));
subplot(2,1,1);
plot(t,f);
ylabel(‘f(t)’);
xlabel(‘Time(s)’);
subplot(2,1,2);
plot(t,y);
ylabel(‘y(t)’);
xlabel(‘Time(s)’);

课后习题

2

syms t a
F = fourier( exp( -a * t ) * Heaviside(t) )

取a为2

syms t
ft = exp(-2t) * heaviside(t);
Fw = fourier( ft );
subplot(2,2,1)
ezplot( ft, [ 0, 4 ] ),grid on;
xlabel( ‘t’ );ylabel( ‘f(t)’ );
subplot(2, 2, 2 )
ezplot( abs( Fw ), [ -12pi 12pi ] ),grid on;
xlabel( ‘\omega’ );ylabel( '|F(j\omega)|’ );
subplot(2, 2, 3 )
ezplot( angle( Fw ), [ -12pi 12pi ] ),grid on;
xlabel( ‘\omega’ );ylabel( ‘\phi(\omega)’ );

3

w = -12 * pi : 0.025 : 12 * pi;
b = [ 1 ]; a = [ 0.08, 0.4, 1 ];
H = freqs( b, a, w );
subplot( 2, 1, 1 );
plot( w, abs( H ) );
grid on;
xlabel( ‘\omega(rad/s)’ );
ylabel( ‘|H(j\omega)|’ );
title( ‘H(jw)幅频特性’ );
subplot(2,1,2);
plot(w,angle (H));grid;
xlabel(’\omega(rad/s)’);
ylabel(’\phi(\omega)’);
title(‘H(jw)的相频特性’);

个人感觉是个低通滤波器 参考链接