天津理工大学中环信息学院1ppt课件_第1页
天津理工大学中环信息学院1ppt课件_第2页
天津理工大学中环信息学院1ppt课件_第3页
天津理工大学中环信息学院1ppt课件_第4页
天津理工大学中环信息学院1ppt课件_第5页
已阅读5页,还剩63页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、天津理工大学中环信息学院天津理工大学中环信息学院电气工程仿真与辅助设计电气工程仿真与辅助设计主要内容:主要内容:1 1、数据构造与语句构造、数据构造与语句构造2 2、根本控制流程构造、根本控制流程构造3 3、M-M-函数设计函数设计4 4、图形可视化、图形可视化5 5、数据处置、数据处置2.1 2.1 数据构造与语句构造数据构造与语句构造第第2 2章章1 1、MATLABMATLAB根本操作根本操作MATLABMATLAB任务界面:任务界面:标题栏、标题栏、菜单栏、菜单栏、工具栏、工具栏、命令窗口、命令窗口、任务空间、任务空间、历史命令窗口历史命令窗口等等命令窗口命令窗口 任务空任务空间间 历

2、史命令历史命令 当出途径当出途径 2 2、编程规范、编程规范 变量命名规那么变量命名规那么 防止运用无意义的字母组合防止运用无意义的字母组合 Temperature Temperature Degree_C/ Degree_F Degree_C/ Degree_F 函数命名规那么函数命名规那么 模块名模块名_ _功能名功能名 Key_GetKey() Key_GetKey() 书写规范书写规范3 3、变量、变量 变量命名规那么:变量命名规那么: 必需以字母开头;必需以字母开头; 可以由字母、数字和下划线混合组成;可以由字母、数字和下划线混合组成; 变量长度应不大于变量长度应不大于3131个;个

3、; 字母区分大小写。字母区分大小写。如:如: 1 1、DegreeDegree、degreedegree 2 2、12Myvar12Myvar、_Myvar _Myvar 正确的变量名如正确的变量名如a a、a1a1、a_1.a_1.在命名的过程中,防止运用无意义的字母组合。在命名的过程中,防止运用无意义的字母组合。 特殊变量特殊变量尽量防止重新赋值尽量防止重新赋值前往变量名称变量名称功能功能变量名称变量名称功能功能ansans用于结果的缺用于结果的缺省变量值省变量值pi pi圆周率圆周率realmaxrealmax最大浮点数值最大浮点数值realminrealmin最小浮点数值最小浮点数值i

4、 i虚数单位虚数单位j j虚数单位虚数单位infinf正无穷大正无穷大NaNNaN不定式不定式(0/0,inf/inf)(0/0,inf/inf)epseps机器浮点运算机器浮点运算误差限误差限narginnargin函数的输入变函数的输入变量数目量数目 回调上一行命令回调上一行命令 回调下一行命令回调下一行命令 左移一个字符左移一个字符 右移一个字符右移一个字符Delet Backspace Delet Backspace 删除字符删除字符Ctrl+C Ctrl+C 中断中断MATLABMATLAB义务义务helphelp协助协助clearclear去除变量去除变量clcclc清屏清屏wha

5、twhat查看当出途径下的查看当出途径下的m m文件文件who, whoswho, whos查看查看workspaceworkspace变量变量load, saveload, save储存、读取储存、读取matmat文件文件在在MATLABMATLAB里共有六种数据类型,每一种里共有六种数据类型,每一种类型可以是一维、二维和多维的。类型可以是一维、二维和多维的。双精度型双精度型数值量常用数值量常用doubledouble, 取值范围取值范围-1.7-1.710308103081.71.71030810308字符串型字符串型单引号扩起来,单引号扩起来,hellohello多维数组多维数组是矩阵的

6、拓展,用是矩阵的拓展,用 括起来。括起来。单元数组单元数组是矩阵的直接扩展,可存储是矩阵的直接扩展,可存储 恣意类型的信息恣意类型的信息,Ai,j,Ai,j类与对象类与对象如传送函数类,如传送函数类,tf(G)tf(G)符号变量符号变量syms asyms a可用于公式推导可用于公式推导 数据的存储与读取 数据存储:save fname A1 A2.Am其中fname为文件名,自动生成fname.mat的数据文件类型。 数据读取:load fname A1 A2.Am 留意:需求存储的数据A1 A2.Am之间用空格隔开, 不能采用逗号隔开举例阐明1 1直接赋值语句直接赋值语句 变量名列表变量名

7、列表= =赋值表达式赋值表达式 等号左边的变量名列表为等号左边的变量名列表为MATLABMATLAB语句的前往语句的前往值,假设一次前往多个结果,那么变量列表用值,假设一次前往多个结果,那么变量列表用 括起来,各变量间用逗号分隔;括起来,各变量间用逗号分隔; 等号右边为表达式,可以是矩阵运算或函数等号右边为表达式,可以是矩阵运算或函数调用,可以由分号调用,可以由分号(;)(;)、逗号、逗号(,)(,)或回车终了。或回车终了。留意:语句结尾有无分号的显示结果区别是什么。留意:语句结尾有无分号的显示结果区别是什么。 假设省略了赋值变量和等号,那么表达式运假设省略了赋值变量和等号,那么表达式运算的结

8、果将赋值保管变量算的结果将赋值保管变量ansans。ansans将存放最后一将存放最后一次无赋值变量语句的结果。次无赋值变量语句的结果。如:如: 2,5,6,7,8,3 2,5,6,7,8,3 运转结果:运转结果:ans=ans= 2 5 6 7 8 3 2 5 6 7 8 32 2函数调用语句函数调用语句 前往变量列表前往变量列表=函数名输入变量列表函数名输入变量列表普通函数名应对应在普通函数名应对应在MATLABMATLAB途径下的一个文件,途径下的一个文件,如,函数名如,函数名my_funmy_fun应该对应在应该对应在my_fun.mmy_fun.m文件;文件;前往变量列表和输入变量列

9、表可以有假设干变量,前往变量列表和输入变量列表可以有假设干变量,变量间用变量间用“,分隔,此外,前往变量还允许用,分隔,此外,前往变量还允许用“空格分隔空格分隔如:如:A_inv=inv(A) A_inv=inv(A) 其中其中invinv求解逆矩阵求解逆矩阵 u,s,v=svd(X) ( u,s,v=svd(X) (其中其中svdsvd对矩阵对矩阵X X进展奇进展奇特值分解,特值分解,s s为对角矩阵,为对角矩阵,X=uX=u* *s s* *v)v) 矩阵的赋值 矩阵用“ 扩起来 “分号表示换行 “,或“空格表示同一行元素分隔如:直接输入 A=1 2, 3; 4 5 6;7, 8 9 冒号

10、操作符 a=0:1:10 a=linspace(0,1,10) a=logspace(1,2,10)矩阵下标按列优先陈列?用两个下标来表示。第i行j列的元素用A(i,j)来表示?用一个下标来表示。对于矩阵按列操作,可用单下标援用A(i) 。?在下标的表达式中运用冒号表示矩阵的一部分。 ?下标援用从序号1开场。矩阵行列删除B=1 2 3;4 5 6; 7 8 9; B(2,: )=B= 1 2 3 7 8 9选择特定行列B=1 2 3;4 5 6; 7 8 9; B(2,: ) B= 4 5 6 例: A= 1 4 7 2 5 8 3 6 9 A(4)+A(8)为了判别矩阵的尺寸,Matlab提

11、供了以下函数: m,n=size(A) %前往矩阵的行数m和列数n M=size(A) %在矩阵M中前往矩阵的行数、 列数 len=length(A) %前往矩阵A行数和列数中的 最大值例子:x=0.3,0.4,0.6,1.5,0.8,1.2x(3)x(1 2 5)b=16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 1b(3,4)b(2 4,1 3)b(2,5,8)b(2;5;8) b(b10) size(b)nrow,ncol=size(b) 单位矩阵eye(n)%生成n维单位阵eye(m,n) %生成mn的单位阵eye(m,n) %生成mn的单位阵eye(size

12、(A)%生成与A矩阵同样大 小的单位矩阵全全1 1矩阵矩阵 ones(n) ones(m,n) ones(m,n) ones(size(A) ones(n) ones(m,n) ones(m,n) ones(size(A)全零矩阵全零矩阵zeros(n) zeros(m,n) zeros(m,n) zeros(size(A)zeros(n) zeros(m,n) zeros(m,n) zeros(size(A)随机矩阵随机矩阵rand(n) rand(m,n) rand(m,n) rand(size(A)rand(n) rand(m,n) rand(m,n) rand(size(A)生成的随机

13、矩阵的各个元素值在生成的随机矩阵的各个元素值在0 0和和1 1之间。之间。对角阵对角阵 diag(n) diag(n) 对角阵对角阵魔术矩阵魔术矩阵 magic(n) magic(n)魔术矩阵的一切行元素和、列元素和、对角线元素和魔术矩阵的一切行元素和、列元素和、对角线元素和均相等均相等转置 B=Ahermit转置,有复数的变共轭复数加减乘 A+B A-B A*B左除 AB 即AX=B的解X=A-1B右除 A/B 即XB=A的解X=AB-1B-1A-1翻转 fliplr(左右翻转)flipud上下翻转 rot90逆时针旋转90乘方 ABA为方阵点运算 A.*B A./B A.B A.B A.

14、A.*B表示矩阵A、B的相应元素之间直接进展乘法。留意:点乘运算A、B矩阵的维数一样逻辑运算& (与) |(或) (非)结果是一个0-1矩阵。当逻辑表达式的值为真时,赋值1,否那么为0 逻辑函数all :当某列的元素都为真时,前往值为1,否那么前往0。最终运算结果为一个0-1行向量any :当向量中至少有一个元素为真时,前往值为1,否那么前往0。最终运算结果为一个0-1行向量 find:用于查找向量中的真元素的下标,前往由一切真元素下标构成的列向量。比较运算在MATLAB里共有六个关系运算符 大于 =大于等于 = 等于 =不等于关系运算符将生成一个0-1矩阵,当运算数相应元素为真时,对

15、应位置上生成1,否那么为0。1、全英文形状输入2、文件类型: 数据文件 *.mat M文件*.m 图片文件 *.fig 模型文件*.mdl3、, ; : ( ) 4、点运算 a= 1 2 3 4b= 1 2 3 4a*b= 7 10 15 22 a.*b= 1 4 9 16 t=1 2 3 4; s=2.t s=2 4 8 162.2 2.2 根本控制流程构造根本控制流程构造第第2 2章章1 1、顺序构造、顺序构造最简单的语句构造,逐条执行语句即可。如:a=1 2 3;4 5 6;b=2 5 8;7 2 0;s=a.*b2 2、循环构造、循环构造MATLAB中的循环语句包括for循环和whil

16、e循环两种类型。根本格式为:for 循环变量=起始值:步长:终止值循环体end 步长的缺省值是1。步长可以在正实数或负实数范围内恣意指定, 对于正数,循环变量的值大于终止值时,循环终了; 对于负数,循环变量的值小于终止值时,循环终了。循环次数确定根本格式为:while 表达式循环体end 其表达式是个逻辑表达式,假设为真,那么执行循环体的内容,执行后再判别表达式能否为真,假设为假那么跳出循环体,向下继续执行,否那么继续执行循环体。 break:从循环体中跳出,并使循环终了循环次数不确定3 3、选择构造、选择构造 转移构造转移构造if-elseif格式假设逻辑表达式的值为真,那么执行语句1,假设

17、逻辑表达式的值为假,那么判别逻辑表达式2的值,假设为真,那么执行语句2,否那么向下执行,假设一切表达式均为假,执行语句el 。if 逻辑表达式1执行语句1elseif 逻辑表达式2执行语句2elseif 逻辑表达式n执行语句nelse执行语句el end循环语句例2另解: s=0; for m=1:1000, s=s+m; if s1000, break; end, end, m3 3、选择构造、选择构造 开关构造开关构造表达式的值和哪种情况的值一样,就执行哪种情况中的语句,假设都不同,那么执行otherwise中的语句。留意:无需像C言语那样在下一个case之前加break语句switch

18、表达式标量或字符串case 值1语句1case 值2语句2case 值n语句notherwise语句owend 2.3 M-2.3 M-函数设计函数设计第第2 2章章1 1、MATLABMATLAB程序设计根本原那么程序设计根本原那么? 后面的内容是程序的注解,要擅长运用注解使程序更具可读性。? 养成在主程序开头用clear指令去除变量的习惯,以消除任务空间中其他变量对程序运转的影响。但留意在子程序中不要用clear。? 参数值要集中放在程序的开场部分,以便维护。要充分利用MATLAB工具箱提供的指令来执行所要进展的运算,在语句运转之后输入分号使其及中间结果不在屏幕上显示,以提高执行速度。?

19、程序尽量模块化,也就是采用主程序调用子程序的方法,将一切子程序合并在一同来执行全部的操作。? 充分利用Debugger来进展程序的调试设置断点、单步执行、延续执行? 设置好MATLAB的任务途径,以便程序运转。2 2、MATLABMATLAB程序的根本组成构造程序的根本组成构造? 阐明? 去除命令:去除workspace中的变量和图形clear,close? 定义变量:包括全局变量的声明及参数值的设定? 逐行执行命令:指MATLAB提供的运算指令或工具箱提供的公用命令? . . .? 控制循环:包含for,if.else,switch,while等语句 ? 逐行执行命令? . . .? 绘图命

20、令:将运算结果绘制出来3 3、MATLABMATLAB程序类型程序类型对于简单的计算可直接在指令窗口中输入指令,但指令数多的时候,那么用源程序文件即M文件设计。M文件:用Matlab言语编写的可以在Matlab环境中运转的程序称为M文件。 M文件以.m为扩展名, M文件分为两种:M-脚本文件命令文件 和M函数文件 ? M-脚本文件:根据用户要求,运用Matlab语 言组成一定功能的Matlab指令集合.运用方法:在MATLAB的提示符下键入该文件名。M-脚本文件运转后,所产生的变量驻留在Matlab任务空间中,只需不用clear命令去除,(或封锁Matlab窗口)这些变量会不断保管在Matla

21、b任务空间中.适用范围:适用于小规模运算例例: :下面语句假设以下面语句假设以.m.m为扩展名存盘为扩展名存盘, ,就构成就构成 M M 脚本文件脚本文件M M 脚本文件:脚本文件: num=2 num=2* *1,2;1,2; den=conv(conv(1,0,1,3),1,2,2); den=conv(conv(1,0,1,3),1,2,2); G1=tf(num,den); G1=tf(num,den); G=ss(G1); G=ss(G1); a,b,c,d=ssdata(G); a,b,c,d=ssdata(G); Ab=a-b Ab=a-b* *c c Bb=b Bb=b Cb=

22、c Cb=c Db=0 Db=0 step(Ab,Bb,Cb,Db) step(Ab,Bb,Cb,Db)? M函数文件:函数文件: M函数文件就是建立一个函数,它可以同函数文件就是建立一个函数,它可以同Matlab的根本函数一样加以运用。的根本函数一样加以运用。 根本构造:根本构造: function 前往变量列表前往变量列表=函数名输入变量列函数名输入变量列表表 注释阐明语句段,由注释阐明语句段,由%引导引导 输入、前往变量格式的检测输入、前往变量格式的检测 函数体语句函数体语句留意:留意:?函数文件的第一行必需包括函数文件的第一行必需包括“function“function这个这个关键字。

23、关键字。?函数文件的文件名必需和函数定义的函数称号函数文件的文件名必需和函数定义的函数称号一样一样, , ?M M函数文件中所定义的变量是内部变量,函数文件中所定义的变量是内部变量, M M函数函数执行完后这些变量随之消逝。执行完后这些变量随之消逝。?输入参数的定义用小括号,假设有多个输输入参数的定义用小括号,假设有多个输入参数那么用逗号分隔;输出参数的定义用中括入参数那么用逗号分隔;输出参数的定义用中括号号,假设有多个输出参数那么用逗号或空格分,假设有多个输出参数那么用逗号或空格分隔。隔。 例:例:function y=mean(x) function y=mean(x) %This is

24、a M function file %This is a M function file m=length(x) m=length(x); y=sum(x)/m y=sum(x)/m ; 将上面程序保管:将上面程序保管:mean1.M mean1.M 构造了一个构造了一个M M函数函数文件,该函数可实现对平均数的求取。随时可调文件,该函数可实现对平均数的求取。随时可调用。用。 如:如:r=1r=1:9999; mean(r) mean(r) 运转得到:运转得到: ans= 50 ans= 501 1函数定义行关键字函数定义行关键字functionfunction functionout1,ou

25、t2,.=filename(in1,in2,.) functionout1,out2,.=filename(in1,in2,.)输入和输出前往的参数个数分别由输入和输出前往的参数个数分别由narginnargin和和nargoutnargout两个两个MATLABMATLAB保管的变量来给出。保管的变量来给出。2 2第一行协助行第一行协助行, ,即即H1H1行行, ,以以% %开头,作为开头,作为lookforlookfor指令搜索指令搜索的行的行3 3函数体阐明及有关注解:以函数体阐明及有关注解:以% %开头,用以阐明函数的作用开头,用以阐明函数的作用及有关内容及有关内容; ;假设不希望显示

26、某段信息,可在它的前面加空行假设不希望显示某段信息,可在它的前面加空行4 4函数体语句:函数体内运用的除前往和输入变量这些在函数体语句:函数体内运用的除前往和输入变量这些在functionfunction语句中直接援用的变量以外的一切变量都是部分变量,即语句中直接援用的变量以外的一切变量都是部分变量,即在该函数前往之后,这些变量会自动在在该函数前往之后,这些变量会自动在MATLABMATLAB的任务空间中去除掉。的任务空间中去除掉。假设希望这些中间变量成为在整个程序中都起作用的变量,那么可假设希望这些中间变量成为在整个程序中都起作用的变量,那么可以将它们设置为全局变量。以将它们设置为全局变量。

27、 2.4 2.4 图形可视化图形可视化第第2 2章章? 图形绘制 )? 图形标注? 图形控制plotyplot(x,y)plot(x1,y1,option1,x2,y2,option2,.)plot控制:控制:功能:在调用函数plot时,可以指定线型,颜 色,和数据点的图标。格式:plotx,y, color_linestyle_marker阐明:参数color_linestyle_marker 为一个字 符串,由颜色、线型、数据点的图标组成。 字符定义字符定义颜色颜色字符定义字符定义颜色颜色y yyellowyellow黄黄m mmagentamagenta洋红洋红c ccyancyan青青

28、r rredred红红g ggreengreen绿绿b bblueblue兰默兰默许许w wwhitewhite白白k kblackblack黑黑字符定义字符定义线型线型字符定义字符定义线型线型- -实线默许实线默许-虚线虚线: :点划线点划线-.-.点连线点连线nonenone无线无线? 图形绘制 ? 图形标注)? 图形控制 text(x,y,字符串)在图形的指定坐标位置(x,y)处,标示单引号括起来的字符串。 gtext(字符串) 利用鼠标在图形的某一位置标示字符串。 title(字符串)在所画图形的最上端显示阐明该图形标题的字符串。 xlabel(字符串),ylabel(字符串设置x,y

29、坐标轴的称号? 图形绘制 ? 图形标注? 图形控制)figure1;figure2;figure(n)翻开不同的图形窗口,以便绘制不同的图形。grid on; grid off 显示/不显示格栅线hold on; hold off 保管/不保管当前的坐标系axisxmin xmax ymin ymax设置坐标轴subplotm,n,k 分割图形显示窗口,m:上下分割个数,n:左右分割个数,k:子图编号semilogx;semilogy 绘制以x或y轴为对数坐标以10为底,y轴(或x轴)为线性坐标的半对数坐标图形。? 预备绘图需求的数据;? 指定绘图的窗口或者区域;? 调用根本绘图命令;? 选择

30、线型、颜色、数据点外形;? 坐标轴控制,包括显示范围、刻度线、比例、网格线;? 标注控制,包括坐标轴称号、标题、相应文本等。 MATLABMATLAB提供了丰富的绘图功能提供了丰富的绘图功能help graph2dhelp graph2d可得到一切画二维图形的命令可得到一切画二维图形的命令help graph3dhelp graph3d可得到一切画三维图形的命令可得到一切画三维图形的命令P26 例2-7 例2-8 例2-92.5 2.5 数据处置数据处置第第2 2章章1 1、矩阵分解、矩阵分解? 三角分解三角分解 ? 特征值分解特征值分解 L,U=lu(A) L,U=lu(A) % %求矩阵求

31、矩阵A A的的LULU分解分解 X Xinv(A) %inv(A) %求矩阵求矩阵A A的逆矩阵的逆矩阵 d=det(A) d=det(A) % %求矩阵求矩阵A A的行列式的行列式 d deigeigA A; V V,DDeigeigA A; V V,DDeigeigA A,nobalancenobalance2 2、多项式的创建、多项式的创建 直接输入法直接输入法例如:例如: A=1,2;3,4 A=1,2;3,4 poly(A) % poly(A) %创建方阵创建方阵A A的特征多项式的特征多项式A=1 2;3 4;poly(A)A=1 2;3 4;poly(A)得到得到ans=1 -5

32、 -2;ans=1 -5 -2;即即x2-5x-2x2-5x-2 poly(a) % poly(a) %假设假设a a为向量为向量bn bn-1b1 b0bn bn-1b1 b0,那么创建那么创建(x-b0 )( x-b1) (x-bn-1 )(x- (x-b0 )( x-b1) (x-bn-1 )(x- bn )bn )生成的多项式的系数向量生成的多项式的系数向量b=1 2 3;poly(b)b=1 2 3;poly(b)得到得到ans=1 -6 11 -6;ans=1 -6 11 -6;即即x3-6x2+11x-6x3-6x2+11x-6多项式常用函数多项式常用函数 roots()%求多项

33、式的根求多项式的根p=poly(A) %求矩阵的特征多项式求矩阵的特征多项式polyval(p,x) %求当多项式求当多项式p的参数为某个的参数为某个特定值特定值x时的多项式的值,时的多项式的值,p为系数向量,为系数向量,polyvalm(p,A) %和和polyval作用一样,输入参作用一样,输入参数值也可为方阵数值也可为方阵A例:例:w=1-32;roots(w)得到得到21即即x2-3x+2=0,得到得到x=1,x=2求求x3-6x2+11x+6,当当x=1时的解时的解a=1-6116;polyval(a,1)卷积和解卷积卷积和解卷积 c=conv(A,B) c=conv(A,B) %

34、%卷积函数卷积函数 ,多项式乘法函数,多项式乘法函数Q,R=deconv (A,B) %Q,R=deconv (A,B) %解卷积函数,多项式除法函数解卷积函数,多项式除法函数导数导数 Pd Pdpolyder(P)polyder(P) % %求多项式求多项式P P的导数的导数 Pm Pmpolyder(P1,P2) polyder(P1,P2) % %求多项式求多项式P1P1和多项式和多项式P2P2乘积的导数乘积的导数 Q,D=polyder(P1,P2) Q,D=polyder(P1,P2) % %求多项式求多项式P1P1除以多项式除以多项式P2P2商的导数商的导数多项式拟合 从最小二乘法

35、的意义上,从最小二乘法的意义上,polyfitpolyfit函数将拟合出函数将拟合出所给数据的多项式系数,其调用格式为:所给数据的多项式系数,其调用格式为:p=polyfit(x,y,n)p=polyfit(x,y,n)其中,其中,x x和和y y为知数据的横坐标和纵坐标向量,为知数据的横坐标和纵坐标向量,n n为多项式的次数。为多项式的次数。多项式插值 多项式插值是指根据给定的有限个样本点,产生另外的估计点多项式插值是指根据给定的有限个样本点,产生另外的估计点以到达数据更为平滑的效果。该技巧在信号处置与图像处置上以到达数据更为平滑的效果。该技巧在信号处置与图像处置上运用广泛。运用广泛。 所用

36、指令有一维的所用指令有一维的interp1interp1、二维的、二维的interp2interp2、三维的、三维的interp3interp3。这些指令分别有不同的方法。这些指令分别有不同的方法methodmethod,设计者可以,设计者可以根据需求选择适当的方法,以满足系统属性的要求。根据需求选择适当的方法,以满足系统属性的要求。Help Help polyfunpolyfun可以得到更详细的内容。可以得到更详细的内容。y=interp1(xs,ys,x,y=interp1(xs,ys,x,methodmethod) ) 在有限样本点向量在有限样本点向量xsxs与与ysys中,插值产生向量

37、中,插值产生向量x x和和y y,所用方法,所用方法定义在定义在methodmethod中,有中,有4 4种选择:种选择:nearestnearest:执行速度最快,输出结果为直角转机:执行速度最快,输出结果为直角转机linearlinear:默许值,在样本点上斜率变化很大:默许值,在样本点上斜率变化很大splinespline:最花时间,但输出结果也最平滑:最花时间,但输出结果也最平滑cubiccubic:最占内存,输出结果与:最占内存,输出结果与splinespline差不多差不多数据处置练习题: 如:如:1 1、多项式、多项式x4-12x3+0 x2+25x+116x4-12x3+0 x

38、2+25x+116求多项式当求多项式当x=5x=5的的值,及多项式值,及多项式=0=0时的根。时的根。P=1 -12 0 25 116P=1 -12 0 25 116;poly2sym(P) poly2sym(P) polyval (P,5) %polyval (P,5) %求此多项式当求此多项式当x=5x=5时的值时的值r=roots(P)%r=roots(P)%求此多项式求此多项式=0=0时的根时的根2 2、a=1 2 3;b=1 2a=1 2 3;b=1 2求求a a与与b b的卷积的卷积c;c;求多项式求多项式c c除以除以多项式多项式a;a;求多项式求多项式a a的导数;求多项式的导数;求多项式a a当当x x为为2 2时的解。时的解。 c=conv(a,b)=1 4 7 6 c=conv(a,b)=1 4 7 6其中,其中,convconv指令可以嵌套运用,如指令可以嵌套运用,如conv(conv(a,b),c)conv(conv(a,b),c) q,r=deconv(c,b) q,r=deconv(c,b)得得q=1 2 3 q=1 2 3 商多项式商多项式 r=0 0 0 r=0 0 0 余多项式余多项式 p

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论