matlab语言复习参考模板

1、一 MATLAB语言的发展PPT语言页：对语言扩展说明：VHDL，UML，C系列语言（C,C+,C#,.NET,JAVA）等；均大量使用面向对象技术。另外，还有专业软件，如PS，周正龙用得很好，忽悠了全国人民说笑而已。 PPT能作页：重点说仿真。什么是仿真？利用模型完成实验研究，我们利用数学模型，建模合仿真离不开MATLAB。举例： 导弹信息处理系统开发，如作用距离指标论证，推理；原理建模，边界条件，作战环境等方面综合考虑，形成最后用于仿真得模型，验证，计算开发实际系统，周期520年；最后强调一点，涉及项目内容不要在上网机上出现，保密是生命线，碰不得。导致，千百人的工作可能白费，巨额经济损失，

2、严重的将危害国家安全。 求解线性系统PPT: 非常方便，一条指令解决。绘图PPT：非常直观，观察函数特性。均有周期性；图2 值大于零，图3调频。二、MATLAB语言的特点及功能2.1 特点编程过程中，自己慢慢体会2.2 功能 8点， PPT开放式结构页， 程序源代码开放，可以修改、编程补充。三、MATLAB安装和启动3.1 MATLAB的安装和工具包选择 必选项内容 3.2 启动MATLAB 界面组成，说明四、MATLAB操作桌面基本用法 第二节课开始4.1 Command Window运行入门 数值、变量及运算在命令窗的表达4.2 界面操作要点 各界面操作简述4.3 工作空间的管理指令数据、

3、变量的查询、保存和定义 问：系统中，比1大的最小数？4.4 Editor/Debugger和脚本编写初步 编程界面说明五、帮助系统及其使用5.1 常用帮助指令 HELP ,LOOKFOR5.2 帮助导航浏览器Help Navigator/Browser4个组成部分，5.3 VIDEO视频演示帮助系统新增加功能，声音，视频讲解六、 MATLAB Notebook简述 1 / 216.1 Notebook的启动6.2 M-book模板的使用 目前，是一种方便的作业工具，文档与运算，画图合为一体， 界面友好性有待加强。七 总结 1 软件运行条件；桌面用法；标点符号用法、命名规则；帮助系统使用；Not

4、ebook使用主体内容：1 概述 12 数值计算 基本内容 2 3 符号计算 14 绘图及可视化 15 程序设计 1 an applewarn: eps a very little number inf nan nargin nariout realmax realminbe care : a./b & b.a is different and never forget '*'warn: *数组乘建立数组a和b，进行乘法运算说明：指令中的“.*”符号表示乘法是在两个数组相同位置上的元素间进行的,把这种乘法称为“数组乘”.“*”符号表示乘法是在两个矩阵间进行，乘法满足代数

5、运算规则。数据显示格式key : ; % key: clf clc clear 变量空间it's a inportant part of matlab.m.fileand we also should know the help window notebook for ferther use.1 在安装MATLAB软件时，哪个组件是必须选择的？假如不“勾选”这个组件，那么就不可能建立MATLAB工作环境。2 数字1.5e2，1.5e3 中的哪个与1500相同吗？3 请指出如下5个变量名中，哪些是合法的？abcd-2xyz_33chana变量ABCDefgh4 在MATLAB环境中，比1

6、大的最小数是多少？5 设 a = -8 , 运行以下三条指令，问运行结果相同吗？为什么？ w1=a(2/3) w2=(a2)(1/3) w3=(a(1/3)26 指令clear, clf, clc各有什么用处？7 以下两种说法对吗？（1）“MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。”（2）MATLAB指令窗中显示的数值有效位数不超过7位。”8 想要在MATLAB中产生二维数组，下面哪些指令能实现目的？ S=1,2,3;4,5,6;7,8;9 S=1 2 3;4 5 6;7 8 9 S=1，2，3；4，5，6；7，8，9 %整个指令在中文状态下输入9 试为 ,z=z1+z2编

7、写一个解题用的M脚本文件？ 1 答： matlab组件2 a1=1.5e2,a2=1.5e3 解：a1 = 150 a2 = 15003 答：xyz_3 , ABCDefgh4 答： 1eps5 w1 = -2.0000 + 3.4641i w2 = 4.0000 w3 = -2.0000 + 3.4641i 答： 不同6 答：clear, 清除内存变量和函数；clf 清除图对象；clc 清除指令窗7 答：（1）错。MATLAB进行数值的表达精度使用双精度进行。指令窗中的数据显示精度可由用户设定，仅仅为了显示简洁才采用较少位数显示。（2）错。指令窗中的数据显示精度可由用户设定，仅仅为了显示简洁

8、才采用较少位数显示。8 答：S=1 2 3;4 5 6;7 8 9能产生二维数组 9 clear; z1=4+3*i; z2=1+2*i; z=z1+z2 a boy数值计算符号计算代数方程的解尽量不要使用inv(a)*b指令，而应采用ab 后者的计算速度比前者快，计算精度高，尤其是矩阵a较大时。另外，除法指令的适应性更强，即使对于非方阵的a，也能给出最小二乘解。傅立叶变换设序列x(n)的长度为M，则x(n)的N（NM）点离散傅立叶变换对定义为clear; close all;N=32;n=0:N-1; k=n;nk=n'*k;xn=exp(j*pi*n/8)WN=exp(-j*2*p

9、i/N);Wnk=WN.nk;Xk=xn*Wnkplot(abs(Xk)试验设计：9 用快速卷积法(FFT变换)计算下面两个序列的卷积，并测试快速卷积和直接卷积的时间 取M20，N15。10 验证频域采样与时域采样的对偶性 试编写MATLAB程序，来验证对频谱X(ejw)进行等间隔采样。（1） 产生一个三角波序列x(n) (2)对M40，计算x(n)的64点DFT,并图示x(n)和X(k)=DFTx(n)，k=0,1,263。(3)对（2）中所得X(k)在0,2*pi上进行32点抽样得X1(k)=X(2k),k=0,1,31(4)求X1(k)的32点IDFT，即x1(n)IDFTX1(k),k

10、=0,1,231(5)画出x1(n)32（x1(n)的周期延拓，周期为32）的波形图，评述x1(n)32与x(n)的关系。并根据频域采样理论加以解释。矩阵创建、运算和数组运算；多项式运算和卷积 数值微积分；代数方程组和微分方程求解函数极值、拟合和插值；随机数据的统计描述 1 正文一、 命令行的基本操作1.1创建矩阵的方法1.2用matlab函数创建矩阵rand(4)-4*4随机矩阵eye(4) 4*4 单位矩阵zeros(4,1) 全部元素都为0的4*1矩阵ones （1,4）全部元素都为1的1*4矩阵1.3. 矩阵的修改数组及其运算小规模数组的直接输入法指 令含 义指 令含 义<小于&

11、gt;=大于等于<=小于等于=等于>大于=不等于标准数组生成函数指令含义指令含义eye产生单位数组（对高维不适用）rand产生均匀分布随机数组 magic产生魔方数组（对高维不适用）randn产生正态分布随机数组ones 产生全1数组zeros产生全“0”数组diag产生对角形数组（对高维不适用）gallery产生特殊的测试矩阵用逻辑函数分 类具 体 描 述数组非0判断all数组A所有元素都不是0，返回1。any数组A不是全0元素，返回1。生成逻辑数组false按指定大小，创建全0逻辑数组。true按指定大小，创建全1逻辑数组。logical创建逻辑数组：1对应输入数组中的非0元素

12、，其余都为0。数据对象判断isempty是否空阵isprime是否质数isfinite是否有限数isreal是否实数isinf是否无穷大isletter是否字母（用于字符串）isnan是否非数isspace是否空格（用于字符串）数据类型判断isa是否指定类别ishandle是否图柄ischar是否字符串islogical是否逻辑类型isglobal是否全局变量isnumeric是否数值类型二、数据的保存与获取产生并保存sin(x)在0,4*pi的波形，然后调用并画图。三、矩阵运算 矩阵的一些特殊操作矩阵的谱分解和矩阵函数函数的数组运算与矩阵运算四、 多项式运算和卷积 1 poly 产生特征多项

13、式系数向量2.roots 求多项式的根3.conv多项式乘运算4.deconv多项式除运算5 多项式微分五、卷积和 fourier分析卷积离散傅立叶变换六 数值微积分 不推荐求数值极限；例子说明；积分例子：积分准确性与使用的算法有关！七、代数方程组求解三种情况各自求解方法及特点。范数，条件数和方程解的精度 ：判断解的可靠性八、微分方程求解多种解法，注意解的格式和步骤。建立函数，然后选择解算指令求解。九、函数极值点两个函数用法十、随机数据的统计描述 交互界面使用和说明十一、拟合与插值 示例和函数介绍1 由指令rand('state',0),A=rand(3,5)生成二维数组A，试

14、求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。clear;rand('state',0),A=rand(3,5) A = 0.9501 0.4860 0.4565 0.4447 0.9218 0.2311 0.8913 0.0185 0.6154 0.7382 0.6068 0.7621 0.8214 0.7919 0.1763 find(A>0.5)Ca=A(find(A>0.5) L=A>0.5islogical(L);X=A(L) L = 1 0 0 0 1 0 1 0 1 1 1 1 1 1 0X = 0.9501 0.6

15、068 0.8913 0.7621 0.8214 0.6154 0.7919 0.9218 0.7382 ans = 1 3 5 6 9 11 12 13 14Ca = 0.9501 0.6068 0.8913 0.7621 0.8214 0.6154 0.7919 0.9218 0.7382 2 已知有理分式，其中，。（1）求该分式的商多项式和余多项式。（2）用程序验算是否成立。%Nx=(3*x3+x)*(x3+0.5),Dx=(x2+2*x-2)*(5*x3+2*x2+1)p1=conv(3 0 1 0,1 0 0 0.5);p2=conv(1 2 -2 ,5 2 0 1);q ,r=de

16、conv(p1,p2);cq='商多项式为'cr='余多项式为'disp(cq,poly2str(q,'s');disp(cr,poly2str(r,'s'); 商多项式为 0.6 s - 1.44余多项式为 -4.4409e-016 s6 + 8.8818e-016 s5 + 21.88 s4 - 5.34 s3 - 5.52 s2 + 4.58 s - 2.88 p3=0.6 -1.44,p4=-4.4409e-016 8.8818e-016 21.88 -5.34 -5.52 4.58 -2.88p=conv(p2,p3)

17、+p4Ps1= poly2str(p,'s') p3 = 0.6000 -1.4400p4 = -0.0000 0.0000 21.8800 -5.3400 -5.5200 4.5800 -2.8800p = 3.0000 0.0000 1.0000 1.5000 0.0000 0.5000 0Ps1 = 3 s6 + 1.3323e-015 s5 + 1 s4 + 1.5 s3 + 8.8818e-016 s2 + 0.5 s Ps2=poly2str(p1,'s') Ps2 = 3 s6 + s4 + 1.5 s3 + 0.5 s 3 a=1:12;利用re

18、shape、rot90、fliplr、 flipud、 diag、tril、triu指令形成新矩阵。a=1:12;b=reshape(a,3,4)c1=rot90(b)c2=flipud(b)c3=tril(b) b = 1 4 7 10 2 5 8 11 3 6 9 12c1 = 10 11 12 7 8 9 4 5 6 1 2 3c2 = 3 6 9 12 2 5 8 11 1 4 7 10c3 = 1 0 0 0 2 5 0 0 3 6 9 0 4 求解二元函数方程组的解clear;x=-2:0.5:2;y=x;X,Y=meshgrid(x,y);F1=sin(X-Y);F2=cos(X

19、+Y);v=-0.2,0,0.2;contour(X,Y,F1,v);hold on,contour(X,Y,F2,v),hold offx0,y0=ginput(2);Fun='sin(x(1)-x(2),cos(x(1)+x(2)'%Fun=inline('sin(x(1)-x(2),cos(x(1)+x(2)','x')xy,f,exit=fsolve(Fun,x0(2),y0(2) Optimization terminated: first-order optimality is less than options.TolFun.xy

20、= 0.7854 0.7854f = 1.0e-016 * 0 0.6123exit = 1 Optimization terminated: first-order optimality is less than options.TolFun.xy = 0.7854 0.7854f = 1.0e-016 * 0 0.6123exit = 1 xy,f,exit=fsolve(Fun,0.8,0.8) Optimization terminated: first-order optimality is less than options.TolFun.xy = 0.7854 0.7854f =

21、 1.0e-015 * 0 -0.1608exit = 1 5 求函数在区间中的最小值点。fun=inline('(sin(5*t)2*exp(0.06*t2)-1.5*t*cos(2*t)+1.8*abs(t+0.5)','t')fmin=fminbnd(fun,-5,5) fun = Inline function: fun(t) = (sin(5*t)2*exp(0.06*t2)-1.5*t*cos(2*t)+1.8*abs(t+0.5)fmin = -1.2850 6 用quad求取的数值积分，并保证积分的绝对精度为。clear;fun=inline(&

22、#39;exp(-abs(x).*abs(sin(x)','x');Iq=quad(fun,-5*pi,1.7*pi,10e-10),Iql=quadl(fun,-5*pi,1.7*pi,1e-9), Iq = 1.0878Iql = 1.0878 7 已知矩阵，运行指令B1=A.(0.5), B2=A(0.5), 可以观察到不同运算方法所得结果不同。（1）请分别写出根据B1, B2恢复原矩阵A的程序。（2）用指令检验所得的两个恢复矩阵是否相等。A=1 2;3 4B1=A.0.5,B2=A0.5 A = 1 2 3 4B1 = 1.0000 1.4142 1.7321

23、2.0000B2 = 0.5537 + 0.4644i 0.8070 - 0.2124i 1.2104 - 0.3186i 1.7641 + 0.1458i RB1=B1.2,RB2=B2*B2 RB1 = 1.0000 2.0000 3.0000 4.0000RB2 = 1.0000 + 0.0000i 2.0000 + 0.0000i 3.0000 - 0.0000i 4.0000 - 0.0000i L1=isequal(RB1,RB2),L1 = 0 8 求的实数解。clear;fun=inline('-0.5+t-10.*exp(-0.2*t).*abs(sin(sin(t)

24、','t'),ft=fzero(fun,0.5) fun = Inline function: fun(t) = -0.5+t-10.*exp(-0.2*t).*abs(sin(sin(t)ft = 2.7341 符号计算1 字符串、元胞和架构数组2 符号对象和符号表达式3 符号表达式和符号函数的操作4 符号微积分5 符号积分变换6 符号线性代数的求解7 符号微分方程的求解8 利用MAPLE的深层符号计算资源运行funtool命令，出现可视化数学分析界be care:自由变量的判断1观察一个数（在此用记述）在以下四条不同指令作用下的异同：a = b = sym( )c

25、= sym( ,'d ' )d = sym( ' ' )%这给出完全准确值在此， 分别代表具体数值 7/3 , pi/3 , pi*3(1/3) ；而异同通过vpa(abs(a-d) , vpa(abs(b-d) , vpa(abs(c-d)等来观察。a1=pi*3(1/3);a=a1,b=sym(a1),c=sym(a1,'d'),d=sym('a1'),ad=vpa(abs(a-d) ,bd= vpa(abs(b-d) , cd=vpa(abs(c-d) a = 4.5310b =5101408179057732*2(-50)

26、c =4.5309606547207899041040946030989d =a1ad =abs(a1-4.5309606547207899041040946030989)bd =abs(a1-4.5309606547207899041040946030989)cd =abs(a1-4.5309606547207899041040946030989) 2说出以下三条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号”对象？3/7+0.1sym(3/7+0.1)vpa(sym(3/7+0.1)class(3/7+0.1) ans =double class(sym(3/7+0.1) ans =sym class(vpa(sym(3/7+0.1) ans =sym 3在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是独立自由变量。sym('sin(w*t)') , sym('a*exp(-X)' ) , sym('z*exp(j*theta)')syms w t,findsym(sin(w*t),1) ans =w syms a X,findsym(a*exp(-X),1) ans =a