Matlab与C语言程序的应用编程接口

matlab实现：常见的离散时间信号

1． 单位抽样序列，或称为离散时间冲激，单位冲激： ? ??=01)(n δ 00≠=n n 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： ???=-01)(k n δ 0≠=n k n 2．单位阶跃序列 ? ??01)(n u 00<≥n n 在MATLAB 中可以利用ones( )函数实现。 );,1(N ones x = 3．正弦序列 )(cos )(0φω+=n A n x 这里， ,,0ωA 和φ都是实数，它们分别称为本正弦信号)(n x 的振幅，角频率和初始相位。 πω200=f 为频率。 4．复正弦序列 n j e n x ω=)( 5．实指数序列 n A n x α=)( 6. 随机序列 长度为N 的随机序列 基本数学函数参考教材P69页以及随后的使用说明。 注意使用行向量，特别是冒号运算符。 举例，长度为N 的实指数序列在MATLAB 中实现： n a x N n .^1 :0=-= 1. 单位采样 长度为N 的单位采样序列u(n)可以通过下面的MATLAB 命令获得：

u=[1 )1,1(-N zeros ]； 延迟M 个采样点的长度为N 的单位采样序列ud(n)(M

VB与MATLAB接口的实现

VB与MATLAB接口的实现: https://www.sodocs.net/doc/7b8371189.html, 实例说明 在本实例中，我们制作一个能够与Matlab进行交互的应用程序。程序运行结果如图78-1所示。 图78-1 运行结果 技术要点 z引用Matlab库 z执行Matlab命令 z结束Matlab 实现过程 ■ 新建项目 打开Visual https://www.sodocs.net/doc/7b8371189.html,，选择“新建项目”，在项目类型窗口中选择“Visual Basic项目”，在模板窗口中选择“Windows应用程序”，在名称域中输入“CnMatlab”，然后选择保存路径。单击“确认”。 ■ 添加控件 向当前窗体添加五个Button按钮，两个Picture控件，一个Label控件，一个Hscroll控件和一个Vscroll控件。单击菜单“项目|添加引用”，选中“Matlab Automation（Version5.3）Type Library”这一项。 注意：本程序只能在安装有Matlab的机器上运行。 ■ 设置属性 将Label控件和Command按钮的Text属性设置为与界面一致。在此不再赘述。 ■ 添加代码 Dim str1 As String ' 显示正弦图 Private Sub Command1_Click(ByVal eventSender As System.Object, ByVal eventArgs As System.EventArgs) Handles Command1.Click Dim matlab As Object matlab = CreateObject("matlab.application") matlab.MinimizeCommandWindow()

MATLABAPI详解

MATLAB、API详解 【例12.1.4-1】有一个绘圆的M脚本文件circle.m如下。希望获得一个MEX绘圆程序。（1）原始的绘圆脚本文件 [circle.m] clf;r=2;t=0:pi/100:2*pi;x=r*exp(i*t); plot(x,'r*');axis('square') （2）对这脚本文件直接编译将因错误而失败 mcc -x circle ??? Error: File "circle" is a Script M-file and cannot be compiled with the current Compiler. Error in ==> H:\MATLAB53\toolbox\compiler\mcc.dll （3）把脚本文件改写成函数文件。 [circle_f.m]: function circle_f(r) clf;t=0:pi/100:2*pi;x=r*exp(i*t); plot(x,'r*');axis('square') （4）再对circle_f.m进行编译，将顺利通过。 mcc -x circle_f %mcc是编译指令，详见12.4节。 （5）运行生成的MEX文件circle_f.dll circle_f(0.5) %调用circle_f绘制一半径为0.5的圆 which circle_f %查询所调用的circle_f的路径全称。

图 12.2.1-1 【Select MATLAB Componets】对话窗的选项局部图 图 12.2.2.1-1 为产生MEX文件所产生的配置屏1

图 12.2.2.1-2 为产生MEX文件所产生的配置屏2 12.1.1.1配置正确性的验证 （1）mex应用程序的验证 cd d:\mywork %把用户目录指定为当前目录 mex my_yprime.c %由my_yprime.c文件生成my_yprime.dll文件my_yprime(1,1:4) %运行my_yprime.dll文件 which my_yprime %获得my_yprime.dll文件的位置信息 ans = 2.0000 8.9685 4.0000 -1.0947 d:\mywork\my_yprime.dll （2）在MATLAB命令窗中验证mcc应用程序 mcc -x my_yprime_m%<1> my_yprime_m(1,1:4) which my_yprime_m ans = 2.0000 8.9685 4.0000 -1.0947 d:\mywork\my_yprime_m.dll （3）在 DOS提示符后验证mex、mcc应用程序

常用信号的MATLAB表示

5 常用信号的MATLAB表示5.1单位冲激函数、单位冲激序列 示例7： t = -5:0.01:5; y = (t==0); subplot(121); plot(t, y, 'r'); n = -5:5; x = (n==0); subplot(122); stem(n, x); 图5 运行结果如图5所示。

程序说明： (1)由n = -5:5得到一个1×11数组n；而在x = (n==0)中，n==0是一个向量运算，即向量n中的每一个元素与0比较是否相等，其比较结果0或1放在x中。这样得到的向量x也是1×11数组，且正好就是单位冲激序列。 (2)在MATLAB中，任何向量x的下标是从1开始的，不能取零或负值，而x(n)中的时间变量n则不此受限制。因此向量x的下标与时间变量n是两个概念，如本例中向量x(n)的下标是从1到11，而时间变量n是从-5到5。所以必须用一个与向量x等长的定位时间变量n，以及向量x，才能完整地表示序列x(n)。在信号的表示和运算中，这一点请务必注意；只有当序列x(n)的时间变量正好是从1开始时，才能省去时间变量n，因为此时向量的下标与时间变量相同。 (3)单位冲激函数的实现方法实际上与单位冲激序列是完全相同的，都是用序列表示。只不过表示连续时间信号的序列中两相邻元素所对应的时间间隔更小，如本例中t的间隔为0.01，而表示离散时间信号的序列中两相邻元素所对应的时间间隔一般为1。 由于单位冲激序列在信号与系统中经常使用，我们专门编制一个函数文件delta.m，在后面的实验部分直接调用该函数即可产生需要的波形。 % delta.m function [x, n] = delta(n1,n2,k) % 产生冲激序列δ(n-k),其中n1<=n<=n 2, n1<=k<=n2

实验1 常见离散信号的MATLAB产生和图形显示

实验1 常见离散信号的MATLAB 产生和图形显示 一、实验目的：加深对常用离散信号的理解 二、实验原理： 1.单位抽样序列：???=01)(n δ 00 ≠=n n 在MATLAB 中可以利用zeros()函数实现。 ; 1)1();,1(==x N zeros x 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n ?δ即：? ??=?01)(k n δ 0≠=n k n 2．单位阶越序列：???0 1 )(n u 00<≥n n 在MATLAB 中可以利用ones()函数实现。(1,)x ones N = 3．正弦序列：)/2sin()(?π+=Fs fn A n x 在MATLAB 中 ) /***2sin(*1 :0fai Fs n f pi A x N n +=?= 4．复正弦序列：n j e n x ?=)( 在MATLAB 中 ) **exp(1 :0n w j x N n =?= 5．指数序列：n a n x =)( 在MATLAB 中 n a x N n .^1:0=?= 三、实验内容： 1、编制程序产生上述5种信号（长度可输入确定），并绘出其图形。 2、讨论正弦序列、复指数序列的性质： （1）绘出信号()zn x n e =，当1126z j π=?+、1126z j π=+、112 z =、62πj z +=、6π j z =时 信号的实部和虚部图；当6 π j z =时信号的周期为多少？

（2）绘出信号() 1.5sin(2*0.1)x n n π=的频率是多少？周期是多少？产生一个数字频率为0.9的正弦序列，并显示该信号，说明其周期。 3、使用帮助功能学习square(方波),sawtooth(锯齿波)和sinc 函数，并绘图。 四、实验要求： 1、预先阅读MATLAB 基础； 2、讨论复指数序列的性质。

MATLAB与FPGA的接口

FPGA器件的开发平台与MATLAB接口仿真 2007-09-03 16:24 FPGA器件的开发平台与MATLAB接口仿真 中南财经政法大学信息学院周巍武汉大学电气工程学院张志杰引言 现场可编程逻辑门阵列FPGA器件的出现是超大规模集成电路技术和计算机辅助设计技术发展的结果。FPGA器件集成度高、体积小，具有通过用户编程实现专门应用功能。 它允许电路设计者利用基于计算机的开发平台，经过设计输入、仿真、测试和校验，直到达到预期的结果。目前使用最多的Quartus II 软件支持几乎所有的EDA工具，并且可以通过命令行或Tcl脚本与第三方EDA工具之间进行无缝连接。但在很多工程设计应用中，由FPGA器件完成的主程序中只完成大量的数学运算，程序调试时以二进制输出的信号可视性差，给设计人员进行仿真、调试带来了很多不便。对于很多工程设计人员来说MATLAB是一种熟悉的具有强大的运算功能和波形仿真、分析功能的软件，如果能将FPGA与MA TLAB接口，就可以快速、准确、直观地对FPGA程序进行校验和仿真，尤其在波形信号处理等工程应用领域具有实际意义。 Quartus II 开发软件 Altera公司的QuartusII软件提供了可编程片上系统（SOPC）设计的一个综合开发环境。Quartus II 开发工具人机界面友好、易于使用、性能优良，并自带编译、仿真功能。QuartusII 软件支持VHDL和Verilog硬件描述语言的设计输入、基于图形的设计输入方式以及集成系统级设计工具。QuartusII软件可以将设计、综合、布局和布线以及系统的验证全部都整合到一个无缝的环境之中，其中也包括和第三方EDA工具的接口。QuartusII设计软件根据设计者需要提供了一个完整的多平台开发环境，它包含整个FPGA和CPLD设计阶段的解决方案。图1说明了QuartusII软件的开发流程。 在实际应用设计中，对程序原理性及可执行性的验证主要集中在程序修改阶段，尤其在

常见连续信号的MATLAB表示

实验名称：常见连续信号的MATLAB 表示 报告人： 姓名班级学号 一、实验目的 1、熟悉常见连续时间信号的意义、特性及波形； 2、学会使用MATLAB 表示连续时间信号的方法； 3、学会使用MATLAB 绘制连续时间信号的波形。 二、实验内容及运行结果 1、运行以上5个例题的程序，保存运行结果。 2、已知信号()t f 的波形如下图所示，试用MATLAB 绘出满足下列要求的信号波形。 <1）()t f -； <2）()2-t f ； <3）()at f <其中a 的值分别为 21= a 和2=a ）； <4）? ?? ??+12 1t f 。 第一题 例题1

程序如下： >> t1=-10:0.5:10。 >> f1=sin(t1>./t1。 >> figure(1> >> plot(t1,f1> >> xlabel('取样间隔p=0.5'>。 >> title('f(t>=Sa(t>=sin(t>/t'>。>> t2=-10:0.1:10。 >> f2=sin(t2>./t2。 >> figure(2> >> plot(t2,f2> >> xlabel('取样间隔p=0.1'>。 >> title('f(t>=Sa(t>=sin(t>/t'>。运行结果如下：

f(t)=Sa(t)=sin(t)/t 取样间隔p＝0.5 f(t)=Sa(t)=sin(t)/t 取样间隔p=0.1例题2 程序如下： >> syms t >> f=sin(t>/t。 >> ezplot(f,[-10,10]>运行结果如下：

sin(t)/t t 例题3： 程序如下： >> t=-1:0.01:4。 >> t0=0。 >> ut=stepfun(t,t0>。>> plot(t,ut> >> axis([-1,4,-0.5,1.5]>运行结果如下:

基于matlab的直接序列扩频通信系统仿真

基于MATLAB的直接序列扩频通信系统仿真 1.实验原理：直接序列扩频（DSSS）是直接利用具有高码率的扩频码系列采用各种调 制方式在发端与扩展信号的频谱，而在收端，用相同的扩频码序去进行解扩，把扩展宽的扩频信号还原成原始的信息。它是一种数字调制方法，具体说，就是将信源与一定的PN码（伪噪声码）进行摸二加。例如说在发射端将"1"用11000100110，而将"0"用00110010110去代替，这个过程就实现了扩频，而在接收机处只要把收到的序列是11000100110就恢复成"1"是00110010110就恢复成"0"，这就是解扩。这样信源速率就被提高了11倍，同时也使处理增益达到10DB以上，从而有效地提高了整机倍噪比。 1.1 直扩系统模型 直接序列扩频系统是将要发送的信息用伪随机码(PN码)扩展到一个很宽的频带上去，在接收端用与发送端相同的伪随机码对接收到的扩频信号进行相关处理，恢复出发送的信号。对干扰信号而言，与伪随机码不相关，在接收端被扩展，使落入信号通频带的干扰信号功率大大降低，从而提高了相关的输出信噪比，达到了抗干扰的目的。直扩系统一般采用频率调制或相位调制的方式来进行数据调制，在码分多址通信中，其调制多采用BPSK、DPSK、QPSK、MPSK等方式，本实验中采取BPSK方式。 直扩系统的组成如图1所示，与信源输出的信号a(t)是码元持续时间为Ta的信息流，伪随机码产生器产生伪随机码c(t)，每个伪随机码的码元宽度为Tc (Tc<

labview与matlab接口的方法

LabVIEW与Matlab接口的方法 The Method of Interfacing Between LabVIEW and Matlab 陈金平 (新疆大学,乌鲁木齐　830008) 0　引言 虚拟仪器技术是计算机技术、现代测控技术和电子仪器技术相互结合、渗透的产物。在虚拟仪器系统中,数据的分析处理、控制、结果输出和用户界面等功能都由软件完成,硬件仅仅是为了解决信号的输入输出,因此,软件是整个仪器系统的核心,从某种意义上可以说:“软件即仪器”。虚拟仪器系统的软件设计可以采用通用的可视化编程语言,如Visual C++、Visual Basic、Delphi等,但更为方便高效的还是专用的虚拟仪器软件开发平台,如美国国家仪器公司(National Instruments,NI)的Lab2 VIEW、LabW indows/C VI,惠普公司的VEE等,而其中首推NI公司的图形化编程语言LabVIEW。 1　LabVIEW的功能及特点 LabVIEW是NI公司推出的一种虚拟仪器软件开发平台,自1986年正式推出,经过短短不到15年的时间,已经发展到以最新板本LabVIEW611为核心,包括控制与仿真、高级数字信号处理、统计过程控制、模糊控制和PID控制等众多附加软件包,运行于W indows NT/98、Linux、M acintosh、Sun和HP-UX等多种平台的工业标准软件开发环境。 LabVIEW在包括航空航天、通信、汽车、半导体和生物医学等众多领域内得到了广泛的应用。其最大的特色是采用编译型图形化编程语言———G语言(G raph2 Pro gramm ing),即用户设计好程序的大体框架后,如同画流程图一般,只需将系统提供的各种图形化功能模块连接起来,就可得到所需的应用软件。LabVIEW中的程序称为VI(virtual instruments),每个VI都由前面板和框图程序以及图标/连接端口三部分组成。 除了具备其它编程语言所提供的常规函数功能外,LabVIEW内部还集成了大量的生成图形界面的模板,如各种表头、旋钮、开关、LE D指示灯、图表等;丰富实用的数值分析、信号处理功能,如FFT变换、各种滤波器、信号发生器等;以及对RS-232、G PI B、VXI、数据采集板卡、网络等多种硬件的设备驱动功能,并免费提供数十家世界知名仪器厂商的几百种源码级仪器驱动,大大方便和简化了用户的设计开发工作。Lab2 VIEW使得过去繁琐、枯燥的软件开发变得简单、方便,尤其适合不熟悉传统文本编程语言(如C、BASIC等)的工程技术人员,被誉为工程师和科学家的语言。 但是,在大型的系统测试和仿真过程中,需要软件进行一些很复杂的数值计算时,LabVIEW的图形化编程语言就显得力不从心,M atlab是一种常用的高效率数学运算工具,它建立在向量、数组和复数矩阵的基础上,使用方便,将它和LabVIEW有机地结合起来会大大减少编程的工作量,提高编程效率。本文通过求解一常微分方程初值问题的例子,介绍了两种编程语言的接口方法。 2　在LabVIEW中调用Matlab语言的方法在测试系统设计和软件开发过程中,数学分析与信号处理是两个不可缺少的重要内容。LabVIEW将数据采集和测试分析中常用的数学和信号分析算法程序集成在一起,提供了先进的数学和信号分析环境,所有的数学分析节点都集中在M athematics子模板中。在此模板中有一M atlab Script节点,利用此节点就可以实现在LabVIEW中对M atlab语言的调用。下面通过具体例子介绍调用方法。举例如下: 用Runge2K utta法计算下列微分方程的解(初值问题): y′=-50y+50x2+2x　,　0≤x≤1 y(0)=1 2.1　编制M文件 启动M atlab610,利用其M文件编辑器编写M文件如下: function y=ff2(x,y) y=-50?y+50?x?x+2?x; 35 LabVIEW与M atlab接口的方法　陈金平

实验一--常见离散信号的MATLAB产生和图形显示

实验一 常见离散信号的MATLAB 产生和图形显示 授课课时：2学时 一、实验目的： （1）熟悉MATLAB 应用环境，常用窗口的功能和使用方法。 （2）掌握MATLAB 在时域内产生常用离散时间信号的方法。 （3）掌握离散信号的基本运算。 （4）掌握简单的绘图命令。 二、实验原理： （一）信号的表示和产生 ① 单位抽样序列 ? ??=01)(n δ 00 ≠=n n 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： ???=-0 1)(k n δ ≠=n k n 参考程序： 例1-1：)2010(()(<<-=n n n x ）δ clear all n1=-10;n2=20;n0=0;%在起点为n1,终点为n2的范围内，于n0处产生冲激。 n=n1:n2;%生成离散信号的时间序列 x=[n==n0];%生成离散信号x(n) stem(n,x);%绘制脉冲杆图

xlabel(' n');ylabel('x(n)');%横坐标和纵坐标的标注说明。 title('Unit Sample Sequence');%图形上方标注图名 axis([-10 20 0 1.2]);%确定横坐标和纵坐标的取值范围 ② 单位阶跃序列 ???=0 1 )(n u 00<≥n n 例1-2：)202((u )(<<-=n n n x ） clear all n1=-2;n2=20;n0=0; n=n1:n2;%生成离散信号的时间序列 x=[n>=n0];%生成离散信号x(n) stem(n,x,'filled'); xlabel('n');ylabel('x(n)'); title('Unit step Sequence'); axis([-2 20 0 1.2]); ③ 正弦序列 )sin()(?+=wn A n x 例1-3：一正弦信号的频率为1HZ ，振幅值幅度A 为1V ，在窗口显示2个周期的信号波形，并对该信号的一个周期进行32点采样获得离散信号并显示该连续信号和离散信号的波形。 参考程序： clear f=1; A=1;nt=2;

Matlab编程以及接口

Matlab编程的初步知识： 1.M脚本文件 Matlab命令集合； 无需输入和输出变量； 与其它脚本文件之间变量透明，即共同使用Matlab的工作空间。 比如 例1： x = (1:10); y =sin(x); plot(x,y) 例2： A = magic (4); B = A'; [C,D = xtimesAB(A,B); 2.函数 同样也是完成具体任务的指令集合，但是用一个名字封装起来，变量对外不透明，需要借助输入变量提供数据，输出变量给出结果。执行完毕后，所用的内存全部释放给Matlab。这样的命令集合体称为函数，封装的名字称为函数名，输入变量和输出变量在函数名前后指定。比如例2的xtimesAB(A,B)定义如下： function [C,D] =xtimesAB(A,B) % % This function can tell the products of e-e and V-C % C = A.*B; D = A*B; end 其中，函数名为xtimesAB，输入变量为A,B；输出变量为C,D。函数名下面%开头的部分为注释内容。在Matlab环境下，可以通过help xtimesAB来显示。 还有一种简单的函数，即inline函数，其特点是随用随定义。比如 1.>> myfun = ‘1+log(r) ‘; 2.>> myfuni=inline(myfun,’r’) 3.>>a=feval(myfuni,10) 4.结果a = 3.3026

我们用得最多的，就是像xtimesAB这样的M函数。与上面这个inline函数对应的M函数为 function y=myfun(r) y=1+log(r); 该函数结尾不含有end，即Matlab不要求必须有end。 使用时，在Matlab命令环境下，直接书写函数名，并给出输入变量，即可以工作。但需要注意的是，Matlab调用的是文件名而不是函数名。但一般两者取一样的名字。因此，关键是文件名。 M函数的流程控制简介如下： I. 循环 5.while a)while condition b)命令语句 c)end 6.for a)for i=1:n b)语句 c)end ii.条件 1.if a)if condition 1 b)语句 c)elseif condition2 d)语句 e)else f)语句 g)end 2.switch a)switch 表达式 b)case case1 c)语句 d)case case2 e)语句 f)otherwise g)语句 h)End iii.跳过：continue iv.跳出：break

第13章 MATLAB外部程序接口技术_习题答案

1 第13章 MATLAB 外部程序接口技术 习题13 一、选择题 1．要在Word 环境下调用MATLAB 的命令，需要调用Word 的（ ）模板。B A ．Normal B ．M-Book C ．自由格式 D ．空白文档 2．在Excel 环境下加载Spreadsheet Link 程序后，会在Excel 窗口的“开始”选项卡中增加一个（ ）命令组。D A ．Excel B ．Spreadsheet C ．Link D ．MATLAB 3．打开一个可读可写的文本文件，其打开方式为 。A A ．rt+ B ．r+ C ．rwt D ．a 4．以下选项中，用于定义指向MAT 文件指针的命令是（ ）。A A ．MA TFile *p; B ．MAT *p; C ．File *p; D ．FIL E *p; 5．关于MATLAB 引擎，下列说法中不正确的是（ ）。C A ．利用MA TLA B 引擎，可以在 C 程序中调用MA TLAB 的函数。 B ．通过MATLAB 引擎，可以提高开发应用程序的效率。 C ．通过MA TLAB 引擎，可以在MA TLAB 中直接调用用C 语言编写的函数 D ．包含MA TLAB 引擎函数的程序的执行效率减低。 二、填空题 1．在Word 与MATLAB 之间进行传递的内容称为 ，由M-Book 文档传向MATLAB 的命令称为 ，M-Book 文档中的MATLAB 命令的执行结果称为 。单元（Cell ），输入单元（Input Cell ），输出单元（Output Cell ） 2．Excel 和MA TLAB 的交互操作，通过 程序来实现。Spreadsheet Link 3．MA TLAB 文件操作的基本步骤是，先 文件，在对文件进行 ，最后 文件。打开，读写，关闭 4．对MAT 文件进行操作的C 程序中，一定要包含 头文件。mat.h 5．MEX 函数在头文件 中得到声明。mex.h 三、应用题 1．在MA TLAB 中创建一个100 × 200的随机矩阵，然后将数据导入到Excel 表格中，在Excel 中调用MATLAB 的差分函数按列进行差分运算。 2．已知)2ln(x y + =π，当x 取-3.0、-2.9、-2.8、…、2.8、2.9、3.0时，求各点的 函数值。要求： （1）将函数值输出到一个数据文件中。

用MATLAB实现常用的离散时间信号及其时域运算

用M A T L A B实现常用的离散时间信号及其时域 运算 WTD standardization office【WTD 5AB- WTDK 08- WTD 2C】

实验四用M A T L A B实现常用的离散时间信号及其时域运算 —— 摘要：在MATLAB中，只能用向量来表示离散时间信号。与连续信号不同，离散时间信号无法用符号运算来表示。用适当的MATLAB语句表示出信号后，就可以利用MATLAB的绘图命令stem来绘出直观的信号波形图，stem是专门用于绘制离散时间信号的。在MATLAB中离散序列的时域运算和变换不能用符号运算来实现，而必须用向量表示的方法，即在MATLAB中离散序列的相加、相乘需表示成两个向量的相加、相乘，因而参加运算的两序列向量必须有相同的维数。 一、实验目的：（1）学习MATLAB语言及其常用指令； （2）学习和掌握用MATLAB语言产生离散时间信号的编程方法； （3）通过编程绘制出离散时间信号的波形，加深理解信号的时域运算。二、实验内容：（1）运用MATLAB的绘图指令绘制离散时间信号； （2）用MATLAB语言实现离散时间信号的时域运算。 三、实验原理：（1）单位阶跃序列和单位样值序列。 离散时间信号只在某些离散的瞬时给出信号的值，因此，它是时间上不连续的序列。单位阶跃序列和单位样值序列在离散时间信号与系统的分析中是两个非常典型的序列，分别记为u(n)和δ(n)。它们的定义分别如下： 1 n≥0 1 n≥0 u(n)= δ(n)= 0 n<0 0 n≠0 若单位阶跃序列的起始点为n0，单位样值序列出现在n0时刻，则表达式分别为：

在MATLAB环境下访问外部函数的共享库文件

在MA TLAB环境下访问外部函数的共享库文件，必须首先把该库文件加载到内存中。一旦加载成功，就能直接在MA TLAB中直接请求关于函数的任何信息。而当不再需要该库时，就应当及时把库文件从内存中卸载以节省内存开销。 加载库 语法：loadlibrary(…shrlib?,?hfile?) 其中shrlib为加载的动态链接库文件名(filename.dll)，hfile为头文件名，它包含函数原型。例如，当加载包含MA TLAB中mx程序的libmx库时，可以使用下列语句。 hfile=[matlabroot?\extern\include\matrix.h?]; loadlibray(…libmx?, hfile) 卸载库 语法：unloadlibrary libmx 使用两个函数可以获取加载库的信息： libfunctions(…libname?) or libfunctions libname libfunctionsview(…libname?) or libfunctionsview libname 这两个函数的不同之处在于显示结果的方式不同，后者是以图形的方式显示在新的窗口中。而前者返回库libmx中有哪些可用的函数。请看示例： libfunctions libmx Methods for class lib.libmx: mxAddField mxGetFieldNumber mxIsLogicalScalarTrue mxArrayToString mxGetImagData mxIsNaN mxCalcSingleSubscript mxGetInf mxIsNumeric mxCalloc mxGetIr mxIsObject mxClearScalarDoubleFlag mxGetJc mxIsOpaque mxCreateCellArray mxGetLogicals mxIsScalarDoubleFlagSet 如果加上命令开头-full，则可以显示函数返回值的细节。 libfunctions libmx -full Methods for class lib.libmx: [mxClassID, MA TLAB array] mxGetClassID(MA TLAB array) [lib.pointer, MA TLAB array] mxGetData(MA TLAB array) [MA TLAB array, voidPtr] mxSetData(MA TLAB array, voidPtr) [lib.pointer, MA TLAB array] mxGetPr(MA TLAB array) [MA TLAB array, doublePtr] mxSetPr(MA TLAB array, doublePtr) uint8 mxIsFinite(double) uint8 mxIsInf(double) 值得注意的是，这两个函数返回值的类型均是MA TLAB的数据类型，虽然函数是利用C语言编写的。调用库函数 一旦库函数被加载到了内存空间，只要指定库名、函数名和变量就可以使用calllib函数调用库中的任何函数了。语法格式： calllib(…libname?,?funcname?,arg1,…,argn) 下列语句显示如何操作： hfile=['C:\MA TLAB7\extern\include\matrix.h'];

Matlab实现M序列的产生及其自相关序列学习资料

M a t l a b实现M序列的产生及其自相关序 列

电子信息工程专业课程设计任务书

1 需求分析 伪随机信号既有随机信号所具有的优良的相关性,又有随机信号所不具备的规律性. 因此,伪随机信号既易于从干扰信号中被识别和分离出来,又可以方便地产生和重复,其相关函数接近白噪声的相关函数, 有随机噪声的优点,又避免了随机噪声的缺点. 伪随机序列具有可确定性、可重复性,易于实现相关接受或匹配接受,故有很好的抗干扰性能. 因此伪随机序列在相关辩识、伪码测距、导航、遥控遥测、扩频通信、多址通信、分离多径、误码测试、线形系统测量、数据加扰、信号同步等方面均有广泛的应用. m序列是伪随机序列中最重要的一种，是最长线性移位寄存器序列，m序列易于实现，具有优良的自相关特性，在直扩通信系统中用于扩展要传递的信号。可以通过移位寄存器,利用MATLAB编程产生m序列。 2 概要设计 m 序列是最长线性反馈移位寄存器序列的简称，m 序列是由带线性反馈的移位寄存器产生的. 由n级串联的移位寄存器和和反馈逻辑线路可组成动态移位寄存器，如果反馈逻辑线路只由模2和构成，则称为线性反馈移位寄存器。 带线性反馈逻辑的移位寄存器设定初始状态后，在时钟触发下，每次移位后各级寄存器会发生变化。其中任何一级寄存器的输出，随着时钟节拍的推移都会产生一个序列，该序列称为移位寄存器序列。 n级线性移位寄存器的如图1所示： 图1 n级线性移位寄存器

图中i C 表示反馈线的两种可能连接方式，i C =1表示连线接通，第n-i 级输出加入反馈中；i C =0表示连接线断开，第n-i 级输出未参加反馈。 因此，一般形式的线性反馈逻辑表达式为 112201 (mod 2)n n n n n i n i i a C a C a C a C a ---==⊕⊕⊕=∑L 将等式左面的n a 移至右面，并将00(1)n n a C a C ==代入上式，则上式可改写为 10 0n i n i C a -==∑ 定义一个与上式相对应的多项式 ()n i i i F x C x ==∑ 其中x 的幂次表示元素的相应位置。式称为线性反馈移位寄存器的特征多项式，特征多项式与输出序列的周期有密切关系.当F(x)满足下列三个条件时，就一定能产生m 序列： (1) F(x)是不可约的，即不能再分解多项式； (2) F(x)可整除1p x +,这里21n p =-; (3) F(x)不能整除1q x +，这里q

MAtlab7.0.4接口完整详细解释并给出多种简易示例

Matlab7.0.4 接口连接说明。PSCAD使用的4.2.1教育版，网络版。 对于大小写的要求，对于系统的变量需要正常的大小写。用户自定义的不需要注意。 切记：每个有MA TLAB的工程都需要设置一下链接，那个勾一定要勾上啊，害了我几天。否则会报错的。如下： 不勾就会报错： 下面的两个设置是一次设置就可以保存的。 1.设置fortran编译器，不能使用GUN编译器，因为它只能用在fortran77中，我用的下面的版本。VF6.6专业版。

在PSCAD中如下设置:

2.设置MA TLAB，对于7.0.4已经可以正常运行。

3.PSCAD系统设置完成后，就可以建case了，开始加模块了。下面是我做的一个简单测试模块。单输入单输出。已成功运行。 testmat1.m :文件如下： function [I1]=testIO(in1)

I1=in1; 模块的制作如下： 首先建立一个模块，这个就自己建了。 上面的三个参数是可以自己改的， 第一个参数是m函数的文件名。 第二个参数是相对路径，即m函数所在的文件夹名，建议仿系统示例目录形式。

第三个参数是不一定需要的。前面两个参数主要是用于调用m函数，表示路径和文件名。 下面是编写的脚本文件，是在DSDYN中。 #STORAGE REAL:2 #LOCAL INTEGER I_CNT STORF(NSTORF) = $INPUT CALL MLAB_INT("%:Dir\$Path", "$Name", "R" , "R") $OUTPUT = STORF(NSTORF+1) NSTORF = NSTORF + 2 联合数组多输入多输出： 首先建立一个INPUT，一个OUTPUT，可以是多维的，下面是一个2维输入输出。注意，这里只有一个变量，使用数组表示，此数组下标是从1开始的。在m文件中也是使用数组。 在main中使用数据组合和分离器。

几种简单连续信号的matlab实现

连续信号及其MATLAB实现（原创） 连续信号又称为模拟信号，其信号存在于整个时间范围内，包括单位冲激信号，单位阶跃信号，斜坡信号，实指数信号，正弦信号，指数调制正弦信号，等等。以下我们将这几种常见的连续信号用MATLAB实现。 1、单位冲激信号 【例1】t=1/A=1/50时，单位脉冲δ(t)的MATLAB实现程序如下： clear all t1=-0.5:0.001:1; K=3; n1=length(t1); u=zeros(1,n1); t2=-0.5:0.001:0; n2=length(t2); u(n2)=K; plot(t1,u) 2、任意函数 【例2】脉冲序列： f(n)=0.4δ(n-2)+0.8δ(n-1)+1.2δ(n)+1.5δ(n+1)+1.0δ(n+2)+0.7δ(n+3)，其MATLAB实现程序如下： clear all t=-2:1:3; N=length(t); x=zeros(1,N); x(1)=0.4; x(2)=0.8; x(3)=1.2; x(4)=1.5; x(5)=1.0; x(6)=0.7; stem(t,x) axis([-2.2 3.2 0 1.7]) 3、单位阶跃函数

【例3】单位阶跃函数的MATLAB实现程序： clear all t=-0.5:0.001:1; t0=0; u=stepfun(t,t0); plot(t,u) axis([-0.6 1 -0.2 1.2]) 4、斜坡函数 【例4】斜坡函数：g(t)=3(t-1)，其MATLAB实现程序如下： clear all t=0:0.001:3; B=3; t0=1; u=stepfun(t,t0); n=length(t); for i=1:n u(i)=B*u(i)*(t(i)-t0); end plot(t,u) axis([-0.2 3.1 -0.2 6.2]) 5、实指数函数 【例5】实指数函数：f(t)=3exp(0.5t)，其实现程序： clear all t=0:0.001:3; B=3; A=3; a=0.5; u=A*exp(a*t); plot(t,u)

wind量化MATLAB接口功能

MATLAB接口主要功能是读取数据并且支持一些基本的函数运算。 1. windmatlab 创建windmatlab数据提取对象； 比如w=windmatlab，即可创建w对象句柄。 help windmatlab，即可查看windmatlab各类接口命令说明。 2. menu 用来显示wind数据向导工具条，也可以指定只显示特定的向导窗体。 比如w.menu(‘wsd’),即只显示wsd向导窗体。 查看menu说明，可以输入help windmatlab/menu 3. isconnected 用来确定创建的wind data object对象是否成功登陆。 比如x=isconnected(w),即判断wind对象w是否已经登陆成功。 4. cancelRequest 可以根据创建的数据请求id，取消数据请求。 比如你订阅了实时数据，返回reqid，那么想取消订阅，则可以通过如下命令行： w.cancelRequest(reqid) 5. wsd 返回选定证券品种的历史序列数据,包括日间的行情数据，基本面数据以及技术数据指标。 [data,codes,fields,times,errorid,reqid] = w.wsd(windCodes,windFields,startTime,endTime,v arargin) 其中，Data返回序列数据结果集。 Codes返回提取数据的WindCode代码。 Fields返回提取的指标名。 Times返回时间序列。 errorID返回错误ID号。 reqID返回请求ID号。

6. wsi 返回指定品种的日内分钟K线数据，包含历史和当天，分钟周期可以指定，技术指标参数可以自定义设置； [data,codes,fields,times,errorid,reqid] = w.wsi(windCodes,windFields,startTime,endTime,va rargin) 7. wst 返回日内盘口买卖十档快照数据和分时成交数据； [data,codes,fields,times,errorid,reqid] = w.wst(windCodes,windFields,startTime,endTime,va rargin) 8. wss 返回指定品种的历史截面数据，比如取沪深300只股票的2012年3季度的净利润财务指标数据；[data,codes,fields,times,errorid,reqid] = w.wss(windCodes,windFields,varargin) 9. wsq 可以返回当天实时指标数据，数据可以一次性请求，也可以通过订阅的方式获取当日实时行情数据； ?一次性请求实时行情数据： [data,codes,fields,times,errorid,reqid] = w.wsq(windCodes,windFields) ?订阅实时行情数据： [data,codes,fields,times,errorid,reqid] = w.wsq(windCodes,windFields,callback) 其中callback为回调函数，用来指定实时指标触发时执行相应的回调函数； 10. wset 返回数据集，包括板块成分、指数成分、ETF申赎成分信息、分级基金明细、融资标的、融券标的、融资融券担保品、回购担保品、停牌股票、复牌股票、分红送转。参数设置为起止日期、板块名称等。[data,codes,fields,times,errorid,reqid]=w.wset(windCodes windFields, varargin) 注意：wset、wpf、tdays、tdaysoffset、tday函数第二个与第四个输出变量没有实际意义，仅仅是为了保持与wsd、wss、wst、wsq函数输出格式的统一性。