MATLAB作业
>>demo
The element type"name"must be terminated by the matching end-tag"".
org.xml.sax.SAXParseException:The element type"name"must be terminated by the matching end-tag"".
at org.apache.xerces.parsers.DOMParser.parse(DOMParser.java:235)
at org.apache.xerces.jaxp.DocumentBuilderImpl.parse(DocumentBuilderImpl.java:201)
at com.mathworks.mde.help.Demos.createNodes(Demos.java:416)
at com.mathworks.mde.help.Demos.access$900(Demos.java:46)
at com.mathworks.mde.help.Demos$LoadDemosThread.run(Demos.java:276)
at com.mathworks.jmi.NativeMatlab.dispatchMTRequests(NativeMatlab.java:250)
Could not parse G:\toolbox\ccslink\ccslink\demos.xml
demo：演示作用
>>A=magic(3)
A=
816
357
492
>>[sum(A),sum(A'),trace(A),trace(rot90(A))]
ans=
1515151515151515
>>t=[0:0.1:2*pi];y=t.*sin(2.*t);plot(t,y)
>>[x,y]=meshgrid(-5:0.1:5);
>>z=(x.^2+y).*exp(-x.^2-(y.^2)/2); >>surf(x,y,z),colorbar
>>who
Your variables are:
A ans t x y z
>>whos
Name Size Bytes Class
A3x372double array ans1x864double array t1x63504double array x101x10181608double array y101x10181608double array z101x10181608double array Grand total is30683elements using245464bytes
who和whos：显示当前工作空间的变量情况
>>which first
'first'not found
which：显示出某个MATLAB文件所在目录
>>workspace
workspace：打开工作空间的浏览器窗口
quit和exit：退出程序
>>path
MATLABPATH
G:\toolbox\matlab\general
G:\toolbox\matlab\ops
G:\toolbox\matlab\lang
G:\toolbox\matlab\elmat
G:\toolbox\matlab\elfun
G:\toolbox\matlab\specfun
G:\toolbox\matlab\datafun
G:\toolbox\matlab\polyfun
G:\toolbox\matlab\funfun
G:\toolbox\matlab\sparfun
G:\toolbox\matlab\scribe
G:\toolbox\matlab\graph2d
G:\toolbox\matlab\graph3d
G:\toolbox\matlab\specgraph
G:\toolbox\matlab\graphics
G:\toolbox\matlab\uitools
G:\toolbox\matlab\strfun
G:\toolbox\matlab\imagesci
G:\toolbox\matlab\iofun
G:\toolbox\matlab\audiovideo
G:\toolbox\matlab\timefun
G:\toolbox\matlab\datatypes
G:\toolbox\matlab\verctrl
G:\toolbox\matlab\codetools
G:\toolbox\matlab\helptools
G:\toolbox\matlab\winfun
G:\toolbox\matlab\demos
G:\toolbox\local
G:\toolbox\matlabxl\matlabxl
G:\toolbox\simulink\simulink
G:\toolbox\simulink\blocks
G:\toolbox\simulink\components
G:\toolbox\simulink\fixedandfloat
G:\toolbox\simulink\fixedandfloat\fxpdemos G:\toolbox\simulink\fixedandfloat\obsolete
G:\toolbox\simulink\simdemos
G:\toolbox\simulink\simdemos\aerospace
G:\toolbox\simulink\simdemos\automotive
G:\toolbox\simulink\simdemos\simfeatures
G:\toolbox\simulink\simdemos\simfeatures\mdlref G:\toolbox\simulink\simdemos\simgeneral
G:\toolbox\simulink\simdemos\simnew
G:\toolbox\simulink\dee
G:\toolbox\shared\dastudio
G:\toolbox\stateflow\stateflow
G:\toolbox\rtw\rtw
G:\toolbox\shared\hds
G:\toolbox\shared\timeseries
G:\toolbox\stateflow\sfdemos
G:\toolbox\rtw\rtwdemos
G:\toolbox\rtw\rtwdemos\rsimdemos
G:\toolbox\rtw\targets\asap2\asap2
G:\toolbox\rtw\targets\asap2\asap2\user
G:\toolbox\rtw\targets\rtwin\rtwin
G:\toolbox\simulink\accelerator
G:\toolbox\rtw\accel
G:\toolbox\aeroblks\aeroblks
G:\toolbox\aeroblks\aerodemos
G:\toolbox\aeroblks\aerodemos\texture
G:\toolbox\bioinfo\bioinfo
G:\toolbox\bioinfo\microarray
G:\toolbox\bioinfo\proteins
G:\toolbox\bioinfo\biomatrices
G:\toolbox\bioinfo\biodemos
G:\toolbox\rtw\targets\common\can\blocks
G:\toolbox\rtw\targets\common\configuration\resource G:\toolbox\rtw\targets\common\tgtcommon
G:\toolbox\rtw\targets\c166\c166
G:\toolbox\rtw\targets\c166\blocks
G:\toolbox\rtw\targets\c166\c166demos
G:\toolbox\ccslink\ccslink
G:\toolbox\ccslink\ccsblks
G:\toolbox\ccslink\ccsdemos
G:\toolbox\cdma\cdma
G:\toolbox\cdma\cdmamasks
G:\toolbox\cdma\cdmamex
G:\toolbox\cdma\cdmademos
G:\toolbox\combuilder\combuilder
G:\toolbox\comm\comm
G:\toolbox\comm\commdemos
G:\toolbox\comm\commdemos\commdocdemos
G:\toolbox\comm\commobsolete
G:\toolbox\commblks\commblks
G:\toolbox\commblks\commmasks
G:\toolbox\commblks\commmex
G:\toolbox\commblks\commblksdemos
G:\toolbox\commblks\commblksobsolete\v2p5
G:\toolbox\commblks\commblksobsolete\v2
G:\toolbox\commblks\commblksobsolete\v1p5
G:\toolbox\compiler
G:\toolbox\control\control
G:\toolbox\control\ctrlguis
G:\toolbox\control\ctrlobsolete
G:\toolbox\control\ctrlutil
G:\toolbox\control\ctrldemos
G:\toolbox\shared\controllib
G:\toolbox\curvefit\curvefit
G:\toolbox\curvefit\cftoolgui
G:\toolbox\shared\optimlib
G:\toolbox\daq\daq
G:\toolbox\daq\daqguis
G:\toolbox\daq\daqdemos
G:\toolbox\database\database
G:\toolbox\database\dbdemos
G:\toolbox\database\vqb
G:\toolbox\datafeed\datafeed
G:\toolbox\datafeed\dfgui
G:\toolbox\dials
G:\toolbox\dspblks\dspblks
G:\toolbox\dspblks\dspmasks
G:\toolbox\dspblks\dspmex
G:\toolbox\dspblks\dspdemos
G:\toolbox\rtw\targets\ecoder
G:\toolbox\rtw\targets\ecoder\ecoderdemos G:\toolbox\rtw\targets\mpt
G:\toolbox\rtw\targets\mpt\mpt
G:\toolbox\rtw\targets\mpt\user_specific G:\toolbox\exlink
G:\toolbox\symbolic\extended
G:\toolbox\filterdesign\filterdesign
G:\toolbox\filterdesign\quantization
G:\toolbox\filterdesign\filtdesdemos
G:\toolbox\finance\finance
G:\toolbox\finance\calendar
G:\toolbox\finance\findemos
G:\toolbox\finance\finsupport
G:\toolbox\finderiv\finderiv
G:\toolbox\finfixed\finfixed
G:\toolbox\fixedpoint\fixedpoint
G:\toolbox\fixedpoint\fidemos
G:\toolbox\fixedpoint\fimex
G:\toolbox\fixpoint
G:\toolbox\ftseries\ftseries
G:\toolbox\ftseries\ftsdemos
G:\toolbox\ftseries\ftsdata
G:\toolbox\ftseries\ftstutorials
G:\toolbox\fuzzy\fuzzy
G:\toolbox\fuzzy\fuzdemos
G:\toolbox\gads
G:\toolbox\gads\gads
G:\toolbox\gads\gadsdemos
G:\toolbox\garch\garch
G:\toolbox\garch\garchdemos
G:\toolbox\rtw\targets\hc12\hc12
G:\toolbox\rtw\targets\hc12\blocks G:\toolbox\rtw\targets\hc12\codewarrior G:\toolbox\rtw\targets\hc12\hc12demos G:\toolbox\hdlfilter\hdlfilter
G:\toolbox\hdlfilter\hdlfiltdemos
G:\toolbox\shared\hdlshared
G:\toolbox\ident\ident
G:\toolbox\ident\idobsolete
G:\toolbox\ident\idguis
G:\toolbox\ident\idutils
G:\toolbox\ident\iddemos
G:\toolbox\ident\idhelp
G:\toolbox\images\images
G:\toolbox\images\imdemos
G:\toolbox\imaq\imaq
G:\toolbox\imaq\imaqdemos
G:\toolbox\instrument\instrument
G:\toolbox\instrument\instrumentdemos G:\toolbox\lmi\lmictrl
G:\toolbox\lmi\lmilab
G:\toolbox\map\map
G:\toolbox\map\mapdemos
G:\toolbox\map\mapdisp
G:\toolbox\map\mapformats
G:\toolbox\map\mapproj
G:\toolbox\mbc\mbc
G:\toolbox\mbc\mbcdata
G:\toolbox\mbc\mbcdesign
G:\toolbox\mbc\mbcexpr
G:\toolbox\mbc\mbcguitools
G:\toolbox\mbc\mbclayouts
G:\toolbox\mbc\mbcmodels
G:\toolbox\mbc\mbcsimulink
G:\toolbox\mbc\mbctools
G:\toolbox\mbc\mbcview
G:\toolbox\physmod\mech\mech
G:\toolbox\physmod\mech\mechdemos
G:\toolbox\physmod\pmimport\pmimport
G:\toolbox\modelsim\modelsim
G:\toolbox\modelsim\modelsimdemos
G:\toolbox\mpc\mpc
G:\toolbox\mpc\mpcdemos
G:\toolbox\mpc\mpcguis
G:\toolbox\mpc\mpcobsolete
G:\toolbox\mpc\mpcutils
G:\toolbox\shared\slcontrollib
G:\toolbox\rtw\targets\mpc555dk
G:\toolbox\rtw\targets\mpc555dk\common\configuration G:\toolbox\rtw\targets\mpc555dk\mpc555demos
G:\toolbox\rtw\targets\mpc555dk\mpc555dk
G:\toolbox\rtw\targets\mpc555dk\pil
G:\toolbox\rtw\targets\mpc555dk\rt\blockset\mfiles G:\toolbox\rtw\targets\mpc555dk\rt\blockset
G:\toolbox\mutools\commands
G:\toolbox\mutools\subs
G:\toolbox\nnet\nnet
G:\toolbox\nnet\nnutils
G:\toolbox\nnet\nncontrol
G:\toolbox\nnet\nndemos
G:\toolbox\nnet\nnobsolete
G:\toolbox\opc\opc
G:\toolbox\opc\opcgui
G:\toolbox\opc\opcdemos
G:\toolbox\optim
G:\toolbox\rtw\targets\osek\osek
G:\toolbox\rtw\targets\osek\osekdemos
G:\toolbox\rtw\targets\osek\blocks
G:\toolbox\rtw\targets\osek\osekworks
G:\toolbox\rtw\targets\osek\proosek
G:\toolbox\pde
G:\toolbox\physmod\powersys\powersys
G:\toolbox\physmod\powersys\powerdemo
G:\toolbox\simulink\reqmgt
G:\toolbox\simulink\reqmgt\rmidemos
G:\toolbox\rf\rf
G:\toolbox\rf\rfdemos
G:\toolbox\rf\rftool
G:\toolbox\rfblks\rfblks
G:\toolbox\rfblks\rfblksmasks
G:\toolbox\rfblks\rfblksmex
G:\toolbox\rfblks\rfblksdemos
G:\toolbox\robust
G:\toolbox\rptgen\rptgen
G:\toolbox\rptgen\rptgendemos
G:\toolbox\rptgenext\rptgenext
G:\toolbox\rptgenext\rptgenextdemos
G:\toolbox\signal\signal
G:\toolbox\signal\sigtools
G:\toolbox\signal\sptoolgui
G:\toolbox\signal\sigdemos
G:\toolbox\slcontrol\slcontrol
G:\toolbox\slcontrol\slctrlguis
G:\toolbox\slcontrol\slctrlutil
G:\toolbox\slcontrol\slctrldemos
G:\toolbox\slestim\slestim
G:\toolbox\slestim\slestmex
G:\toolbox\slestim\slestdemos
G:\toolbox\slestim\slestguis
G:\toolbox\sloptim\sloptim
G:\toolbox\sloptim\sloptguis
G:\toolbox\sloptim\sloptdemos
G:\toolbox\sloptim\sloptobsolete
G:\toolbox\simulink\slvnv
G:\toolbox\simulink\simcoverage
G:\toolbox\simulink\simcoverage\simcovdemos G:\toolbox\splines
G:\toolbox\stats
G:\toolbox\symbolic
G:\toolbox\rtw\targets\tic2000\tic2000
G:\toolbox\rtw\targets\tic2000\tic2000blks
G:\toolbox\rtw\targets\tic2000\tic2000demos G:\toolbox\shared\etargets\etargets
G:\toolbox\rtw\targets\tic6000\tic6000
G:\toolbox\rtw\targets\tic6000\tic6000blks
G:\toolbox\rtw\targets\tic6000\tic6000demos G:\toolbox\shared\etargets\rtdxblks
G:\toolbox\vr\vr
G:\toolbox\vr\vrdemos
G:\toolbox\wavelet\wavelet
G:\toolbox\wavelet\wavedemo
G:\toolbox\webserver\webserver
G:\toolbox\webserver\wsdemos
G:\toolbox\rtw\targets\xpc\xpc
G:\toolbox\rtw\targets\xpc\target\build\xpcblocks
G:\toolbox\rtw\targets\xpc\xpcdemos
G:\toolbox\rtw\targets\xpc\xpc\xpcmngr
G:\toolbox\rtw\targets\xpc\target\kernel\embedded
G:\work
path：显示MATLAB包含的路径
>>cd
G:\work
cd：显示当前路径
>>addpath a
>In path at115
In addpath at95
addpath：将一个新的字符串表示的子目录名添加到MATLAB包含的路径下rmpath：从当前的MATLAB包含的某个路径下清除某个字符串表示的子目录>>format long
>>pi
ans=
3.14159265358979
format long：小数点后精确到7位（单精度）或14位（双精度）有效数字>>format short e
>>pi
ans=
3.1416e+000
format short e：5位科学记数表示
>>format long e
>>pi
ans=
3.141592653589793e+000
format long e：16位科学记数表示
>>format short g
>>pi
ans=
3.1416
format short g：从format short和format short e中自动选择最佳计数方式>>format long g
>>pi
ans=
3.14159265358979
format long g：从format long和format long e中自动选择最佳计数方式>>format hex
>>pi
ans=
400921fb54442d18
format hex：十六进制表示
>>format rat
>>pi
ans=
355/113
format rat：近似有理数表示
clc：清除命令窗口
cla：清除当前坐标系下的内容
clf：清除图形窗口
close：关闭所有的命令窗口
copyfile：复制文件
>>pwd
ans=
G:\work
pwd：显示当前的工作目录名称
>>mkdir a
>>dir
...a
mkdir：建立目录
dir：获取当前目录并保存在列表中
>>delete a
Warning:'G:\work\a'is a https://www.360docs.net/doc/713824121.html,e rmdir to delete directories. delete：删除文件
MATLAB编程作业
《Matlab 编程训练》 作业 专 业 学生姓名 班级 学 号 指导教师 完成日期
实训一 MATLAB 语言介绍和数值计算 1．先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存变量。 12 2sin851z e =+ . 2. 已知 1234413134787,2033657327A B --???? ????==???? ????-???? ，求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） A+6*B:
A-B+I: (2)A*B和A.*B A*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] c=A*B 结果： A.*B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=A.*B 结果：
(3)A^3和A.^3 A^3程序： A=[12 34 -4;34 7 87;3 65 7] E=A^3 结果： A．^3程序： A=[12 34 -4;34 7 87;3 65 7] C=A.^3 (4)A/B及B\A A/B程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] C=A/B 结果：
B\A程序： A=[12 34 -4;34 7 87;3 65 7] B=[1 3 -1;2 0 3;3 -2 7] D=B\A 结果： (5)将矩阵C=B\A的右下角2*2子矩阵赋给D, 并（3）保存变量（mat文件）程序： A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; C=B*inv(A); D=C(2:3,2:3) 结果：
(完整版)MATLAB常用函数大全
一、MATLAB常用的基本数学函数 abs(x)：纯量的绝对值或向量的长度 angle(z)：复数z的相角(Phase angle) sqrt(x)：开平方 real(z)：复数z的实部 imag(z)：复数z的虚部 conj(z)：复数z的共轭复数 round(x)：四舍五入至最近整数 fix(x)：无论正负，舍去小数至最近整数 floor(x)：地板函数，即舍去正小数至最近整数ceil(x)：天花板函数，即加入正小数至最近整数rat(x)：将实数x化为分数表示 rats(x)：将实数x化为多项分数展开 sign(x)：符号函数(Signum function)。 当x<0时，sign(x)=-1； 当x=0时，sign(x)=0; 当x>0时，sign(x)=1。 rem(x,y)：求x除以y的馀数 gcd(x,y)：整数x和y的最大公因数 lcm(x,y)：整数x和y的最小公倍数 exp(x)：自然指数 pow2(x)：2的指数 log(x)：以e为底的对数，即自然对数或 log2(x)：以2为底的对数 log10(x)：以10为底的对数 二、MATLAB常用的三角函数 sin(x)：正弦函数 cos(x)：余弦函数
tan(x)：正切函数 asin(x)：反正弦函数 acos(x)：反馀弦函数 atan(x)：反正切函数 atan2(x,y)：四象限的反正切函数 sinh(x)：超越正弦函数 cosh(x)：超越馀弦函数 tanh(x)：超越正切函数 asinh(x)：反超越正弦函数 acosh(x)：反超越馀弦函数 atanh(x)：反超越正切函数 三、适用於向量的常用函数有： min(x): 向量x的元素的最小值 max(x): 向量x的元素的最大值 mean(x): 向量x的元素的平均值 median(x): 向量x的元素的中位数 std(x): 向量x的元素的标准差 diff(x): 向量x的相邻元素的差 sort(x): 对向量x的元素进行排序（Sorting）length(x): 向量x的元素个数 norm(x): 向量x的欧氏（Euclidean）长度sum(x): 向量x的元素总和 prod(x): 向量x的元素总乘积 cumsum(x): 向量x的累计元素总和cumprod(x): 向量x的累计元素总乘积 dot(x, y): 向量x和y的内积 cross(x, y): 向量x和y的外积 四、MATLAB的永久常数
matlab期末大作业
电气学科大类 Modern Control Systems Analysis and Design Using Matlab and Simulink Title: Automobile Velocity Control Name: 巫宇智 Student ID: U200811997 Class：电气0811
电气0811 巫宇智 Catalogue Preface (3) The Design Introduction (4) Relative Knowledge (5) Design and Analyze (6) Compare and Conclusion (19) After design (20) Appendix (22) Reference (22)
Automobile Velocity Control 1.Preface: With the high pace of human civilization development, the car has been a common tools for people. However, some problems also arise in such tendency. Among many problems, the velocity control seems to a significant challenge. In a automated highway system, using the velocity control system to maintain the speed of the car can effectively reduce the potential danger of driving a car and also will bring much convenience to drivers. This article aims at the discussion about velocity control system and the compensator to ameliorate the preference of the plant, thus meets the complicated demands from people. The discussion is based on the simulation of MATLAB. Key word: PI controller, root locus
MATLAB基础训练作业(含答案)
实验一 MATLAB 工作环境熟悉及简单命令的执行 一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。 二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单 （Menu ）、工具条 （Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。 三、实验步骤： 1、启动MATLAB ，熟悉MATLAB 的桌面。 2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。 （1）（365-52?2-70）÷3 （2）>>area=pi*2.5^2 （3）已知x=3，y=4，在MATLAB 中求z ： ()2 3 2 y x y x z -= （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。 m1=????? ???? ???115 14 4 12679810115 133216 执行以下命令 >>m1( 2 , 3 ) >>m1( 11 ) >>m1( : , 3 ) >>m1( 2 : 3 , 1 : 3 ) >>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) （5）执行命令>>help abs 查看函数abs 的用法及用途，计算abs( 3 + 4i ) （6）执行命令 >>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y) （6）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。 四、思考题 1、以下变量名是否合法？为什么？ （1）x2 （2）3col （3）_row （4）for 2、求以下变量的值，并在MATLAB 中验证。 （1）a = 1 : 2 : 5 ; （2）b = [ a' a' a' ] ; （3）c = a + b ( 2 , : )
MATLAB程序设计作业
Matlab程序设计 班级 姓名 学号
《MATLAB程序设计》作业 1、考虑如下x-y 一组实验数据： x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot的原始数据、一次拟合曲线和三次拟合曲线，给出MATLAB代码和运行结果。 代码如下： x=[1,2,3,4,5,6,7,8,9,10]; y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]; plot(x,y); title('原始数据'); p=polyfit(x,y,1); q=polyval(p,x); figure,plot(x,q); title('一次拟合'); p=polyfit(x,y,2); q=polyval(p,x); figure,plot(x,q); title('二次拟合'); 运行结果如下：
1 2 3 4 5 6 7 8 9 10 12 3 4 5 6 7 8 原始数据 123 456789 102 2.5 3 3.54 4.5 55.56 6.57一次拟合 123456789 101 2 3 4 5 6 7 二次拟合 2、在[0，3π]区间，绘制y=sin(x)曲线（要求消去负半波，即(π，2π)区间内的函数值置零），求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 代码如下： clear clc x=(0:0.01:3*pi); y=sin(x); plot(x,y); y1=(y>=0).*y; figure,plot(x,y1);
matlab 常用函数汇总
matlab 常用函数汇总 编程2008-07-10 21:45:20 阅读46 评论0 字号：大中小订阅matlab常用函数 图形注释 Title 图形标题 Xlabel X轴标记 Ylabel Y轴标记 Text 文本注释 Gtext 用鼠标放置文本 Grid 网格线 MATLAB编程语言 Function 增加新的函数 Eval 执行由MA TLAB表达式构成的字串 Feval 执行由字串指定的函数 Global 定义全局变量 程序控制流 If 条件执行语句 Else 与if命令配合使用 Elseif 与if命令配合使用 End For,while和if语句的结束 For 重复执行指定次数（循环） While 重复执行不定次数（循环） Break 终止循环的执行 Return 返回引用的函数 Error 显示信息并终止函数的执行 交互输入 Input 提示用户输入 Keyboard 像底稿文件一样使用键盘输入 Menu 产生由用户输入选择的菜单 Pause 等待用户响应 Uimenu 建立用户界面菜单 Uicontrol 建立用户界面控制 一般字符串函数 Strings MATLAB中有关字符串函数的说明 Abs 变字符串为数值 Setstr 变数值为字符串 Isstr 当变量为字符串时其值为真 Blanks 空串 Deblank 删除尾部的空串 Str2mat 从各个字符串中形成文本矩阵 Eval 执行由MA TLAB表达式组成的串 字符串比较 Strcmp , , , 比较字符串 Findstr 在一字符串中查找另一个子串
Upper 变字符串为大写 Lower 变字符串为小写 Isletter 当变量为字母时，其值为真 Isspace 当变量为空白字符时，其值为真 字符串与数值之间变换 Num2str 变数值为字符串 Int2str 变整数为字符串 Str2num 变字符串为数值 Sprintf 变数值为格式控制下的字符串 Sscanf 变字符串为格式控制下的数值 十进制与十六进制数之间变换 Hex2num 变十六进制为IEEE标准下的浮点数Hex2dec 变十六制数为十进制数 Dec2hex 变十进制数为十六进制数 建模 Append 追加系统动态特性 Augstate 变量状态作为输出 Blkbuild 从方框图中构造状态空间系统Cloop 系统的闭环 Connect 方框图建模 Conv 两个多项式的卷积 Destim 从增益矩阵中形成离散状态估计器Dreg 从增益矩阵中形成离散控制器和估计器Drmodel 产生随机离散模型 Estim 从增益矩阵中形成连续状态估计器Feedback 反馈系统连接 Ord2 产生二阶系统的A、B、C、D Pade 时延的Pade近似 Parallel 并行系统连接 Reg 从增益矩阵中形成连续控制器和估计器Rmodel 产生随机连续模型 Series 串行系统连接 Ssdelete 从模型中删除输入、输出或状态ssselect 从大系统中选择子系统 模型变换 C2d 变连续系统为离散系统 C2dm 利用指定方法变连续为离散系统 C2dt 带一延时变连续为离散系统 D2c 变离散为连续系统 D2cm 利用指定方法变离散为连续系统 Poly 变根值表示为多项式表示 Residue 部分分式展开 Ss2tf 变状态空间表示为传递函数表示 Ss2zp 变状态空间表示为零极点表示
春MATLAB仿真期末大作业
MATLAB仿真 期末大作业 姓名：班级：学号：指导教师：
2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。
仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。
MATLAB第一章作业答案
第一章 M A T L A B 概况与基本操作 1.选择题（每题2分，共20分）： (1)最初的MATLAB 核心程序是采用D 语言编写的。 (2)即将于2011年9月发布的MATLAB 新版本的编号为C 。 2011Ra 2011Rb R2011b R2011a (3)在默认设置中，MATLAB 中的注释语句显示的颜色是B 。 A.黑色 B.绿色 C.红色 D.蓝色 (4)如果要以科学计数法显示15位有效数字，使用的命令是A 。 long e long long g long d (5)在命令窗口新建变量a 、b ，如果只查看变量a 的详细信息，使用的命令为A 。 a a (6)如果要清除工作空间的所有变量，使用的命令为C 。 all C.两者都可 D.两者都不可 (7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上B 。 A.冒号 B.分号 C.空格 D.逗号 (8)如果要重新执行以前输入的命令，可以使用D 键。 A.下箭头↓ B.右箭头→ C.左箭头← D.上箭头↑ (9)如果要查询函数det 的功能和用法，并显示在命令窗口，应使用命令C 。 D.三者均可 (10)如果要启动Notebook 文档，下列D 操作是可行的。 A.在命令窗口输入notebook 命令 B.在命令窗口输入notebook filename 命令 C.在Word 中启动M-book 文档 D.三者均可 2.填空题（每空1分，共20分）： (1)MATLAB 是matrix 和laboratory 两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler 和Jack Little 。 (2)在MATLAB 的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为褐红色，注释显示的字体为绿色，错误信息显示的字体为红色。 (3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format loose ，各行之间不添加空行的命令为format compact 。 (4)在MATLAB 中，各种标点符号的作用是不同的。例如，空格的作用是分隔数组每行各个元素，逗号的作用是分隔数组每行各个元素或函数的各个输入参数，分号的作用是作为不显示命令结果的命令行的结尾或分隔数组各列，冒号的作用是生成一维数组或表示数组全部元素，百分号的作用是引导一行注释，…的作用是连接相邻两行，感叹号的作用是调用操作系统命令。 3.程序设计题（每题10分，共40分） (1)以25m/s 的初速度向正上方投球（g=s 2 ），计算到达最高点的时间tp 以及球从出发点到最高点的距离hp 。 解：根据物理学知识，物体上抛运动的速度与经过的时间之间的关系为0p p v v gt =-，因此所需要的时间为0p p v v t g -= 。而到达最高点时的速度0p v =，因此可根据此公式求出tp ： v0=25;g=;vp=0; tp=(v0-vp)/g tp =
matlab程序设计第三章课后习题答案
1. p138 第6题在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。 >> t=0:0.01:pi; >> x1=t; >> y1=2*x1-0.5; >> x2=sin(3*t).*cos(t); >> y2=sin(3*t).*sin(t); >> plot(x1,y1,'r-',x2,y2,'g-') >> axis([-1,2,-1.5,1]) >> hold on >> s=solve('y=2*x-0.5','x=sin(3*t)*cos(t)','y=sin(3*t)*sin(t)'); >> plot(double(s.x),double(s.y),'*'); 截图：
p366 第4题绘制极坐标曲线，并分析对曲线形状的影响。 function [ output_args ] = Untitled2( input_args ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here theta=0:0.01:2*pi; a=input('请输入a的值：'); b=input('请输入b的值：'); n=input('请输入n的值：'); rho=a*sin(b+n*theta); polar(theta,rho,'k'); end 下面以a=1，b=1，n=1的极坐标图形为基础来分析a、b、n的影响。
对a的值进行改变：对比发现a只影响半径值的整倍变化 对b的值进行改变：对比发现b的值使这个圆转换了一定的角度
对n的值进行改变：对比发现当n>=2时有如下规律 1、当n为整数时，图形变为2n个花瓣状的图形 2、当n为奇数时，图形变为n个花瓣状的图形 分别让n为2、3、4、5
matlab 常用函数(1)
A axis() axis([xmin xmax ymin ymax]) sets the limits for the x- and y-axis of the current axes. axis([xmin xmax ymin ymax zmin zmax cmin cmax]) sets the x-, y-, and z-axis limits and the color scaling limits (see caxis) of the current axes. axis equal sets the aspect ratio so that the data units are the same in every direction. The aspect ratio of the x-, y-, and z-axis is adjusted automatically according to the range of data units in the x, y, and z directions C clf Clear current figure window G grid off/on The grid function turns the current axes' grid lines on and off. H hold on/off ●The hold function determines whether new graphics objects are added to the graph or replace objects in the graph. ●hold on retains the current plot and certain axes properties so that subsequent graphing commands add to the existing graph. ●hold off resets axes properties to their defaults before drawing new plots. hold off is the default
期末大作业报告
期末大作业报告 课程名称：数字图像处理 设计题目：车牌识别 学院：信息工程与自动化学院 专业：计算机科学与技术 年级：xxxxx 学生姓名：xxxxxxx（学号xxxxxxxxxxxxx） 指导教师：xxxx 日期：20XX.6.10 教务处制 车牌识别 摘要:数字图像处理技术是20世纪60年代发展起来的一门新兴学科，随着图像处理理论和方法的进一步完善，使得数字图像处理技术在各个领域得到了广泛应用，并显示出广阔的应用前景。MATLAB既是一种直观、高效的计算机语言,同时又是一个科学计算平台。它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。根据它提供的500多个数学和工程函数，工程技术人员和科学工作者可以在它的集成环境中交互或编程以完成各自的计算。MATLAB中集成了功能强大的图像处理工具箱。由于MA TLAB语言的语法特征与C语言极为相似，而且更加简单，更加符合科技人员对数学表达式的书写格式，而且这种语言可移植性好、可扩展性强，再加上其中有丰富的图像处理函数，所以MA TLAB在图像处理的应用中具有很大的优势。车牌识别技术是智能交通系统的重要组成部分，在近年来得到了很大的发展。本文从预处理、边缘检测、车牌定位、字符分割、字符识别五个方面，具体介绍了车牌自动识别的原理。并用MATLAB软件编程来实现每一个部分,最后识别出汽车牌照。 关键词:车牌识别、数字图像处理、MATLAB
一、设计原理 车辆牌照识别系统的基本工作原理为：将摄像头拍摄到的包含车辆牌照的图像通过视频卡输入到计算机中进行预处理，再由检索模块对牌照进行搜索、检测、定位，并分割出包含牌照字符的矩形区域，然后对牌照字符进行二值化并将其分割为单个字符，然后输入JPEG或BMP 格式的数字，输出则为车牌号码的数字。牌照自动识别是一项利用车辆的动态视频或静态图像进行牌照号码、牌照颜色自动识别的模式识别技术。其硬件基础一般包括触发设备、摄像设备、照明设备、图像采集设备、识别车牌号码的处理机等，其软件核心包括车牌定位算法、车牌字符分割算法和光学字符识别算法等。某些牌照识别系统还具有通过视频图像判断车辆驶入视野的功能称之为视频车辆检测。一个完整的牌照识别系统应包括车辆检测、图像采集、牌照识别等几部分。当车辆检测部分检测到车辆到达时触发图像采集单元，采集当前的视频图像。牌照识别单元对图像进行处理，定位出牌照位置，再将牌照中的字符分割出来进行识别，然后组成牌照号码输出。 二、设计步骤 1. 提出总体设计方案: （1）车牌图像预处理方法 因为车牌图像都是在室外拍摄的，所以不可避免地会受到光照、气候等因素的影响，而且拍摄者的手部抖动与车辆的移动会造成图像的模糊。要去除这些干扰就得先对车牌图像进行预处理。由于当前数码相机的像素较高，原始图像的数据一般比较大，输入的彩色图像包含大量颜色信息，会占用较多的存储空间，且处理时也会降低系统的执行速度。因此对图像进行识别等处理时，常将彩色图像转换为灰度图像，以加快处理速度。对图像进行灰度化处理后常用的方法是图像二值化、去除背景图像、增强处理、边缘检测、滤波等处理等。
MATLAB作业4参考答案
MATLAB 作业四参考答案 1、 用2sin(103)y t =+在（0，3）区间内生成一组较稀疏的数据，并用一维数据插值的方法对给出的数据进行曲线拟合，并将结果与理论曲线相比较。 【求解】类似于上面的例子，可以用几乎一致的语句得出样本数据和插值效果。 >> t=0:0.2:3; y=sin(10*t.^2+3); plot(t,y,'o') ezplot('sin(10*t^2+3)',[0,3]); hold on x1=0:0.001:3; y1=interp1(t,y,x1,'spline'); plot(x1,y1) 由于曲线本身变换太大，所以在目前选定的样本点下是不可能得出理想插值效果的，因为样 本数据提供的信息量不够。为了得到好的插值效果，必须增大样本数据的信息量，对本例来 说，必须在快变化区域减小样本点的步长。 >> hold off t=[0:0.1:1,1.1:0.04:3]; y=sin(10*t.^2+3); plot(t,y,'o') ezplot('sin(10*t^2+3)',[0,3]); hold on x1=0:0.001:3; y1=interp1(t,y,x1,'spline'); plot(x1,y1) 2、 用2422 3 1(,)sin()3x y f x y e xy x y x y --= ++原型函数生成一组网络数据或随机数据，分别拟合出曲面，并和原曲面进行比较。 【求解】由下面的语句可以直接生成一组网格数据，用下面语句还可以还绘制出给定样本点是三维表面图。 >> [x,y]=meshgrid(0.2:0.2:2); z=exp(-x.^2-y.^4).*sin(x.*y.^2+x.^2.*y)./(3*x.^3+y); surf(x,y,z) 选择新的密集网格，则可以通过二元插值得出插值曲面。对比插值结果和新网格下的函数值精确解，则可以绘制出绝对插值误差曲面。由插值结果可见精度是令人满意的。 >> [x1,y1]=meshgrid(0.2:0.02:2); z1=interp2(x,y,z,x1,y1,'spline'); surf(x1,y1,z1) >> z0=exp(-x1.^2-y1.^4).*sin(x1.*y1.^2+x1.^2.*y1)./(3*x1.^3+y1); surf(x1,y1,abs(z1-z0)) 现在假设已知的样本点不是网格形式分布的，而是随机分布的，则可以用下面语句生成样本点，得出分布的二维、三维示意图。 >> x=0.2+1.8*rand(400,1); y=0.2+1.8*rand(400,1); % 仍生成(0.2,2) 区间的均匀分布随机数 z=exp(-x.^2-y.^4).*sin(x.*y.^2+x.^2.*y)./(3*x.^3+y); plot(x,y,'x') figure, plot3(x,y,z,'x')
实验二--MATLAB程序的设计(含实验报告)
实验二 MATLAB 程序设计 一、 实验目的 1.掌握利用if 语句实现选择结构的方法。 2.掌握利用switch 语句实现多分支选择结构的方法。 3.掌握利用for 语句实现循环结构的方法。 4.掌握利用while 语句实现循环结构的方法。 5.掌握MATLAB 函数的编写及调试方法。 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 M 文件的编写： 启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正 三、 实验容 1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。并输入几组典型值加以检验。 （提示：提示输入使用input 函数） 2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。 要求：（1）用switch 语句实现。 （2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 （提示：注意单元矩阵的用法） 3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1
《科学计算与MATLAB》期末大作业
杭州电子科技大学信息工程学院《科学计算与MATLAB》期末大作业
给出程序、图、作业分析，程序需加注释。 1. 试编写名为fun.m 的MATLAB 函数，用以计算下述的值： ?? ? ??-<->=t t n t t t n t f 的)4/sin()(si 对所有)4/sin(其他情况)sin(的)4/sin()(si 对所有)4/sin()(ππππ 绘制t 关于函数f(t)的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。 function y=fun()%定义函数 % t=-6*pi:pi/10:6*pi; %定义变量范围 y = (sin(pi/4)).*(sin(t)>sin(pi/4))+(sin(-pi/4)).*(sin(t)
2.解以下线性方程组 ??? ??=+=++=--3 530 42231 321321x x x x x x x x A=[2 -1 -1;1 1 4;3 0 5];%输入矩阵 B=[2;0;3]; %输入矩阵 X = A\B %计算结果 3.已知矩阵? ? ??? ???? ???=44434241 3433323124232221 14131211A 求： (1)A(2:3,2:3) (2)A(:,1:2) (3)A(2:3,[1,3]) (4)[A,[ones(2,2)；eye(2)]]
A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];%输入矩阵A(2:3,2:3) %输出矩阵 A(:,1:2) %输出矩阵 A(2:3,[1,3]) %输出矩阵 [A,[ones(2,2);eye(2)]] %输出矩阵
Matlab编程与应用习题和一些参考答案
Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设???? ??????------=81272956313841A ，??????????-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ?? ????++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 （mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g)
MATLAB常用函数
数字信号处理与MATLAB 实现 1． n1=[ns:nf]； x1=[zeros(1,n0-ns),1,zeros (1,nf-n0)]； %单位抽样序列的产生 2. subplot(2,2,4) 画2行2列的第4个图 3. stem(n,x) %输出离散序列，(plot 连续) 4. 编写子程序可调用 4.1 单位抽样序列)(0n n -δ生成函数impseq.m [x,m]=impseq(n0,ns,nf); %序列的起点为ns ，终点为nf ，在n=n0点处生成一个单位脉冲 n=[-5:5];x1=3*impseq(2,-5,5)-impseq(-4,-5,5) x1 = 0 -1 0 0 0 0 0 3 0 0 0 n=[-5:5];x1=3*impseq(2,-4,5)-impseq(-4,-5,4) %起点到终点长度要一致 x1 = 0 -1 0 0 0 0 3 0 0 0 4.2 单位阶跃序列)(0n n u -生成函数stepseq.m [x,n]=stepseq(no,ns,nf) %序列的起点为ns ，终点为nf ，在n=n0点处生成一个单位阶跃 4.3 两个信号相加的生成函数sigadd.m [y,n]=sigadd(x1,n1,x2,n2) 4.4 两个信号相乘的生成函数sigmult.m [y,n]=sigmult(x1,n1,x2,n2) 4.5 序列移位y(n)=x(n-n0)的生成函数sigshift.m [y,n]=sigshift(x,m,n0) 4.6 序列翻褶y(n)=x(-n)的生成函数sigfold.m [y,n]=sigfold(x,n) 4.7 evenodd.m 函数可以将任一给定的序列x(n)分解为xe(n)和xo(n)两部分 [xe,xo,m]=evenodd(x,n) 4.8 序列从负值开始的卷积conv_m, conv 默认从0开始 function [y,ny]=conv_m(x,nx,h,nh) 有{x(n):nx1≤n ≤nx2},{h(n):nh1≤n ≤nh2}, 卷积结果序列为 {y(n):nx1+nh1≤n ≤nx2+nh2} 例. 设1132)(-++=z z z X ,1225342)(-+++=z z z z X ,求)()()(21z X z X z Y += 程序： x1=[1,2,3];n1=-1:1; x2=[2,4,3,5];n2=-2:1; [y,n]=conv_m(x1,n1,x2,n2)
MATLAB期末大作业模板
MATLAB应用技术 期末大作业 专业： 姓名： 学号： 分数
一、在一个图形窗口中以子图形式同时绘制正弦、余弦、正切、余切曲线。请写下完整代码，展示图形结果。（请标注题图和坐标轴，用不同颜色和不同线型分别绘制以上曲线）。（15分） 二、某公司员工的工资计算方法如下。 （1）工作时数超过120小时者，超过部分加发15%。 （2）工作时数低于60小时者，扣发700元。 （3）其余按每小时84元发。 根据员工的工时数，计算应发工资。请写下完整的程序代码，并任意输入一工时数（使用input 函数），将结果展示（使用disp 函数）利用该代码进行计算工资，请写下计算结果。（15分） 三、编写一个函数文件，使其能够产生如下的分段函数： ?? ? ??≥<<≤-=66225.0,25.05.15.0)(x x x x x x f 请编写完整的函数文件（保存函数文件名为hanshu.m ），并编写脚本文件代码，任意输入x 值（使用input 函数），在脚本文件中调用函数文件求)(x f ，展示结果（使用disp 函数），请写下计算结果。（15分） 四、将5个学生的6门功课的成绩存入矩阵P 中，进行如下处理： （1）分别求每门课的最高分、最低分及相应学生的序号。 （2）分别求每门课的平均分和标准差。 （3）5门课总分的最高分、最低分及相应学生序号。 （4）将5门课总分按从大到小顺序存入score 中，相应学生序号存入num 。 请将各小题的运行代码完整写下来，并写下运行结果。（20分） 五、请利用所学的MATLAB 知识，自主设计一个图形用户界面，请完整记录它的设计过程，需提供文字、代码和图片，以充分说明设计的图形用户界面可实现
MATLAB作业答案
题4.1 控制系统结构如图4.A 所示： 图4.A 习题 （1） 利用MA TLAB 对以上单位负反馈控制系统建立传递函数模型； （2） 将第一问中求得的传递函数模型转化为零极点增益形式和状态空间形式。 解：（1）MA TLAB 程序代码如下： z=[-1];p=[-1,-3];k=2 [num,den]=zp2tf(z,p,k) g_tf=tf(num,den) 程序运行结果如下： Transfer function: 2 s + 2 ------------- s^2 + 4 s + 3 由计算结果可知，系统的传递函数模型为： G （s ）= 3 4222 +++s s s (2) MA TLAB 程序代码如下： num=[2,2];den=[1,4,3] [z,p,k]=tf2zp(num,den) [A,B,C,D]=tf2ss(num,den) g_zpk=zpk(z,p,k) g_ss=ss(A,B,C,D) 程序运行结果如下： Zero/pole/gain: 2 (s+1) ----------- (s+3) (s+1) a = x1 x2 x1 -4 -3 x2 1 0 b = u1 x1 1 x2 0
c = x1 x2 y1 2 2 d = u1 y1 0 由计算结果可知，系统的零极点增益模和状态空间模型分别为： G （s ）=)1)(3()1(2+++s s s 和[]? ????=??????+??????-=x y u x x 2 2010 13- 4. 题5.1 设单位负反馈控制系统的开环传递函数为G （s ）= ) 177(2 ++s s s K 。 （1） 试绘制K=10、100时闭环系统的阶跃响应曲线，并计算稳态误差、上升时间、超调 量和过渡过程时间； （2） 绘制K=1000时闭环系统的阶跃响应曲线，与K=10、100所得结果相比较，分析增 益系数与系统稳定性的关系； （3） 利用roots 命令，确定使系统稳定时K 的取值范围。 解：（1） 题 6.1 已知单位负反馈控制系统的前向传递函数分别为G （s ）= ) 4)(2()1(2 +++s s s s K 、 ) 164)(1() 1(2 ++-+s s s s s K 和 ) 15)(7)(5)(3() 8(2 +++++s s s s s s K ，试利用MA TLAB 分别绘制各系 统的根轨迹图。 解：（1）G （s ）= ) 4)(2()1(2 +++s s s s K MA TLAB 程序代码如下： num=[1,1];den=conv([1,0],conv([1,0],[1,6,8])); sys=tf(num,den)
matlab程序设计作业
Matlab程序设计作业 姓名: 学号: 专业:
? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线，给出 原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10
7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合 三次拟合
x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在［0, 3n区间，绘制y二Sin(x)曲线(要求消去负半波，即(n 2n)区间内的函数值置零)，求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >>