基于Matlab的脉冲编码调制（PCM）系统设计与仿真

1、PAGE 武汉理工大学学科基础课群综合设计PAGE 课程设计任务书学生姓名： 专业班级： 指导教师： 工作单位： 题 目: 脉冲编码调制（PCM）的实现 初始条件：1、MATLAB软件；2、脉冲编码调制相关知识。要求完成的主要任务: 1、任务实现脉冲编码调制（PCM）技术的三个过程：采样、量化与编码。2、要求用仿真软件对其进行验证，使其满足以下要求：（1）模拟信号的最高频率限制在4KHZ以内；（2）分别实现64级电平的均匀量化和压缩率的非均匀量化；（3）按照13折线A律特性编成8位码。时间安排：第1，2天：分析题目，方案设计；第3，4，5天：软件设计；第6，7天：系统仿真；第8天：答辩，完成设

2、计说明书。指导教师签名： 年 月 日系主任（或责任教师）签名： 年 月 日目 录 TOC o 1-3 h z u HYPERLINK l _Toc251342766 摘 要 PAGEREF _Toc251342766 h I HYPERLINK l _Toc251342767 Abstract PAGEREF _Toc251342767 h II HYPERLINK l _Toc251342768 1 绪论 PAGEREF _Toc251342768 h 1 HYPERLINK l _Toc251342769 2 MATLAB简介 PAGEREF _Toc251342769 h 2 HYPER

3、LINK l _Toc251342770 2.1 MATLAB软件简介 PAGEREF _Toc251342770 h 2 HYPERLINK l _Toc251342771 2.2 MATLAB程序设计方法 PAGEREF _Toc251342771 h 2 HYPERLINK l _Toc251342773 3 PCM脉冲编码原理4 HYPERLINK l _Toc251342774 3.1 模拟信号的抽样及频谱分析4 HYPERLINK l _Toc251342775 3.1.1 信号的采样4 HYPERLINK l _Toc251342776 3.1.2 抽样定理 PAGEREF _T

4、oc251342776 h 4 HYPERLINK l _Toc251342777 3.1.3 采样信号的频谱分析5 HYPERLINK l _Toc251342779 3.2 量化 PAGEREF _Toc251342779 h 5 HYPERLINK l _Toc251342780 3.2.1 量化的定义 PAGEREF _Toc251342780 h 5 HYPERLINK l _Toc251342781 3.2.2 量化的分类 PAGEREF _Toc251342781 h 6 HYPERLINK l _Toc251342783 3.2.3 MATLAB的A律13折线量化 PAGERE

5、F _Toc251342783 h 12 HYPERLINK l _Toc251342784 3.3 PCM编码 PAGEREF _Toc251342784 h 12 HYPERLINK l _Toc251342785 3.3.1 编码的定义 PAGEREF _Toc251342785 h 12 HYPERLINK l _Toc251342786 3.3.2 码型的选择 PAGEREF _Toc251342786 h 13 HYPERLINK l _Toc251342787 3.3.3 PCM脉冲编码的原理 PAGEREF _Toc251342787 h 13 HYPERLINK l _Toc

6、251342794 4 PCM的MATLAB实现15 HYPERLINK l _Toc251342794 4.1 PCM抽样的MATLAB实现15 HYPERLINK l _Toc251342794 4.2 PCM量化的MATLAB实现18 HYPERLINK l _Toc251342794 4.2.1 PCM均匀量化的MATLAB实现18 HYPERLINK l _Toc251342794 4.2.2 PCM A律非均匀量化的MATLAB实现20 HYPERLINK l _Toc251342794 4.3 PCM A律13折线编码的MATLAB实现22 HYPERLINK l _Toc251

7、342788 5 结果分析及总结25 HYPERLINK l _Toc251342790 参考文献26PAGE II摘 要本设计结合PCM的抽样、量化、编码原理，利用MATLAB软件编程和绘图功能，完成了对脉冲编码调制（PCM）系统的建模与仿真分析。课题中主要分为三部分对脉冲编码调制（PCM）系统原理进行建模与仿真分析，分别为采样、量化和编码原理的建模仿真。同时仿真分析了采样与欠采样的波形、均匀量化与A律13折线非均匀量化的量化性能及其差异。通过对脉冲编码调制（PCM）系统原理的仿真分析，设计者对PCM原理及性能有了更深刻的认识，并进一步掌握MATLAB软件的使用。 关 键 词：脉冲编码调制（

8、PCM） 均匀与非均匀量化MATLAB仿真AbstractIn this design, combination the Simulink emulatation function and the S- functions spread function of MATLAB software, have completed the systematic emulatation and modeling for pulse code modulation( PCM). In this design，divide into 3 parts mainly, emulate to build moul

9、d and emulate analysis for the principle of pulse code modulation( PCM) systematic. They are modeling and emulatation of sampling, quantizing and ecoding. At the same time, emulate to analyse the waveform of sampling and owe sampling , the quantizing error of uniform quantizing and nonuniform quanti

10、zing. Through this design,the designer has a more profound understanding of PCM principles and performance , and further master the use of MATLAB software.Keywords: Pulse coding modulation ( PCM)uniform and non-uniform quantitative MATLAB simulationPAGE 421 绪论数字通信作为一种新型的通信手段，早在20世纪30年代就已经提出。在1937年，英

11、国人里费（A.H.Reeves）提出了脉冲编码调制（PCM）方式。从此揭开了近代数字传输的序幕。PCM系统的优点是：抗干扰性强；失真小；传输特性稳定，远距离再生中继时噪声不累积，而且可以采用有效编码、纠错编码和保密编码来提高通信系统的有效性、可靠性和保密性。另外，由于PCM可以把各种消息（声音、图像、数据等等）都变换成数字信号进行传输，因此可以实现传输和交换一体化的综合通信方式，而且还可以实现数据传输与数据处理一体化的综合信息处理。故它能较好地适应信息化社会对通信的要求。PCM的缺点是传输带宽宽、系统较复杂。但是，随着数字技术的飞跃发展这些缺点也不重要。因此，PCM是一种极有发展前途的通信方式

12、。2 MATLAB简介2.1 MATLAB软件简介MATLAB和 HYPERLINK /view/30608.htm t _blank Mathematica、 HYPERLINK /view/127864.htm t _blank Maple并称为三大数学软件。它在数学类科技应用软件中在 HYPERLINK /view/920695.htm t _blank 数值计算方面首屈一指。MATLAB可以进行 HYPERLINK /view/10337.htm t _blank 矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、

13、 HYPERLINK /view/14662.htm t _blank 图像处理、 HYPERLINK /view/1345304.htm t _blank 信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对 HYPERLINK /view/10075.htm t _blank C， HYPERLINK /view/36402.htm

14、t _blank FORTRAN， HYPERLINK /view/824.htm t _blank C+ ， HYPERLINK /view/29.htm t _blank JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的

15、问题。其具有以下特点：友好的工作平台和编程环境；简单易用的程序语言；强大的科学计算机数据处理能力；出色的图形处理功能；应用广泛的模块集合工具箱；实用的程序接口和发布平台；应用软件开发（包括用户界面）。2.2 MATLAB程序设计方法MATLAB有两种工作方式：一种是交互式的命令行工作方式；另一种是M文件的程序工作方式。在前一种工作方式下，MATLAB被当做一种高级数学演算纸和图形表现器来使用，MATLAB提供了一套完整的而易于使用的编程语言，为用户提供了二次开发的工具，下面主要介绍MATLAB控制语句和程序设计的基本方法。用MATLAB语言编写的程序，称为M文件。M文件有两类：命令文件和函数文

16、件。两者区别在于：命令文件没有输入参数，也不返回输出参数；而函数文件可以输入参数，也可以返回输出参数。命令文件对MATLAB工作空间的变量进行操作，而且函数文件中定义的变量为局部变量，当函数文件执行完毕时，这些变量被清除。M文件可以使用任何编辑程序建立和编辑，而一般常用的是使用MATLAB提供的M文件窗口。首先从MATLAB命令窗口的File菜单中选择New菜单项，在选择M-file命令，将得到的M文件窗口。在M文件窗口输入M文件的内容，输入完毕后，选择此窗口File菜单的save as命令，将会得到save as 对话框。在对话框的File 框中输入文件名，再选择OK按钮即完成新的M文件的建

17、立。然后在从MATLAB 命令窗口的File 菜单中选择Open对话框，则屏幕出现Open对话框，在Open对话框中的File Name 框中输入文件名，或从右边的directories框中打开这个M文件。在M文件所在的目录，再从File Name 下面的列表框中选中这个文件，然后按OK按钮即打开这个M文件。在M文件窗口可以对打开的M文件进行编辑修改。在编辑完成后，选择File菜单中的Save命令可以把这个编辑过的M文件报存下来。当用户要运行的命令较多或需要反复运行多条命令时，直接从键盘逐渐输入命令显得比较麻烦，而命令文件则可以较好地解决这一问题。我们可以将需要运行的命令编辑到一个命令文件中，

18、然后再MATLAB命令窗口输入该命令文件的名字，就会顺序执行命令文件中的命令。3 PCM脉冲编码原理3.1 模拟信号的抽样及频谱分析3.1.1 信号的采样离散时间信号通常是有连续时间信号经周期采样得到的。完成采样功能的器件称为采样器，下图所示为采样器的示意图。图中Xa(t)表示模拟信号，Xa(nt)表示采样信号，T为采样周期，n=0，1，2，。一般可以把采样器视为一个每隔T秒闭合一次的电子开关S。在理想情况下，开关闭合时间满足T。实际采样过程可视为脉冲调幅过程，Xa(t)为调制信号，被调脉冲载波p(t)是周期为T、脉宽为的周期脉冲串。当0时的理想采样情况是实际采样的一种科学的、本质的抽象，同时

19、可使数学推导得到简化。下面主要讨论理想采样。图3.1 采样器示意图及波形图3.1.2 抽样定理抽样也称取样、采样，是把时间连续的模拟信号变换为时间离散信号的过程。抽样定理是指：一个频带限制在（0，fH）内的时间连续信号m(t)，如果以T1/2fH秒的间隔对它进行等间隔抽样，则m(t)将被所得到的抽样值完全确定。这意味着，若m(t)的频谱在某一角频率H上为零，则m(t)中的全部信息完全包含在其间隔不大于1/2fH秒的均匀抽样序列里。换句话说，在信号最高频率分量的每一个周期内起码应抽样两次。根据抽样脉冲的特性，抽样分为理想抽样、自然抽样（亦称曲顶取样）、瞬时抽样（亦称平顶抽样）；根据被抽样信号的性

20、质，抽样又分为低通抽样和带通抽样。虽然抽样种类很多，但是间隔一定时间，抽样连续信号的样值，把信号从时间上离散，这是各种抽样共同的作用，抽样是模拟信号数字化及时分多路的理论基础。我们考察一个频带限制在(0,fH)赫的信号m(t)。假定将信号m(t)和周期性冲击函数(t)相乘，如图所示，乘积函数便是均匀间隔为T秒的冲激序列，这些冲激的强度等于相应瞬时上的m(t)值，它表示对函数m(t)的抽样。我们用ms(t)表示此已抽样的函数，即有ms(t)=m(t)(t)上述关系如下图所示。 图3.2 抽样示意图3.1.3 采样信号的频谱分析频谱分析自然要使用快速傅里叶变换FFT了，对应的命令即 fft ，简单

21、使用方法为：Y=fft(b,N)，其中b即是采样数据，N为fft数据采样个数。一般不指定N，即简化为Y=fft(b)。Y即为FFT变换后得到的结果，与b的元素数相等，为复数。以频率为横坐标，Y数组每个元素的幅值为纵坐标，画图即得数据b的幅频特性；以频率为横坐标，Y数组每个元素的角度为纵坐标，画图即得数据b的相频特性。对于现实中的情况，采样频率fs一般都是由采样仪器决定的，即fs为一个给定的常数；另一方面，为了获得一定精度的频谱，对频率分辨率F有一个人为的规定，一般要求F100秒；由采样时间ts和采样频率fs即可决定采样数据量，即采样总点数N=fs*ts。这就从理论上对采样时间ts和采样总点数N

22、提出了要求，以保证频谱分析的精准度。3.2 量化3.2.1 量化的定义模拟信号进行抽样以后，其抽样值还是随信号幅度连续变化的，即抽样值m(kT)可以取无穷多个可能值，如果用N个二进制数值信号来代表该样值的大小，以便利用数字传输系统来传输该样值的信息，那么N个二进制信号只能同M=2N个电平样值相对应，而不能同无穷多个电平值相对应。这样一来，抽样值必须被划分成M个离散电平，此电平被称作量化电平。或者说，采用量化抽样值的方法才能够利用数字传输系统来实现抽样值信息的传输。利用预先规定的有限个电平来表示模拟抽样值的过程称为量化。抽样是把一个时间连续信号变换成时间离散的信号，而量化则是将取值连续的抽样变换

23、成取值离散的抽样。通常，量化器的输入是随机模拟信号。可以用适当速率对此随机信号m(t)进行抽样，并按照预先规定，将抽样值m(kT)变换成M个电平q1，q2，qM之一,有mq(kTs)=qi，若mi-1m(kTs)V,h(i)=V;end if f(i)=-V,h(i)=-V;end flag=0; for j=2:L/2+1 if(flag=0) if(f(i)p(j) h(i)=p(j-1); flag=1; end; end; end; for j=L/2+2:L+1 if(flag=0) if(f(i)=0 if(x(i)=-t) y(i)=-(a*-x(i)/(1+log(a); el

24、se y(i)=-(1+log(a*-x(i)/(1+log(a); end endend仿真结果：图4.5 A律量化波形4.3 PCM A律13折线编码的MATLAB实现PCM均匀量化的MATLAB程序设计按如下步骤进行：（1）确定输入模拟信号；（2）根据给均匀量化的原理确定非均匀量化的算法程序；（3）将上述编码的十进制数转化成8位二进制数。PCM抽样的MATLAB实现源程序如下：function a_13code()t=0:0.000025:0.00025;y=sin(8000*pi*t)z=line13(y)c=pcmcode(z)function y=line13(x)x=x/max(

25、x);z=sign(x);x=abs(x);for i=1:length(x) if(x(i)=0)&(x(i)=1/64)&(x(i)=1/32)&(x(i)=1/16)&(x(i)=1/8)&(x(i)=1/4)&(x(i)=1/2)&(x(i) a_13codey = Columns 1 through 7 0 0.5878 0.9511 0.9511 0.5878 0.0000 -0.5878 Columns 8 through 11 -0.9511 -0.9511 -0.5878 -0.0000z = Columns 1 through 7 0 0.9045 1.0000 1.000

26、0 0.9045 0.0000 -0.9045 Columns 8 through 11 -1.0000 -1.0000 -0.9045 -0.0000c = 1 0 0 0 0 0 0 0 0 1 1 1 0 0 1 1 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 0 1 1 1 0 0 1 1 0 0 0 0 0 0 0 0 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 0 0 0 0 0 0 05结果分析及总结根据仿真的波形图和输出地量化、编码值可以得到以下结论：当抽样频率大于或等于

27、输入连续信号的频率2倍时，就可以无失真恢复原始信号；当不满足上述条件时就会出现频率混叠失真，不能恢复原始信号。均匀量化输出波形图清晰地显示处均匀量化的特征，每个量阶都是均匀分布的，每个间隔都是相等的。由于量化级数是64，所以从图中看到的结果不是那么明显，和输入波形相比几乎没什么变化。将A律非均匀量化的结果和A律13折线近似量化进行比较，两者压缩特性很接近。13折线输出的码组序列也符合要求。通过本次课程设计，我较系统地掌握有关PCM脉冲编码调制的设计思想和设计方法，主要对MATLAB的仿真方法，开发环境等有了一定的了解并对其进行测试和加以应用的知识得到学习。掌握了用程序对信号进行分析的基本方法，

28、并画出波形图。以前对PCM编码的方法只是在理论上，经过这次课设，加深了对PCM编码的基本原理理解，并对其在实际中的应用有了一定了解。通过这次课设我认识到在以后的学习中，不仅要有刻苦钻研的精神，还要有创新精神，对自己感兴趣的一定要用心去学。而在本次课程设计在刚开始时，由于对MATLAB的应用不太熟，觉得做起来有些棘手，当发现可以用软件仿真来实现PCM的编码过程后，我便去图书馆查阅了相关书籍，进一步熟悉了MATLAB编程方法和PCM的MATLAB实现原理，这是我最终顺利完成PCM系统设计的前提。随着设计的完成，我也逐渐掌握了PCM编码的工作原理及PCM系统的工作过程，通过应用软件仿真来实现各种通信

29、系统的设计，进一步地，可以完成硬件上的实现，以增强动手能力和学业技能。总体来说，这次课程设计使我受益匪浅。在摸索该如何设计PCM系统使之实现所需功能的过程中，不仅体验到了动手的乐趣，而且培养了我的设计思维，增加了实际操作能力。在让我体会到设计艰辛的同时，更让我体会到了成功的喜悦。今后，我要更加严格要求自己，主动寻找通信专业方面进行软件仿真及硬件设计的机会，通过实践逐步提升自己的专业技术水平，为今后发展奠定坚实的基础。参考文献1 樊昌信 曹丽娜，通信原理第六版，国防工业出版社，20072 周开利，邓春晖主编 MATLAB基础及其应用教程，北京大学出版社，20073 董振海,精通MATLAB 7

30、编程与数据库应用，电子工业出版社,20074 陈怀琛，MATLAB在电子信息课程中的应用（第二版），电子工业出版社，20065 张德丰主编，MATLAB通信工程仿真，机械工业出版社，2010本科生课程设计成绩评定表姓 名性 别男专业班级题 目： 脉冲编码调制（PCM）的实现答辩或质疑记录：1.你是怎样实现PCM系统的设计？有没有其他方法？答：本次设计我利用MATLAB编程实现了PCM，根据PCM采样、量化、编码原理，将各模块功能写成子函数形式，在主函数中对其进行调用。其他方法有：用MATLAB的simulink仿真，用电路仿真软件如Multisim进行仿真等。 2.问：非均匀量化与均匀量化相比

31、的优点是什么？ 答：首先，当输入量化器的信号具有非均匀分布的概率密度时，非均匀量化器的输出端可以得到较高的平均信号量化噪声功率比；其次，非均匀量化时，量化噪声功率的均方根基本上与信号抽样值成比例。因此量化噪声对大、小信号的影响大致相同，即改善了小信号时的信号量噪比。 3.问：对于PCM的采样你是怎么进行设计验证的？ 答：根据设计要求，模拟信号的频率要在4KHz以内，我选取200/2Hz。然后将正常抽样与有失真的抽样进行对比。根据抽样定理，当满足fs2fh时，才能不失真的还原出原始信号。因此，我选用200Hz作为正常抽样频率，再用100Hz作为失真的抽样频率。对比两者输出波形便可得到符合抽样定理

32、的结论。成绩评定依据：最终评定成绩（以优、良、中、及格、不及格评定）指导教师签字：_ 年 月 日附录资料：MATLAB的30个方法1 内部常数pi 圆周率 exp(1)自然对数的底数ei 或j 虚数单位Inf或 inf 无穷大 2 数学运算符a+b 加法a-b减法a*b矩阵乘法a.*b数组乘法a/b矩阵右除ab矩阵左除a./b数组右除a.b数组左除ab 矩阵乘方a.b数组乘方-a负号 共轭转置.一般转置3 关系运算符=等于大于=大于或等于=不等于4 常用内部数学函数 指数函数exp(x)以e为底数对数函数log(x)自然对数，即以e为底数的对数log10(x)常用对数，即以10为底数的对数lo

33、g2(x)以2为底数的x的对数开方函数sqrt(x)表示x的算术平方根绝对值函数abs(x)表示实数的绝对值以及复数的模三角函数（自变量的单位为弧度）sin(x)正弦函数cos(x)余弦函数tan(x)正切函数cot(x)余切函数sec(x)正割函数csc(x)余割函数反三角函数 asin(x)反正弦函数acos(x)反余弦函数atan(x)反正切函数acot(x)反余切函数asec(x)反正割函数acsc(x)反余割函数双曲函数 sinh(x)双曲正弦函数cosh(x)双曲余弦函数tanh(x)双曲正切函数coth(x)双曲余切函数sech(x)双曲正割函数csch(x)双曲余割函数反双曲函

34、数 asinh(x)反双曲正弦函数acosh(x)反双曲余弦函数atanh(x)反双曲正切函数acoth(x)反双曲余切函数asech(x)反双曲正割函数acsch(x)反双曲余割函数求角度函数atan2(y,x)以坐标原点为顶点，x轴正半轴为始边，从原点到点（x，y）的射线为终边的角，其单位为弧度，范围为（ ， 数论函数gcd(a,b)两个整数的最大公约数lcm(a,b)两个整数的最小公倍数排列组合函数factorial(n)阶乘函数，表示n的阶乘 复数函数 real(z)实部函数imag(z)虚部函数abs(z)求复数z的模angle(z)求复数z的辐角，其范围是（ ， conj(z)求复

35、数z的共轭复数求整函数与截尾函数ceil(x)表示大于或等于实数x的最小整数floor(x)表示小于或等于实数x的最大整数round(x)最接近x的整数最大、最小函数max(a，b，c，)求最大数min(a，b，c，)求最小数符号函数 sign(x)5 自定义函数-调用时：“返回值列=M文件名（参数列）”function 返回变量=函数名（输入变量） 注释说明语句段（此部分可有可无）函数体语句 6进行函数的复合运算compose(f,g) 返回值为f(g(y)compose(f,g,z) 返回值为f(g(z)compose(f,g,x,.z) 返回值为f(g(z)compose(f,g,x,y

36、,z) 返回值为f(g(z)7 因式分解syms 表达式中包含的变量 factor(表达式) 8 代数式展开syms 表达式中包含的变量 expand(表达式)9 合并同类项syms 表达式中包含的变量 collect(表达式，指定的变量)10 进行数学式化简syms 表达式中包含的变量 simplify(表达式)11 进行变量替换syms 表达式和代换式中包含的所有变量 subs(表达式，要替换的变量或式子，代换式)12 进行数学式的转换调用Maple中数学式的转换命令，调用格式如下：maple(Maple的数学式转换命令) 即：maple(convert(表达式，form)将表达式转换成f

37、orm的表示方式 maple(convert(表达式，form, x) 指定变量为x，将依赖于变量x的函数转换成form的表示方式（此指令仅对form为exp与sincos的转换式有用） 13 解方程solve(方程，变元) 注：方程的等号用普通的等号： = 14 解不等式调用maple中解不等式的命令即可，调用形式如下： maple(maple中解不等式的命令)具体说，包括以下五种：maple( solve（不等式）) maple( solve（不等式，变元） ) maple( solve（不等式，变元） ) maple( solve（不等式，变元） ) maple( solve（不等式，变

38、元） )15 解不等式组调用maple中解不等式组的命令即可，调用形式如下： maple(maple中解不等式组的命令) 即：maple( solve（不等式组，变元组） )16 画图方法：先产生横坐标的取值和相应的纵坐标的取值，然后执行命令： plot(x,y) 方法2：fplot(f(x)，xmin,xmax) fplot(f(x)，xmin,xmax,ymin,ymax) 方法3：ezplot(f(x) ezplot(f(x) ,xmin,xmax) ezplot(f(x) ,xmin,xmax,ymin,ymax) 17 求极限（1）极限：syms x limit(f(x), x, a) （2）单侧极限：左极限：syms x limit(f(x), x, a,left)右极限：syms x limit(f(x), x, a,right) 18 求导数diff(f(x) diff(f(x),x) 或者：syms x diff(f(x) syms x diff(f(x), x) 19 求高阶导数 diff(f(x)，n) diff(f(x),x，n)或者：sy