MATLAB实验上机所用word(数值计算篇)
第三章数值计算3.1LU分解和恰定方程组的解Matlab矩阵的分解形式主要有：三角分解、正交化、特征值分解。
3.1.1LU三角分解、行列式和逆1、LU分解是利用Gaussian列主元消去法进行的。
为了保证主元消去策略的实施，一般来说，必须对被分解矩阵实施行置换因此有：LU=PA 式中L为主对角元为1的下三角矩阵，U是上三角矩阵，P是由0或1组成的行置换矩阵：［L UP］=lu(A)2、A-1=U-1L-1P Matlab根据以上算法编制了相应的指令：det(a) 求矩阵a的行列式inv(a) 求矩阵a的逆矩阵3.1.2恰定方程组的解在求解方程式，尽量不要使用inv(a)*b指令，而应采用a\b，后者计算速度比前者快，精确度高。
【例3.1.2-1】“求逆”法和“左除”法解恰定方程的性能对比（1）randn('state',0);A=gallery('randsvd',100,2e13,2); ％产生条件数为2e13的100阶随机矩阵x=ones(100,1);b=A*x;（2）ticxi=inv(A)*b;ti=toc（3）tic;xd=A\b;td=toc;3.1.3矩阵特征值和矩阵函数矩阵A与向量x相乘,即表示矩阵对向量的变换(transformation)，一般说来，向量在变换的作用下将发生旋转（rotation），反射(reflection)和放大缩小。
但对于任何一个矩阵来说，中存在那么一些特殊的向量，再对其变化的作用下，向量的方向不变，而仅长短发生变化。
这种向量就是所谓的特征向量（eigenvector），它满足方程 Ax=λx3.1.3.1特征值和特征向量的求取d=eig(a) 仅仅计算a的特征值[v d]=eig(a) 计算矩阵a的特征向量阵v和特征对角阵d，使av=vd成立。
【例3.1.3.1-1】简单实阵的特征值问题。
A=[1,-3;2,2/3];[V,D]=eig(A)V =0.7746 0.77460.0430 - 0.6310i 0.0430 + 0.6310iD =0.8333 + 2.4438i 00 0.8333 - 2.4438i3.2数据分析3.2.1 基本统计函数指令note: 1、median(x) 当x为向量时，先把x元素有小到大排列，然后在新排成数组中；n为奇取（n+1）/2个元素；n为偶数时取n/2与(n/2+1)元素的平均值，作为总位数。
2、输入为向量，则运算对整个向量进行的。
若输入是数组，那么指令运算时按列进行的。
结果为一行向量。
3、Max(a,b)取a b矩阵中加大的数.clear all; a=reshape(1:9,3,3);b=reshape(1:16,4,4);median(a),median(b) ans =2 5 8ans =2.5000 6.5000 10.5000 14.5000clear all;a=reshape(1:16,4,4);b=max(a),c=max(max(a)),b =4 8 12 16c =16[x y]=max(a)x =4 8 12 16y =4 4 4 4【例】hist(histogram)指令的使用示例。
randn('state',1),rand('state',31)x=randn(1000,1);y=rand(1000,1);%图 3.2-13.3函数的数值导数(1) 数值差分和导数、偏导数dx=diff(x) ％求X相邻行元素间的一阶差分。
dx=diff(x,n) ％求X相邻元素间的n阶差分。
dx=diff(x,n,dim) ％在dim指定的维上，求X相邻元素间的n阶差分。
note：1、diff是基于前向差分概念设计的，即dx/dt=(x(t+h)-x(t))/h2、数值导数的求取应尽量避免。
3、一般需先通过多项式拟合，或通过样条拟合，然后再从元数据拟合函数求导。
4、利用多项式求导。
【例3.3-1】y=sin(x)的导数为cos(x)，利用数值求导及多项式拟合求导，并进行比较。
x=0:pi/10:2*pi; y=sin(x);plot(x,y);hold on;y1=diff(y);plot(x(2:end),y1);plot(x,cos(x),':r');p=polyfit(x,y,5);p1=polyder(p);y2=polyval(p1,x);p（2）数值梯度[fx fy]=gradient(F,h)note:1、数值梯度使用的场合与数值导数相同2、MATLAB约定：数组的“行”数据点沿x轴取得，“列”数据点沿y轴取的。
3、h为步长。
（3）方向导数的可视化quiver(x,y,u,v,scale) 在（x,y）二维平面点上，画（u,v）方向的箭头。
scale 表示箭头的长度，默认值为1。
【例3.3-2】用一个简单矩阵表现diff和gradient指令计算方式。
F=[1,2,3;4,5,6;7,8,9]Dx=diff(F)Dx_2=diff(F,1,2)[FX,FY]=gradient(F)[FX_2,FY_2]=gradient(F,0.5)【例 3.3-3】研究偶极子的电势和电场强度。
设在),(b a 处有电荷q +，在),(b a --处有电荷q -。
那么在电荷所在平面上任何一点的电势和场强分别为)11(4),(0-+-=r r q y x V πε，V E -∇= 。
其中2222)()(,)()(b y a x r b y a x r +++=-+-=-+。
9010941⋅==πεk 。
又设电荷6102-⋅=q ，5.1=a ，5.1-=b 。
clear;clf;q=2e-6;k=9e9;a=1.5;b=-1.5;x=-6:0.6:6;y=x;[X,Y]=meshgrid(x,y);rp=sqrt((X-a).^2+(Y-b).^2);rn=sqrt((X+a).^2+(Y+b).^2); V=q*k*(1./rp-1./rn); [ex,ey]=gradient(-V);ae=sqrt(ex.^2+ey.^2);ex=ex./ae;ey=ey./ae; cv=linspace(min(min(V)),max(max(V)),49); contourf(X,Y,V,cv,'k-') axis('square')title('\fontname{隶书}\fontsize{22}偶极子的场'),hold on quiver(X,Y,ex,ey,0.7) plot(a,b,'wo',a,b,'w+')plot(-a,-b,'wo',-a,-b,'w-')xlabel('x');ylabel('y'),hold off3.4 函数的零点和极点很多情况下，需要求函数数组子矩阵的零点和极点3.4.1 多项式的根出于计算考虑，Matlab 不对多项式直接求根，而是通过求他的伴随矩阵的特征值进行，对于阶数不超过20的多项式，这种处理方法被认为是最适当可靠的计算方法。
指令：roots(p) 求多项式P 的根，其中p 为多项式系数的行向量。
【例 3.4.1-1】利用多项式计算y=sin(x)，在零到2＊pi 区间内5阶多项式拟合，并计算此多项式的的零点。
clear all;x=0:pi/10:2*pi;y=sin(x); p=polyfit(x,y,5); roots(p)clear all;p=[1 -5 6]; poly2str(p,'x'),roots(p)一元函数的零点与多项式不同，任意函数f(x)=0可能有零点，也可能没有零点，可能有一个零点，也可能有多个零点，因此，很难说出一个通用解法。
一般来说，零点的数值计算过程是，先猜测一个初始零点或该令点所在的区间，然后通过一些计算，是猜测值不断精确化，或是猜测区间不断收缩，直到达到预先指定的精度，终止计算。
指令: z=fzero(fun,xo)note:1、该函数只能求去一元连续函数穿越横轴的零点，而不会确定连续函数曲线，那种子触及横轴而不穿越横轴的零点。
例如abs(sin(x))。
2、fun 为字符串或内联函数（inline function ）4、x0时表示零点初始猜测，此点应该靠近零点5、无法判断x0时，可利用plot指令。
Example 1Calculate π by finding the zero of the sine function near 3.x = fzero(@sin,3)x =3.1416Example 2To find the zero of cosine between 1 and 2x = fzero(@cos,[1 2])x =1.5708Note that cos(1) and cos(2) differ in sign.Example 3To find a zero of the function f(x) = x3–2x–5, write an anonymous function f:f = @(x)x.^3-2*x-5;Then find the zero near 2:z = fzero(f,2)z =2.0946Because this function is a polynomial, the statement roots([1 0 -2 -5]) finds the same real zero, and a complex conjugate pair of zeros.2.0946-1.0473 + 1.1359i-1.0473 - 1.1359iIf fun is parameterized, you can use anonymous functions to capture the problem-dependent parameters. For example, suppose you want to minimize the objective function myfun defined by the following function file:function f = myfun(x,a)f = cos(a*x);Example 4求sin(x)靠近pi的零点z=fzero('sin(x)',3)Example 5求sin(x)/x的主瓣上的零点clear all;x=-2*pi:pi/10:2*pi;y=sin(x+eps)./(x+eps);figuregrid ;holdplot(x,y);hold off[a b]=ginput(2);fs=inline('sin(x)/x','x');z1=fzero(fs,a(1));z2=fzero(fs,a(2));3.4.2函数极值点许多科学研究和工程计算问题口可以归结为一个极值问题,如能量最小,时间最短,最佳拟合等等。
(完整word版)含答案《MATLAB实用教程》
第二章 MATLAB 语言及应用实验项目实验一 MATLAB 数值计算三、实验内容与步骤1.创建矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a（1（2）用（3）用（42.矩阵的运算（1）利用矩阵除法解线性方程组。
⎪⎪⎩⎪⎪⎨⎧=+++=-+-=+++=+-12224732258232432143214321421x x x x x x x x x x x x x x x 将方程表示为AX=B ，计算X=A\B 。
（2）利用矩阵的基本运算求解矩阵方程。
已知矩阵A 和B 满足关系式A -1BA=6A+BA ，计算矩阵B 。
其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=7/10004/10003/1A ，Ps: format rata=[1/3 0 0;0 1/4 0;0 0 1/7];b=inv(a)*inv(inv(a)-eye(3))*6*a（3）计算矩阵的特征值和特征向量。
已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=1104152021X ，计算其特征值和特征向量。
（4）Page:322利用数学函数进行矩阵运算。
已知传递函数G(s)=1/(2s+1),计算幅频特性Lw=-20lg(1)2(2w )和相频特性Fw=-arctan(2w),w 的范围为[0.01，10]，按对数均匀分布。
3.多项式的运算（1）多项式的运算。
已知表达式G(x)=(x-4)(x+5)(x 2-6x+9),展开多项式形式，并计算当x 在[0，20]内变化时G(x)的值，计算出G(x)=0的根。
Page 324（2）多项式的拟合与插值。
将多项式G(x)=x 4-5x 3-17x 2+129x-180，当x 在[0，20]多项式的值上下加上随机数的偏差构成y1，对y1进行拟合。
对G(x)和y1分别进行插值，计算在5.5处的值。
Page 325 四、思考练习题1.使用logspace 函数创建0～4π的行向量，有20个元素，查看其元素分布情况。
Ps: logspace(log10(0),log10(4*pi),20) (2) sort(c,2) %顺序排列 3.1多项式1）f(x)=2x 2+3x+5x+8用向量表示该多项式，并计算f(10)值. 2）根据多项式的根[-0.5 -3+4i -3-4i]创建多项式。
matlab数值计算实验报告
matlab数值计算实验报告Matlab数值计算实验报告引言：Matlab是一种广泛应用于科学与工程领域的高级计算机语言和环境，它提供了丰富的函数库和工具箱，方便用户进行数值计算、数据分析和可视化等任务。
本实验报告将介绍我在使用Matlab进行数值计算实验中的一些经验和心得体会。
一、数值计算方法数值计算方法是一种利用数值近似来解决实际问题的方法，它在科学和工程领域具有广泛的应用。
在Matlab中，我们可以利用内置的函数和工具箱来实现各种数值计算方法，例如插值、数值积分、数值微分等。
二、插值方法插值是一种通过已知数据点来推测未知数据点的方法。
在Matlab中，我们可以使用interp1函数来进行插值计算。
例如，我们可以通过已知的一些离散数据点，利用interp1函数来估计其他位置的数值。
这在信号处理、图像处理等领域具有重要的应用。
三、数值积分数值积分是一种通过分割曲线或曲面来近似计算其面积或体积的方法。
在Matlab中，我们可以使用quad函数来进行数值积分计算。
例如，我们可以通过quad函数来计算某个函数在给定区间上的积分值。
这在概率统计、物理学等领域具有广泛的应用。
四、数值微分数值微分是一种通过数值逼近来计算函数导数的方法。
在Matlab中，我们可以使用diff函数来进行数值微分计算。
例如，我们可以通过diff函数来计算某个函数在给定点上的导数值。
这在优化算法、控制系统等领域具有重要的应用。
五、数值求解数值求解是一种通过数值近似来计算方程或方程组的根的方法。
在Matlab中，我们可以使用fsolve函数来进行数值求解计算。
例如，我们可以通过fsolve函数来求解某个非线性方程的根。
这在工程计算、金融分析等领域具有广泛的应用。
六、实验应用在本次实验中，我使用Matlab进行了一些数值计算的应用实验。
例如，我利用插值方法来估计某个信号在给定位置的数值，利用数值积分方法来计算某个曲线下的面积，利用数值微分方法来计算某个函数在给定点的导数值，以及利用数值求解方法来求解某个方程的根。
MATLAB上机实验实验报告
MATLAB上机实验一一、实验目的初步熟悉MATLAB 工作环境，熟悉命令窗口，学会使用帮助窗口查找帮助信息。
命令窗口二、实验内容(1) 熟悉MATLAB 平台的工作环境。
(2) 熟悉MATLAB 的5 个工作窗口。
(3) MATLAB 的优先搜索顺序。
三、实验步骤1. 熟悉MATLAB 的5 个基本窗口①Command Window (命令窗口)②Workspace (工作空间窗口)—③Command History (命令历史记录窗口)④Current Directory (当前目录窗口)⑤Help Window (帮助窗口)(1) 命令窗口(Command Window)。
在命令窗口中依次输入以下命令：>>x=1>> y=[1 2 34 5 67 8 9];>> z1=[1:10],z2=[1:2:5];>> w=linspace(1,10,10);>> t1=ones(3),t2=ones(1,3),t3=ones(3,1)>> t4=ones(3),t4=eye(4)x =1z1 =1 2 3 4 5 6 7 8 9 10 t1 =1 1 11 1 11 1t2 =1 1 1t3 =111t4 =1 1 11 1 11 1 1t4 =1 0 0 00 1 0 00 0 1 00 0 0 1思考题：①变量如何声明，变量名须遵守什么规则、是否区分大小写。
答：（1）变量声明1.局部变量每个函数都有自己的局部变量，这些变量只能在定义它的函数内部使用。
当函数运行时，局部变量保存在函数的工作空间中，一旦函数退出，这些局部变量将不复存在。
脚本（没有输入输出参数，由一系列MATLAB命令组成的M文件）没有单独的工作空间，只能共享调用者的工作空间。
当从命令行调用，脚本变量存在基本工作空间中；当从函数调用，脚本变量存在函数空间中。
2.全局变量在函数或基本工作空间内，用global声明的变量为全局变量。
MATLAB数值计算.doc
MATLAB数值计算第5章 MATLAB数值计算 MATLAB的数学计算,数值计算，符号计算其中符号计算是指使用未定义的符号变量进行运算，而数值计算不允许使用未定义的变量。
注:数值计算有误差。
数值计算受到计算机所保留的有效位数的限制，因此每一次运算都会有一定的误差，重复的多次数值计算就可能会造成很大的累积误差。
5.1 矩阵和数组MATLAB最基本也是最重要的功能就是进行实数或复数矩阵的运算。
矩阵和数组的概念在MATLAB的运算中，经常要使用标量、向量、矩阵和数组，这几个名称的定义如下: , 标量:是指1×1的矩阵，即为只含一个数的矩阵。
, 向量:是指1×n或n×1的矩阵，即只有一行或者一列的矩阵。
, 矩阵:是一个矩形的数组，即二维数组，其中向量和标量都是矩阵的特例，0×0矩阵为空矩阵([])。
, 数组:是指n维的数组，为矩阵的延伸，其中矩阵和向量都是数组的特例。
5.1.1矩阵输入(前已讲“向量的输入”，相同～) (1) 矩阵元素应用方括号([])括住;(2) 每行内的元素间用逗号或空格隔开;(3) 行与行之间用分号或回车键隔开;(4) 元素可以是数值或表达式。
1. 通过显式元素列表输入矩阵c=[1 2;3 4;5 3*2] % [ ]表示构成矩阵,分号分隔行,空格分隔元素c =1 23 45 6用回车键代替分号分隔行:c=[1 23 45 6]1 23 45 612. 通过语句生成矩阵(1) 使用from:step:to方式生成向量(2) 使用linspace和logspace函数生成向量3. 由矩阵生成函数产生特殊矩阵MATLAB提供了很多能够产生特殊矩阵的函数，各函数的功能如表5.2所示。
表5.2 矩阵生成函数例子函数名功能输入结果zeros(2,3) zeros(m,n) 产生m×nans = 的全0矩 0 0 0 阵 0 0 0ones(2,3) ones(m,n) 产生m×nans = 的全1矩1 1 1 阵 1 1 1rand(2,3) rand(m,n) 产生均匀ans = 分布的随 0.9501 0.6068 0.8913 机矩阵，元 0.2311 0.4860 0.7621 素取值范围0.0,1.0。
实验五+MATLAB数值计算(含实验报告)
实验五 MATLAB 数值计算一、实验目的1．掌握求数值导数和数值积分的方法。
2．掌握代数方程数值求解的方法。
3．掌握常微分方程数值求解的方法。
二、实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。
设计提示1．参考本节主要内容，学习并理解相关函数的含义及调用方法。
三、实验内容1．线性系统方程：分别使用左除（\）和求逆（inv ）求解下面系统方程的解：⎪⎩⎪⎨⎧=+=+=++377251463c b b a c b a2. 数值积分：使用quad 和trapz 求解⎰-503/dx xe x 的数值积分，并与其解析解9243/5+--e 相比较；3. 请完成教材P154页中实验指导环节的实验内容第2题4. 请完成教材P155页中思考练习的第3题（1），并绘制解在该求解区间（即[0,5]）上的图像；。
5、请完成教材P164页实验指导环节的实验内容第5题。
（提示：该函数的符号导数，可以通过函数diff 求得。
首先定义符号变表达式，如求sin(x)的一阶符号导数，可以先定义f=’sin(x)’;df=diff(f);可求得df=cos(x)。
其中df 即为函数f 的一阶符号导数）。
四、实验报告要求（包含预习报告要求和最终报告要求）1.实验名称2.实验目的3.实验设备及条件4.实验内容及要求5.实验程序设计指程序代码。
6.实验结果及结果分析实验结果要求必须客观,现象。
结果分析是对实验结果的理论评判。
7.实验中出现的问题及解决方法8. 思考题的回答五、实验报告的提交方式Word文档，命名方式:实验号_你的学号_姓名例如本次实验：实验一_000000001_张三.doc（信息101提交报告邮箱）：E_mail: *******************(网络工程101提交作业邮箱)：E_mail: *******************（注意网络班的M是大写的)下一次课前提交,过期不收!六、参考文献参考教材和Matlab帮助文件。
MATLAB实验上机所用word(字符串及其运算篇)
第四章字符串、元胞和构架数组与符号计算MATLAB 6.x版的内建数据类型（Built-in data type）就有5种以上，此外还有许多其他专门设计的类（Class），如符号类、内联函数类、控制工具包中的线性时不变模型类、神经网络类等。
就程序设计而言，MATLAB 6.x版采用了面向对象编程技术。
数据和编程的改变使用户能更简捷而自然地解决复杂的计算问题（如符号计算问题、多变量控制系统问题、神经网络问题）。
本章内容根据MATLAB6.5编写，但绝大部分内容适用于其他MATLAB6.x版本。
4.1字符串数组4.1.1字符串入门字符变量的创建方式是：在指令窗中，先把待建的字符方在“单引号对中，再按’Enter‘。
注意单引号对必须在英文状态输入。
【例3.1.1-1】先请读者实际操作本例，以体会数值量与字符串的区别。
cleara=12345.6789class(a)a_s=size(a)b='S'class(b)b_s=size(b)4.1.2串数组的属性和标识【例3.1.2-1】本例演示：串的基本属性、标识和简单操作。
（1）创建一个由19个字符组成的串，这19个字符必须被放在但引号对内。
a='This is an example.'（2）串数组的大小在以上赋值后，变量就是一个串数组。
该串的每个字符（英文字母、空格和标点都是平等的）占据一个元素位。
串数组的大小可以通过下面指令获得。
size(a)（3）串数组的元素标识在一维串数组中，matlab按自左至右的次序用自然数数码标识字符的位置。
a14=a(1:4) ％提出一个字数串ra=a(end:-1:1) ％字符串的倒排（4）串数组的ASCII码字符串的存储是用ASCII码实现的，指令double、abs都可以用来获取串数组对应的ASCII码数值数组，指令char可以把ASCII数组变为串数组。
ca=double(a)char(ca)note: 1、在char把数字转换为字符时，负数将导致“警告”；非整数部分将被截尾。
matlab上机实验报告
《MATLAB 语言及其应用实验》上机实验报告1姓名： 学号：请各位同学将命令和结果由MATLAB 中复制（或截屏）后保存至word 文档中。
1. 计算矩阵 与 的和、乘积、各元素的乘积以及Kronecker 积。
2.对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
3.对矩阵912563827⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦进行特征值分解，给出特征值与特征向量。
4.已知向量[4＋5i, 6, 7-1.5j, 5+9j]，以其为对角线元素生成对角矩阵。
5. 生成含有8个元素的高斯向量a，要求元素的均值、方差分别为5、4，再生成含有7个元素的随机向量，要求元素服从上的均匀分布。
求两者的卷积。
6. 求解多项式x3-7x2+2x+40的根。
7. 计算多项式乘法(x2+2x+2)(x2+5x+4)。
8. 对下式进行部分分式展开：9. 已知y=sin(x)，x从0到2 ，间隔 x=0.02 ，求y的最大值、最小值、均值和标准差。
10. 用符号函数法求解方程at2+b*t+c=0。
11. 用符号计算验证三角等式：sin( 1)cos( 2)-cos( 1)sin( 2) =sin( 1- 2)12. 因式分解：13. ⎥⎥⎦⎤⎢⎢⎣⎡=)sin()log(12x x e x x a f ax ，用符号微分求df/dx ，并用符号定积分来验证。
14. 求解非线性方程组1211223sin 2cos x x x x x x -=⎧⎨+=⎩。
>> syms x1 x2>> [x1,x2]=solve(x1-3*x2-sin(x1),2*x1+x2-cos(x2))x1 =0.49662797440907460178544085171994x2 =0.0067214622395756734146654770697884。
《MATLAB与数值分析》第一次上机实验报告
电子科技大学电子工程学院标准实验报告（实验）课程名称MATLAB与数值分析学生姓名：李培睿学号：2013020904026指导教师：程建一、实验名称《MATLAB与数值分析》第一次上机实验二、实验目的1. 熟练掌握矩阵的生成、加、减、乘、除、转置、行列式、逆、范数等运算操作。
（用.m文件和Matlab函数编写一个对给定矩阵进行运算操作的程序）2. 熟练掌握算术符号操作和基本运算操作，包括矩阵合并、向量合并、符号转换、展开符号表达式、符号因式分解、符号表达式的化简、代数方程的符号解析解、特征多项式、函数的反函数、函数计算器、微积分、常微分方程的符号解、符号函数的画图等。
（用.m文件编写进行符号因式分解和函数求反的程序）3. 掌握Matlab函数的编写规范。
4、掌握Matlab常用的绘图处理操作，包括：基本平面图、图形注释命令、三维曲线和面的填充、三维等高线等。
（用.m文件编写在一个图形窗口上绘制正弦和余弦函数的图形，并给出充分的图形注释）5. 熟练操作MATLAB软件平台，能利用M文件完成MATLAB的程序设计。
三、实验内容1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。
并以x，y为坐标显示图像x(n+1) = a*x(n)-b*(y(n)-x(n)^2);y(n+1) = b*x(n)+a*(y(n)-x(n)^2)2. 编程实现奥运5环图，允许用户输入环的直径。
3. 实现对输入任意长度向量元素的冒泡排序的升序排列。
不允许使用sort 函数。
四、实验数据及结果分析题目一：①在Editor窗口编写函数代码如下：并将编写的函数文件用“draw.m”储存在指定地址；②在Command窗口输入如下命令：③得到图形结果如下：题目二：①在Editor窗口编写函数代码如下：并将编写的函数文件用“circle.m”储存在指定地址；②再次在Editor窗口编写代码：并将编写的函数文件用“Olympic.m”储存在指定地址；③在Command窗口输入如下指令（半径可任意输入）：④按回车执行，将在图形窗口获得五环旗：题目三：①在Editor窗口编写函数代码如下：并用.将编写的函数文件用“qipaofa.m”储存在指定地址；②在Command窗口输入一组乱序数值，则可以得到升序排序结果如下：五、总结及心得体会1.要熟悉MATLAB编译软件的使用方法，明白有关语法，语句的基本用法，才可以在编写程序的时候游刃有余，不至于寸步难行。
(完整word)Matlab实验报告
实验一:Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境.2．学习使用图形函数计算器命令funtool及其环境。
二、实验内容熟悉Matlab操作环境,认识命令窗口、内存工作区窗口、历史命令窗口;学会使用format 命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息;学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1)y=sin（x)；（2） y=(1+x)^3*（2-x）；求下列函数的符号积分(1)y=cos(x）;(2)y=1/（1+x^2）；（3）y=1/sqrt（1—x^2)；(4）y=(x1)/(x+1）/(x+2)求反函数(1)y=（x-1)/(2＊x+3）; （2） y=exp(x）；（3） y=log（x+sqrt（1+x^2));代数式的化简(1）(x+1)*（x-1）*(x-2）/（x-3）/（x—4)；(2）sin（x）^2+cos(x）^2；（3）x+sin(x）+2*x—3*cos(x）+4＊x*sin（x）;2．函数与参数的运算操作。
从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1）^2(2）y2=（x+2)^2(3) y3=2＊x^2 (4） y4=x^2+2 (5） y5=x^4 (6） y6=x^2/2 3．两个函数之间的操作求和(1)sin（x)+cos（x) (2） 1+x+x^2+x^3+x^4+x^5乘积(1)exp(—x）＊sin（x） (2） sin(x)＊x商(1)sin(x）/cos（x); （2) x/（1+x^2); （3) 1/(x—1)/(x—2）; 求复合函数(1)y=exp（u) u=sin（x） (2） y=sqrt（u) u=1+exp(x^2）(3） y=sin（u） u=asin（x） (4) y=sinh(u) u=-x实验二：MATLAB基本操作与用法一、实验目的1.掌握用MATLAB命令窗口进行简单数学运算。
《MATLAB及应用》实验报告1 MATLAB数值计算
电气工程学院实验报告实验项目名称MATLAB数值计算所属课程名称MATLAB及应用实验类型上机实验实验日期2014-03- 5指导教师XXXXXX班级XXXXXXXXX学号XXXXXXXXXX 姓名XXXXXXXX成绩一、实验名称MATLAB数值计算二、实验目的（1）掌握MATLAB变量的使用（2）掌握MATLAB数组的创建（3）掌握MATLAB数组和矩阵的运算（4）熟悉MATLAB多项式的运用三、实验原理1. 矩阵分析矩阵转置：单引号（’）矩阵的旋转：rot90(A,k)，功能是将矩阵A旋转90度的k倍，缺省值是1 矩阵的左右翻转：fliplr(A)矩阵的上下翻转：flipud(A)矩阵的逆：inv(A)，与A^(-1)等价矩阵的行列式：det(A)矩阵的秩：rank(A)矩阵的迹：trace(A)将矩阵化为最简式：rref(A)矩阵的特征值与特征向量：(1) E = eig(A)；矩阵A的所有特征值构成向量E；(2) [V,D]=eig(A)；A的所有特征值构成对角阵D，A的特征向量构成V的列向量；2. 多项式多项式的建立：若多的项的全部根构成的向量为X，则以X为根的多项式为poly(X)多项式的根：roots(p)计算以向量p为系数的多项式的根，包括重根，复根多项式求值：polyval(p,x)，p是多项式的系数，x可以是一个数也可以是一个矩阵多项式求拟合次数：polyfit(x,y,n)，x可以是一个数也可以是一个矩阵，y 是x对应的数或矩阵多项式的四则运算：(1)P1+P2；(2)P1-P2；(3)conv(P1,P2), (4)deconv(P1,P2)四、实验内容1. 已知矩阵11 12 13 1421 22 23 24A= 31 32 33 344142 43 44(1) A(:,1) (2) A(2,:) (3) A(:,2:3)(4) A(2:3,2:3) (5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:) (9) ones(2,2)(10) eye(2) (11) [A,[ones(2,2);eye(2)]] (12) diag(A)(13) diag(A,1) (14) diag(A,-1) (15) diag(A,2)2. 已知2354127354217623A-⎡⎤⎢⎥--⎢⎥=⎢⎥-⎢⎥--⎣⎦.求（1）A的逆；（2）A的行列式；（3）A的迹；（4）A的所有特征向量和特征值。