华理-蒋志强-matlab_第二次实验报告
(最新版)MATLAB实验报告
![(最新版)MATLAB实验报告](https://img.taocdn.com/s3/m/315fed3afe00bed5b9f3f90f76c66137ee064ff5.png)
(最新版)MATLAB实验报告实验一典型环节的MATLAB仿真一、实验目的1.熟悉MATLAB桌面和命令窗口,初步了解SIMULINK功能模块的使用方法。
2.通过观察典型环节在单位阶跃信号作用下的动态特性,加深对各典型环节响应曲线的理解。
3.定性了解各参数变化对典型环节动态特性的影响。
二、SIMULINK的使用MATLAB中SIMULINK是一个用来对动态系统进行建模、仿真和分析的软件包。
利用SIMULINK功能模块可以快速的建立控制系统的模型,进行仿真和调试。
1.运行MATLAB软件,在命令窗口栏“>>”提示符下键入simulink命令,按Enter 键或在工具栏单击按钮,即可进入如图1-1所示的SIMULINK仿真环境下。
2.选择File菜单下New下的Model命令,新建一个simulink 仿真环境常规模板。
3.在simulink仿真环境下,创建所需要的系统。
以图1-2所示的系统为例,说明基本设计步骤如下:1)进入线性系统模块库,构建传递函数。
点击simulink下的“Continuous”,再将右边窗口中“Transfer Fen”的图标用左键拖至新建的“untitled”窗口。
2)改变模块参数。
在simulink仿真环境“untitled”窗口中双击该图标,即可改变传递函数。
其中方括号内的数字分别为传递函数的分子、分母各次幂由高到低的系数,数字之间用空格隔开;设置完成后,选择OK,即完成该模块的设置。
3)建立其它传递函数模块。
按照上述方法,在不同的simulink 的模块库中,建立系统所需的传递函数模块。
例:比例环节用“Math”右边窗口“Gain”的图标。
4)选取阶跃信号输入函数。
用鼠标点击simulink下的“Source”,将右边窗口中“Step”图标用左键拖至新建的“untitled”窗口,形成一个阶跃函数输入模块。
5)选择输出方式。
用鼠标点击simulink下的“Sinks”,就进入输出方式模块库,通常选用“Scope”的示波器图标,将其用左键拖至新建的“untitled”窗口。
Matlab第二次实验报告
![Matlab第二次实验报告](https://img.taocdn.com/s3/m/3f5ac193ec3a87c24028c46b.png)
Matlab第二次实验报告实验目的:1、了解plot函数和subplot函数的基本用法和matlab绘图的基本原理。
2、了解图形的属性设置。
比如画图的颜色,画图采用的线性标识符等。
二:实验基本知识1:1. 单窗口单曲线绘图;2. 单窗口多曲线绘图;3. 单窗口多曲线分图绘图;4. 多窗口绘图;5.可任意设置颜色与线型;6.图形加注功能;7.fplot——绘制函数图函数;8.ezplot——符号函数的简易绘图函数9:subplot—将画图区域分块函数。
实验内容:将高数课本后的18个图画出上机练习程序第1——6图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1);x=-50:1:50;y=a.*x.^3;subplot(3,2,1);plot(x,y,'r');title('y=a*x^3');xlabel('X');ylabel('Y');x=0:1:50;y=sqrt(a.*x.^3);subplot(3,2,2);plot(x,y,x,-y,'r');title('y^2=a*x^3');xlabel('X');ylabel('Y');x=-3:0.1:3;y=exp(-x.^2);subplot(3,2,3);plot(x,y,'r');title('y=e^x^2');xlabel('X');ylabel('Y');x=-3:0.1:3;y=8*a^3./(x.^2+4*a^2);subplot(3,2,4);plot(x,y,'r');title('y=8*a^3./(x.^2+4*a^2)'); xlabel('X');ylabel('Y');x=0:0.01:5;y=sqrt(x.^3./(a*2.-x)); subplot(3,2,5);plot(x,y,x,-y,'r');title('y^2*(2a-x)=x^3'); xlabel('X');ylabel('Y');%t=-1:0.01:5;%x=3*a.*t./(1+t.^3);%y=3*a*t.^2./(1+t.^3);subplot(3,2,6);ezplot('x.^3+y.^3-3*3*x.*y'); %plot(x,y);title('x^3+y^3-3axy=0'); xlabel('X');ylabel('Y');输入a=3:显示结果第7——10图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0:0.1:2*pi;x=a*cos(t).^3;y=a*sin(t).^3;subplot(2,2,1);plot(x,y,'b')title('x^1.5+y^1.5=a^1.5'); xlabel('X');ylabel('Y');t=-2*pi:0.1:2*pi;x=a.*(t-sin(t));y=a.*(1-cos(t));subplot(2,2,2);plot(x,y,'b')title('°ÚÏß');xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*(1-cos(t));subplot(2,2,3);plot(p.*cos(t),p.*sin(t),'b'); title('p=a(1-cos£¨t£©)'); xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*t;subplot(2,2,4);plot(p.*cos(t),p.*sin(t),'b'); title('p=at');xlabel('X');ylabel('Y');输入a=1,第11题图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º') figure(1);t=-2*pi:pi/100:2*pi;p=exp(a.*t);plot(p.*cos(t),p.*sin(t)); title('p=e^at');xlabel('X');ylabel('Y');输入a=0.1,显示结果:clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0.2*pi:pi/100:100*pi;p=a./t;plot(p.*cos(t),p.*sin(t));title('p=at');xlabel('X');ylabel('Y');输入a=1,显示结果:第13——14题图:clear,clca=input('ÇëÊäÈëaµÄÖµ£º');figure(1);t=0:pi/100:2*pi;p=(a^2.*sin(2*t)).^0.5;subplot(1,2,1);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*sin(2t)');xlabel('X');p=(a^2.*cos(2*t)).^0.5;subplot(1,2,2);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*cos(2t)');xlabel('X');ylabel('Y');输入a=1,显示结果:第15-16题图:clear;clc;a=input('ÇëÊäÈëa:');figure(1)t=0:pi/200:pi;p=a.*cos(3*t);subplot(1,2,1);plot(p.*cos(t),p.*sin(t),'r--');title('p=a*cos(3t)');xlabel('X');p=a.*sin(3*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'r--'); title('p=a*sin(3t)');xlabel('X');ylabel('Y');输入a=1,显示结果:第17——18题图:clear;clc;a=input('ÇëÊäÈëa:');figure(1);t=-pi:pi/200:pi;subplot(1,2,1);p=a.*sin(2*t);plot(p.*cos(t),p.*sin(t),'k-.'); title('p=a*sin(2t)');xlabel('X');ylabel('Y');p=a.*cos(2*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'k-.');title('p=a*cos(2t)');xlabel('X');ylabel('Y');输入a=1,显示结果:实验心得:这次实验主要是学习matlab软件的二维绘图功能:matlab软件还有丰富的图形修饰功能,如改变线条的形式和颜色;除此之外其还可以以多种形式如在一个图形中显示多个函数图形,一个figure中显示多个坐标系,还有一个程序中就可以用多个figure一起来表示不同的函数……在画图过程中你可能因为角度的区间取值不一样,画出来的图形可能跟书本上给出的图形有很大的差异,但没关系,只要你思路是对的,区间可以自己慢慢改动知道跟书本上给出的图形想接近。
MATLAB实验报告2
![MATLAB实验报告2](https://img.taocdn.com/s3/m/3449f9c308a1284ac85043a7.png)
实验二用MATLAB实现线性系统的时域分析[实验目的]1.研究线性系统在典型输入信号作用下的暂态响应;2.熟悉线性系统的暂态性能指标;3.研究二阶系统重要参数阻尼比ξ对系统动态性能的影响;4.熟悉在MATLAB下判断系统稳定性的方法;5.熟悉在MATLAB下求取稳态误差的方法。
[实验原理]MATLAB中有两类用于求解系统时域响应的方法。
一、用MATLAB函数(命令)进行暂态响应分析1 求取线性连续系统的单位阶跃响应的函数——step基本格式为:step(sys)、step(num,den)、step(A,B,C,D)、step(sys,t) 、step(sys1,sys2,…,t) 、y=step(sys,t)、[y,t]=step(sys)、[y,t,x]=step(sys)其中模型对象的类型如下:sys = tf(num,den) 多项式模型、sys = zpk(z,p,k) 零点极点模型、sys = ss(a,b,c,d) 状态空间模型参数无t,表示时间向量t的范围自动设定。
参数有t,表示给定时间向量t,应该有初值,时间增量,末值,如t=0:0.01:2。
前5种函数可以绘出阶跃响应曲线;后3种函数不绘阶跃响应曲线,而是返回响应变量y,时间向量t,以及状态变量x。
2 求取线性连续系统的单位脉冲响应的函数——impulse基本格式为:impulse(sys)、impulse(num,den) 、impulse (sys,tf)、impulse (sys,t) 、impulse (sys1,sys2,…,t) 、y=impulse(sys,t) 、[y,t]=impulse(sys) 、[y,t,x]=impulse (sys)3 求取线性连续系统的单位斜坡响应MATLAB没有直接求系统斜坡响应的功能函数。
在求取控制系统的斜坡响应时,通常用阶跃响应函数step()求取传递函数为G (s)/s的系统的阶跃响应,则其结果就是原系统G (s)的斜坡响应。
MATLAB实验报告(2)
![MATLAB实验报告(2)](https://img.taocdn.com/s3/m/9487e8fc5122aaea998fcc22bcd126fff7055d06.png)
MATLAB实验报告(2)广西大学实验报告纸姓名:徐律1002100405电气工程学院自动化102班2012年12月27日【实验名称】MATLAB的图形绘制【实验要求】1、学习MATLAB图形绘制的基本方法;2、熟悉和了解MATLAB图形绘制程序编辑的基本指令;3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面,并添加图形的各种标注;4、掌握plot、subplot的指令格式和语法。
【实验基本知识】1、plot(x,y)------绘制由x,y所确定的曲线;2、多组变量绘图:plot(x1, y1, 选项1, x2, y2, 选项2, ……);3、双Y轴绘图:plotyy()函数;4、图形窗口的分割;5、图形编辑窗口的使用。
【实验内容】【1】二维曲线绘图基本指令演示。
本例运作后,再试验plot(t), plot(Y), plot(Y,t) ,以观察产生图形的不同。
t=(0:pi/50:2*pi)';k=0.4:0.1:1;Y=cos(t)*k;plot(t,Y)plot(t,Y); plot(t);plot(Y);plot(Y,t); %图形如下图所示【2】用图形表示连续调制波形Y=sin(t)sin(9t)及其包络线。
t=(0:pi/100:pi)';y1=sin(t)*[1,-1];y2=sin(t).*sin(9*t);t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y2,'b',t3,y3,'bo')axis([0,pi,-1,1])【3】在一个图形窗口绘制正弦和余弦曲线,要求给图形加标题“正弦和余弦曲线”,X轴Y轴分别标注为“时间t”和“正弦、余弦”,在图形的某个位置标注“sin(t)”“cos(t)”,并加图例,显示网格,坐标为正方形坐标系。
实验程序:实验结果:Sin(t)Cos(t)【4】绘制向量x=[1 3 0.5 2.5 2]的饼形图,并把3对应的部分分离出来。
MATLAB课后上机实验报告
![MATLAB课后上机实验报告](https://img.taocdn.com/s3/m/0bb1431cac02de80d4d8d15abe23482fb5da0216.png)
《MATLAB程序设计与应用》实验指导书实验一 matlab 集成环境使用与运算基础1,先求下列表达式的值,然后显示matlab 工作空间的使用情况并保存全部变量。
(1)0122sin851z e =+程序:.>> z1=2*sin(85*pi/180)/(1+exp(2)) 结果: z1 =0.2375(2)22121(0.4552i z In x x +⎡⎤=+=⎢⎥-⎣⎦其中 程序:>> x=[2,1+2*i;-0.45,5];>> z2=0.5*log(x+sqrt(1+x*x)) 结果: z2 =0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i(3)0.3,9.2,8.2,...,8.2,9.2,0.3,23.0)3.0sin(23.03.03---=+++-=-a aIn a e e z a a 提示:利用冒号表达式生成a 向量,求各点函数值时用点乘运算。
程序:>> a=-3.0:0.1:30;>> z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin((a+0.3)*pi/180)+log((0.3+a)/2) 结果: z3 =1.0e+003 *Columns 1 through 40.0003 + 0.0031i 0.0003 + 0.0031i 0.0003 + 0.0031i 0.0002 + 0.0031iColumns 5 through 80.0002 + 0.0031i 0.0001 + 0.0031i 0.0001 + 0.0031i 0.0000 + 0.0031i Columns 9 through 12-0.0000 + 0.0031i -0.0001 + 0.0031i -0.0001 + 0.0031i -0.0002 + 0.0031i Columns 13 through 16-0.0003 + 0.0031i -0.0003 + 0.0031i -0.0004 + 0.0031i -0.0005 + 0.0031i Columns 17 through 20-0.0006 + 0.0031i -0.0007 + 0.0031i -0.0008 + 0.0031i -0.0009 + 0.0031i Columns 21 through 24-0.0010 + 0.0031i -0.0012 + 0.0031i -0.0014 + 0.0031i -0.0016 + 0.0031i Columns 25 through 28-0.0019 + 0.0031i -0.0023 + 0.0031i -0.0030 + 0.0031i -0.0370 Columns 29 through 32-0.0030 -0.0023 -0.0019 -0.0016 Columns 33 through 36-0.0014 -0.0012 -0.0010 -0.0009 Columns 37 through 40-0.0008 -0.0007 -0.0006 -0.0005 Columns 41 through 44-0.0004 -0.0003 -0.0003 -0.0002 Columns 45 through 48-0.0001 -0.0001 -0.0000 0.0000 Columns 49 through 520.0001 0.0001 0.0002 0.0002 Columns 53 through 560.0003 0.0003 0.0003 0.0004 Columns 57 through 600.0004 0.0005 0.0005 0.0005 Columns 61 through 640.0006 0.0006 0.0006 0.0007 Columns 65 through 680.0007 0.0007 0.0008 0.0008 Columns 69 through 720.0008 0.0008 0.0009 0.0009 Columns 73 through 760.0009 0.0010 0.0010 0.0010 Columns 77 through 800.0011 0.0011 0.0011 0.0011 Columns 81 through 840.0012 0.0012 0.0012 0.0013 Columns 85 through 880.0013 0.0013 0.0013 0.0014 Columns 89 through 920.0014 0.0014 0.0015 0.0015 Columns 93 through 960.0015 0.0016 0.0016 0.0016 Columns 97 through 1000.0017 0.0017 0.0017 0.0018 Columns 101 through 1040.0018 0.0018 0.0019 0.0019 Columns 105 through 1080.0020 0.0020 0.0020 0.0021 Columns 109 through 1120.0021 0.0022 0.0022 0.0023 Columns 113 through 1160.0023 0.0024 0.0024 0.0025 Columns 117 through 1200.0025 0.0026 0.0026 0.0027 Columns 121 through 1240.0027 0.0028 0.0029 0.0029 Columns 125 through 1280.0030 0.0031 0.0031 0.0032 Columns 129 through 1320.0033 0.0034 0.0034 0.0035 Columns 133 through 1360.0036 0.0037 0.0038 0.0039 Columns 137 through 1400.0040 0.0041 0.0042 0.0043 Columns 141 through 1440.0044 0.0045 0.0046 0.0047 Columns 145 through 1480.0049 0.0050 0.0051 0.0053 Columns 149 through 1520.0054 0.0056 0.0057 0.0059 Columns 153 through 1560.0060 0.0062 0.0064 0.0066 Columns 157 through 1600.0068 0.0069 0.0071 0.0074 Columns 161 through 1640.0076 0.0078 0.0080 0.0083 Columns 165 through 1680.0085 0.0088 0.0090 0.0093 Columns 169 through 1720.0096 0.0099 0.0102 0.0105 Columns 173 through 1760.0108 0.0112 0.0115 0.0119 Columns 177 through 1800.0123 0.0127 0.0131 0.0135 Columns 181 through 1840.0139 0.0144 0.0148 0.0153 Columns 185 through 1880.0158 0.0163 0.0168 0.0174 Columns 189 through 1920.0180 0.0185 0.0191 0.0198 Columns 193 through 1960.0204 0.0211 0.0218 0.0225 Columns 197 through 2000.0233 0.0241 0.0249 0.0257 Columns 201 through 2040.0265 0.0274 0.0284 0.0293 Columns 205 through 2080.0303 0.0313 0.0324 0.0335 Columns 209 through 2120.0346 0.0358 0.0370 0.0382 Columns 213 through 2160.0395 0.0409 0.0423 0.0437 Columns 217 through 2200.0452 0.0467 0.0483 0.0500 Columns 221 through 2240.0517 0.0534 0.0552 0.0571 Columns 225 through 2280.0591 0.0611 0.0632 0.0654 Columns 229 through 2320.0676 0.0699 0.0723 0.0748 Columns 233 through 2360.0773 0.0800 0.0827 0.0856 Columns 237 through 2400.0885 0.0915 0.0947 0.0979 Columns 241 through 2440.1013 0.1047 0.1083 0.1121 Columns 245 through 2480.1159 0.1199 0.1240 0.1282 Columns 249 through 2520.1326 0.1372 0.1419 0.1467 Columns 253 through 2560.1518 0.1570 0.1624 0.1679 Columns 257 through 2600.1737 0.1796 0.1858 0.1921 Columns 261 through 2640.1987 0.2055 0.2125 0.2198 Columns 265 through 2680.2273 0.2351 0.2431 0.2514 Columns 269 through 2720.2600 0.2689 0.2781 0.2876 Columns 273 through 2760.2974 0.3076 0.3180 0.3289 Columns 277 through 2800.3401 0.3517 0.3637 0.3761 Columns 281 through 2840.3889 0.4021 0.4158 0.4299 Columns 285 through 2880.4446 0.4597 0.4753 0.4915 Columns 289 through 2920.5082 0.5254 0.5433 0.5617 Columns 293 through 2960.5807 0.6004 0.6208 0.6418 Columns 297 through 3000.6636 0.6861 0.7093 0.7333 Columns 301 through 3040.7581 0.7838 0.8103 0.8377 Columns 305 through 3080.8660 0.8952 0.9254 0.9567 Columns 309 through 3120.9890 1.0223 1.0568 1.0924 Columns 313 through 3161.1292 1.1673 1.2066 1.2472Columns 317 through 3201.2892 1.3326 1.3774 1.4237Columns 321 through 3241.4715 1.5210 1.5721 1.6249Columns 325 through 3281.6794 1.7357 1.7940 1.8541Columns 329 through 3311.9163 1.98052.0468(4)⎪⎩⎪⎨⎧=<≤<≤<≤+--=5.2:5.0:0,322110,121,2224t t t t t t t t z 其中提示:用逻辑表达式求分段函数值。
Matlab实验报告
![Matlab实验报告](https://img.taocdn.com/s3/m/457d0bc74028915f804dc248.png)
实验一Matlab语言编程一实验目的1.熟悉Matlab语言及其环境,掌握编程方法。
2.要求认真听取实验指导老师的讲解及其演示。
二具体实验内容、步骤、要求1.运行交互式学习软件,学习Matlab语言。
2.在Matlab的命令窗口下输入INTRO然后,根据显示出来的幻灯片右面的按钮进行操作,可按START→NEXT→NEXT 按钮一步一步的运行,观察。
3.自编程序并完成上机编程,调试,运行,存盘。
(1)利用Matlab运行给中矩阵运算,例如:A=则运算结果:A(:, 1)ans =11213141A(2,:)ans =21 22 23 24A(1:2,2:3)ans =12 1322 23A(:,1:2)ans =11 1221 2231 3241 42ones(2,2)ans =1 11 1eye(2)ans =1 00 1(2)绘制数学函数模型绘制函数在命令窗口输入:>> t=[0:0.01:8];>> y=1-2*exp(1)-t.*sin(t);>> plot(t,y)结果如下:4.理解函数文件与命令文件的区别,并自编函数文件并调用。
实验二数值积分算法练习及函数调用一实验目的理解数值积分的算法,熟练掌握Matlab函数的调用。
方程(2-2)为,x(0)=1,t=[0 3]Euler方法求解为:>> t=[0:0.1:3];x0=1;h1=0.1;n=length(t);u=x0;uu(1)=u;for i=2:ndu=u-t(i-1)^2;u=du*h1+u;uu(i)=u;end>> plot(t,uu)其图为:而ode23画出的图为:二者比较可知ode23曲线没有Euler图那么的圆滑。
式2-1为du/dt=u-2t/u, u(0)=1,t=[0 1]其ode23算法为:fun21.m函数:function xdot=fun21(t,x)xdot = x-2*t/x;fxyz21.m函数:t0=0;tf=1;tol=1e-6;x0=1;trace=1;[t,x]=ode23('fun21',t0,tf,x0,tol,trace);plot(t,x)其图为:其Euler算法的图为:通过比较得知ode23算法的精度要高于Euler算法。
MATLAB上机实验报告
![MATLAB上机实验报告](https://img.taocdn.com/s3/m/3bcdb3e66294dd88d0d26b2e.png)
Matlab上机实验报告实验二读入MATLAB下自带图像pout.tif1)利用亮度变换函数,调整图像亮度。
a)调整范围设定[0 1],[1 0],观测显示效果;b)调整范围设定[0.5 0.75],[1 0],观测显示效果。
解:a, I=imread('pout.tif');colormap;imshow(I);j=imadjust(I,[0 1],[1 0],1.5);figure;subimage(j);b, >> I=imread('pout.tif');colormap;imshow(I);j=imadjust(I,[0.5 0.75],[1 0],1.5);figure;subimage(j);2)利用对比度拉伸函数,压缩高值灰度(c值自行设定)。
解:I=imread('pout.tif');colormap;subplot(1,2,1);imshow(I);xlabel('a)原始图像');J=double(I);J=100*log(J+1);I=uint8(J);subplot(1,2,2);subimage(J);xlabel('b)非线性变换');3)利用直方图函数,生成并绘制图像直方图。
解:I=imread('pout.tif');subplot(1,2,1);imshow(I);title('原始图像');subplot(1,2,2);imhist(I);4)利用直方图修正函数,生成均衡化后的图像直方图(n值自定设定)。
解:I=imread('pout.tif');figure(1);subplot(1,2,1);imshow(I);xlabel('a)原始图像');J=histeq(I);figure(1);subplot(1,2,2);imshow(J);xlabel('b)直方图均衡');figure(2);imhist(I,100);figure(3);imhist(J,100);实验三1.运行例3、4,显示并分析输出结果,说明逆滤波和维纳滤波的区别。
MATLAB实验报告
![MATLAB实验报告](https://img.taocdn.com/s3/m/a9a321dc0408763231126edb6f1aff00bed57077.png)
MATLAB实验报告实验⼀ MATLAB 环境的熟悉与基本运算⼀、实验⽬的及要求1.熟悉MATLAB 的开发环境;2.掌握MATLAB 的⼀些常⽤命令;3.掌握矩阵、变量、表达式的输⼊⽅法及各种基本运算。
⼆、实验内容1、熟悉MATLAB 的开发环境: ① MATLAB 的各种窗⼝:命令窗⼝、命令历史窗⼝、⼯作空间窗⼝、当前路径窗⼝。
②路径的设置:建⽴⾃⼰的⽂件夹,加⼊到MATLAB 路径中,并保存。
? 设置当前路径,以⽅便⽂件管理。
2、学习使⽤clc 、clear,了解其功能与作⽤。
3、矩阵运算:已知:A=[1 2;3 4]; B=[5 5;7 8]; 求:A*B 、A 、*B,并⽐较结果。
4、使⽤冒号选出指定元素:已知:A=[1 2 3;4 5 6;7 8 9]; 求:A 中第3列前2个元素;A 中所有列第2,3⾏的元素; 5、在MATLAB 的命令窗⼝计算: 1))2sin(π2) 5.4)4.05589(÷?+ 6、关系及逻辑运算1)已知:a=[5:1:15]; b=[1 2 8 8 7 10 12 11 13 14 15],求: y=a==b ,并分析结果 2)已知:X=[0 1;1 0]; Y=[0 0;1 0],求: x&y+x>y ,并分析结果 7、⽂件操作1)将0到1000的所有整数,写⼊到D 盘下的data 、txt ⽂件 2)读⼊D 盘下的data 、txt ⽂件,并赋给变量num 8、符号运算1)对表达式f=x 3-1 进⾏因式分解2)对表达式f=(2x 2*(x+3)-10)*t ,分别将⾃变量x 与t 的同类项合并 3)求3(1)xdz z +?三、实验报告要求完成实验内容的3、4、5、6、7、8,写出相应的程序、结果实验⼆ MATLAB 语⾔的程序设计⼀、实验⽬的1、熟悉 MATLAB 程序编辑与设计环境2、掌握各种编程语句语法规则及程序设计⽅法3、函数⽂件的编写与设计4、了解与熟悉变量传递与赋值⼆、实验内容1.编写程序,计算1+3+5+7+…+(2n+1)的值(⽤input 语句输⼊n 值)。
MatLAB第二次实验
![MatLAB第二次实验](https://img.taocdn.com/s3/m/3bb0a73610661ed9ad51f39a.png)
自动控制理论实 验 报 告姓 名 班 级 成 绩_______________学 号 日 期 合作者________________________________________________________________________________________________一、实验题目:利用MATLAB 求系统的的瞬态性能指标 二、实验目的:1、熟悉并掌握MatlAB 软件在工程计算中的应用。
2、求系统的的瞬态性能指标。
三、实验仪器:电脑、MatlAB 软件; 四、实验内容用impulse 函数求单位脉冲响应,step 求求单位阶跃响应。
程序:%a_1.mt=[0:0.01:0.8] %num_G=[50];tao=0;den_G=[0.05 1+50*tao 50];G1=tf(num_G ,den_G);tao=0.0125;den_G=[0.05 50*tao tao 50];G2=tf(num_G ,den_G); tao=0.025;den_G=[0.05 1+50*tao 50];G3=tf(num_G ,den_G); %[y1,T]=impulse(G1,t);[y1a,T]=step(G1,t); [y2,T]=impulse(G2,t);[y2a,T]=step(G2,t); [y3,T]=impulse(G3,t);[y3a,T]=step(G3,t); %subplot(1,2,1)plot(T,y1,'r-',T,y2,'g-',T,y3,'b-')legend('tao=0','tao=0.0125','tao=0.025') xlabel('t(sec)') ylabel('x(t)'); grid on;subplot(1,2,2)plot(T,y1a,'r-',T,y2a,'g-',T,y3a,'b-')legend('tao=0','tao=0.0125','tao=0.025') grid on;xlabel('t(sec)') ylabel('x(t)');二、用lsim函数求任意输入响应程序:%3_2.mclc;t=[0:0.01:1];u=sin(2*pi*t);%tao=0.025;num_G=[50];den_G=[0.05 1+50*tao 50];G=tf(num_G,den_G);%y=lsim(G,u,t);%plot(t,u,'r-',t,y,'g-',t,u','b',t,y','r',t,u'-y,'b-','linewidth',1)legend('u(t)','xo(t)','u(t)','xo(t)','e(t)=u-y')grid;xlabel('t(sec)')ylabel('x(t)');三、利用MATLAB求系统的的瞬态性能指标程序:close all;%t=0:0.001:1;%yss=1;dta=0.02;%num_G=[50];tao=0;den_G=[0.05 1+50*tao 50];G1=tf(num_G,den_G);tao=0.0125;den_G=[0.05 1+50*tao 50];G2=tf(num_G,den_G);tao=0.025;den_G=[0.05 1+50*tao 50];G3=tf(num_G,den_G);y1=step(G1,t);y2=step(G2,t);y3=step(G3,t);%r=1;while y1(r)<yss;r=r+1;endtr1=(r-1)*0.001;%[ymax,tp]=max(y1);tp1=(tp-1)*0.001;%mp1=(ymax-yss)/yss;%s=1001;while y1(s)>1-dta & y1(s)<1+dta;s=s-1;endts1=(s-1)*0.001;%τ=0.0125r=1;while y2(r)<yss;r=r+1;endtr2=(r-1)*0.001;[ymax,tp]=max(y2);tp2=(tp-1)*0.001;mp2=(ymax-yss)/yss;s=1001;while y2(s)>1-dta &y2(s)<1+dta;s=s-1;endts2=(s-1)*0.001;%τ=0.025r=1;while y3(r)<yss;r=r+1;endtr3=(r-1)*0.001;[ymax,tp]=max(y3);tp3=(tp-1)*0.001;mp3=(ymax-yss)/yss;s=1001;while y3(s)>1-dta & y3(s)<1+dta;s=s-1;endts3=(s-1)*0.001;%[tr1 tp1 mp1 ts1;tr2 tp2 mp2 ts2;tr3 tp3 mp3 ts3]结果:ans =0.0640 0.1050 0.3509 0.3530 -------系统10.0780 0.1160 0.1523 0.2500 -------系统20.1070 0.1410 0.0415 0.1880 -------系统3任课教师:。
MATLAB实验报告(2)
![MATLAB实验报告(2)](https://img.taocdn.com/s3/m/c8a4a4a268dc5022aaea998fcc22bcd126ff4221.png)
MATLAB实验报告(2)仲恺农业⼯程学院实验报告纸_⾃动化学院_(院、系)_⼯业⾃动化_专业_144_班_Matlab仿真控制实践课程实验⼀MATLAB绘图基础⼀、实验⽬的了解MATLAB常⽤命令和常见的内建函数使⽤。
熟悉矩阵基本运算以及点运算。
掌握MATLAB绘图的基本操作:向量初始化、向量基本运算、绘图命令plot,plot3,mesh,surf 使⽤、绘制多个图形的⽅法。
⼆、实验内容建⽴并执⾏M⽂件multi_plot.m,使之画出如图的曲线。
三、实验⽅法四、实验要求1.分析给出的MA TLAB参考程序,理解MA TLAB程序设计的思维⽅法及其结构。
2.添加或更改程序中的指令和参数,预想其效果并验证,并对各语句做出详细注释。
对不熟悉的指令可通过HELP查看帮助⽂件了解其使⽤⽅法。
达到熟悉MA TLAB画图操作的⽬的。
3.总结MATLAB中常⽤指令的作⽤及其调⽤格式。
五、实验思考1、实现同时画出多图还有其它⽅法,请思考怎样实现,并给出⼀种实现⽅法。
(参考程序如下)2、思考三维曲线(plot3)与曲⾯(mesh, surf)的⽤法,(1)绘制参数⽅程233,)3cos(,)3sin()(t z e t t y e t t t x t t ===--的三维曲线;(2)绘制⼆元函数xyy x ex x y x f z ----==22)2(),(2,在XOY 平⾯内选择⼀个区域(-3:0.1:3,-2:0.1:2),然后绘制出其三维表⾯图形。
(以下给出PLOT3和SURF 的⽰例)绘制题⽬要求曲⾯:%绘制⼆元函数,在XOY平⾯内选择⼀个区域(-3:0.1:3,-2:0.1:2)仲恺农业⼯程学院实验报告纸_⾃动化学院_(院、系)_⼯业⾃动化_专业_144_班_Matlab 仿真控制实践课程实验⼆:基于Simulink 的控制系统仿真实验⽬的1.掌握MATLAB 软件的Simulink 平台的基本操作; 2.能够利⽤Simulink 平台研究PID 控制器对系统的影响;实验原理PID (⽐例-积分-微分)控制器是⽬前在实际⼯程中应⽤最为⼴泛的⼀种控制策略。
MATLAB第二次上机实验报告
![MATLAB第二次上机实验报告](https://img.taocdn.com/s3/m/88646e88ec3a87c24028c473.png)
电子科技大学电子工程学院标准实验报告(实验)课程名称MATLAB与数值分析学生姓名:学号:指导教师:一、实验名称实验二 线性方程组求解和函数的数值逼近二、实验目的通过上机实验,使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。
实验涉及的核心知识点:病态方程求解、矩阵分解和方程组求解、Lagrange 插值。
实验重点与难点:算法设计和MATLAB 编程三、实验内容1. 对高阶多项式()()()()()2011220k p x x x x x k ==---=-∏编程求下面方程的解()190p x x ε+=并绘图演示方程的解与扰动量ε的关系。
2. 对220n =,生成对应的Hilbert 矩阵,计算矩阵的条件数;通过先确定解获得常向量b 的方法,确定方程组()n H x b =最后,用矩阵分解方法求解方程组,并分析计算结果。
3. 对函数()21125f x x =+ []1,1x ∈-的Chebyshev 点()()21cos 21k k x n π⎛⎫-= ⎪ ⎪+⎝⎭,1,2,,1k n =+编程进行Lagrange 插值,并分析插值结果。
四、实验数据及结果分析1. 对高阶多项式()()()()()2011220k p x x x x x k ==---=-∏编程求下面方程的解()190p x x ε+=并绘图演示方程的解与扰动量ε的关系。
p=[1,-1]; for i=2:20 n=[1,-i];p=conv(p,n); % 求多项式乘积 endm=zeros(1,21); % m 的最高次幂为20,有21项 hold on x=1:20;d=[-1,0,0.1,0.5,1]; for i=1:5delt=d(i); m(2)=delt;y=(roots(p+m))'; % 求多项式的根 plot(x,y,'-o','color',[i/5,i/20,i/10]); endtitle('方程p(x)=0的解与扰动量delt 的关系')legend('delt=-1','delt=0','delt=0.1','delt=0.5','delt=1')24681012141618200102030405060方程p(x)=0的解与扰动量delt 的关系delt=-1delt=0delt=0.1delt=0.5delt=12.对220n =,生成对应的Hilbert 矩阵,计算矩阵的条件数;通过先确定解获得常向量b 的方法,确定方程组()n H x b =最后,用矩阵分解方法求解方程组,并分析计算结果。
华理-蒋志强-matlab_第二次实验报告
![华理-蒋志强-matlab_第二次实验报告](https://img.taocdn.com/s3/m/6b6dc638bcd126fff7050bf9.png)
华东理工大学商学院数值分析实验报告姓名:李彬班级:金融103学号:10102469实验名称:MATLAB编程基础和拉格朗日插值指导教师:蒋志强实验地址:奉贤校区信息楼319室实验时间:2012/10/1113:30-15:10【实验步骤】问题1.已知x=[012 2.53 3.54],y=[2.50.50.5 1.5 1.5 1.1240],编程求拉格朗日插值函数Pn(2.5),Pn(4.3),并画出拉格朗日插值函数Pn(x)在x ∈[04]上的图像。
(1)定义拉格朗日插值函数,将其保存在LB_ex_1.m中,具体实现程序如下:function P=LB_ex_1(a,x,y)%计算拉格朗日插值多项式,x为自变量,y为因变量,a为插值点m=length(x);n=length(y);if m~=nerror('x和y的长度需一致')endP=zeros(size(a));for j=1:nL=1;for i=1:nif i~=jL=L.*(a-x(i))/(x(j)-x(i));endendP=P+L*y(j);endend(2)具体编程实现如下:>>x=[0122.533.54];>>y=[2.50.50.51.51.51.1240];>>a=0:0.1:4;>>P=LB_ex_1(a,x,y);>>plot(a,P,x,y,'o');>>xlabel('x');>>ylabel('P');问题2.课本109页习题3编程实现如下:>>x=0:1:6;>>y=[2,2.1592,3.1687,5.4332,9.1411,14.406,21.303];>>plot(x,y);>>xlabel('x');ylabel('y');问题3.课本110页习题5(1)仍然使用拉格朗日公式进行插值,使用LB_ex_1.m,取x0为0,x1为1,记为x=[01],相应e x记为y,y=exp(x).插值点记为a,插值结果记为P,误差记为E,E=exp(a)-P(2)实现步骤如下:>>x=[01];>>y=exp(x);>>a=-1:0.2:2;>>P=LB_ex_1(a,x,y);>>E=exp(a)-P;>>plot(a,E);>>hold on>>plot([-12],[00],'--'); >>xlabel('x');ylabel('E');【分析讨论】本次实验中,我。
MATLABA第二次上机实验
![MATLABA第二次上机实验](https://img.taocdn.com/s3/m/25682f44fe4733687e21aa86.png)
实验二MATLAB矩阵分析与处理2.产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P,且求其行列式的值Hh和Hp以及他们的条件数Th和Tp,判断哪个矩阵性能更好。
为什么?程序代码及运算结果如下:> H=hilb(5)H =1.0000 0.5000 0.3333 0.2500 0.20000.5000 0.3333 0.2500 0.2000 0.16670.3333 0.2500 0.2000 0.1667 0.14290.2500 0.2000 0.1667 0.1429 0.12500.2000 0.1667 0.1429 0.1250 0.1111>> P=pascal(5)P =1 1 1 1 11 2 3 4 51 3 6 10 151 4 10 20 351 5 15 35 70>> Hh=det(H)Hh =3.7493e-012>> Hp=det(P)Hp =1>> Th=cond(H)Th =4.7661e+005>> Tp=cond(P)Tp =8.5175e+003由运算结果知Tp的值越接近1,故矩阵P的性能越好。
5.下面是一个线性方程组:A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;](1)求方程的解。
(2)将方程右边向量元素b3改为0.53,在求解,并比较b3的变化和解的相对变化。
(3)计算系数矩阵A的条件数并分析结果。
程序代码及结果如下:>> A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6;]A =0.5000 0.3333 0.25000.3333 0.2500 0.20000.2500 0.2000 0.1667>> b=[0.95;0.67;0.52;]b =0.95000.67000.5200>> x=inv(A)*bx =1.20000.60000.6000>> b1=[0.95;0.67;0.53]b1 =0.95000.67000.5300>> x=inv(A)*b1x =3.0000-6.60006.6000>> B=cond(A)B =1.3533e+003由B的值知,B离1较远,故矩阵x的稳定性较差,改变b3的值为0.53,虽b3的值只增加了0.01,所得结果却发生的很大变化,这都是矩阵A的稳点性所决定的。
MATLAB实验报告——运用MATLAB求解和分析线性时不变系统资料
![MATLAB实验报告——运用MATLAB求解和分析线性时不变系统资料](https://img.taocdn.com/s3/m/fee882507fd5360cba1adb6e.png)
MATLAB实验报告课程名称 MATLAB程序设计实验日期 2015 年 05 月 11 日学生姓名学号班级实验名称运用MATLAB求解和分析线性时不变系统实验仪器MATLAB7.1 Windows XP实验目的1.熟悉线性时不变LTI系统在典型激励信号下的响应及其特征。
2.掌握线性LTI系统单位冲激响应的求解方法。
3.熟悉MATLAB的相关函数的调用格式和作用。
4.会用MATLAB对系统进行时域分析。
实验基本原理1.Impulse函数功能:计算并画出系统的冲激响应。
调用格式: Impulse(sys):其中sys 可以是利用命令tf,zpk或ss建立的系统函数。
Impulse(sys,t)计算并画出系统在向量t定义的时间内的冲激响应。
Y=impulse(sys,t):保存系统的输出值。
2.Step函数功能:计算并画出系统的阶跃响应曲线。
调用格式: Step(sys):其中sys可以是利用命令tf,zpk,或ss 建立的系统。
Step(sys,t):计算并画出系统在向量t定义的时间内的阶跃响应。
3.Lsim函数功能:计算并画出系统在任意输入下的零状态响应。
调用格式: Lsim(sys,x,t):其中sys可以是利用命令tf,zpk或ss建立的系统函数,x是系统的输入,t定义的是时间范围。
Lsim(sys,x,t,zi):计算出系统在任意输入和零状态下的全响应,sys必须是状态空间形式的系统函数,zi是系统的原始状态。
4.roots函数功能:计算齐次多项式的根。
调用格式: r=roots(b):计算多项式b的根,r为多项式的根。
5.impz函数功能: 求离散系统单位脉冲响应,并绘制其时域波形。
调用格式:impz(b ,a) :以默认方式绘出向量a , b 定义的离散系统的单位脉冲响应的离散时域波形.impz(b ,a ,n) :绘出由向量a , b定义的离散系统在0—n (n必须为整数)离散时间范围内的单位序列响应的时域波形.impz(b ,a ,n1:n2) : 绘出由向量a , b定义的离散系统在n1—n2(n1 , n2必须为整数,且n1<n2)离散时间范围内的单位序列响应的时域波形。
[IT计算机]MATLAB实验报告
![[IT计算机]MATLAB实验报告](https://img.taocdn.com/s3/m/1286c93d4a35eefdc8d376eeaeaad1f34693113e.png)
MATLAB实验报告姓名:学号:班级:专业:指导老师:实验一MATLAB运算基础一、实验目的:1、熟悉启动和退出MATLAB的方法。
2、熟悉MATLAB命令窗口的组成。
3、掌握建立矩阵的方法、4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。
二、实验内容:1.1、>> t=0:0.5:2.5t =0 0.5000 1.0000 1.5000 2.0000 2.5000>>y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*((t>=1)&(t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2)&(t <3))y =0 0.2500 0 1.2500 1.0000 2.2500分析:该实验用逻辑表达式求出了分段函数的值。
1.2、>> p=rem([100:999],21)==0;>> sum(p)ans =43>> p1=rem([100:999],21),A=find(p1==0),length(A)ans =43分析:该实验求出了[100,999]之间能被21整除的数的个数为43.1.3、>> ch='ABcdefGHd',k=find(ch>='A'&ch<='Z'),ch(k)=[]ch =ABcdefGHdk =1 2 7 8ch =cdefd分析:该实验建立了一个字符串向量,删除了其中的大写字母。
1.4、A=[1 2 3;4 5 6;7 8 9],[i,j]=find(A>=5),for n=1:length(i)m(n)=A(i(n),j(n))endm>> mm =7 5 8 6 9分析:该实验求出了矩阵A中大于或等于5的元素为7 5 8 6 9.1.5、>>a11=input('all='),a12=input('a12='),a21=input('a21='),a22=input('a22='),A=[a 11,a12;a21,a22],DA=det(A),IN=inv(A),EA=eig(A)a11=23a12=12a21=11a22=32A =23 1211 32DA =604IN =0.0530 -0.0199-0.0182 0.0381EA =15.161039.8390分析:该实验求出了矩阵A中的行列式的值、逆和特征值。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
华东理工大学商学院
数值分析实验报告
姓名:李彬
班级:金融103
学号:10102469
实验名称:MATLAB编程基础和拉格朗日插值
指导教师:蒋志强
实验地址:奉贤校区信息楼319室
实验时间:2012/10/1113:30-15:10
【实验步骤】
问题1.已知x=[012 2.53 3.54],y=[2.50.50.5 1.5 1.5 1.1240],编程求拉格朗日插值函数Pn(2.5),Pn(4.3),并画出拉格朗日插值函数Pn(x)在x ∈[04]上的图像。
(1)定义拉格朗日插值函数,将其保存在LB_ex_1.m中,具体实现程序如下:
function P=LB_ex_1(a,x,y)
%计算拉格朗日插值多项式,x为自变量,y为因变量,a为插值点
m=length(x);
n=length(y);
if m~=n
error('x和y的长度需一致')
end
P=zeros(size(a));
for j=1:n
L=1;
for i=1:n
if i~=j
L=L.*(a-x(i))/(x(j)-x(i));
end
end
P=P+L*y(j);
end
end
(2)具体编程实现如下:
>>x=[0122.533.54];
>>y=[2.50.50.51.51.51.1240];
>>a=0:0.1:4;
>>P=LB_ex_1(a,x,y);
>>plot(a,P,x,y,'o');
>>xlabel('x');
>>ylabel('P');
问题2.课本109页习题3
编程实现如下:
>>x=0:1:6;
>>y=[2,2.1592,3.1687,5.4332,9.1411,14.406,21.303];
>>plot(x,y);
>>xlabel('x');ylabel('y');
问题3.课本110页习题5
(1)仍然使用拉格朗日公式进行插值,使用LB_ex_1.m,取x0为0,x1为1,记为x=[01],相应e x记为y,y=exp(x).插值点记为a,插值结果记为P,误差记为E,E=exp(a)-P
(2)实现步骤如下:
>>x=[01];
>>y=exp(x);
>>a=-1:0.2:2;
>>P=LB_ex_1(a,x,y);
>>E=exp(a)-P;
>>plot(a,E);
>>hold on
>>plot([-12],[00],'--'); >>xlabel('x');ylabel('E');
【分析讨论】
本次实验中,我。