基于MATLAB信号与系统中抽样定理的仿真-(最终版)

1、PAGE - 分类号 编号烟 台 大 学 文 经 学 院毕 业 论 文（设 计）基于MATLAB信号与系统中抽样定理的仿真Signal and System Based on MATLAB simulation sampling theorem系 别：电子信息与计算机科学系 专 业： 通信技术 班 级： 姓 名： 学 号： 指导老师： （讲师） 2013年 6 月 1 日烟台大学文经学院- 基于MATLAB信号与系统中抽样定理的仿真姓 名： 导 师： 2013年 6 月 1 日烟台大学文经学院- 烟台大学文经学院毕业论文（设计）任务书系（部）：电子信息与计算机科学系姓名学号毕业届别专业通信技术

2、毕业论文（设计）题目基于MATLAB的信号与系统中抽样定理的仿真指导教师学历硕士研究生职称讲师所学专业物理电子学具体要求(主要内容、基本要求、主要参考资料等)：主要内容：基于MATLAB的信号与系统中抽样定理的仿真，利用MATLAB在数字信号处理中的基本应用，并会对结果用所学知识进行分析。 基本要求：掌握MATLAB的基本操作，掌握基于MATLAB的通信系统的设计与实现的基本工作原理，理解系统中各信号抽样仿真的原理。主要参考资料： 1 楼顺天.基于MATLAB的系统分析与设计信号处理M.西安:西安电子科技大学出版社2 邹理和.数字信号处理M.北京:国防工业出版社,1988.3941滞后，这就是

3、离散系统最常用零阶保持器的主要原因之一。进度安排：2013年3月5日前，确定选题及指导教师2013年3月5日至3月31日，进行毕业设计调研，完成大概设计2013年4月1日至4月20日，进行毕业设计，写论文2013年4月20日至4月25日，对内容和机构进行第一遍修改2013年5月1日前，进行第二遍修改2013年5月10日6月1日 设计作品验收、论文装订、毕业答辩指导教师（签字）： 年 月 日系部、院（系）意见： 系（部）主任或教学院长签字： 年 月 日备注：PAGE 摘要 与连续时间信号相比，离散时间信号的处理更加灵活、方便，应用更加广泛。在实际应用中，通常先将连续时间信号转换成离散时间信号，前

4、提是抽样过程必须满足抽样定理。抽样定理在连续时间信号与离散时间信号之间架起了一座桥梁。本文利用MATLAB对抽样过程进行仿真，并对抽样定理的条件与结论进行验证和讨论，通过研究可以得到符合实际应用的结论，同时加深对抽样定理的理解和把握。关键词 MATLAB; 抽样定理；仿真 Abstract: The processing for discretetime signal is more flexible，convenient，and more widely used，compared with continuoustime signa1In practice，continuoustime sig

5、nal is usually changed into discrete-time signal firstly，provided that the sampling process meets the sampling theoremSampling theorem sets up a bridge between continuous-time signal and discretetime signa1But must the restrictions and conclusion of the sampling theorem coincide with the actual appl

6、ication? In this paper，sampling process is simulated by using MATLAB ,and the conditions and conclusion of samplingtheorem is validated and discussion，through which shows that the restrictions and the conclusion of sampling theorem agrees with the practice；meanwhile，understanding and grasp of sampli

7、ng theorem is deepenedKey words：MATLAB；sample theory; simulation目 录 TOC o 1-3 h z u HYPERLINK l _Toc201119338 第1章 绪论 PAGEREF _Toc201119338 h 1 HYPERLINK l _Toc201119339 1.1 MATLAB语言的特点 PAGEREF _Toc201119339 h 1 HYPERLINK l _Toc201119340 1.2 MATLAB产品主要的应用领域 PAGEREF _Toc201119340 h 1 HYPERLINK l _Toc2

8、01119341 1.3 抽样定理简介1 HYPERLINK l _Toc201119342 第2章 抽样信号 PAGEREF _Toc201119342 h 4 HYPERLINK l _Toc201119343 2.1 抽样信号原理 PAGEREF _Toc201119343 h 4 HYPERLINK l _Toc201119344 2.2 模拟信号算法 PAGEREF _Toc201119344 h 5 HYPERLINK l _Toc201119345 2.2.1 模拟信号频率计算. PAGEREF _Toc201119345 h 5 HYPERLINK l _Toc20111934

9、6 2.2.2 采样信号频率计算. PAGEREF _Toc201119346 h 5 HYPERLINK l _Toc201119347 2.2.3 模拟信号实现 PAGEREF _Toc201119347 h 6 HYPERLINK l _Toc201119348 2.3 本章小结 PAGEREF _Toc201119348 h 8 HYPERLINK l _Toc201119363 第3章 MATLAB程序仿真 PAGEREF _Toc201119363 h 9 HYPERLINK l _Toc201119364 3.1仿真目的9 HYPERLINK l _Toc201119365 3.

10、2仿真原理说明及设计内容9 HYPERLINK l _Toc201119365 3.3设计内容10 HYPERLINK l _Toc201119365 3.4仿真设计实现：信号的产生和频域分析10 HYPERLINK l _Toc201119370 结论17 HYPERLINK l _Toc201119371 致谢18 HYPERLINK l _Toc201119372 参考文献 PAGEREF _Toc201119372 h 19 HYPERLINK l _Toc201119373 附录A PAGEREF _Toc201119373 h 20 HYPERLINK l _Toc20111937

11、4 附录B PAGEREF _Toc201119374 h 28PAGE V 烟台大学文经学院毕业论文（设计）PAGE 29 第1章 绪论MATLAB（Matrix Laboratory）是1984年美国Math Works公司产品，MATLAB的推出得到了各个领域专家学者的广泛关注，并越来越多的应用到我们的学习生活中来，是目前通信工程上最广泛应用的软件之一。Matrix Laboratory意为“矩阵实验室”，最初的MATLAB只是一个数学计算工具。但现在的MATLAB已经远不仅仅是一个“矩阵实验室”，它已经成为一个集概念设计、算法开发、建模仿真，实时实现于一体的集成环境，它拥有许多衍生子集

12、工具。MATLAB的基本数据单位是 HYPERLINK /view/10337.htm t _blank 矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完全相同的事情简捷得多.在新的版本中也加入了对 HYPERLINK /view/10075.htm t _blank C, HYPERLINK /view/36402.htm t _blank FORTRAN, HYPERLINK /view/824.htm t _blank c+ , HYPERLINK /view/29.htm t _blank JAVA的支持.可以直接调用,用

13、户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用。1.1MATLAB语言的特点（1）语言简洁紧凑, 语法限制不严，程序设计自由度大，可移植性好。（2）运算符、库函数丰富。（3）强大的数值（矩阵）运算功能。（4）界面友好、编程效率高。（5）图形功能强大。1.2MATLAB产品主要的应用领域（1）仿真和建模（2）实时仿真（3）自动控制（4）信号处理与通信（5）数据分析/科学计算（6）算法开发（7）图形和可视化法（8）独立应用开发（9）其他领域 1.3抽样定理简介定义所谓抽样。就是对时间连续的信号隔一定的时间间隔T抽取一个瞬时幅度值（样值），抽样是由抽样门完成的。 在一个频带限

14、制在（0，）内的时间连续信号，如果以的时间间隔对它进行抽样，那么根据这些抽样值就能完全恢复原信号。或者说，如果一个连续信号的频谱中最高频率不超过，当抽样频率时，抽样后的信号就包含原连续的全部信息。意义抽样定理指出，由样值序列无失真恢复原信号的条件是 ，为了满足抽样定理，要求模拟信号的频谱限制在0之内（为模拟信号的最高频率）。为此，在抽样之前，先设置一个前置 HYPERLINK /view/204522.htm t _blank 低通滤波器，将模拟信号的带宽限制在以下，如果前置低通滤波器特性不良或者抽样频率过低都会产生折叠噪声。例如，话音信号的最高频率限制在3400HZ，这时满足抽样定理的最低的

15、抽样频率应为=6800HZ，为了留有一定的防卫带，CCITT规定话音信号的抽样率=8000HZ，这样就留出了8000-6800=1200HZ作为滤波器的防卫带。应当指出，抽样频率不是越高越好，太高时，将会降低信道的利用率（因为随着升高，数据传输速率也增大，则数字信号的带宽变宽，导致信道利用率降低。）所以只要能满足,并有一定频带的防卫带即可。以上讨论的抽样定理实际上是对低通信号的情况而言的，设模拟信号的 HYPERLINK /view/476354.htm t _blank 频率范围为，带宽B=-.如果B，则称之为 HYPERLINK /view/1892207.htm t _blank 带通信

16、号，载波12路群信号（频率范围为60108KHZ）就属于带通型信号。对于低通型信号来讲，应满足的条件，而对于带通型信号，如果仍然按照这个抽样，虽然能满足样值频谱不产生重叠的要求，但是无疑太高了（因为带通信号的高），将降低信道频宽的利用率，这是不可取的。适用范围抽样定理在实际应用中应注意在抽样前后模拟信号进行滤波，把高于二分之一抽样频率的频率滤掉。这是抽样中必不可少的步骤。抽样定理抽样定理：设时间连续信号，其最高截止频率为 ，如果用时间间隔为T的开关信号对进行抽样时，则就可被样值信号唯一地表示。定理分类时域抽样定理频带为F的连续信号可用一系列离散的采样值,.来表示,只要这些采样点的时间间隔t1/

17、2F，便可根据各采样值完全恢复原来的信号。 这是时域采样定理的一种表述方式。时域采样定理的另一种表述方式是：当时间信号函数的最高频率分量为时,的值可由一系列采样间隔小于或等于的采样值来确定,即采样点的重复频率。图为模拟信号和采样样本的示意图。时域采样定理是采样误差理论、随机变量采样理论和多变量采样理论的基础。1频域抽样定理对于时间上受限制的连续信号（即当tT时,=0,这里T=是信号的持续时间），若其频谱为,则可在频域上用一系列离散的采样值 采样值来表示,只要这些采样点的频率间隔 频率间隔。用MATLAB软件来进行抽样定理的实现，可以很容易的实现，能够达到实验室所不能达到的理想化（器件各种参数）

18、，用软件实现还可以节省大量的人力及物力。第2章 抽样信号2.1 抽样信号原理设 是连续时间信号, 的傅立叶变换为： (2.1)设 为周期冲激脉冲信号, 为采样周期, (2.2)以表示采样输出, 则： (2.3)根据傅立叶变换性质,采样信号 的傅立叶变换为 (2.4)其中 P(j) 为 的傅立叶变换, 因此, (2.5)从 (2.5) 式我们又可得到 (2.6)采样信号的傅立叶变换的另一种表示形式为 (2.7)上式说明采样后的信号频谱是采样之前信号频谱以整数倍采样频率的平移然后叠加形成的。式(2.7)用于离散时间信号的傅立叶变换的计算。是连续的模拟频率,。 2.2 模拟信号算法2.2.1 模拟信

19、号频率计算.在计算机中模拟信号是无法存储的, 但当模拟信号 为有限 (-L t P ， L, P 是正数) 时,以充分小的时间取 的值，就可得到一根平滑的曲线，用于近似分析。本文中, 选用双指数函数， L= P 。此时, 可用一数组 表示, 数组的大小是从- L 的点起到L 点之间间隔为$ t 的点数，并用 m 表示各个点。当 用 表示时, （2.8）设 ,则 则为矩阵 Y 的转置。考虑到频率范围较大, 需将化分成一定间隔的频率点。设的频率范围为 (-,) ，取频率间隔为 ，K为常数，根据实际情况来选择，则正频率点为0，2, , 。负频率点与正频率点对称，共有 2K+1 个频率点。由此y()

20、中后一个数值表示频率点.根据上述原理，式(2.8)变为 ，其中=,，t=-L,-L+，P-,P,W=k()，k=-K,-K+1，K-1,K,为t矩阵的转置.2.2.2 采样信号频率计算.设采样周期为,则采样信号 =, 其中=L.根据,其中n=-N,-(N-1),0,(N-1),NW=k()k=-K,-K+1,K-1,KK为常数,为n矩阵的转置。这里以频率f为横轴,不同于文献2。. 2.2.3模拟信号实现1)模拟信号及其傅立叶变换图形2.取,用MATLAB编程画出的和的傅立叶变换图形如图1所示.程序中利用了0,并取L=P=5ms,时间间隔为0.05ms。（a）模拟信号； （b）连续傅立叶变换图2

21、.1 用MATLAB编程画出的及当时,。为了与后面的采样信号对比,取=2P(6000),K=2000，共有4001个频率点，,傅立叶变换，以f为横轴。2)=0.2ms时,的采样信号及其傅立叶变换图形.从图1可以认为,模拟信号的最高频率=2kHz,现=0.2ms,则=1=5000Hz2,采样频率大于模拟信号最高频率的2倍,满足抽样定理.其它条件L=P=5ms,=2P(6000),K=2000.采样信号及其傅立叶变换如图2.2所示。（a）采样信号（b）离散时间傅立叶变换图2.2 =0.2ms,采样信号和离散时间傅立叶变换3)=0.5ms时,的采样信号及其傅立叶变换图形.因为=0.5ms,则=1=2

22、000Hz2,采样频率小于模拟信号最高频率的2倍,不满足抽样定理.其它条件L=P=5ms,=2P(6000),K=2000.采样信号x2及其傅立叶变换X2如图2.3所示(a)采样信号 (b)离散时间傅立叶变换图2.3 =0.5ms,采样信号和离散时间傅立叶变换 2.3本章小结图2.12.3利用了相同的坐标。图2.2以5kHz的频率采样,图2.3以2kHz的频率采样,采样后的信号频谱是采样之前信号频谱以整数倍采样频率的平移、叠加形成;当采样频率大于模拟信号最高频率2倍时,采样信号的频谱不互相重叠,模拟信号能从采样信号中不失真地还原;当采样频率小于模拟信号最高频率2倍时,采样信号的频谱互相重叠,模

23、拟信号不能从采样信号中还原.这完全符合抽样定理3.图形是根据原理用MATLAB编程实现的，实现起来非常简单，而且结论非常可靠明确。第3章 MATLAB程序仿真3.1仿真目的1)熟悉抽样定理、信号的抽样过程；2)通过实验观察欠采样时信号频谱的混叠现象；3)掌握抽样前后信号的频谱的变化，加深对抽样定理的理解；4)掌握抽样频率的确定方法。3.2仿真原理说明及设计内容低通抽样定理：一个频带限制在（0，）赫内的时间连续信号，如果以秒的时间间隔对它进行等间隔（均匀抽样，则将被所得到的抽样值完全确定。此定理告诉我们：若的频谱在某一角频率上以上为零，则中的全部信息完全包含在其间隔不大于秒的均匀抽样序列里。抽样

24、速率（每秒钟的抽样点数）应不小于，否则，若抽样速率，则会产生失真，这种失真叫混叠失真。3.3设计内容产生一个连续的时间连续信号，并对其进行频谱分析，绘制时域波形图和频域波形图。对产生的连续信号进行抽样，并绘制抽样后的时域波形图，和频域波形图。（3）改变抽样频率，分别对原始连续信号抽样，绘制抽样后的时域和频域波形，最后对得到的波形进行分析。从而验证抽样定理。3.4仿真设计实现：信号的产生和频域分析 用MATLAB产生一个连续的信号，；根据抽样定理，在MATLAB中编写源程序代码，画出原信号时域波形和频域波形，再分别用不同的频率的抽样脉冲对其进行抽样，在MATLAB中实现不同频率抽样时，时域和频域

25、波形的效果对比，验证抽样定理。原始信号的时域波形和频域波形的源程序代码如下：t0=10;%定义时间长度ts=0.001; % 抽样周期fs=1/ts;df=0.5;% 频率的分辨率t=-t0/2:ts:t0/2;%定义时间序列x=sin(200*t);m=x./(200*t);w=t0/(2*ts)+1;m(w)=1;%定义在t=0时刻的值为1m=m.*m; m=50.*m;%定义函数sinc（200t）subplot(2,1,1); plot(t,m);xlabel(时间); title(原信号的时域波形)axis(-0.15,0.15,-1,50);M,mn,dfy=fftseq(m,ts

26、,df);%傅里叶变换，程序在后面M=M/fs;f=0:dfy:dfy*length(mn)-dfy-fs/2;%定义频率序列subplot(2,1,2); plot(f,abs(fftshift(M);xlabel(频率);axis(-500,500,0,1);title(原信号的频域波形); 原始信号在抽样频率为200HZ和400HZ的抽样脉冲下抽样波形(此时的满足抽样定理)，其时域和频域波形如下：%200hz抽样t0=10;%定义时间长度ts1=0.005; %充足抽样周期fs1=1/ts1;df=0.5;% 频率的分辨率t1=-t0/2:ts1:t0/2;%定义时间序列x1=sin(2

27、00*t1);m1=x1./(200*t1);w1=t0/(2*ts1)+1;m1(w1)=1;%定义在t=0时刻的值为1m1=m1.*m1; m1=50.*m1;%定义函数sinc（200t）subplot(2,1,1);stem(t1,m1);xlabel(时间); title(抽样频率2fhfs时时域波形)axis(-0.15,0.15,-1,50);M1,mn1,dfy1=fftseq(m1,ts1,df);%傅里叶变换，程序在后面M1=M1/fs1;N1=M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1;f1=-7*dfy1*length(mn1):dfy

28、1:6*dfy1*length(mn1)-dfy1-fs1/2;%定义频率序列subplot(2,1,2);plot(f1,abs(fftshift(N1)xlabel(频率);axis(-500,500,0,1);title(抽样频率fh1/2fs时信号的频域波形); （3）原始信号在抽样频率为100HZ的抽样脉冲下的抽样波形（此时不满足抽样定理），其时域和频域波形如下所示：%100hz抽样t0=10;%定义时间长度ts1=0.01; %欠抽样周期fs1=1/ts1;df=0.5;% 频率的分辨率t1=-t0/2:ts1:t0/2;%定义时间序列x1=sin(200*t1);m1=x1./(

29、200*t1);w1=t0/(2*ts1)+1;m1(w1)=1;%定义在t=0时刻的值为1m1=m1.*m1; m1=50.*m1;%定义函数sinc（200t）subplot(2,1,1);stem(t1,m1);xlabel(时间); title(抽样不足时时域波形)axis(-0.15,0.15,-1,50);M1,mn1,dfy1=fftseq(m1,ts1,df);%傅里叶变换，程序在后面M1=M1/fs1;N1=M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1;f1=-7*dfy1*length(mn1):dfy1:6*dfy1*length(mn1)

30、-dfy1-fs1/2;%定义频率序列subplot(2,1,2);plot(f1,abs(fftshift(N1)xlabel(频率);axis(-500,500,0,1);title(抽样不足（fh2fs）时信号的频域波形);（4）因为在MATLAB的库函数中没有傅里叶变换函数，而在分析抽样定理时需要观察频域的抽样波形，需要用到傅里叶变换，故编写了子函数fftseq（），实现其频域变换，供上述程序调用其代码如下：function M,m,df=fftseq(m,ts,df)% M,m,df=fftseq(m,ts,df)fs=1/ts;if nargin = 2n1=0;else n1=f

31、s/df;endn2=length(m);n=2(max(nextpow2(n1),nextpow2(n2);M=fft(m,n);m=m,zeros(1,n-n2);df=fs/n;end 结论通过MATLAB的编程实现对连续信号进行抽样的仿真实验，进一步加深了对我们队抽样定理的理解。首先通过以上仿真的波形图，我们可以很直观的看到，原始信号的频率约为64HZ，当抽样频率为原始信号频率的2倍及其以上（实验采用的是200HZ和400HZ）时，在抽样信号的频域图上，可以看出信号的频谱的相邻的周期内的波形不会发生混叠现象（即抽样只是实现了信号频谱在频域内周期的复制和搬移而已），通过一个合适的低通滤波

32、器，就可以得到和原频谱一样的的波形，从而能够无失真的重建原始信号。而当抽样频率小于原始信号频谱的2倍时，从仿真的频域波形可以看出相邻频谱之间发生混叠从而与原来信号的频谱波形不一样。这样的波形通过低通滤波器也不能恢复原来的波形。从而验证了抽样定理的正确性。 致谢大学三年学习时光已经接近尾声，在此我想对我的母校，我的父母、亲人们，我的老师和同学们表达我由衷的谢意。感谢我的家人对我大学三年学习的默默支持；感谢我的母校烟台大学文经学院给了我在大学三年深造的机会，让我能继续学习和提高；感谢烟台大学文经学院的老师和同学们三年来的关心和鼓励。老师们课堂上的激情洋溢，课堂下的谆谆教诲；同学们在学习中的认真热情

33、，生活上的热心主动，所有这些都让我的三年充满了感动。 这次毕业论文设计我得到老师和同学的帮助，其中我的论文指导老师刘美娟老师对我的关心和支持尤为重要。每次遇到难题，我最先做的就是向刘老师寻求帮助，而刘老师每次不管忙或闲，总会抽空来找我面谈，然后一起商量解决的办法。刘老师平日里工作繁多，但我做毕业设计的每个阶段，从选题到查阅资料，论文提纲的确定，中期论文的修改，后期论文格式调整等各个环节中都给予了我悉心的指导。这几个月以来，刘老师不仅在学业上给我以精心指导，同时还在思想给我以无微不至的关怀，在此谨向刘老师致以诚挚的谢意和崇高的敬意。同时，本篇毕业论文的写作也得到了同学的热情帮助。感谢在整个毕业设

34、计期间和我密切合作的同学，和曾经在各个方面给予过我帮助的伙伴们，在此，我再一次真诚地向帮助过我的老师和同学表示感谢！ 参考文献1 AlnVO,AlanSW.Signals&systemsM.NewJersey:Prentice-HallInt,Inc,1997.5145222 楼顺天.基于MATLAB的系统分析与设计信号处理M.西安:西安电子科技大学出版社,1999.1361473 邹理和.数字信号处理M.北京:国防工业出版社,1988.3941滞后，这就是离散系统最常用零阶保持器的主要原因之一。4 苏金明,阮沈勇.MATLAB6.1实用指南1M2.北京:电子工业出版社,2002.5 姚东.M

35、ATLAB命令大全1M2.北京:人民邮电出版社,2000.6 郑阿奇主编,曹弋,赵阳编著.MATB实用教程1M2.北京:电子工业出版社,2004.7现代通信原理与技术第二版，张辉 曹丽娜 主编。8 MATLAB 电子仿真与应用教程第二版，王华 李友军 刘建存 主编。9MATLAB 7.0/Simulink 6.0 应用实例仿真与高效算法开发 黄永安 李文成 高小科 主编。附录A主程序清单function w(action) if nargin1, action=initialize;end; if strcmp(action,initialize), oldFigNumber=watchon;

36、 figNumber=figure( . Name,抽样定理及其恢复演示实验 NumberTitle,off, . Visible,off, . DoubleBuffer,on, . BackingStore,off, . Colormap,); axes( . Units,normalized, . Position,0.05 0.40 0.75 0.90, . Visible,off, . NextPlot,add); text(0,0)奈奎斯特抽样定理 HorizontalAlignment,center,. FontWeight,bold,. FontSize,10,. Color,0

37、,0,1); axis(-1 1 -1 1); axes( . Units,normalized, . Position,0.05 0.30 0.75 0.90, . Visible,off, . NextPlot,add); text(0,0,要想抽样后能够不失真地还原出原信号 HorizontalAlignment,center,. FontWeight,bold,. FontSize,10,. Color,0,0,1); axis(-1 1 -1 1); axes( . Units,normalized, . Position,0.05 0.20 0.75 0.90, . Visible

38、,off, . NextPlot,add); text(0,0,抽样频率必须大于两倍信号谱的最高频率。,. HorizontalAlignment,center,. FontWeight,bold,. FontSize,10,. Color,0,0,1); axis(-1 1 -1 1); %= % 建立信息窗口 top=0.30; left=0.40; right=0.70; bottom=0.05; labelHt=0.03; spacing=0.005; % 首先建立信息窗的框架 frmBorder=0.02; frmPos=left-frmBorder bottom-frmBorder

39、 . (right-left)+2*frmBorder (top-bottom)+2*frmBorder; uicontrol( . Style,frame, . Units,normalized, . Position,frmPos, . BackgroundColor,0.50 0.50 0.50); % 然后建立文本标签 labelPos=left top-labelHt (right-left) labelHt; uicontrol( . Style,text, . Units,normalized, . Position,labelPos, . BackgroundColor,0.5

40、0 0.50 0.50, . ForegroundColor,1 1 1, . String,信息窗); % 然后建立可编文本区域 mcwPos=left bottom (right-left) top-bottom-labelHt-spacing; mcwHndl=uicontrol( . Style,edit, . HorizontalAlignment,left, . Units,normalized, . Max,10, . BackgroundColor,1 1 1, . Position,mcwPos, . Callback,w(quickeval); % 存储此句柄以便将来使用

41、set(gcf,UserData,mcwHndl); %= % 关于所有按钮的信息 labelColor=0.8 0.8 0.8; top=0.95; bottom=0.05; left=0.75; yInitLabelPos=0.90; left=0.75; labelWid=0.21; labelHt=0.05; btnWid=0.21; btnHt=0.05; % 同一命令下标签和按钮之间的空隙 btnOffset=0.003; %某一按钮和下一命令标签之间的空隙 spacing=0.05; %= % 控制台框架 frmBorder=0.02; yPos=0.05-frmBorder;

42、frmPos=left-frmBorder yPos btnWid+2*frmBorder 0.9+2*frmBorder; h=uicontrol( . Style,frame, . Units,normalized, . Position,frmPos, . BackgroundColor,0.50 0.50 0.50); %= % 模拟信号及采样频率命令按钮 btnNumber=1; yLabelPos=top-(btnNumber-1)*(btnHt+labelHt+spacing); labelStr=模拟信号及采样频率 ; labelList=|矩形波fs2000Hz|矩形波fs1

43、500Hz|矩形波fs800Hz|正弦波fs1800Hz|正弦波fs1200Hz|正弦波fs600Hz|锯齿波fs2500Hz|锯齿波fs2000Hz|锯齿波fs800Hz; cmdList=str2mat( . % 按下窗体右端的按钮可以看到抽样定理与恢复演示实验,. square1(2000),%模拟信号表达式为y(t)=2*u(t-0.05)-u(t)-u(t-0.0625),最高频率为750赫兹,fs2000赫兹即大于模拟信号最高频率的2倍时采样序列的频谱不混叠,恢复信号不失真,. square1(1500),%模拟信号表达式为y(t)=2*u(t-0.05)-u(t)-u(t-0.0

44、625),最高频率为750赫兹,fs1500赫兹即等于模拟信号最高频率的2倍时采样序列的频谱不混叠,恢复信号不失真,. square1(800),%模拟信号表达式为y(t)=2*u(t-0.05)-u(t)-u(t-0.0625),最高频率为750赫兹,fs800赫兹即小于模拟信号最高频率的2倍时采样序列的频谱混叠,恢复信号失真, . sine(1800),%模拟信号表达式为 y(t)=200*exp(-50*sqrt(2)*pi*t)*sin(-50*sqrt(2)*pi*t)最高频率为600赫兹,fs1800赫兹即大于模拟信号最高频率的2倍时采样序列的频谱不混叠,恢复信号不失真,. sin

45、e(1200),%模拟信号表达式为y(t)=200*exp(-50*sqrt(2)*pi*t)*sin(-50*sqrt(2)*pi*t)最高频率为600赫兹,fs1200赫兹即等于模拟信号最高频率的2倍时采样序列的频谱不混叠,恢复信号不失真,. sine(600),%模拟信号表达式为y(t)=200*exp(-50*sqrt(2)*pi*t)*sin(-50*sqrt(2)*pi*t)最高频率为600赫兹,fs600赫兹即小于模拟信号最高频率的2倍时采样序列的频谱混叠,恢复信号失真,. sawtooth1(2500),%模拟信号表达式为y(t)=20*(t-0.05)最高频率为1000赫兹,

46、fs2500赫兹大于模拟信号最高频率的2倍时采样序列的频谱不混叠,恢复信号不失真 sawtooth1(2000),%模拟信号表达式为y(t)=20*(t-0.05)最高频率为1000赫兹,fs2000赫兹即等于模拟信号最高频率的2倍时采样序列的频谱不混叠,恢复信号不失真 ,. sawtooth1(800),%模拟信号表达式为y(t)=20*(t-0.05)最高频率为1000赫兹,fs800赫兹即小于模拟信号最高频率的2倍时采样序列的频谱混叠,恢复信号失真); callbackStr=w eval; % 标签信息 labelPos=left yLabelPos-2*labelHt labelWi

47、d labelHt; uicontrol( . Style,text, . Units,normalized, . Position,labelPos, . BackgroundColor,labelColor, . HorizontalAlignment,left, . String,labelStr); % 选项按钮信息 btnPos=left yLabelPos-2*labelHt-btnHt-btnOffset btnWid btnHt; hndl1=uicontrol( . Style,popup, . Units,normalized, . Position,btnPos, . S

48、tring,labelList, . Callback,callbackStr, . UserData,cmdList); %= uicontrol( . Style,pushbutton, . Units,normalized, . Position,left bottom+2*btnHt+spacing btnWid 1.5*btnHt, . String,关于, . Callback,w(info); % 关闭按钮 uicontrol( . Style,pushbutton, . Units,normalized, . Position,left bottom btnWid 1.5*bt

49、nHt, . String,关闭, . Callback,close(gcf); % 显示图形 hndlList=mcwHndl hndl1 ; watchoff(oldFigNumber); set(figNumber,Visible,on, . UserData,hndlList); w(eval);elseif strcmp(action,info), helpwin(mfilename)elseif strcmp(action,quickeval), % 执行小命令窗内的命令 hndlList=get(gcf,UserData); mcwHndl=hndlList(1); evalmc

50、w(mcwHndl);elseif strcmp(action,eval) hndlList=get(gcf,UserData); n=length(hndlList); mcwHndl=hndlList(1); evalStr=str2mat(); for count=2:n, newStrList=get(hndlList(count),UserData); newStrVal=get(hndlList(count),Value); newStr=deblank(newStrList(newStrVal,:); if isempty(newStr), evalStr=str2mat(eva

51、lStr,newStr); end; end; set(mcwHndl,String,evalStr); evalmcw(mcwHndl);end; % if strcmp(action, .附录BMATLAB is Matrix Laboratories (Matrix Laboratory) means. In addition to a superior numerical ability, it also provides a professional level of symbols, word processing, visualization and modeling and s

52、imulation of real-time control functions.MATLAB is the basic unit of data matrix, and the mathematical expression of the directive, commonly used in the form of projects are very similar, they use MATLAB to solving problems than with C, FORTRAN languages such as End of the same things much more simp

53、le in the new Version also joined the C, FORTRAN, c + +, JAVA support can directly call, the user can also prepare themselves to the utility into MATLAB functions in their own convenience after the call, in addition many of the fans are prepared MATLAB A number of classic, users can directly downloa

54、d can be used, very convenient.MATLAB is based on the matrix, but because of his openness, and mathwork also absorbed as the advantages of software such as maple, MATLAB become a powerful mathematical software Popular MATLAB 6.5/7.0 including the hundreds of functions within the main pack 30 and sev

55、eral tool kits (Toolbox). Kits can be divided into functional tool kits and tool kits disciplines. Functional tool kits used to expand the MATLAB The symbols, visual modeling and simulation, word processing and real-time control functions. Disciplinary tool kit is relatively strong professional tool

56、 kits and tool kits control, signal processing tool kit, communications tool kits and so are such.Open so that users MATLAB widely welcomed. In addition to internal functions, all MATLAB main package of documents and various tool kits are readable may amend the document, the source of users to modif

57、y or add your own programming tools for the new structure Package.The advantages and characteristics of MATLAB(1) friendly work platform and programming environment MATLAB from a series of tools. These tools user-friendly MATLAB functions and paper, which is used in tool-graphical user interface. MA

58、TLAB, including desktop and command window, history command window, the editor and debugger, the path for users to browse and search assistance, work space, the paper browser. With the commercialization of MATLAB and the continuing escalation of the software itself, MATLAB user interface is more ref

59、ined, more close to the standard Windows interface, a more interactive, more simple operation. And a new version of the MATLAB provide a complete on-line enquiries, help system, greatly facilitate the use of the user. Simple programming environment provides a more complete debugging systems, procedu

60、res do not have to direct the compiler can run, but also to timely report to the mistakes and errors cause analysis. (2) easy-to-use programming language MATLAB an advanced matrix / array language, it contains control sentence, function, data structure, input and output and object-oriented programmi