MATLAB习题

Matlab

楼宇091 黄娟 20090271

1.将图形窗口分成两个窗格，分别绘制出函数：

在[0,3]区间上的曲线，并利用axis 调整轴刻

度纵坐标刻度，使y1在[0,12]区间上，y2在

[-2,1.5]区间上。

>> x=0:0.1:3;

>> subplot(2,1,1)

>> y1=2*x+5;

>> plot(x,y1)

>> subplot(2,1,2)

>> y2=x.^2-3*x+1;

>> plot(x,y2)

>> axis([0,3,-2,1.5])

00.51 1.52 2.5346

00.51 1.52 2.53

-2-1

135

2221+-=+=x x y x y

>> performance=input('请输入百分制成绩: '); >> if (performance>100) | (performance<0) error

end

>> switch floor(performance/10)

case{0,1,2,3,4,5}

disp('E');

case{6}

disp('D');

case{7}

disp('C');

case{8}

disp('B');

case{9}

disp('A');

case{10}

disp('A+');

otherwise

disp('成绩不合理');

‘请输入百分数制成绩：80’

end

matlab练习题及答案

第1章MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？答：MATLAB语言有一下特点： 1、起点高。 2、人机界面适合科技人员。 3、强大而简易的作图功能。 4、智能化程度高。 5、功能丰富，可扩展性强。 1.2 MATLAB系统由那些部分组成？答：MATLAB系统主要由开发环境、MATLAB数学函数库、MMATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。

1.5 如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器。（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到

MATLAB自适应滤波去噪

《MATLAB自适应滤波去噪》课程设计报告 1.课程设计目的此次课程设计目的是为了让我们学会使用MATLAB进行计算机仿真，使用自适应滤波法设计一个语音去噪声电路。培养我们的电路设计思路及其算法，明白理论与实践相结合的重要性，培养了我们的实际操作能力以及锻炼我们对实际问题的分析与解决的能力。 2.课程设计内容 2.1 LMS自适应算法原理自适应过程一般采用典型LMS自适应算法，但当滤波器的输入信号为有色随机过程时，特别是当输入信号为高度相关时，这种算法收敛速度要下降许多，这主要是因为输入信号的自相关矩阵特征值的分散程度加剧将导致算法收敛性能的恶化和稳态误差的增大。此时若采用变换域算法可以增加算法收敛速度。变换域算法的基本思想是：先对输入信号进行一次正交变换以去除或衰减其相关性，然后将变换后的信号加到自适应滤波器以实现滤波处理，从而改善相关矩阵的条件数。因为离散傅立叶变换DFT 本身具有近似正交性，加之有FFT快速算法，故频域分块LMS FBLMS 算法被广泛应用。 FBLMS算法本质上是以频域来实现时域分块LMS算法的，即将时域数据分组构成N个点的数据块，且在每块上滤波权系数保持不变。其原理框图如图2所示。FBLMS 算法在频域内可以用数字信号处理中的重叠保留法来实现，其计算量比时域法大为减少，也可以用重叠相加法来计算，但这种算法比重叠保留法需要较大的计算量。块数据的任何重叠比例都是可行的，但以50％的重叠计算效率为最高。对FBLMS算法和典

型LMS算法的运算量做了比较，并从理论上讨论了两个算法中乘法部分的运算量。本文从实际工程出发，详细分析了两个算法中乘法和加法的总运算量，其结果为：复杂度之比＝FBLMS实数乘加次数／LMS实数乘加次数＝（25Nlog2N＋2N－4）／[2N(2N－1)] 采用ADSP的C语言来实现FBLMS算法的程序如下： for(i＝0;i＜＝30;i＋＋) {for(j＝0;j＜＝n－1;j＋＋) {in[j]＝input[i×N＋j;] rfft(in,tin,nf,wfft,wst,n); rfft(w,tw,wf,wfft,wst,n); cvecvmlt(inf,wf,inw,n); ifft(inw,t,O,wfft,wst,n); for(j＝0,j＜＝N－1;j＋＋) {y[i×N＋j]＝O[N＋j]．re; e[i×N＋j]＝re fere[i×N＋j]－y[i×N＋j]; temp[N＋j]＝e[i×N＋j;}

MATLAB-课后习题集

第2章 MATLAB 概论 1、与其他计算机语言相比较，MA TLAB 语言突出的特点是什么？ 2、 MATLAB 系统由那些部分组成？ 3、安装MATLAB 时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 4、 MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 5、如何启动M 文件编辑/调试器？ 6、存储在工作空间中的数组能编辑吗？如何操作？ 7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 8、如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 9、在MA TLAB 中有几种获得帮助的途径？第3章 MATLAB 数值运算 3.1 在MA TLAB 中如何建立矩阵?? ? ? ??194375，并将其赋予变量a ？ 3.2 有几种建立矩阵的方法？各有什么优点？ ①直接输入法，如a=[2 5 7 3]，优点是输入方法方便简捷； ②通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改； ③由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵； ④通过数据文件建立，该方法可以调用由其他软件产生数据。 3.3 在进行算术运算时，数组运算和矩阵运算各有什么要求？进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 3.4 数组运算和矩阵运算的运算符有什么区别？在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 3.5 计算矩阵??????????897473535与???? ? ?????638976242之和。

matlab课后习题解答第二章

第2章符号运算习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗不能从显示形式判断数据类型，而必须依靠class指令。〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗理解自由符号变量的确认规则。〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1)

ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。〖目的〗理解subexpr 指令。〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8（1）通过符号计算求t t y sin )(=的导数 dt dy 。（2）然后根据此结果，求- =0t dt dy 和2 π = t dt dy 。〖目的〗 diff, limit 指令的应用。如何理解运行结果。〖解答〗 syms t

matlab习题及答案2

MATLAB 基本运算 1．在MATLAB 中如何建立矩阵?? ?? ??194375，并将其赋予变量a ？>>a=[573;491] 2．在进行算术运算时，数组运算和矩阵运算各有什么要求？进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 3．数组运算和矩阵运算的运算符有什么区别？在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 4．计算矩阵??????????897473535与??????????638976242之和。>>a=[535;374;798]; >>b=[242;679;836]; >>a+b ans = 7 779 1413151214 5．计算??????=572396a 与?? ????=864142b 的数组乘积。>>a=[693;275]; >>b=[241;468]; >>a.*b ans = 12 36384240 6．“左除”与“右除”有什么区别？在通常情况下，左除x=a\b 是a*x=b 的解，右除x=b/a 是x*a=b 的解，一般情况下，a\b ≠b/a 。

7．对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。>>A=[492;764;357]; >>B=[372628]’; >>X=A\B X = -0.5118 4.0427 1.3318 8．已知：???? ??????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 >>a=[123;456;789]; >>a.^2 ans = 1 4916 253649 6481 >>a^2 ans = 30 364266 81961021261509．[]7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。相当于a=[11011]。 10．在sin(x )运算中，x 是角度还是弧度？在sin(x)运算中，x 是弧度，MATLAB 规定所有的三角函数运算都是按弧度进行运算。

自适应滤波LMS与RLS地matlab实现

MATLAB 仿真实现LMS 和RLS 算法的二阶AR 模型及仿真结果分析一、题目概述：二阶AR 模型如图1a 所示，可以如下差分方程表示： )()()2()1()()(21n d n v n x a n x a n v n x +=----= （1）图1a 其中，v(n)是均值为0、方差为0.965的高斯白噪声序列。，为描述性参数， .95,0,195.021=-=a a 设x(-1)=x(-2)=0，权值， =0.04①推导最优滤波权值（理论分析一下）。②按此参数设置，由计算机仿真模拟权值收敛曲线并画出，改变步长在此模拟权值变化规律。③对仿真结果进行说明。④应用RLS 算法再次模拟最优滤波权值。解答思路：（1）高斯白噪声用normrnd 函数产生均值为0、方差为0.965的正态分布随机1*N 矩阵来实现。随后的产生的信号用题目中的二阶AR 模型根据公式（1）产生，激励源是之前产生的高斯白噪声。（2）信号长度N 取为2000点，用以观察滤波器权值变化从而估计滤波器系数，得到其收敛值。（3）仿真时分别仿真了单次LMS 算法和RLS 算法下的收敛性能以及100次取平均后的LMS 和RLS 算法的收敛性能，以便更好的比较观察二者的特性。（4）在用不同的分别取3个不同的值仿真LMS 算法时，值分别取为

0.001,0.003,0.006；用3个不同的λ值仿真RLS 算法时λ值分别取为1,0.98,0.94，从而分析不同步长因子、不同遗忘因子对相应算法收敛效果的影响。二、算法简介 1．自适应算法的基本原理自适应算法的基本信号关系如下图所示：图 1b 自适应滤波器框图输入信号x(n)通过参数可调的数字滤波器后产生输出信号y(n)，将其与参考信号d(n)进行比较，形成误差信号e(n)。e(n)通过某种自适应算法对滤波器参数进行调整，最终是e(n)的均方值最小。当误差信号e(n)的均方误差达到最小的时候，可以证明信号y(n)是信号d(n)的最佳估计。 2. LMS 算法简介 LMS 算法采用平方误差最小的原则代替最小均方误差最小的原则，信号基本关系如下： 1 ()()() ()()()(1)()2()()(0,1,2,....1) N i i i y n w n x n i e n d n y n w n w n e n x n i i N μ-=-=-+=+-=-∑ （2）写成矩阵型式为： ()()()()()() (1)()2()() T y n W n X n e n d n y n W n W n e n X n μ==-+=+ （3）式（3）中，W(n) 为n 时刻自适应滤波器的权值， 011()[(),(),....()]T N W n w n w n w n -=，N 为自适应滤波器的阶数，本设计中取为2000；X( n) 为n 时刻自适应滤波器的参考输入矢量，由最近N 个信号采样值构成，()[(),(1),....(1)]T X n x n x n x n N =--+；d ( n) 是期望的输出值；e ( n) 为自适应滤波器的输出误差调节信号(简称失调信号) ；μ是控制自适应速度与稳定性的增益常数，又叫收敛因子或步长因子。

matlab课后习题及答案详解

第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？

matlab基础练习题带答案

Matlab 基础练习题常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（）（A ）_num （B ）num_ （C ）num- （D ）-num 2、在MA TLAB 中，要给出一个复数z 的模，应该使用（）函数。（A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、下面属于MATLAB 的预定义特殊变量的是？（）（A ）eps （B ）none （C ）zero （D ）exp 4、判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放一个中文字符‘啊’需要占用2个字节。（错，都是2个字节） 5、判断：MA TLAB 中，i 和j （对） 6、判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（错，后面还有很多位小数） 7、在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字节。 9、在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。数组 1、在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序的逆序排列输出，应该使用下面的（）指令。（A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（）。（A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都赋值为0，应该在指令窗中输入（）（A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时，若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（）（A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应该是L1 =___ 1 1 1 1 1 0 0 0 0___。

MATLAB习题及答案 (入门)

习题 1，计算?? ?? ??=572396a 与??? ???=864142b 的数组乘积。 a.*b ans = 12 36 3 8 42 40 2，对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。 inv(a)*b ans = -0.5118 4.0427 1.3318 3，已知：??? ? ? ?????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 a.*a ans = 1 4 9 16 25 36 49 64 81 >> a^2 ans = 30 36 42 66 81 96 102 126 150

4，角度[]604530=x ，求x 的正弦、余弦、正切和余切。(应用sin,cos,tan.cot) >> sin(x./pi) ans = -0.1242 0.9826 0.2465 >> cos(x./pi) ans = -0.9923 -0.1857 0.9692 >> tan(x./pi) ans = 0.1252 -5.2915 0.2543 >> cot(x./pi) ans = 7.9894 -0.1890 3.9321 5，将矩阵?? ????=7524a 、??????=3817b 和??? ???=2695c 组合成两个新矩阵：（1）组合成一个4?3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ?? ??? ???? ???237 912685 574 [reshape(a,4,1),reshape(b,4,1),reshape(c,4,1)] ans = 4 7 5 5 8 6 2 1 9 7 3 2 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即

自适应波束形成与Matlab程序代码注解

1.均匀线阵方向图（1）matlab 程序 clc; clear all; close all; imag=sqrt(-1); element_num=32;%阵元数为8 d_lamda=1/2;%阵元间距d与波长lamda的关系 theta=linspace(-pi/2,pi/2,200); theta0=0;%来波方向 w=exp(imag*2*pi*d_lamda*sin(theta0)*[0:element_num-1]'); for j=1:length(theta) a=exp(imag*2*pi*d_lamda*sin(theta(j))*[0:element_num-1]'); p(j)=w'*a; end patternmag=abs(p); patternmagnorm=patternmag/max(max(patternmag)); patterndB=20*log10(patternmag); patterndBnorm=20*log10(patternmagnorm); figure(1) plot(theta*180/pi,patternmag); grid on; xlabel('theta/radian') ylabel('amplitude/dB') title([num2str(element_num) '阵元均匀线阵方向图','来波方向为' num2str(theta0*180/pi) '度']); hold on; figure(2) plot(theta,patterndBnorm,'r'); grid on; xlabel('theta/radian') ylabel('amplitude/dB') title([num2str(element_num) '阵元均匀线阵方向图','来波方向为' num2str(theta0*180/pi) '度']); axis([-1.5 1.5 -50 0]);

自适应载波同步及其Matlab仿真

成都理工大学工程技术学院本科毕业论文自适应载波同步及其Matlab仿真作者姓名：专业名称：指导老师：年月日

摘要自适应滤波算法的研究是现在社会自适应信号处理中最为活跃的研究课题之一。找寻收敛速度快，计算简单，数值稳定性好的自适应滤波算法是研究人员不断努力追求的目标。本设计在论述自适应滤波基本原理的基础上，说明了几种当前几种典型的自适应滤波算法和应用。并对这几种典型自适应滤波算法的性能特点进行简单的比较，给出了算法性能的综合评价。载波同步是无线通信接收机的主要功能之一，其对通信系统质量的提高至关重要。随着新算法涌现和芯片处理速度的提高，不同的解决方案不断的提出。自适应载波同步是一种依据自适应算法的同步方法，内容新颖。本课题在介绍自适应算法和载波同步问题的基础上，详细讨论了平方差分环路法和锁相环路法，具体包括代价函数、代价函数的导数、迭代公式和原理图等，并在论文的第三部分给出了这两种方法的Matlab仿真。仿真结果验证了这两种方法在跟踪载波相位方面是满足要求的，且收敛速度较快。关键词：自适应滤波载波同步平方差分环路锁相环路法

Abstact The research of adaptive filtering algorithm is one of the most activity tasks, the goal that researchers want to pursue is to find an adaptive filtering algorithm that converge fast and compute simplely. Based on the basis adaptive filtering principle, this paper introduces several typical adaptive algorithms and applications, then compares those algorithm's characters and gives the orithm performance evaluation. Carrier synchronization is one of the main functions of Wireless communications receiver，it is essential for the improvement in the quality of the communication system. With the emergence of new algorithms and the speed improvement of chip processing, different solutions is proposed continuously. Adaptive carrier synchronization is a synchronization method based on adaptive algorithms, and its content is innovative. Based on the introducing of adaptive algorithm and carrier synchronization, this issue has a detailed discussion of the square difference method and the PLL loop method, including its cost function, cost function derivative, iterative formula and schematic, etc. And the third part of the paper gives two methods of Matlab simulation.Simulation results show the two methods with tracking the carrier phase is to meet the requirements, and convergence speedly. Keywords:adaptive filter, carrier synchronization, differential circle square , phase-locked loop method

matlab习题及答案

2. 用MATLAB 语句输入矩阵A 和B 3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来，赋给B 矩阵，用magic(8)A =命令生成A 矩阵，用上述命令检验一下结果是不是正确。 4．用数值方法可以求出∑=++++++==63 63622284212i i S Λ，试不采用循环的形式求出和式的数值解。由于数值方法是采用double 形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。 5．选择合适的步距绘制出下面的图形。（1）)/1sin(t ，其中)1,1(-∈t ；（2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 6. 试绘制出二元函数2 2 2 2 )1(1)1(1),(y x y x y x f z +++ +-= =的三维图和三视图 7. 试求出如下极限。（1）x x x x 1)93(lim +∞ →；（2）1 1lim 0-+→→xy xy y x ；（3）2 2)()cos(1lim 2 2 220 0y x y x e y x y x +→→++- 8. 已知参数方程? ??-==t t t y t x sin cos cos ln ，试求出x y d d 和3 /2 2d d π=t x y 9. 假设?-=xy t t e y x f 0 d ),(2 ，试求2222 22y f y x f x f y x ??+???-?? 10. 试求出下面的极限。（1）??????-++-+-+-∞→1)2(1 161141121lim 2222n n Λ；（2）)131211( lim 2 222π πππn n n n n n n ++++++++∞ →Λ 11. 试求出以下的曲线积分。（1）?+l s y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=， )20(π≤≤t 。

matlab课后答案完整版

ones表示1矩阵 zeros表示0矩阵 ones(4)表示4x4的1矩阵 zeros(4)表示4x4的0矩阵 zeros(4,5)表示4x5的矩阵 eye(10,10)表示10x10的单位矩阵rand(4,5)表示4x5的伴随矩阵 det(a)表示计算a的行列式 inv(a)表示计算a的逆矩阵 Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩 [v,d]=eig(a)对角矩阵 b=a’表示求a矩阵的转置矩阵 sqrt表示求平方根 exp表示自然指数函数 log自然对数函数 abs绝对值第一章一、5（1） b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7]; >> e=b*c e = 5271 11574 -11336 664 1978 3112 （2）a=50:1:100 二、1 、x=-74; y=-27; z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)))+pi) z = 2、a=::; >> b=exp*a).*sin(a+ 3、x=[2 4; 5]; y=log(x+sqrt(1+x.^2))/2 y =4、a*b表示a矩阵和b矩阵相乘 a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列 A(m,:)表示取a矩阵第m行的全部元素 A(:,n)表示取a矩阵的第n列全部元素 A./B表示a矩阵除以b矩阵的对应元素， B.\A等价于A./B A.^B表示两个矩阵对应元素进行乘方运算A.^2表示a中的每个元素的平方 A^2表示A*A 例：x=[1,2,3]; y=[4,5,6]; z=x.^y z= 1 3 2 729 指数可以是标量（如y=2）.底数也可以是标量（如x=2） 5、a=1+2i; >> b=3+4i; >> c=exp((pi*i)/6) c = + d=c+a*b/(a+b) d = + 第二章二、4、（1） y=0;k=0; >> while y<3 k=k+1; y=y+1/(2*k-1); end >> display([k-1,y-1/(2*k-1)]) ans = 第三章二1（1）

Matlab习题

习题 1 1. 执行下列指令，观察其运算结果, 理解其意义： (1) [1 2;3 4]+10-2i (2) [1 2; 3 4].*[0.1 0.2; 0.3 0.4] (3) [1 2; 3 4].\[20 10;9 2] (4) [1 2; 3 4].^2 (5) exp([1 2; 3 4]) (6)log([1 10 100]) (7)prod([1 2;3 4]) (8)[a,b]=min([10 20;30 40]) (9)abs([1 2;3 4]-pi) (10) [1 2;3 4]>=[4,3;2 1] (11)find([10 20;30 40]>=[40,30;20 10]) (12) [a,b]=find([10 20;30 40]>=[40,30;20 10]) (提示：a 为行号，b 为列号) (13) all([1 2;3 4]>1) (14) any([1 2;3 4]>1) (15) linspace(3,4,5) (16) A=[1 2;3 4];A(:,2) 2. 执行下列指令，观察其运算结果、变量类型和字节数，理解其意义： (1) clear; a=1,b=num2str(a),c=a>0, a= =b, a= =c, b= =c (2) clear; fun='abs(x)',x=-2,eval(fun),double(fun) 3. 本金K 以每年n 次，每次p %的增值率(n 与p 的乘积为每年增值额的百分比)增加，当增加到rK 时所花费的时间为 ) 01.01ln(ln p n r T += (单位：年) 用MA TLAB 表达式写出该公式并用下列数据计算：r =2, p =0.5, n =12. 4．已知函数f (x )=x 4 -2x 在(-2, 2)内有两个根。取步长h =0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。（提示：求近似根等价于求函数绝对值的最小值点） ? 5. (1) 用z=magic(10)得到10阶魔方矩阵； (2) 求z 的各列元素之和； (3) 求z 的对角线元素之和(提示：先用diag(z)提取z 的对角线)； (4) 将z 的第二列除以3;

用matlab实现自适应图像阈值分割最大类方差法代码

%用matlab实现自适应图像阈值分割最大类方差法代码clear; warning off; SE = strel('diamond',4); BW1 = imread('cameraman.tif'); BW2 = imerode(BW1,SE); BW3 = imdilate(BW2,SE); BW4 = BW1-BW3; %rgb转灰度 if isrgb(BW4)==1 I_gray=rgb2gray(BW4); else I_gray=BW4; end figure,imshow(I_gray); I_double=double(I_gray);%转化为双精度 [wid,len]=size(I_gray); colorlevel=256; %灰度级 hist=zeros(colorlevel,1);%直方图 %threshold=128; %初始阈值 %计算直方图 for i=1:wid for j=1:len m=I_gray(i,j)+1; hist(m)=hist(m)+1; end end hist=hist/(wid*len);%直方图归一化 miuT=0; for m=1:colorlevel miuT=miuT+(m-1)*hist(m); end xigmaB2=0; for mindex=1:colorlevel threshold=mindex-1; omega1=0; omega2=0; for m=1:threshold-1 omega1=omega1+hist(m); end omega2=1-omega1; miu1=0; miu2=0; for m=1:colorlevel if m

MATLAB习题参考答案

第1章 MATLAB简介 1、MA TLAB的主要特点有：①语言简洁，编程效率高。②人机界面友善，交互性好。③绘图功能强大，便于数据可视化。④学科众多、领域广泛的MATLAB工具箱。⑤源程序的开放性。 MATLAB的典型应用领域有：①自动控制②汽车③电子④仪器仪表⑤生物医学⑥信号处理⑦通信等。 2、填空题 ⑴命令窗口、命令历史窗口、当前目录窗口 ⑵查阅、保存、编辑 ⑶清除图形窗、清除命令窗口中显示内容、清除MATLAB工作空间中保存的变量。 3、如果想查看某一变量具体内容或者对其修改操作，可以在工作空间中双击该变量名称，可以打开数组编辑器，在数组编辑器中可以查看变量的具体内容，也可以对其修改。如果想删除MATLAB内存中的变量，可以在工作空间中选中该变量，然后利用工作空间窗口的菜单命令或工具条中的快捷图标进行删除。 4、1+2+3+4+5+... （+ 后面可以直接跟...，也可以在+和...中加一个空格。） 6+7+8+9 1+2+3+4+5 ... （5后面必须跟一个空格，不能直接跟...，否则报错，这在预置一个+6+7+8+9 大数组时很重要。）第2章矩阵与数值数组 1、填空题： ⑴非数、无穷大、机器零阈值，浮点数相对精度,eps= 2.2204e-016。 ⑵全下标、单下标。 2、阅读程序题：（本题主要考察数组的寻访、赋值和简单运算，提示：带；的语句不显示结果） ⑴ans = 2 3 7 Sa = 10 20 30 A = 1 20 5 30 9 10 4 6 8 10 ⑵ A = 1 3 5 7 9 11 13 15

2 4 6 8 10 12 14 16 ans = 1 5 9 13 2 6 10 14 3 7 11 15 4 8 12 16 A = 0 0 5 7 0 0 13 15 2 4 0 0 10 12 0 0 ⑶ ans = -1 -4 6 4 ans = 3 0 5 -2 ans = 3 6 9 12 ans = 3 6 9 12 3、A=magic(4)； L=A<10 L = 0 1 1 0 1 0 0 1 1 1 1 0 1 0 0 1 4、省略。第3章字符串、元胞和构架数组 1. 直接创建法：S1=['Where there is life,' 'there is hope. '] %注意第2行要加入6个空格使其与第一行字符数相等，否则报错。 S2=str2mat('Where there is life,','there is hope.') S3=strvcat('Where there is life,','there is hope.') %后两种方法则不用考虑两行字符的数目 2. 填空题：A（2，3）；A{2，3} 3. 阅读程序题： ⑴subch = ABc12

基于MATLAB_Simulink机器人鲁棒自适应控制系统仿真研究

第18卷第7期系统仿真学报?V ol. 18 No. 7 2006年7月Journal of System Simulation July, 2006 基于MATLAB/Simulink机器人鲁棒自适应控制系统仿真研究高道祥，薛定宇（东北大学教育部暨辽宁省流程工业综合自动化重点实验室，沈阳 110004）摘要：介绍了一种在MATLAB/Simulink环境下进行机器人鲁棒自适应控制系统仿真的方法，利用Matlab软件强大的数值运算功能，将系统模型用Matlab语言编写成M-Function（或S-Function）文件，通过User-Defined-Function模块嵌入到Simulink仿真环境中，可以充分发挥Simulink模块实时的动态仿真功能，简化仿真模型的设计，修改和调整。基于M-Function建立机器人系统模型的方法可以推广到其他复杂控制系统的建模，SimMechanics在建立多自由度连杆机器人受控对象仿真模型时，简单可靠。关键词：机器人；Matlab/Simulink；SimMechanics；仿真；鲁棒自适应控制中图分类号：TP391.9 文献标识码：A 文章编号：1004-731X(2006) 07-2022-04 Simulation Research of Robust Adaptive Control System for Robotic Manipulators Based on MATLAB/Simulink GAO Dao-xiang, XUE Ding-yü (Key Laboratory of Process Industry Automation, Ministry of Education, Northeast University, Shenyang 110004, China) Abstract: A simulation method of robust adaptive control was proposed for the robotic manipulator system. The method took the advantage of the powerful computing function of Matlab to programme M-function (or S-Function) for the system model by Matlab language and embedded it to the Simulink by User-Defined-Function module. The real time dynamic simulating function of Simulink would be exerted adequately and the design, modification and adjust of the system model could be greatly simplified. The method of constructing manipulator control system model based on M-Function could be generalized to the other complicated control system and SimMechanics would make the n-links manipulator model conveniently and credibly. Key words: robotic manipulator; Matlab/Simulink; SimMechanics; simulation; robust adaptive control 引言一个新的控制算法在付诸使用之前，无论从经济原因还是技术角度，都需要经过仿真阶段来测试控制系统的性能和缺陷。尤其对复杂系统控制的研究，虽然仿真并不能说明控制算法是绝对合理与可靠的，但随着仿真技术的发展，仿真的确是系统设计必不可少的中间步骤。 Matlab/Simulink以其强大的数学运算能力，方便实用的绘图功能及语言的高度集成性成为系统仿真和自动控制领域首选的计算机辅助设计工具。Simulink可以将可视化的模块很方便地组成系统模型的仿真框图，对于使用普通Simulink模块不易搭建的复杂控制系统，用Matlab语言编写M-Function或S-Function文件，通过User-Defined-Function 模块嵌入到Simulink仿真环境中，大大扩充了Simulink的功能。对于机器人这类的复杂控制系统，通过Simulink可以很方便的建立其仿真模型。机器人控制系统仿真模型中不易采用普通Simulink模块搭建的部分是控制器模型和受控对象—机器人模型，可以采用Matlab语言编写M-Function实现控制器和机器人模型。收稿日期：2005-05-09 修回日期：2005-08-02 作者简介：高道祥（1972-），男，山东蓬莱人，博士生，研究方向为机器人鲁棒自适应控制。薛定宇（1963-），男，辽宁沈阳人，教授，博导，研究方向为控制系统CAD，机器人控制。另外，由于SimMechanics提供了机构的仿真模块集，对于n自由度的连杆机器人，也可以采用SimMechanics模块进行组建。 1 n连杆机器人的仿真模型如果不考虑摩擦力等外界干扰的作用，机器人的动力学方程可由下式描述， τ= + +) ( ) , ( ) (q G q q q C q q M (1) 式中，q q q , ,是1 × n向量，表示各个关节的位置，速度，加速度。) (q M是n n×阶对称正定的惯量矩阵。q q q C ) , (是1 × n向量，表示离心力和哥氏力项。) (q G是1 × n向量，表示重力项。τ表示外界输入的控制力矩。由式(1)可以看出，机器人的动力学模型是一个高度复杂，强耦合的非线性时变方程，尤其模型的复杂程度随着自由度的增加呈指数上升趋势。因此，在用Matlab/Simulink进行机器人控制系统的仿真研究时，需要寻求一种简单可靠却行之有效的方法建立机器人控制系统仿真模型。采用M-Function定制的Simulink模块与普通模块一样具有输入和输出向量，控制器与机器人受控对象的仿真模型函数可用如下关系式描述， ) , , , , (0q q q q q f d d d = τ(2) ) , , (1τ q q f q =(3)