通信系统基于（15,11）汉明编码的matlab仿真 clear m=4; n=2^m-1; k=11; msg1=randint(1,6000*k,[0,1])%产生信息序列 code0=vec2mat(msg1,k); code1=encode(msg1,n,k,'hamming/binary');%进行汉明纠错编码 code2=reshape(code1,90,1000)'; code3=zeros(1000,90);%设置零矩阵，以便储存交织后的序列 for i=1:1000 temp=code2(i,:); temp1=reshape(temp,15,6);%按15*6交织 code3(i,:)=reshape(temp1',1,90); end [row,cloums]=size(code3); code4=code3'; [row1,cloums1]=size(code4); code5=reshape(code4,6,15000)'; %将交织后的序列转换为6行，15000列矩阵msg2=zeros(15000,6); for j=2:6, msg2(:,j)=xor(code5(:,j),code5(:,(j-1))); end msg2(:,1)=code5(:,1);%此时得到的msg2为gray映射后的信号序列 code8=msg2';%转置，便于比特符号转换 A=[32,16,8,4,2,1]; %为比特符合转换所设的序列 msg3=A*code8;%生成符号序列 msg4=qammod(msg3,64);%将符号序列进行64QAM调制 [row2,cloums2]=size(msg4); dB=0:1:20; for k=1:length(dB), snr=10.^(dB(k)./10); %信噪比 sgma=sqrt(63./(6*snr));%标准差 b1=real(msg4)+sgma*randn(row2,cloums2);%分路叠加噪声 b2=imag(msg4)+sgma*randn(row2,cloums2);%分路叠加噪声 rx=complex(b1,b2); %————————量化判决——————————% for m=1:row2, for n=1:cloums2, if ((b1(m,n)<-6)) b11(m,n)=-7;
系统仿真 § 4.1控制系统的数学模型 1、传递函数模型（tranfer function） 2、零极点增益模型（zero-pole-gain） 3、状态空间模型（state-space） 4、动态结构图（Simulink结构图） 一、传递函数模型(transfer fcn-----tf) 1、传递函数模型的形式 传函定义：在零初始条件下，系统输出量的拉氏变换C（S）与输入量的拉氏变换R（S）之比。 C(S) b1S m+b2S m-1+…+b m G（S）=----------- =- -------------------------------- R(S) a1S n + a2S n-1 +…+ a n num(S) = ------------ den(S) 2、在MATLAB命令中的输入形式 在MATLAB环境中，可直接用分子分母多项式系数构成的两个向量num、den表示系统： num = [b1, b2, ..., b m]; den = [a1, a2, ..., a n]; 注：1）将系统的分子分母多项式的系数按降幂的方式以向量的形式输入两个变量，中间缺项的用0补齐，不能遗漏。 2）num、den是任意两个变量名，用户可以用其他任意的变量名来输入系数向量。 3）当系统种含有几个传函时，输入MATLAB命令状态下可用n1,d1;n2,d2…….。 4）给变量num，den赋值时用的是方括号；方括号内每个系数分隔开用空格或逗号；num,den方括号间用的是分号。 3、函数命令tf( ) 在MATLAB中，用函数命令tf( )来建立控制系统的传函模型，或者将零极点增益模型、状态空间模型转换为传函模型。 tf( )函数命令的调用格式为： 圆括号中的逗号不能用空格来代替 sys = tf ( num, den ) [G= tf ( num, den )]
VB与MATLAB接口的实现: https://www.doczj.com/doc/4f14292444.html, 实例说明 在本实例中，我们制作一个能够与Matlab进行交互的应用程序。程序运行结果如图78-1所示。 图78-1 运行结果 技术要点 z引用Matlab库 z执行Matlab命令 z结束Matlab 实现过程 ■ 新建项目 打开Visual https://www.doczj.com/doc/4f14292444.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()
学习中心/函授站_ 成都学习中心 姓名赵洪学号7020140122093 西安电子科技大学网络与继续教育学院 2015学年上学期 《MATLAB与系统仿真》期末考试试题 （综合大作业） 考试说明： 1、大作业于2015年4月3日公布，2015年5月9日前在线提交； 2、考试必须独立完成，如发现抄袭、雷同、拷贝均按零分计。 3、程序设计题（三（8，10））要求写出完整的程序代码，并在matlab软件环境调试并运行通过，连同运行结果一并附上。 一、填空题（1? ×25=25?） 1、Matlab的全称为矩阵实验室。 2、在Matlab编辑器中运行程序的快捷键是：F5 。 3、Matlab的工作界面主要由以下五个部分组成，它们分别是：菜单栏、 工具栏、当前工作目录窗口、工作空间管理窗口和命令窗口。 4、在Matlab中inf表示：无穷大；clc表示：清空命令窗口中的显示内容；more表示：在命令窗口中控制其后每页的显示内容行数；who表示：查阅Matlad内存变量名；whos表示：列出当前工作空间所有变量。 5、在Matlab命令窗口中运行命令Simulink 可以打开Simulink模块库浏览器窗口。 6、求矩阵行列式的函数：det ；求矩阵特征值和特征向量的函数eig 。 7、Matlab预定义变量ans表示：没有指定输出变量名；eps表示：系统精度 ；nargin表示：函数输入参数的个数。 8、Matlab提供了两种方法进行程序分析和优化，分别为：通过Profiler工具优化和通过tic和toc函数进行优化。 9、建立结构数组或转换结构数组的函数为：struct ； 实现Fourier变换在Matlab中的对应函数为：fourier() ；Laplace变换的函数：Laplace() 。
补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样
课程设计报告 题目PID控制器应用 课程名称控制系统仿真院部名称机电工程学院专业 班级 学生姓名 学号 课程设计地点 课程设计学时 指导教师 金陵科技学院教务处制成绩
一、课程设计应达到的目的 应用所学的自动控制基本知识与工程设计方法，结合生产实际，确定系统的性能指标与实现方案，进行控制系统的初步设计。 应用计算机仿真技术，通过在MATLAB软件上建立控制系统的数学模型，对控制系统进行性能仿真研究，掌握系统参数对系统性能的影响。 二、课程设计题目及要求 1．单回路控制系统的设计及仿真。 2．串级控制系统的设计及仿真。 3．反馈前馈控制系统的设计及仿真。 4．采用Smith 补偿器克服纯滞后的控制系统的设计及仿真。 三、课程设计的内容与步骤 (1)．单回路控制系统的设计及仿真。 （a）已知被控对象传函W(s) = 1 / (s2 +20s + 1)。 （b）画出单回路控制系统的方框图。 （c）用MatLab的Simulink画出该系统。 （d）选PID调节器的参数使系统的控制性能较好，并画出相应的单位阶约响应
曲线。注明所用PID调节器公式。PID调节器公式Wc（s）=50(5s+1)/(3s+1) 给定值为单位阶跃响应幅值为3。 有积分作用单回路控制系统 无积分作用单回路控制系统
大比例作用单回路控制系统 （e）修改调节器的参数，观察系统的稳定性或单位阶约响应曲线，理解控制器参数对系统的稳定性及控制性能的影响？ 答：由上图分别可以看出无积分作用和大比例积分作用下的系数响应曲线，这两个PID调节的响应曲线均不如前面的理想。增大比例系数将加快系统的响应，但是过大的比例系数会使系统有比较大的超调，并产生振荡，使稳定性变坏；增大积分时间有利于减小超调，减小振荡，使系统的稳定性增加，但是系统静差消除时间变长，加入微分环节，有利于加快系统的响应速度，使系统超调量减小，稳定性增加。 (2)．串级控制系统的设计及仿真。 （a）已知主被控对象传函W 01(s) = 1 / (100s + 1)，副被控对象传函W 02 (s) = 1 / (10s + 1)，副环干扰通道传函W d (s) = 1/(s2 +20s + 1)。 （b）画出串级控制系统方框图及相同控制对象下的单回路控制系统的方框图。（c）用MatLab的Simulink画出上述两系统。
一、实验目的 1.利用MATLAB 实现循环卷积。 2.比较循环卷积与线性卷积的区别。 二、实验条件 PC 机，MATLAB7.0 三、实验内容 1）循环卷积的定义：两个序列的N 点循环卷积定义为： )0()()()]()([1 0N n m n x m h n x n h N k N N <≤-=?∑-= 利用MATLAB 实现两个序列的循环卷积可以分三个步骤完成： (1)初始化：确定循环点数N ，测量输入2个序列的长度。 (2)循环右移函数：将序列x(n)循环右移，一共移N 次（N 为循环卷积的循环次数），最后将每次循环成的新序列组成一个矩阵V 。 (3)相乘：将x(n)移位后组成的矩阵V 与第二个序列h(n)对应相乘，即得循环卷积结果。程序如下： 程序一： clear;close all ; N=10; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1); stem(xxn1,x1); subplot(3,1,2); stem(xxn2,x2); x11=fft(x1,N);
x12=fft(x2,N); y11=x11.*x12; y1=ifft(y11,N); subplot(3,1,3); n=0:length(y1)-1; stem(n,y1,'.'); title('循环卷积的结果'); xlabel('n');ylabel('y1(n)'); 运行后所得图形如下： 观察所得的循环卷积结果发现并没有呈现周期性的序列，因此将程序做下列改变。程序二： clear;close all; N=40; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; x2=[x2,x2,x2,x2]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1);
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 仿真大作业附题目+ 完整报告内容
设二阶控制系统如图1所示，其中开环传递函数 ) 1(10 )2()(2+=+=s s s s s G n n ξωω 图1 图2 图3 要求： 1、分别用如图2和图3所示的测速反馈控制和比例微分控制两种方式改善系统的性能，如果要求改善后系统的阻尼比ξ =0.707，则和 分别取多少？ 解： 由)1(10 )2()(2 += +=s s s s s G n n ξωω得10 21,10,102===ξωωn t K d T
对于测速反馈控制，其开环传递函数为：) 2()s (2 2n t n n K s s G ωξωω++=； 闭环传递函数为：2 2 2)2 (2)(n n n t n s K s s ωωωξωφ+++= ； 所以当n t K ωξ2 1+=0.707时，347.02)707.0(t =÷?-=n K ωξ； 对于比例微分控制，其开环传递函数为：)2()1()(2 n n d s s s T s G ξωω++=； 闭环传递函数为：) )2 1(2)1()(2 22 n n n d n d s T s s T s ωωωξωφ++++=； 所以当n d T ωξ2 1 +=0.707时，347.02)707.0(=÷?-=n d T ωξ； 2、请用MATLAB 分别画出第1小题中的3个系统对单位阶跃输入的响应图； 解： ①图一的闭环传递函数为： 2 22 2)(n n n s s s ωξωωφ++=，10 21 ,10n ==ξω Matlab 代码如下： clc clear wn=sqrt(10); zeta=1/(2*sqrt(10)); t=0:0.1:12; Gs=tf(wn^2,[1,2*zeta*wn,wn^2]); step(Gs,t)
matlab 基本语句 1.循环语句for for i=s1:s3：s2 循环语句组 end 解释：首先给i赋值s1；然后，判断i是否介于s1与s2之间；如果是，则执行循环语句组,i=i+s3（否则，退出循环.）；执行完毕后，继续下一次循环。 例：求1到100的和，可以编程如下： sum=0 for i=1:1:100 sum=sum+i end 这个程序也可以用while语句编程。 注：for循环可以通过break语句结束整个for循环. 2.循环语句while 例：sum=0;i=1; while(i<=100) sum=sum+i;i=i+1; end 3.if语句 if(条件) 语句 end if(条件) 语句 else 语句 end if(条件) 语句 elseif 语句 end 4.关系表达式：
=,>,<,>=,<=,==(精确等于) 5.逻辑表达式：|(或),&(且) 6.[n,m]=size(A)(A为矩阵) 这样可以得到矩阵A的行和列数 n=length(A)，可以得到向量A的分量个数；如果是矩阵，则得到矩阵A的行与列数这两个数字中的最大值。 7.！后面接Dos命令可以调用运行一个dos程序。 8.常见函数： poly():为求矩阵的特征多项式的函数，得到的为特征多项式的各个系数。如 a=[1,0,0;0,2,0;0,0,3]，则poly(a)＝1 -6 11 -6。相当于poly(a)＝1入^3+(-6)入^2+11入+(-6)。 compan():可以求矩阵的伴随矩阵. sin()等三角函数。 MATLAB在数学建模中的应用(3) 一、程序设计概述 MATLAB所提供的程序设计语言是一种被称为第四代编程语言的高级程序设计语言，其程序简洁，可读性很强，容易调试。同时，MATLAB的编程效率比C/C++语言要高得多。 MATLAB编程环境有很多。常用的有： 1.命令窗口 2.word窗口 3.M-文件编辑器，这是最好的编程环境。 M-文件的扩展名为“.m”。M-文件的格式分为两种： ①λ M-脚本文件，也可称为“命令文件”。 ② M-函数文件。这是matlab程序设计的主流。λ 保存后的文件可以随时调用。 二、MATLAB程序结构 按照现代程序设计的观点，任何算法功能都可以通过三种基本程序结构来实现，这三种结构是：顺序结构、选择结构和循环结构。其中顺序结构是最基本的结构，它依照语句的自然顺序逐条地执行程序的各条语句。如果要根据输入数据的实际情况进行逻辑判断，对不同的结果进行不同的处理，可以使用选择结构。如果需要反复执行某些程序段落，可以使用循环结构。 1 顺序结构 顺序结构是由两个程序模块串接构成。一个程序模块是完成一项独立功能的逻辑单元，它可以是一段程序、一个函数，或者是一条语句。 看图可知，在顺序结构中，这两个程序模块是顺序执行的，即先执行<程序
第七章系统仿真的MATLAB实现 由于计算机技术的高速发展，我们可以借助计算机完成系统的数字仿真。综前所述，数字仿真实质上是根据被研究的真实系统的模型，利用计算机进行实验研究的一种方法。仿真的主要过程是：建立模型、仿真运行和分析研究仿真结果。仿真运行就是借助一定的算法，获得系统的有关信息。 MATLAB是一种面向科学与工程计算的高级语言，它集科学计算、自动控制、信号处理、神经网络和图像处理等学科的处理功能于一体，具有极高的编程效率。MATLAB是一个高度集成的系统，MATLAB提供的Simulink是一个用来对动态系统进行建模、仿真和分析的软件包，它支持线性和非线性系统，能够在连续时间域、离散时间域或者两者的混合时间域里进行建模，它同样支持具有多种采样速率的系统。在过去几年里，Simulink已经成为数学和工业应用中对动态系统进行建模时使用得最为广泛的软件包。 MATLAB仿真有两种途径：（1）MATLAB可以在SIMULINK窗口上进行面向系统结构方框图的系统仿真；（2）用户可以在MATLAB的COMMAND窗口下，用运行m文件，调用指令和各种用于系统仿真的函数，进行系统仿真。这两种方式可解决任意复杂系统的动态仿真问题，前者编辑灵活，而后者直观性强，实现可视化编辑。 下面介绍在MATLAB上实现几类基本仿真。 7.1 计算机仿真的步骤 在学习计算机仿真以前，让我们先总结一下计算机仿真的步骤。 计算机仿真，概括地说是一个“建模—实验—分析”的过程，即仿真不单纯是对模型的实验，还包括从建模到实验再到分析的全过程。因此进行一次完整的计算机仿真应包括以下步骤：
（1）列举并列项目 每一项研究都应从说明问题开始，问题由决策者提供或由熟悉问题的分析者提供。 （2）设置目标及完整的项目计划 目标表示仿真要回答的问题、系统方案的说明。项目计划包括人数、研究费用以及每一阶段工作所需时间。 （3）建立模型和收集数据 模型和实际系统没有必要一一对应，模型只需描述实际系统的本质或者描述系统中所研究部分的本质。因此，最好从简单的模型开始，然后进一步建立更复杂的模型。 （4）编制程序和验证 利用数学公式、逻辑公式和算法等来表示实际系统的内部状态和输入/输出的关系。建模者必须决定是采用通用语言如MATLAB、FORTRAN、C还是专用仿真语言来编制程序。在本教材中，我们选择的是MATLAB和其动态仿真工具Simulink。 （5）确认 确认指确定模型是否精确地代表实际系统。它不是一次完成，而是比较模型和实际系统特性的差异，不断对模型进行校正的迭代过程。 （6）实验设计 确定仿真的方案、初始化周期的长度、仿真运行的长度以及每次运行的重复次数。 （7）生产性运行和分析 通常用于估计被仿真系统设计的性能量度。利用理论定性分析、经验定性分析或系统历史数据定量分析来检验模型的正确性，利用灵敏度分析等手段来检验模型的稳定性。 （8）文件清单和报表结果 （9）实现
仿真报告 课程名称：自动化技术导论 报告题目：MATLAB/simulink系统仿真分析 班级 姓名 学号 xxxxxx自动化学院 2016年4月 软件版本：MATLAB R2010b MATLAB强处理能力 MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而且经过了各种优化和容错处理。在通常情况下，可以用它来代替底层编程语言，如C和C++ 。在计算要求相同的情况下，使用MATLAB的编程工作量会大大减少。MATLAB的这些函数集包括从最简单最基本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数。函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等。 MATLAB图形处理 MATLAB自产生之日起就具有方便的数据可视化功能，以将向量和矩阵用图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用于科学计算和工程绘图。新版本的MATLAB对整个图形处理功能作了很大的改进和完善，使它不仅在一般数据可视化软件都具有的功能（例如二维曲线和三维曲面的绘制和处理等）方面更加完善，而且对于一些其他软件所没有的功能（例如图形的光照处理、色度处理以及四维数据的表现等），MATLAB同样表现了出色的处理能力。同时对一些特殊的可视化要求，例如图形对话等，MATLAB也有相应的功能函数，保证了用户不同层次的要求。另外新版本的MATLAB还着重在图形用户界面（GUI）的制作上作了很大的改善，对这方面有特殊要求的用户也可以得到满足。 MATLAB对许多专门的领域都开发了功能强大的模块集和工具箱。一般来说，它们都是由特定领域的专家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。领域，诸如数据采集、数据库接口、概率统计、样条拟合、优化算法、偏微分方程求解、神经网络、小波分析、信号处理、图像处理、系统辨识、控制系统设计、LMI控制、鲁棒控制、模型预测、模糊逻辑、金融分析、地图工具、非线性控制设计、实时快速原型及半物理仿真、嵌入式系统开发、定点仿真、DSP与通讯、电力系统仿真等，都在工具箱（Toolbox）家族中有了自己的一席之地。MATLAB程序接口
循环结构：for语句 格式： for 循环变量=表达式1：表达式2：表达式3 循环体 end 【注】：表达式1：循环变量初值， 表达式2：步长，为1时，可省略； 表达式3：循环变量终值。 或： for循环变量=矩阵表达式 循环体 end 【注】：执行过程是依次将矩阵的各列元素赋给循环变量，然后执行循环体语句，直至各列元素处理完毕。 2 while语句： 格式： while（条件） 循环体 end 【注】：条件成立时，执行循环体 3
break语句&& continue语句： break：破坏，破坏循环，终止循环的进行，跳出循环，程序将执行循环语句的下一语句。 continue：继续，循环继续，程序将跳过循环体中剩下的语句，继续下一次循环。 4 循环的嵌套—多重循环结构 5 选择结构：if-else语句 格式： if 表达式 程序模块 end 或 if 表达式 程序模块1 else 程序模块2 end 6 switch语句： 格式：
switch 表达式 case 数值1 程序模块1 case 数值2 程序模块2 case 数值3 程序模块3 ...... otherwise 程序模块n end 执行过程：首先计算表达式的值， 然后将其结果与每一个case后面的数值依次进行比较， 如果相等，则执行该case的程序模块； 如果都不相等，则执行otherwise模块中的语句。 switch语句可以替代多分支的if语句，而且switch语句简洁明了，可读性更好。 7 matlab中一些基本知识： END 注意事项 for循环可以通过break语句结束整个for循环
simulink仿真 -1<ξ<0 >> step(tf(4^2,[1,2*(-0.5)*4,4^2])) ξ<-1 >> step(tf(4^2,[1,2*(-1.5)*4,4^2])) ξ＝0 >> step(tf(4^2,[1,2*0*4,4^2])) 0<ξ<1 >> figure >> step(tf(4^2,[1,2*0.1*4,4^2])) >> step(tf(4^2,[1,2*0.2*4,4^2])) >> step(tf(4^2,[1,2*0.3*4,4^2])) >> step(tf(4^2,[1,2*0.4*4,4^2])) >> step(tf(4^2,[1,2*0.5*4,4^2])) >> step(tf(4^2,[1,2*0.6*4,4^2])) >> step(tf(4^2,[1,2*0.7*4,4^2])) >> step(tf(4^2,[1,2*0.8*4,4^2])) >> step(tf(4^2,[1,2*0.9*4,4^2])) ωn不变，ζ减小
ξ=1 >> figure >> step(tf(4^2,[1,2*1*4,4^2])) ξ>1 >> hold on >> step(tf(4^2,[1,2*2.0*4,4^2])) >> step(tf(4^2,[1,2*4.0*4,4^2])) >> step(tf(4^2,[1,2*8.0*4,4^2])) ωn不变，ζ减小 ξ＝0.5，改变ωn时的情况： >> figure >> step(tf(1^2,[1,2*0.5*1,1^2])) （ωn＝1）
>> hold on >> step(tf(2^2,[1,2*0.5*2,2^2])) （ωn＝2）>> step(tf(4^2,[1,2*0.5*4,4^2])) （ωn＝4）>> step(tf(8^2,[1,2*0.5*8,8^2])) （ωn＝8） ζ不变，ωn增大 曲线拟合程序 >> figure >> x=[0,0.2,0.4,0.6,0.8,1.0,1.2]; >> y=[1.135,1.135,1.216,1.351,1.534,1.737,2.0,]; >> plot(x,y,'.') >> hold on >> x1=[0:0.1:1.2]; >> y1=1+0.6*x1+0.2*x1.^2; >> plot(x1,y1) >> y1=1+0.7*x1; >> plot(x1,y1)
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接口的方法 陈金平
while语法： while expression statements end 说明：while expression, statements, end 计算一个表达式，并在该表达式为true 时在一个循环中重复执行一组语句。表达式的结果非空并且仅包含非零元素（逻辑值或实数值）时，该表达式为true。否则，表达式为false。 示例代码如下： function [sum] = summation(ratio, head, top) sum = 0; while (head <= top) sum = sum + ratio ^ head; head = head + 1; end end
假设ratio = 2，head = 0，top = 63 matlab循环语句for怎么用？ matlab中for语句使用方法和应用实例 for循环语句 1、一般格式为： for x（循环变量）= array（数组） commands（执行的循环代码） end 2、array可以是一个数字，也可以是数组，例如输入：for a=5 for a=1:5 for a=1:1:5（以1为步长到5）
只不过在a=1:5和a=1:1:5时，会显示之间的结果，a=5时只显示最后结果。a的变动就是第一次循环a=1，第二次循环a=2，第三次循环a=3，第四次循环a=4，第五次循环a=5。 3、(commands)就是命令，其中的命令行可以很多很多，最常见的就是调用上面说的a变动比如： for a=10 %循环10次 s=a+1 %循环语句 end %结束 上面的a不需要再指定，a的变化就是上面说的先是1，然后是2，3，…那么s就是先s=1+1=2，然后再s=2+1=3，s=3+1=4，…，s=10+1=11，循环结束，就是一个连续加s的指令，最后=11。 4、for语句可以嵌套的，和C一样 for a=5 %第一个for循环
创新实践报告 报告题目: 基于matlab的通信系统仿真学院名称: 信息工程学院 姓名: 班级学号: 指导老师: 二O一四年十月十五日
一、引言 现代社会发展要求通信系统功能越来越强,性能越来越高,构成越来越复杂;另一方面,要求通信系统技术研究与产品开发缩短周期,降低成本,提高水平。这样尖锐对立的两个方面的要求,只有通过使用强大的计算机辅助分析设计技术与工具才能实现。在这种迫切的需求之下,MA TLAB应运而生。它使得通信系统仿真的设计与分析过程变得相对直观与便捷,由此也使得通信系统仿真技术得到了更快的发展。通信系统仿真贯穿着通信系统工程设计的全过程,对通信系统的发展起着举足轻重的作用。通信系统仿真具有广泛的适应性与极好的灵活性,有助于我们更好地研究通信系统性能。通信系统仿真的基本步骤如下图所示: 二、仿真分析与测试 (1)随机信号的生成 利用Matlab中自带的函数randsrc来产生0、1等概分布的随机信号。源代码如下所示: global N N=300; global p
p=0、5; source=randsrc(1,N,[1,0;p,1-p]); (2)信道编译码 1、卷积码的原理 卷积码(convolutional code)就是由伊利亚斯(p 、Elias)发明的一种非分组码。在前向纠错系统中,卷积码在实际应用中的性能优于分组码,并且运算较简单。 卷积码在编码时将k 比特的信息段编成n 个比特的码组,监督码元不仅与当前的k 比特信息段有关,而且还同前面m=(N-1)个信息段有关。 通常将N 称为编码约束长度,将nN 称为编码约束长度。一般来说,卷积码中k 与n 的值就是比较小的整数。将卷积码记作(n,k,N)。卷积码的编码流程如下所示。 可以瞧出:输出的数据位V1,V2与寄存器D0,D1,D2,D3之间的关系。根据模2加运算特点可以得知奇数个1模2运算后结果仍就是1,偶数个1模2运算后结果就是0。 2、译码原理 卷积码译码方法主要有两类:代数译码与概率译码。代数译码主要根据码本身的代数特性进行译码,而信道的统计特性并没有考虑在内。目前,代数译码的主要代表就是大数逻辑解码。该译码方法对于约束长度较短的卷积码有较好的效果,并且设备较简单。概率译码,又称最大似然译码,就是基于信道的统计特性与卷积 码的特点进行计算。在现代通信系统中,维特比译码就是目前使用最广泛的概率 译码方法。 02 1V D D =⊕01232V D D D D =⊕⊕⊕
控制系统仿真 [教学目的] 掌握数字仿真基本原理 控制系统的数学模型建立 掌握控制系统分析 [教学内容] 一、控制系统的数学模型 sys=tf(num,den)%多项式模型，num为分子多项式的系数向量，den为分母多项式的系%数向量，函数tf（）创建一个TF模型对象。 sys=zpk(z,p,k)%z为系统的零点向量，p为系统的极点向量，k为增益值，函数zpk（）创建一个ZPK模型对象。 （一）控制系统的参数模型 1、TF模型 传递函数 num=[b m b m-1b m-2…b1b0] den=[a m a m-1a m-2…a1a0] sys=tf(num,den) 【例1】系统的传递函数为。 >>num=[01124448]; >>den=[11686176105]; >>sys=tf(num,den); >>sys Transfer function: s^3+12s^2+44s+48 ------------------------------------- s^4+16s^3+86s^2+176s+105 >>get(sys) >>set(sys) >>set(sys,'num',[212])
>>sys Transfer function: 2s^2+s+2 ------------------------------------- s^4+16s^3+86s^2+176s+105 【例2】系统的传递函数为。 >>num=conv([20],[11]); >>num num= 2020 >>den=conv([100],conv([12],[1610])); >>sys=tf(num,den) Transfer function: 20s+20 ------------------------------- s^5+8s^4+22s^3+20s^2 【例3】系统的开环传递函数为，写出单位负反馈时闭环传递函数的TF模型。>>numo=conv([5],[11]); >>deno=conv([100],[13]); >>syso=tf(numo,deno); >>sysc=feedback(syso,1) Transfer function: 5s+5 ---------------------- s^3+3s^2+5s+5 【例4】反馈系统的结构图为： R
MATLAB编程如何不用for循环--以DFT变换矩阵为例 缘起：大家都知道MATLAB中用for循环编写的代码执行起来效率不高，如何用矩阵和向量的运算提升效率对每一个用matlab的人来说都是很有必要的，但是此项功夫高手一般不愿意给初学者讲，此功夫是高手和低手的分水岭，高手们更是拿此功夫在初学者面前炫耀。本人当初怀着很恭敬的心向高手请教，高手笑笑说这要我自己编。出于让后来人受益，帮助和我一样无助的求知者。本人今天话了一天时间将此问题研究下，并且将代码毫无保留的公布出来。希望大家能够受益，阿弥陀佛！ 上述W矩阵的第一列代表直流成分，第二列到最后一列是信号的交流成分，可以看出倍频关系！我以前不知道DFT 可以通过矩阵表示。注意matlab中dftmtx实现上述W矩阵的时候没有用1/sqrt(N) 进行归一化！可以通过dftmtx(2)验证，没有1/sqrt(2)。 例1 DFT matrix 是Hermitian的 Nfft = 8; xn=rand(1,Nfft); y=dftmtx(Nfft)*xn.' %结果是个列向量y=dftmtx(Nfft)*x.' 和y=fft(x,Nfft)是等价的 y = xn*dftmtx(Nfft) %结果是个行向量 y=fft(xn,Nfft) %结果是个行向量 %dft变换公式，n代表时域采样点，k代表频域采样点 Y(k)=sum(x.*exp(-j*2*pi*n*k)) 相应的，dftmtx(Nfft)产生的矩阵中，第k行，n列元素=exp(-j*2*pi*k*n/Nfft)，与x.'相乘正好对应fft变换后的每个频点值。
例1：双重循环求DFT N = 8; x=rand(1,N); for k=0:N-1 sum=0; % 注意每个X(k)的值不应该受上次计算的影响 for m=0:N-1 w=exp(-j*2*pi/N.*m.*k); %DFT matrix 的每一行的元素是不同的 sum=sum+x(m+1).*w; % 这个循环实际上是计算DFT matrix的每一行与信号x的内积 end X(k+1)=sum; %Matlab下标从1开始 end % 注意X 是个行向量，是个数组，我以前不知道X(k)的循环赋值的结果是个行向量 y=fft(x,N) %验证 例2：单循环求DFT %% x 是行向量的Version，结果X也是行向量 N = 8; x=rand(1,N); m=0:N-1; sum=0; for k=0:N-1 %核心思想是内积运算：x躺着，后面的必须站着 X(k+1)= x*exp(-j*2*pi/N.*m'.*k) ; % 以前这样的表达物理含义是不明确的--> X(k+1)= x*exp(-j*2*pi/N.*m.*k).' ; end y=fft(x,N) %验证 %% x 是列向量的Version，结果X是行向量，fft(x)的结果是列向量 N = 8; x=rand(1,N)'; m=0:N-1; sum=0; for k=0:N-1 % 核心思想是内积运算：x躺着，后面的必须站着 X(k+1)= exp(-j*2*pi/N.*m.*k)*x ; % 以前这样的表达物理含义是不明确的--> X(k+1)= x*exp(-j*2*pi/N.*m.*k).' ; end y=fft(x,N) %验证结果y是个列向量，X是y的转置。因为X按下标赋值的结果是个行向量！ 例3：不用循环求DFT %% x 是行向量的Version，结果X也是行向量 N = 8; xn=rand(1,N); %一次去掉2个循环，不要试图一次去掉一个 n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n'*k; %看看如何不用for循环来实现，向量的外积N-by-1 multiply 1-by-N WNnk=WN.^nk; %WNnk就是dftmtx(8)，dftmtx原来是这样构造的 Xk=xn*WNnk; % 究竟nk这个向量外径是n'*k还是k'*n要看内积在哪一维执行