第八周 MATLAB在通信系统中的应用.ppt

一、周期信号频谱的MATLAB实现 频谱Fn一般为复数，可分别利用abs和angle函数获得其幅频特性和相频特性。 其调用格式分别为 x=abs(Fn) y=angle(Fn) 周期信号的频谱Fn 为离散信号，可以用stem 例1 试用MATLAB画出图示周期三角波信号的频谱。 解：周期信号的频谱为 画三角波信号频谱的MATLAB程序 N=8; n1= -N:-1; %计算n=-N到-1的Fourier系数 c1= -4*j*sin(n1*pi/2)/pi^2./n1.^2; c0=0; %计算n=0时的Fourier系数 n2=1:N; %计算n=1到N的Fourier系数 c2= -4*j*sin(n2*pi/2)/pi^2./n2.^2; cn=[c1 c0 c2]; n= -N:N; subplot(2,1,1); stem(n,abs(cn));ylabel('Cn的幅度'); subplot(2,1,2); stem(n,angle(cn)); ylabel('Cn的相位');xlabel('\omega/\omega0'); 程序运行结果 利用MATLAB进行系统频域分析 连续系统频响特性的计算 周期信号通过系统的响应 离散系统频响特性的计算 一、连续系统频响特性的计算 b 分子多项式系数 a 分母多项式系数 w 需计算的H(jw)的抽样点 (数组w中少需包含两个w的抽样点)。 计算频响的MATLAB函数 H=freqs(b,a,w) 一、连续系统频响特性的计算 例1 三阶归一化的Butterworth低通滤波器的系统函数为 w=linspace(0,5,200); b=[1];a=[1 2 2 1]; h=freqs(b,a,w); subplot(2,1,1); plot(w,abs(h)); subplot(2,1,2); plot(w,angle(h)); 试画出|H(jw)| 和?(w)。 一、连续系统频响特性的计算 三阶Butterworth低通滤波器的幅度响应和相位响应 二、离散系统频率响应的计算 计算频率响应的MATLAB函数 b 分子的系数 a 分母系数 w 抽样的频率点(至少2点)， w在0~2p之间 幅频特性: abs, 相频特性: angle h = freqz(b,a,w) b=[1]; w=linspace(0,2*pi,512); h2=freqz(b,a2,w); plot(w/pi,abs(h1),w/pi,abs(h2),':'); 二、离散系统频响特性的计算 a1=[1 -0.9]; a2=[1 0.9]; h1=freqz(b,a1,w); legend('\alpha=0.9','\alpha=-0.9'); 二、离散系统频响特性的计算 利用MATLAB进行连续系统的复频域分析 部分分式展开的MATLAB实现 H(s)的零极点与系统特性的MATLAB计算 一、部分分式展开的MATLAB实现 [r,p,k]=residue(num,den) num,den分别为F(s)分子多项式和分母多项式的系数向量。 r为部分分式的系数，p为极点，k为多项式 二、H(s)的零极点与系统特性的MATLAB计算 计算多项式根roots的函数可用于计算H(s)的零极点。 r=roots(N) %计算多项式N的根 H(s)零极点分布图可用pzmap函数画出，调用形式为 pzmap(sys) 表示画出sys所描述系统的零极点图。 例1 用部分分式展开法求F(s)的反变换。 format rat %将结果数据以分数的形式输出 num=[1 2]; den=[1 4 3 0]; [r,p]=residue(num,den) 运行结果为 r = -1/6 ,-1/2 ,2/3 p = -3 ,-1 ,0 故F(s)可展开为 例2 用部分分式展开法求F(s)的反变换。 num=[2 3 0 5]; den=conv([1 1],[1 1 2]); %将因子相乘的形式转换成多项式的形式 [r,p,k]=residue(num,den) magr=abs(r) %求r的模 angr=angle(r) %求r的相角 例2 用部分分式展开法求F(s)的反变换。 运行结果为 r =-2.0000 + 1.1339i, -2.0000 - 1.1339i, 3.0000 p =-0.5000 + 1.3229i, -0.5000 - 1