MATLAB仿真实例

1、一、实验目的1二、实验题目1三、实验内容13.1 傅里叶变换与傅里叶反变换13.2 题目一：正弦信号波形及频谱23.2.1仿真原理及思路23.2.2程序流程图33.2.3仿真程序及运行结果33.2.4实验结果分析53.3 题目二：单极性归零(RZ)波形及其功率谱53.3.1 仿真原理及思路53.3.2 程序流程图63.3.3 仿真程序及运行结果63.3.4 实验结果分析83.4 题目三：升余弦滚降波形的眼图及其功率谱83.4.1 仿真原理及思路83.4.2 程序流程图83.4.3 仿真程序及运行结果83.4.4 实验结果分析：103.5 题目四：完成PCM编码及解码的仿真113.5.1 仿真原

2、理及思路113.5.2 程序流程图123.5.3 仿真程序及运行结果123.5.4 实验结果分析153.6附加题一：最佳基带系统的PeEbNo曲线，升余弦滚降系数a=0.5，取样值的偏差是Ts/4163.6.1 仿真原理及思路163.6.2 程序流程图163.6.3 仿真程序及运行结果163.6.4 实验结果分析183.7附加题二：试作出PeEb/No曲线。升余弦滚降系数a=0.5,取样时间无偏差，但信道是多径信道，C(f)=|1-0.5-j2ft|,t=Ts/2183.7.1 仿真原理及思路183.7.2 程序流程图193.7.3 仿真程序及运行结果193.7.4 实验结果分析21四、实验心

3、得21一、实验目的 学会MATLAB软件的最基本运用。MATLAB是一种很实用的数学软件，它易学易用。MATLAB对于许多的通信仿真类问题来说是很合适的。 了解计算机仿真的基本原理及方法，知道怎样通过仿真的方法去研究通信问题。 加深对通信原理课程有关内容的理解。二、实验题目必做题：1) 正弦信号波形及频谱；2) 单极性归零(RZ)波形及其功率谱，占空比为50%；3) 升余弦滚降波形的眼图及其功率谱。滚降系数为0.5。发送码元取值为0、2。 选做题：完成PCM编码及解码的仿真。 附加题：1) 最佳基带系统的PeEbNo曲线，升余弦滚降系数a=0.5，取样值的偏差是Ts/4；2) 试作出PeEb/

4、No曲线。升余弦滚降系数a=0.5,取样时间无偏差，但信道是多径信道，C(f)=|1-0.5-j2f|t,t=Ts/2。三、实验内容3.1 傅里叶变换与傅里叶反变换对于确定信号f(t)，其傅里叶变换为：F(f)=Lf(t)e-2jnftdtF(f)傅里叶反变换为：f(t)=LF(f)e2ff0在通信原理仿真中，傅里叶变换与傅里叶反变换会经常用到，我们可以利用MATLAB的快速傅里叶变换函数fft与快速傅里叶反变换函数ifft编写傅里叶变换子程序与傅里叶反变换子程序。其程序代码如下：傅里叶变换子程序：%傅里叶变换子程序functionX=t2f(x)globaldtdfNtfT%X=t2f(x)

5、%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同，并为2的整幂。%本函数需要一个全局变量dt(时域取样间隔)H=fft(x);X=H(N/2+1:N),H(1:N/2).*dt;end傅里叶反变换子程序：%傅里叶反变换子程序functionx=f2t(X)globaldtdftfTN%x=f2t(X)%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同并为2的整幂%本函数需要一个全局变量dt(时域取样间隔)X=X(N/2+1:N),X(1:N/2);x=ifft(X)/dt;%x=tmp(N/2+1:N),tmp(1:N/2);end3.2 题目一：正弦信号波形及频谱3.2.1仿真

6、原理及思路一般来说，任意信号s(t)是定义在时间区间(-0+)上的连续函数，但所有计算机的CPU都只能按指令周期离散运行，同时计算机也不能处理(-a)这样一个时间段。为此一TTT我们把s(t)按区间-T,T截短为s(t)，再对s(t)按时间间隔At均匀取样得到个样_22TTAt值。仿真时我们用这个样值集合来表示信号s(t)。显然At反映了仿真系统对信号波形的分辨率，At越小则仿真的精确度越高。据通信原理所学，信号被取样以后的频谱是频率的周期函数，其重复周期是丄。如果信号的最高频率为f，那么必须有f丄才能保证不AtHH2At发生混叠失真。我们称B=丄为仿真系统的系统带宽。如果我们的仿真程序中设定

7、的采s2At样间隔是At，那么我们不能用此仿真程序来研究带宽大于B的信号或系统。s此外，信号s(t)的频谱S(f)通常来说也是定义在频率区间(-0+8)上的连续函数，所以仿真频域特性时，我们也必须把S(f)截短并取样。考虑到系统带宽是B，我们把频域的截短s区间设计为-B,B,然后再按间隔Af均匀取样得到2Bs个样值。纣反映了仿真系统在频域上的分辨率。频域离散的信号对应到时域是一个周期信号，其周期为丄。类似前面的Af分析，如果我们的仿真程序中设定的采样间隔是A，那么我们就不能仿真截短时间超过占;的信号。考虑到时域截短时间为T,我们把频域的取样间隔设计为Af=丄。这样一来，T时域的总取样点数及频域

8、的总取样点数都相等，为N=丄。要提高仿真的精度，就必须AtAf降低时域取样间隔At及频域取样间隔Af，也就是要加大总取样点数N。这说明仿真的精度与仿真系统的运算量直接有关。为了处理上的方便，我们今后规定采样点数N为2的整幂。首先，设定正弦信号的采样点数为2k，时域采样间隔为dt=0.01，频域采样间隔为df=1/(N*dt)，所以定义一个时域t的N维矢量t=linspace(-T/2,T/2,N)，定义一个频域f的N维矢量f=linspace(-Bs,Bs,N)；其次，定义一个余弦函数s=cos(2/3*pi*f*t)，其中f=1,并对其进行傅里叶变换及傅00里叶反变换；最后，画出该余弦函数波

9、形图与频谱图。3.2.2 程序流程图产生余弦信号及频谱流程图3.2.3仿真程序及运行结果仿真程序：%实验一：正弦信号波形及其频谱closeallclearallglobaldtdfNtfT%全局变量N=2人14;%采样点数dt=0.01;%时域采样间隔df=1/(N*dt);%频域采样间隔T=N*dt;%截短时间Bs=N*df/2;%系统带宽t=linspace(-T/2,T/2,N);f=linspace(-Bs,Bs,N);s=sin(2/3*pi*t);S=t2f(s);a=f2t(S);figure(l)set(l,'Position',10,350,600,200)%

10、设定窗口位置及大小figure(2)set(2,'Position',10,50,600,200)%设定窗口位置及大小figure(l)as=abs(S);%求模plot(f,as)axis(-2,+2,1.1*min(as),1.1*max(as)xlabel('f(MHz)')ylabel('Ps(f)')gridonfigure(2)plot(t,a)axis(-5,5,1.1*min(a),1.1*max(a);xlabel('t')ylabel('s(t)')gridon程序运行结果：取样点数=2k,k=

11、14，得到如下波形：t324实验结果分析由上图可以看出，余弦函数的频谱为两个脉冲信号的叠加，两个脉冲信号分别在11f二+-和f二一-处。333.3题目二：单极性归零(RZ)波形及其功率谱3.3.1仿真原理及思路P二卜P(f)df二lim十J+Tf(t)2dt功率信号f(t)的平均功率：fTts2T_t。该信号的双边功率谱密度P(f)为：P(f)二其中S(f)是s(t)截短后的傅氏变换，|S(f)|2是s(t)的能TT8T，TTT量谱，是S(t)在截短时间内的功率谱。对于仿真系统，若x是时域取样值矢量，X是对应的傅氏变换，那么X的功率谱便为矢量P=(X.*conj(X)/T。对于采用归零(RZ)

12、及不归零(NRZ)矩形脉冲波形的数字信号，可以用简单的方法信号矢量s。设a是码元矢量，N是总取样点数，M是总码元数，L是每个码元内的点数，那么NRZ信号可这样获得：s=zeros(1,N);forii=1:L,s(ii+0:M-1*L)=a;end对于，若Rt是要求的占空比，dt是仿真系统的时域采样间隔，则RZ信号的产生方法是：s=zeros(1,N);forii=1:Rt/dt,s(ii+0:M-1*L)=a;end首先，利用rand函数产生一个尺寸为1*M的矩阵，其元素按均匀分布随机取值于区间0,1,并用round函数对其四舍五入，得到一个随机产生的0,1序列；其次，利用for循环产生码元

13、长度为L点，码元为0,1，且占空比为50%的单极性归零码，并画出其波形图；最后，计算该RZ的功率谱密度，并画出其波形。3.3.2 程序流程图产生0,1分布的随机序列卜RZ波形及其功率谱密度仿真流程图3.3.3仿真程序及运行结果仿真程序：%实验二：画出单极性归零码及其功率谱closeallclearallglobaldttdfNN=2人14;%采样点数L=64;%每码元的采样点数M=N/L;%码元数Rb=2;%码速率为2Mb/sTs=1/Rb;%码元间隔dt=Ts/L;%时域采样间隔Rt=0.5;%占空比df=1/(N*dt);%频域采样间隔T=N*dt;%截短时间Bs=N*df/2;%系统带宽

14、t=linspace(-T/2,T/2,N);%时域横坐标f=linspace(-Bs,Bs,N);%频域横坐标EP=zeros(1,N);forjj=1:100a=round(rand(1,M);%产生M个取值0，1等概的随机码s=zeros(1,N);%产生一个N个元素的零序列forii=1:Rt*Ts/dts(ii+0:M-1*L)=a;%产生单极性归零码endQ=t2f(s);%付氏变换P=Q.*conj(Q)/T;%P为单极性归零码的功率EP=(EP*(jj-1)+P)/jj;%累计平均endaa=30+10*log10(EP+eps);%加eps以避免除以零figure(1)set

15、(l,'Position',10,350,600,200)%设定窗口位置及大小figure(2)set(2,'Position',10,50,600,200)%设定窗口位置及大小figure(l)plot(f,aa,'g')xlabel('f(MHZ)')ylabel('Ps(f)')axis(-15,+15,1.1*min(aa),1.1*max(aa)gridonfigure(2)plot(t,s,'r')xlabel('t(ms)')ylabel('s(t)(V)

16、9;)axis(-10,+10,1.1*min(s),1.1*max(s)gridon程序结果：输入取样点数=2k,k=14，得到如下波形：3.3.4实验结果分析由上图结果可以看单极性归零码序列的功率谱密度不仅含有离散的直流分量及连续谱,而且还包含离散的时钟分量和奇次谐波分量，功率谱主瓣宽度为码元速率2RbMb/s。根据2数字PAM信号功率谱密度公式P(f)=巧G(f)|2+m£G(-)5(f-)得到：功率谱sTt1T2ttTssmss分为两个部分，第一部分是连续谱，形状取决于G*f);第二部分是离散线谱，相邻线谱频率间隔为1/Ts。若序列的均值ma为零，则第二部分为零，即离散线谱消

17、失，单极性码的均值不为零，故都存在直流分量。3.4题目三：升余弦滚降波形的眼图及其功率谱3.4.1仿真原理及思路眼图是数字信号在示波器上重复扫描得到的显示图形。若示波器的扫描范围是Na个码元，那么画眼图的方法是：tt=0:dt:Na*L*dt;holdonforii=1:Na*L:N-Na*Lplot(tt，s(ii+1:Na*L);end首先，产生一个sinc函数g二沁竺2，并定义升余弦滤波器的传递函数为1兀t/Tsg二C0如“丁)，得到升余弦脉冲为g二Sin(兀“2co如”-)；其次，对升余弦脉2l-(2a兀t/T)2兀t/Tl-(2a兀t/T)2sss冲进行傅里叶变换，画出升余弦信号的频

18、谱图；最后，再对所得函数进行傅里叶反变换并对结果取实数部分，并利用for循环，画出Na个码元宽度的眼图。3.4.2 程序流程图产生升余弦信号卜傅氏变换计算功率谱卜f傅氏逆变换卜升余弦滚降波形的眼图及其功率谱仿真流程图3.4.3仿真程序及运行结果仿真程序：%实验三：升余弦滚降波形的眼图及其功率谱closeallclearall%采样点数%每码元的采样点数%码元数%码速率是2Mb/s%码元间隔%时域采样间隔%频域采样间隔%截短时间%系统带宽%示波器扫描宽度为4个码元globaldttdfNN=2人14;L=32;M=N/L;Rb=2;Ts=1/Rb;dt=Ts/L;df=1/(N*dt);T=N*

19、dt;Bs=N*df/2;Na=4;alpha=input('滚降系数=0.5');t=-T/2+dt/2:dt:T/2;%时域横坐标f=-Bs+df/2:df:Bs;%频域横坐标g1=sin(pi*t/Ts)./(pi*t/Ts);g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).人2);g=g1.*g2;%升余弦脉冲波形G=t2f(g);figure(1)set(1,'Position',10,350,600,200)%设定窗口位置及大小figure(2)set(2,'Position',10,50,600,2

20、00)%设定窗口位置及大小holdongridonaxis(-3,+3,-50,50)xlabel('tinus')ylabel('s(t)inV')EP=zeros(size(f)+eps;forii=1:50a=sign(randn(1,M);a=a+1;imp=zeros(1,N);%产生冲激序列imp(L/2:L:N)=a/dt;S=t2f(imp).*G;%升余弦信号的傅氏变换s=f2t(t2f(imp).*G);%升余弦信号的时域波形s=real(s);P=S.*conj(S)/T;%升余弦信号的功率谱EP=(EP*(ii-1)+P+eps)/ii;

21、endfigure(2)tt=O:dt:Na*L*dt;forjj=l:Na*L:N-Na*Lplot(tt,s(jj:jj+Na*L);axis(0,+2,-1,3)gridonend%作眼图figure(1)plot(f,30+10*log10(EP),'g');gridonaxis(-3,+3,-50,50)xlabel('f(MHz)')ylabel('Ps(f)(dBm/MHz)')程序运行结果：输入滚降系数=0.5，得到如下波形:-1012f(MHz)N工运EBP)%>.EH.、武50=60204oi01t2t4ilte2tin

22、us3.4.4实验结果分析：由上图可以清晰地看出Q=0.5时的眼图，定性的判断此时系统的误码严重程度，可以得到最佳抽样时刻约为0.22us、0.72us、1.22us、1.72us。该程序还可以输入不同的a值得到不同滚降系数时的眼图，a越接近1,信号成形的波形越好，眼图的质量也越好，而在一定码元速率下随着a值的增加，信号占用的带宽就越大,频带利用率下降。3.5题目四：完成PCM编码及解码的仿真3.5.1仿真原理及思路在PCM中，对模拟信号进行抽样、量化，将量化的信号电平值转化为对应的二进制码组的过程称为编码，其逆过程称为译码或解码。从理论上看，任何一个可逆的二进制码组均可用于PCM。但是目前最

23、常见的二进制码有三类：二进制自然码（NBC）、折叠二进制码组（FBC）、格雷二进制码（RBC）。在PCM中实际采用的是折叠二进制码。Ira图3.5.1A律13折线由表3-1可见，如果把16个量化级分成两部分：07的8个量化级对于于负极性样值脉冲，815的8个量化级对应于正极性样值脉冲。自然二进制码就是一般的十进制正整数的二进制表示，在16个量化级中：24=16，采用4位码元表示为：b二23、b二22、二21、123b二2o的有无组合来构成。比如第11个量化级可表示为411=23+0+21+2o=8+0+2+1其对应的码组可表示为：1011,其余依次类推。本程序中采用自然码的编码方式。表3.5.

24、1自然码、折叠码、格雷码样值脉冲极性量化级自然二进制码b1b2b3b4折叠二进制码b1b2b3b4格雷码b1b2b3b4f二8kHz,对每个抽样值进行A律或者卩律对数压缩非均匀量化及非线性编码，每个样值s用八位二进制代码表示，这样，每路标准话路的比特率为64kbps。表3.5.2中给出的是信号正值的编码，负值编码是对称的，其绝对值与此表相同。整个信号动态范围共分13个段落，各段落的量化间隔都不同，并且有2的倍数关系。每个段落内位均匀分层量化，共16层。每个样值用8比特来表示，即bbbbbbbb。这8比12345678特分为三部分：b为极性码，0代表负值，1代表正值。bb称为段落码，表示段落的1

25、24号码，其值为07,代表8个段落。bbbb表示每个段落内均匀分层的位置，其值为015,5678代表一段落内的16个均匀量化间隔。在PCM解码时，根据八比特码确定某段落内均匀分层的位置，然后去其量化间隔的中间值作为量化电平。本程序首先产生一个正弦信号，并对其进行采样量化，生成一个幅值矩阵；然后利用编码子函数对此矩阵中的每个元素按照A律13折线编码规则编码，并产生一个输出码组矩阵；最后利用解码子函数对输出码组矩阵解码，并画出编码前与解码后的波形图。3.5.2 程序流程图A律十三折线编码译码流程图3.5.3仿真程序及运行结果仿真程序：%实验四：PCM编码及解码仿真clearall;closeall

26、;globaldttdfNt=0:0.01:10;x=sin(t);v=max(x);xx=x/v;%normalizesxx=floor(xx*4096);y=pcm_encode(sxx);yy=pcm_decode(y,v);drawnowfigure(1)set(l,'Position',10,350,600,200)%设定窗口位置及大小plot(t,x)title('samplesequence');figure(2)set(2,'Position',10,50,600,200)%设定窗口位置及大小plot(t,yy)title(

27、9;pcmdecodesequence');编码子程序：functionout=pcm_encode(x)%xencodetopcmcoden=length(x);%-4096<x<4096fori=1:nifx(i)>0out(i,1)=1;elseout(i,1)=0;endifabs(x(i)>=0&abs(x(i)<32out(i,2)=0;out(i,3)=0;out(i,4)=0;step=2;st=0;elseif32<=abs(x(i)&abs(x(i)<64out(i,2)=0;out(i,3)=0;out(i

28、,4)=1;step=2;st=32;elseif64<=abs(x(i)&abs(x(i)<128out(i,2)=0;out(i,3)=1;out(i,4)=0;step=4;st=64;elseif128<=abs(x(i)&abs(x(i)<256out(i,2)=0;out(i,3)=1;out(i,4)=1;step=8;st=128;elseif256<=abs(x(i)&abs(x(i)<512out(i,2)=1;out(i,3)=0;out(i,4)=0;step=16;st=256;elseif512<=a

29、bs(x(i)&abs(x(i)<1024out(i,2)=1;out(i,3)=0;out(i,4)=1;step=32;st=512;elseif1024<=abs(x(i)&abs(x(i)<2048out(i,2)=1;out(i,3)=1;out(i,4)=0;step=64;st=1024;elseif2048<=abs(x(i)&abs(x(i)<4096out(i,2)=1;out(i,3)=1;out(i,4)=1;step=128;st=2048;endifabs(x(i)=4096out(i,2:8)=1111111;

30、elsetmp=floor(abs(x(i)-st)/step);t=dec2bin(tmp,4)-48;%函数dec2bin输出的是ASCII字符串，48对应0out(i,5:8)=t(1:4);endendout=reshape(out',1,8*n);解码子程序：functionout=pcm_decode(in,v)%decodetheinputpcmcode%in:inputthepcmcode8bitssample%v:quantizedleveln=length(in);in=reshape(in',8,n/8)'slot(1)=0;slot(2)=32;

31、slot(3)=64;slot(4)=128;slot(5)=256;slot(6)=512;slot(7)=1024;slot(8)=2048;step(1)=2;step(2)=2;step(3)=4;step(4)=8;step(5)=16;step(6)=32;step(7)=64;step(8)=128;fori=1:n/8ss=2*in(i,1)-1;tmp=in(i,2)*4+in(i,3)*2+in(i,4)+1;st=slot(tmp);dt=(in(i,5)*8+in(i,6)*4+in(i,7)*2+in(i,8)*step(tmp)+0.5*step(tmp);out(

32、i)=ss*(st+dt)/4096*v;end程序结果：3.5.4实验结果分析从上图可以看出，该程序能将正弦信号正确采样量化编码解码。正弦信号经A律13折线编码解码后，两条曲线基本一致。此图可以清晰地看出当信号幅度较小时，解码后的波形是平滑的曲线，而当信号幅度较大时（比如波峰和波谷附近）可以看到细小的锯齿状波形，此现象说明了A率13折线编码当输入信号小时，段落小，量化级间隔小；当输入信号大时，段落大，量化级间隔大的特点，可以有效减低小信号的量化误差。3.6附加题一：最佳基带系统的PeEbNo曲线，升余弦滚降系数a=0.5，取样值的偏差是Ts/43.6.1仿真原理及思路最佳基带系统框图如下所示

33、白噪声n0/2图3.6.1最佳基带系统首先，产生一个随机基带信号，并计算出信道参数；其次，将此基带信号依次通过发送滤波器、传输信道、接收滤波器，得到一个接收信号的矩阵；然后，对比发送矩阵和接收矩阵，当出现不同码元时误码增加1；最后，画出基带系统的PeEbNo曲线。3.6.2 程序流程图最佳基带系统PeEbNo曲线设计流程图3.6.3仿真程序及运行结果仿真程序：closeallclearallglobaldttfdfNT%全局变量N=2人14;%采样点数L=8;%每码元的采样点数M=N/L;%码元数Rb=2;%码速率是2Mb/sTs=1/Rb;%码元间隔dt=Ts/L;%时域采样间隔df=1/(

34、N*dt);%频域采样间隔T=N*dt;%截短时间Bs=N*df/2;%系统带宽alpha=0.5;%滚降系数=0.5t=linspace(-T/2,T/2,N);%时域横坐标f=linspace(-Bs,Bs,N)+eps;%频域横坐标figure(1)set(1,'Position',10,50,300,200)%设定窗口位置及大小hr1=sin(pi*t/Ts)./(pi*t/Ts);hr2=cos(alpha*pi*t/Ts)./(l-(2*alpha*t/Ts).人2);hr=hr1.*hr2;HR=abs(T2F(hr);GT=sqrt(HR);GR=GT;forl

35、oopl=l:20Eb_N0(loopl)=(loopl-l)%分贝值变为真值eb_n0(loop1)=10人(Eb_N0(loop1)/10);Eb=l;n0=Eb/eb_n0(loop1);sita=n0*Bs;n_err=0;%信道噪声谱密度%噪声功率%误码计数forloop2=1:5a=sign(randn(1,M);%发送码元imp=zeros(1,N);imp(L/2:L:N)=a/dt;IMP=T2F(imp);n_ch=sqrt(sita)*randn(size(t);%信道噪声nr=real(F2T(T2F(n_ch).*GR);sr=real(f2t(IMP.*HR)+nr

36、;y=sr(L/4:L:N);%以L/4为起点，步长L,取样点Naa=sign(y);%接收码元n_err=n_err+length(find(aa=a);endPe(loop1)=n_err/(M*loop2);%误码率Pe'=n_err/(M*loop2)figure(1)set(1,'Position',10,10,600,400)%设定窗口位置及大小semilogy(Eb_N0,Pe,'g');eb_nO=1O.人(Eb_N0/10);%还原为真值holdonsemilogy(Eb_N0,0.5*erfc(sqrt(eb_n0);axis(0,9

37、,le-4,l)xlabel('Eb/NO')ylabel('Pe')end实验结果：误码率曲线：绿色线为实际曲线，蓝色为理论误码率曲线Eb/NO3.6.4实验结果分析由上图可以看处，当取样时间偏差为Ts/4的时候，误码率明显高于无偏差的理论误码率.原因可以从信号的眼图看出，没有偏差时的取样点在眼睛睁开最在处，判决效果最佳,而偏差Ts/4之后取样信号的幅度下降，比较容易受噪声干扰,误码率上升。3.7附加题二：试作出PeEb/No曲线。升余弦滚降系数a=0.5，取样时间无偏差，但信道是多径信道，C(f)=l1-0.5-j2ft|,t=Ts/23.7.1仿真原理及思

38、路仿真原理同最“佳基带系统的PeEbNo曲线”仿真原理，只是信道不同。3.7.2程序流程图最佳基带系统PeEbNo曲线设计流程图3.7.3仿真程序及运行结果仿真程序：%多径信道传输的PeEb/No曲线。升余弦滚降系数a=0.5,取样时间无偏差。closeallclearallglobaldttfdfNTN=2人14;%采样点数L=32;%每码元的采样点数M=N/L;%码元数Rb=2;%码速率是2Mb/sTs=1/Rb;%码元间隔dt=Ts/L;%时域采样间隔df=1/(N*dt);%频域采样间隔T=N*dt;%截短时间Bs=N*df/2;%系统带宽alpha=0.5;%滚降系数t=linspace(-T/2,T/2,N);%时域横坐标f=linspace(-Bs,Bs,N)+eps;%频域横坐标figure(1)set(1,