2023年南京邮电大学matlab软件

南京邮电大学通信学院软件课程设计实验报告模块名称：＿__MATLAB软件设计专业班级:＿_通信工程_17班姓名：_＿__张步涛______学号：____1０0０1０１8实验日期：2023年６月１7—２８日实验报告日期:202３年7月1日一、规定练习的实验部分在时间区间[0,10]中,绘出曲线。程序：t=0：０.１:10;y＝1-exp((-0.5)*ｔ）.*ｃoｓ(2*t）;pｌoｔ(t,y,'r-');sｈg结果：写出生成如图E2-1所示波形的MATＬAB脚本M文献。图中虚线为正弦波，规定它的负半波被置零,且在处被削顶。程序：t=liｎspacｅ(０,3*pi,500);ｙ=sin(t);ａ=ｓｉn(pｉ/3）;z＝（y>=0).*y；z=(y>=a)．＊a+(ｙ<a).*z;pｌot（t，y,＇:ｒ＇）;holdoｎ；plot(t,ｚ,'b-')xlａｂel(＇ｔ＇)，ylabｅl('z=f(t)'),titｌe（'逐段解析函数')legｅnd('ｙ=sin(t)',＇ｚ＝f(t)',４),ｈｏldoff结果:令，运营[U,S,Ｖ]=sｖｄ(A);toｌ=S（1,１)*3*ｅpｓ；然后回答以下问题：sｕm(diag(S)＞toｌ)-ｒank(Ａ)是多少？S(1,１)-norm（A)=0是多少?sqｒｔ(sｕｍ(ｄｉag（S＊Ｓ）)）－noｒｍ（A,'ｆro')的结果是什么？S（1,1）/Ｓ(3,3)-cond(A)的结果是什么?S（1，1）＊S(２，2）*S(3,3)－deｔ（Ａ)<tol的运营结果是什么？V(:,1）'*nulｌ(A）得到什么结果?abs(A＊null（A))<tol得到什么结果?Ｕ（:,1:2）＝=oｒtｈ(A)的运营结果是什么?程序:clear;ｃlc;ｄisp('设A＝[1,2,３；4,５,6;7，8,9］，得＇)A=[1，2,3；4,5,６;7,８,9]%一个3*3矩阵ｄisｐ('设[U,S，V]=svｄ(A）,得')［U,S,Ｖ］=svd(A)%奇异值分解(SVD)ＵSＶ:Ｕm＊ｍ酉矩阵Vn＊n酉矩阵S对角阵dｉsｐ('设tol=S(1,1)*3＊eps，得')ｔｏl=S(1，１)＊3＊epsｄisp('设ａ=suｍ(dｉag(S)>ｔol）－rank(A),求a')ａ=sum(diag（S)>tol)-rank(Ａ)dｉsｐ('设b＝S（1,１)-ｎoｒm(A）,求b')ｂ＝Ｓ(1,1)-norm(Ａ)disp('设c=sqrt(ｓum(diag（S＊S)))－ｎoｒm(A,frｏ),求c'）c＝sｑrt(sｕm(diag(Ｓ*S)))-noｒm(A,'fro')ｄisp('设d=S(1,1)/S(3,3)-cｏｎd(A），求ｄ＇)ｄ=S(1,１)/S(3,3)-ｃｏｎd(A)disp('设e=S（1,１）*Ｓ(２,2)＊S（3，3）-deｔ(A)<tol,求e')e＝Ｓ（1,1）＊Ｓ（2,２)*S(3,3)-dｅｔ（Ａ)<toldisp('设f=V(:,1)*null(A)，求f＇)ｆ=Ｖ（:,1）＇*ｎulｌ(A)dｉsp（＇设ｇ=abs(Ａ*null(A))<tol,求ｇ')g=ａbs(Ａ*nuｌl(Ａ))<toｌｄisp（'设h=U(：，1：２)＝=orth(Ａ),求ｈ')h=U（:,１:2)=＝orth(A)结果：设A=[1,2,3;4,5,6；7，8，9]，得Ａ=1２345６789设［U,Ｓ,V］=svｄ（A)，得Ｕ=-0.２1480.88720.4082－0．5２０60．２496-0．816５－0．８263-0．38790.４0８2S=1６．848１０0０1.068400０0.000０Ｖ=-０．47９7-0.776７-０.408２-0.5724-０.075７0．8１６5－0.６６５１0．６２5３-0．40８2设ｔoｌ=S(1,１)*3*eps,得tol=1.1２23e-０14设ａ=suｍ(diag（S)>tｏｌ）－raｎk（Ａ),求aa＝0设b＝S(１,1)－norm(A），求bb＝0设c=ｓqrｔ(ｓum（diag(S*Ｓ))）-ｎorｍ(A，fro），求cc=3.5５２7ｅ-015设d=Ｓ(1,１)/S(3,3）-cｏｎd(A）,求dd＝-8设ｅ＝S(1,1)＊Ｓ(2,２)*S(3，3)-det(A)<tol，求ee＝１设f＝V（:,1)*nuｌl(A）,求ff=0设ｇ=abｓ（Ａ*null(A)）<toｌ,求gｇ=111设h=Ｕ（:，１：2)＝=ortｈ(Ａ），求ｈh=111111结果分析：由上面的结果得出如下的题目答案:（1)sum(diag(S)>tol)－rａnk(A)的结果是0；（２)Ｓ(１,１）-noｒm（A)=０的结果是０；（3)ｓqrt(ｓum(diag(S*S)))-norm(A,'fｒo＇)的结果是3.5５27e-015；（4)S(1，1）/S（3,3）－cond(A）的结果是-８;(５)Ｓ（1,1）*Ｓ(2,2）*S(3，3)-ｄet(A)<tｏl的运营结果是1；（6)V(：，１)'*null(A)的结果是0;(７）abs（A*null(Ａ))<toｌ的结果是111;(8)U(:，１：2)＝＝ｏｒth（A)的运营结果是111１11求积分，。程序：ｔ=0：ｐi／1００:2＊pi;ｆ=abs(siｎ(cos(t)))；Ｆ=cumｓuｍ(f）＊ｐi/1０0;pｌｏｔ（t，F)；xlabel（'x轴＇）;yｌabｅl('y轴＇）；ｔｉtｌe('描绘y＝ａbs（sｉn（ｃｏｓ(t)))在(0，x)间积分曲线［ｘ∈(0,2π)]')gridｏｎ;结果：求方程的解。程序：clｅar;cｌc;ｄisp(＇计算方程组x＾2+y^2＝1x*y＝2的根ｘy')［x,y]＝solve('x^2+y^２=1','ｘ*ｙ=2')结果：计算方程组x^２+ｙ^2＝1x*y=2的根xyx=[-1/2*(1/2*5^(1/2）+1/2＊i＊３＾（1/2))^3+１/4*5^(1/２)+1/4*i*3＾（1／2)][-1/２＊(１/２*5^（１／2）-1/2*i*3^(１/2))^3+1/4*５^(1/２）-1/４*i*3^(1/2)][-1/２*(－1／2*5^(1/2)＋１/２＊ｉ＊3^（1／２))^3-1/4*5＾(1/２)＋1/4*i*3^(1/2）][-1/2*(-1／２＊5^(1/2)-1/2*i*３^（１/２）)＾3-１/4*5^(1／2)-1/4＊ｉ＊3^(1/2)]y=[１／2＊5^（1/２)+１／2＊i＊3^(1/２)]［１/2*５^(1/２)-1／2＊i*3^（1/２）][－1/2*5^(1／2)+１/2＊i*3^(1/2)][-1/2＊5^（１/2)-1/2*i＊3^(1/2)]在某激励条件下，二阶系统归一化响应可表达为，其中为阻尼系数，,。请用不同的颜色或线型，在同一张图上,绘制取值下系统在区间内的响应曲线,并规定用和对它们相应的两条曲线进行醒目的文字标志。程序:b=0.2:0.2:２.0；t=0:0.１：18;color＝['r'，＇g.','bo','cx＇,'m+','y＊'，'kv',＇r-－','ｇ:','b-．'];forｉ=1：10p（i）=sqrt(1-b(i)．*b(i));q（i)=atａｎ（sqrt(1-ｂ（i).*b(i）).／b(i));ｙ=1-1./p(i).＊ｅxp(－b(i）.＊t).＊sｉn(p(ｉ)．*t+q(ｉ）)；ifi＝＝1plｏt(ｔ，y，colｏr(ｉ))text(t(30),y(3０）,'\ｌeftarrｏwξ＝０.2','FontSize'，18)ｈoldonelｓｅiｆi==１0ｐｌoｔ（t,y,color(i))texｔ（t（50)，ｙ(50）,'\ｌefｔarrowξ=2.0＇,'ＦontSize',18)hｏldｏｎｅlseplｏｔ(t,ｙ,color(i))ｈｏldｏneｎdendend结果:构建一个简朴的全波整流模型,并用示波器分两路同时观测原信号和整流后的信号波形。规定：信源模块被重新命名为“输入正弦波”；信宿模块被重新命名为“示波器”；连接到信宿的信号线上分别标注“原信号”和“整流信号”；在模型窗中添加注释文本。(1)、原理图：各模块参量设立:［示波器]ParａｍetｅｒsNuｍbeｒｏfaxes：2其他默认[其余模块]默认选项结果(2)、原理图：各模块参数设定［Gain]Ｇain:-1其他默认[其余模块］默认选项结果运用SIMＵLINK及其标准模块设计一个低通滤波器，从受噪声干扰的多频率混合信号中获取１0Hｚ的信号。在此,而各模块参量设立：［SineWave]Freqｕencｙ(ｒad／sec）:１0其他默认[SｉｎｅWavｅ1]Fｒｅquｅｎcy(rａd/ｓec)：１０0Phaｓe(rａｄ）:pi／2其他默认[RandｏmNumbｅｒ]Vａｒianｃｅ：0.2Samplｅtimｅ:10０其他默认[Sｕｍ］Ｌiｓtofsigns:|＋++其他默认[ＡnalｏgＦilterＤeｓign］Passbanｄｅdgefrｅquency(rad／s)：１0其他默认[其余模块］默认选项结果：已知系统的状态方程为：，其中,请构建该系统的仿真模型,并用XYＧｒａpｈ模块观测相轨迹。原理图:参数配置；[Sum]Liｓtｏfsigns:|--+其他默认 [integｒator]Initialｃoｎditioｎ:0．2ﻫ[inteｇrator1] ﻩInｉｔiａｌｃondiｔｉon:0.2结果:1０.ＣD７４HＣ00仿真（数字电路)实现方法:这里就是四个与非门的实现,运用LogｉcalOperaｔor模块做出四个与非门，还要再和Vcｃ和Gnd相与,再输出。然后用Ｓubsysteｍ封装成子系统即可;原理图:（1)、封装前的原理图(2)、分装后的测试原理图结果：分析：前两路分别为产生“1１00”和“0１01”的脉冲发生器的波形,后四路波形则分别为四个与非门的输出；可以看到，四路输出信号与输入信号均符合与非逻辑的关系,符合7400的真值表,设计对的!１1．2FSK信号的生成与接受滤波规定:产生二进制［0,１]随机信号波,信号波特率为1０00Ｂ对此信号进行２FSＫ调制,f1＝1500Ｈz,f２=3０50Ｈz加入高斯白噪声在接受端进行带通滤波抽样速率为20ＫHz.观测滤波前后信号的波形和频谱。实现方法：运用BernｏulliBinaｒyGenerator模块产生随机信号；运用SinｅＷave模块产生f1=1500Ｈz，f2=305０Hz的载波，结合Sｗitcｈ模块进行2FSＫ调制;运用AWGNChaｎnｅｌ模块构建高斯白噪声信道，其参数设立为:信噪比:１0dB,输入信号功率:1W;运用DiｇｉtａlFiｌtｅｒＤesign设计带通滤波器，分别为1500Hz和3050Hz滤波然后运用相干解调，在将两个载波相加,用符号函数判决出原信号。原理图:(ｚｏushｕＳimulink１1．mdl、zｏushｕSimulink1１b．mdl)参数设计：ﻩ［SineWave]: ﻩFrequency:1５00Hzﻩﻩ其他：默认ﻩ[SineＷａｖe1]: ﻩＦｒｅquｅncy:３050Hｚ 其他:默认ﻩ[ＢｅｒnouｌliBinaryGenerａtｏｒ]ﻩ Sampletiｍe:1/1000ﻩ 其他:默认ﻩ[Swｉtch]:ﻩﻩThreshｏiｄ:１/2ﻩ 其他:默认 ［ＡＷGＮChanｎeｌ]ﻩ SNＲ:10dｂ ﻩPower：1W 其他：默认 [ＤigitalＦilterDｅsign] Ｂａnｄpass；ﻩ Fs:20ｋHz ﻩFｓtop1：0.8KHzﻩ Fｐaｓs１:１.2KHz: Fｓtop２:2ＫＨzﻩ Ｆpass：２.2.5KＨz:ﻩ 其他:默认ﻩ[DｉgiｔalFiｌterDesign1］ﻩﻩBanｄpａss; ﻩＦs:2０ｋHzﻩ Ｆstoｐ1：2.5KＨz Ｆｐass２:2.8KHz: Fｓｔoｐ2:3.5KＨzﻩﻩFpasｓ2:4KHz: ﻩ其他：默认调制的原理图：调制加解调的原理图：结果:频谱图：调制、解调的波形图：结果分析:原信号通过２FＳＫ调制,通过噪声信道，再进行滤波、相干解调，最后又恢复出原信号，结果符合题目规定，系统设计对的!滤波器等相关设备的参数也设立恰当!实验成功！创建一个简朴的离散多速率系统:单位阶跃信号通过具有不同速率的采样后分别用作两个离散传递函数的输入。这两个离散传递函数有相同的有理分式，但采样时间和时间偏置二元对分别设为[1,0．1]和［0.7，0]。规定：观测这两个离散传递函数的输出有什么不同;用不同的颜色标帜不同采样速率系统。（通信系统)原理图：结果二、模拟数字电路仿真实验数字逻辑电路基础熟悉常用逻辑单元的特性,学会运用基本逻辑单元(与、或、非、异或、Ｒ－S触发器、Ｄ触发器、J-K触发器等),修改参数、增减输入/输出端。组合逻辑电路仿真设计二／四线译码器实现方法:运用LogicOperaｔor模块实现非门和与非门;运用PｕｌseGｅｎeｒａｔoｒ模块产生010１和０011序列码，用来检查设计出的子系统。原理图：封装前原理图:结果及分析分析：当!EN端置0时，此时,随着第二、三两路信号输入的不同,作为输出端的Y0,Y1，Y2,Y3依次输出低电平，符合真值表规定,证明设计的子系统工作正常。分析：当!EN端置1时,子系统没有被使能。此时,作为输出端的Y0,Ｙ１,Y2,Y3没有随着输入A0，A1作出相应变化,而是一直输出高电平,证明设计的子系统工作正常。设计对的！设计四选一数据选择器能仿真测试，并设计成子系统元件实现方法:运用的模块同上;地址信号用PulsｅGeneｒaｔor产生；Ｄ0到D3也由PulseGeｎerａｔor模块产生，为了区分各路信号,将其信号周期依次设立为0.０6２5s,0．1２5s,0.25s和0．５秒。原理图结果及分析:前两路分别为输入信号，组合产生0０,01，10，11的地址信号；第四到第七这四路信号为D０到D1，它们的频率各不相同;最后一路为输出信号;可以看出,随着输入的地址信号的变化，输入一次输出该地址信号相应的信号。子系统设计成功!！时序逻辑电路仿真设计四位二进制计数器（带置位和清零)实现方法：ﻫ运用JK触发器和逻辑门构成之。运用比较系数的方法可得电路的驱动方程为检查电路将P,T,LD,ＣR置为1,然后观测示波器输出。原理图:封装前的原理图：封装后：结果:从图中可以看出，脉冲在上升沿有效，计数器按00０0~１1１1方式进行计数。当Ｐ、T、L(置数）、CR(清零）有一个为低电平时的输出为全零。P、T为０时，计数器不能工作；L为0置数，D3~D0全零，输出也全零。ＣＲ为零时清零，所以输出也全零。三、数字信号解决仿真实验(一)、运用Kaiser窗函数设计的FIR低通滤波器进行数字滤波运用Ｋａｉser窗函数,设计具有如下指标的ＦＩＲ低通滤波器:ｆs=2０KHz,ｆpass＝4kHz,ｆstop＝5KHz,Aｐass=０.1ｄB，Ａstoｐ=80dＢ。归一化低通滤波器的重要性能指标有：绝对指标:[0，wｐ]为滤波器的通带，δｐ为可允许的通带波动；[ｗs,π］为滤波器的阻带,δs为阻带波动。相对指标:Ap为用ｄB表达的通带波动;As为用dB表达的阻带波动相对指标与绝对指标的关系为：Ａp＝-2０ｌg（（1-δｐ)/(1+δｐ))As=-20lg（δｓ)设计基本思绪:一方面选择一个符合规定的抱负滤波器(这里是抱负低通滤波器)；由于抱负滤波器的冲击响应是非因果的且无限长,为了能用ＦIR滤波器实现，必须用适当的窗函数来截取,从而得到线性相位和因果的FＩＲ滤波器。而滤波器的特性与窗函数的长度（也就决定了FIR滤波器的长度)密切相关。 一个截止频率为wｃ的抱负低通滤波器的单位冲击响应为：ｈ(n）＝sｉn［wc(n-a)]/[π(n-a）］a为采样延迟,相应的频谱特性为:Ｈ(w）=exp(-jａw),当|w|<wc时。用一个窗函数w(ｎ)去截取后，得到:hd(n）=ｈ（n).w（n）窗函数满足在［０,M］内关于a对称，在其他区域为0。在频域,FIR滤波器的频率响应由抱负滤波器的频响和窗函数的频响卷积得到。对给定的阻带衰减，Kａiser窗提供了最大的主瓣宽度，从而提供最陡的过渡带。窗函数的表达式为：w(n)=I0｛βsｑｕａｒe(1-(1-2n/M)2)}其中I0为零阶Beｓsｅl函数。Ｋａiseｒ窗的优点在于它可以通过改变参数β和Ｍ来改变阻带衰减和过渡带宽。在MATLAＢ的ＳimLinｋ中,有相关的模块供调用。题目：１生成２KHｚ和8KHz的混合信号，使该信号通过上述滤波器，观测滤波前后的波形变化和频谱分布情况。(基本规定）实现方法：一方面用两个SｉneWａvｅ模块分别产生２kHz和8kHz的正弦信号,之后用add模块将两者混合。用ＤｉｇitalFiｌｔeｒDesiｇn设计低通滤波器,具体设立如图所示：原理图：、结果、图为滤波之后和之前的信号频谱，可以观测到,滤波之后,处在通带之外的8kＨz的正弦分量消失了。上图是示波器的波形,其中第一路为２kHz和８kＨｚ的混合信号，第二路为滤波之后的信号，可见低通滤波器符合规定！2在改信号中加入高斯白噪声，观测滤波前后的波形变化和频谱分布情况。（基本规定）实现方法:在上题基础上增长AWGN信道，加入高斯白噪声。原理图：、结果：、分析：图为滤波之后和之前的信号频谱，可以观测到,滤波之后，处在通带之外的８kHz的正弦分量消失了。分析:上图是加噪声后示波器的波形，其中第一路为2kHｚ和8kHz和噪声的混合信号,第二路为滤波之后的信号（尚有低频的噪声)!3、将上述滤波器改为中心频率为１0０KHz的带通滤波器，信号源为带外的90ＫＨｚ和带内的100ＫHz的混合正弦信号，重做1,２问。实现方法:运用ＤigitalFｉlｔｅｒDesign设计带通滤波器，具体设立为:Fs=1000kＨｚ,Fsｔop1＝９9kHz,Fｐaｓｓ1＝99.５kHｚ,Ｆpasｓ２=100.5kHz,Fstｏp2=10１kHz。原理图:、结果:无噪声的频谱分析：无噪声的波形分析:加入噪声的频谱分析:加入噪声的波形分析：（二)、DＳB-SC信号的生成与解调规定:用离散(ＤSP）的方法生成DＳB信号载波频率为15０ＫHz,音频为500Hz和202３Ｈz的混合音。加入高斯白噪声带通滤波，分别用正交解调和检波的方法实现解调。抽样速率为600Hz。观测各点的信号波形和频谱。Ａ.正交解调法实现方法本题目规定用离散的方法生成DＳB信号，因此我们用编写M文献方式而不是用Siｍulink来实现。生成正弦音频信号时，对时间t用向量表达，可以实现离散化。DSB调制与解调分别使用amoｄ,ademｏd函数。加入高斯白噪声用awgn函数。带通滤波用ｃｈeｂ1oｒd和ｆftfilt函数实现。画频谱用ffｔ函数实现。程序代码（zｏｕsｈuＣode３00２.m）：ｔ＝０：1/(600*1000-1)：0．0１;m１=sin(500*2*pｉ*t);％5０0Ｈz正弦信号m2=sin(２*pi*2023*ｔ);％２0２3Hz正弦信号m=m1＋ｍ２;%混合信号figｕre(1)；ｓubplｏｔ(3,1,1);ｐlot(t,m）;titlｅ（＇原始信号')；z=amｏd(m，1．5*1０^5,6*10^5);%调制ｓubplot(3,1,２）,ｐlｏt（t，z）;ｔitle('DSB-SC信号图');y=awgn(ｚ,５);%加入高斯白噪声sｕｂｐlｏt（3，1,３),pｌot(t,y）;tiｔle（'加入噪声后的DSB-SC图，ＳＮＲ=５');ｆｉgure(2);f=(0：lengｔh（fft(z）)-1）*（600*100０-1)/lengtｈ（ｆft(ｚ)）;subplot(3,１，１);ploｔ（f,abs(ｆft(z)));ｔitle('调制后的频谱图＇）;％带通滤波器参数As=100;Ap=１;fs=６＊1０^5;%采样频率fｐ1=400；fｐ2=1.6*10＾5；%通频带fｓ1=３00;ｆｓ2＝1.7＊10^5；%截止带wp=[fp1ｆｐ２]/（fs/２）；ws=[fs１fs2]/(fｓ/2);[n，wn]=cｈeb1orｄ(wp,ｗs,Ap，Aｓ);ｂ=fir1(n,wn,ｃｈebwin(n＋1,Aｓ))；ｘ＝fｆｔfｉlt(b，y);z１=aｄｅｍｏｄ(x,１.5＊10^５,6*１0＾５)；%滤波法解调ｓubｐlot(３，1,2);ｐlot（t,z１)；title（'解调后的波形')f=（0:lengtｈ(fft(z１))-1)*（６00*10００-１）/lｅngth(fft(z1)）;suｂｐloｔ（３,1,3）；pｌot（f，ａbs(ffｔ（ｚ1))）;axiｓ([0６０0０03500]）;title('解调后的频谱图')；结果及分析:加入高斯白噪声后对源信号干扰比较大，解调出来的波形失真较严重,但解调后的频谱图是对的的。B.检波法实现方法:此方法与正交解调方法类似,只但是在解调时加入了直流信号,以便取出包络,取出包络后再减去直流信号得到源信号。取包络用取绝对值函数ａbｓ实现。程序:(zoｕｓhuCodｅ３00202.ｍ）Ａ=3;%包络解调加入直流信号t=０：1/(600*１000－1)：0.01；m1=sin(500*2*ｐi*t);％500Hz正弦信号m２＝ｓin(2＊pｉ*202３*t);%2023Hz正弦信号ｍ=m1＋m２+A；%包络解调加入直流信号fiｇｕrｅ(1);sｕbpｌot（3，1，1）;%混合信号plot(ｔ,m）;titｌe('和直流信号叠加在一起的原始信号');ｚ＝amod（m，1．5*10^5，６*10^5,'ａmｄsb-ｓc');%调制subpｌot(3,１，2),ｐlot（t,z);tｉtｌe('DＳB-ＳＣ信号图');ｙ=ａwｇｎ（ｚ,5)；%加入高斯白噪声ｓubpｌot(３,１,３）,plot(t,ｙ);tiｔle('加入噪声后的DSＢ-SＣ图，SＮR=5')；fiｇuｒe(2）;f=(０:lｅngｔh(fft(ｚ))-1)*（6０0＊1００0－1）/lengtｈ(fft(z)）;ｓｕｂplot(３,1,１）；plot(f,abs(ｆfｔ(z）));titｌe('调制后的频谱图')；%带通滤波器参数Aｓ=１0０;Ａp=１;ｆs＝6＊1０^5;%采样频率fp１=400；fp2=1.6*10^５;%通频带fｓ1=3００；fs2＝１.7*１0^5;%截止带ｗp=［fp1fp2]／（ｆs/2);ws＝[ｆs1ｆｓ2]/(fｓ/２)；［n,wn]=cheb1oｒd（wp,ws,Aｐ,As);b=fｉr1(ｎ,wn,cｈｅbwｉn(n+１,As))；%ｆiｒ滤波器ｚ1=ｐi/2*ｆftfiｌt(b,abs(y));%找出包络并进行缩放ｓuｂｐlｏｔ(3,1,２)；plot(ｔ,z1);titｌe('包络图')z1=z1-A;%去除直流信号f＝（0:lengｔh(fft(z1))－1)*（600