MATLAB模拟试题

一、填空题

1. MA TLAB 命令窗口中可用 clear 命令清除工作区中的变量；用 clc 命令清除命令窗口中的内容。

2. MA TLAB 中的运算符包括数学运算符、 、 。

3 MA TLAB 中的M 文件有 命令文件 和 函数文件 。

4. MA TLAB 中的函数类型有嵌套函数、 、 、私有函数和匿名函数。

5. MA TLAB 中用系统预定义的特殊变量 ans 表示函数输入参数个数。

6. MA TLAB 中用系统预定义的特殊变量 inf 表示无穷大。

7. MA TLAB 中用 符号定义单元数组。

8. MA TLAB 中用 命令打开图形保持功能。

9 MA TLAB 中如果一个变量名为123，该变量名是否合法 否。

10. MA TLAB 中用 函数完成矩阵的求逆运算，用 函数求矩阵的行列式。

二、阅读下述MA TLAB 脚本程序

A=[3 5 7 9;0 2 4 1;6 4 4 8];

B=ones(3,3);

b=A([1 2],2:3);

B(4,4)=5;

A([1,2],:)=[];

C=[A;B(:,1:3)'];

D=reshape(C,2,8);

写出b 、A 、B 、C 、D 的结果。

三、已知矩阵A 为4×4的正态分布的随机矩阵，编写脚本程序实现矩阵B, B 矩阵为对角矩阵，主对角线元素为A 矩阵的主对角线元

素加1，并将矩阵B 存入文件result.mat 文件中。

四、已知多项式532)(3++=x x x a 和13)(2

4++=x x

x b ，编写脚本程序计算 1） 多项式a(x)+b(x)

2） a(x)*b(x)

2） 多项式a(x)/b(x)的部分分式展开。 五、用MA TLAB 编写函数，计算()?????<≤-<≤=其他

021210x f x x

x x

，函数的输入参数为x ，输出参数为f ，函数名为fun 。 六、编写MATLAB 程序，以0.2为步长绘制正弦曲线()t πsin 和余弦曲线()t πcos 在[0,2π]之间的图形，分别位于图形窗口的二个子图中，正弦曲线以红色虚线绘制，并且用星号标注每个节点，余弦曲线以红色实线绘制，用圆标注每个节点。

———————————————————————————————————————

一、填空题

1. MA TLAB 中用 函数查找数组中的非零元素，用 函数对数组进行排序。

2. MA TLAB 中的M 文件包括 和 。

3. MA TLAB 中用系统预定义的特殊变量 表示浮点数的精度。

4. MA TLAB 中8位有符号整数用 定义。

5. MA TLAB 中小括号、转置运算符和冒号运算符，具有最高优先级的为 。

6. MA TLAB 中用系统预定义的特殊变量 表示函数输出参数个数。

7. MA TLAB 中用 符号表示其后为程序注释。

8. 已知矩阵A,则矩阵A 的共轭转置可表示为 。

9. MA TLAB 中用 命令显示坐标轴下的网格线。

10. MA TLAB 中如果一个函数名为test ，该函数名是否合法 （填是或否）。

11. MA TLAB 中用 创建一个空数组。

12. MA TLAB 命令窗口中可用 命令查看当前路径；用 命令查看搜索路径。

二、阅读下述MA TLAB 脚本程序

A=[1 2 3 4;5 6 7 8;9 0 0 0];

a=diag(A);

b=A(:,[1 3 4]);

A(:,3)=[];

B=repmat(A,1,2);

c=length(B);

写出a 、b 、c 、A 、B 的结果。

三、已知矩阵A 为4×4的均匀分布的随机矩阵，编写脚本程序实现矩阵B, B 矩阵由A 矩阵的第1、3列，第1、2行元素构成，并将矩阵B 存入文件result.mat 文件中。

四、已知矩阵X = [3 2 0; -5 0 7; 0 0 1]，编写脚本程序查找元素值大于2小于9的元素向量及其下标向量，并对矩阵X 按列维方向升序排序。

五、已知多项式12)(23++=x

x x a 和1)(2

+=x x b ，编写脚本程序计算： 1） 多项式a(x)+b(x)

2） 多项式a(x)*b(x)。

3） 多项式a(x)的根。

六、用MA TLAB 编写函数，当函数输入为一个参数时，输出为输入的平方；当函数输入为2个参数时，输出为输入之和。函数定义为：

function out=test （in1，in2）

七、编写MATLAB 程序，计算Fibonacci 数列中的第一个大于10000的元素及n 值。Fibonacci 定义为：

()???>-+-≤=2)2()1(21n f n n f n f n

八、编写MA TLAB 程序，以0.2为步长绘制均值为0，标准差分别为1、2和3的正态分布概率密度函数在[-3,3]之间的图形。

MATLAB数学实验第二版答案(胡良剑)

数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点 >> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans = 1.2500

数学实验练习题(MATLAB)

注意:在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上）. 第一次练习题 1.求解下列各题： 1)30sin lim x mx mx x ->- 2)(4)cos ,1000.0=x mx y e y 求 3）21/2 0mx e dx ?(求近似值，可以先用inline 定义被积函数，然后用quad 命令) 4）4 224x dx m x +? 5 0x =展开(最高次幂为8). 2.对矩阵21102041A m -?? ?= ? ?-?? ，分别求逆矩阵，特征值，特征向量，行列式，并求矩阵,P D （D 是对角矩阵），使得1A PDP -=。 3. 已知2 1(),()2f x e x μσ=--分别在下列条件下画出)(x f 的图形： (1)/600m σ=,μ分别为0,1,1-(在同一坐标系上作图); (2)0μ=,σ分别为1,2,4,/100m (在同一坐标系上作图). 4.画 （1）sin 020cos 02100x u t t y u t u t z m ??=≤≤?=?≤≤??=?

(2) sin()03,03z mxy x y =≤≤≤≤ （3）sin()(/100cos )02cos()(/100cos )02sin x t m u t y t m u u z u π π=+?≤≤?=+?≤≤?=? 的图（第4题只要写出程序）. 5．对于方程50.10200 m x x --=，先画出左边的函数在合适的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会. 第二次练习题 判断迭代收敛速度的程序 x0=1;stopc=1;eps=10^(-8);a=1;c=1;b=2*c;d=a;k=0; f=inline('(a*x+b)/(c*x+d)'); kmax=100; while stopc>eps&k

MATLAB数学实验练习题

MATLAB数学实验练习题

————————————————————————————————作者：————————————————————————————————日期：

“MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) > ezplot('exp(x)-3*x^2') >> grid on ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2）1sin 02 x x - =至 少三个根

>> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2sin cos 0x x x -= 所有根

>> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3,x,0) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x)

MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案

“”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;

>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =

0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1

6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5

8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =

matlab数学实验练习题

Matlab 数学实验 实验一 插值与拟合 实验内容： 预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。 （1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,c o s 10 ≤≤-=x x y （4）、22),exp(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(t e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面 一组t ，V 数据确定0V 和τ。 实验二 常微分方程数值解试验 实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。 实验内容：

实验三地图问题 1.下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量： 以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，这样就得到了表中的数据（单位mm）。 根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土 的近似面积，并与它的精确值41288km2比较。

matlab图像处理综合实验实验报告

《数字图像处理》 实验报告 学院： 专业： 班级： 姓名： 学号： 实验一 实验名称：图像增强 实验目的：1.熟悉图像在Matlab下的读入，输出及显示； 2.熟悉直方图均衡化； 3.熟悉图像的线性指数等； 4.熟悉图像的算术运算及几何变换. 实验仪器：计算机，Matlab软件 实验原理： 图像增强是为了使受到噪声等污染图像在视觉感知或某种准则下尽量的恢复到原始图像的水平之外，还需要有目的性地加强图像中的某些信息而抑制另一些信息，以便更好地利用图像。图像增强分频域处理和空间域处理，这里主要用空间域的方法进行增强。空间域的增强主要有：灰度变换和图像的空间滤波。 图像的直方图实际上就是图像的各像素点强度概率密度分布图，是一幅图像所有像素集合的最基本统计规律，均衡化是指在每个灰度级上都有相同的像素点过程。 实验内容如下： I=imread('E:\cs.jpg');%读取图像 subplot(2,2,1),imshow(I),title('源图像') J=rgb2gray(I)%灰度处理 subplot(2,2,2),imshow(J) %输出图像 title('灰度图像') %在原始图像中加标题 subplot(2,2,3),imhist(J) %输出原图直方图

title('原始图像直方图') I=imread('E:\cs.jpg');%读取图像 subplot(1,2,1),imshow(I); subplot(2,2,1),imshow(I),title('源图像') J=rgb2gray(I)%灰度处理 subplot(2,2,2),imshow(J),title('灰度变换后图像') J1=log(1+double(J)); subplot(2,2,3),imshow(J1,[]),title('对数变换后') 指数运算: I=imread('E:\dog.jpg'); f=double(I); g=(2^2*(f-1))-1 f=uint8(f); g=uint8(g); subplot(1,2,1);subimage(f),title('变换一') 00100200 源图像灰度变换后图像对数变换后

浅析Matlab数学实验报告

数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)

cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);

南邮MATLAB数学实验答案(全)

第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).

MATLAB上机实验练习题答案

数学建模 MATLAB上机实验练习题 1、给出一个系数矩阵A[234;541;132]，U=[123]，求出线性方程组的一个精确解。 2、给出两组数据x=[00.30.81.11.62.3]’y=[0.820.720.630.600.550.50]’，我们可以简单的 认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。 3、解线性方程 4、通过测量得到一组数据： 5、已知一组测量值 6、从某一个过程中通过测量得到： 分别采用多项式和指数函数进行曲线拟合。 7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近 插值法，双线性插值法，二重三次方插值法） 8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。分别使用不同的 颜色，线形和标识符。 9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额

占总销售额的比例，并分离第三种产品的切片。 X=19.322.151.6 34.270.382.4 61.482.990.8 50.554.959.1 29.436.347.0 10、对应时间矢量t，测得一组矢量y t00.30.8 1.1 1.6 2.3 y0.50.82 1.14 1.25 1.35 1.40 采用一个带有线性参数的指数函数进行拟合，y=a0+a1e-t+a2te-t，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。 11、请创建如图所示的结构数组（9分） 姓名编号指标 江明顺071023身高：176，体重：82 于越忠060134身高：168，体重：74 邓拓050839身高：182，体重：77 12、创建如图所示的元胞数组。（9分） 13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。 年份19901992199419961998200020022004200620082010 产量（万吨）75.99591.972105.711123.203131.669150.697179.323203.212226.505249.633256.344 14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。 time=[2.55.07.510.013.017.020.030.040.050.060.070.0] res=[0.290.560.771.051.361.522.002.272.813.053.253.56]

matlab 数学实验 迭代 _ 蛛网图(免积分)

数学实验—实验报告（免积分） 一、实验项目：Matlab实验三—迭代 二、实验目的和要求 a.熟悉MATLAB软件的用户环境，掌握其一般目的命令和MATLAB数组操作与 运算函数； b.掌握MATLAB软件的绘图命令，能够熟练应用循环和选择结构实现各种循环 选择功能； c.借助MATLAB软件的绘图功能，对函数的特性进行探讨，广泛联想，大胆猜 想，发现进而证实其中的规律。 三、实验内容 问题一：将方程53 x x x +-+=改写成各种等价的形式进行迭代 5210 观察迭代是否收敛，并给出解释。 问题二：迭代以下函数，分析其收敛性。 4 f(x)=x-a 使用线性连接图、蛛网图或费根鲍姆图对参数a进行讨论和观察，会得到什么结论？ 问题一： （1）画图 x1=-6:0.01:6; x2=-3:0.01:3; x3=-1:0.01:1; x4=-0.8:0.01:-0.75; y1=x1.^5 +5*x1.^3-2*x1+1; y2=x2.^5 +5*x2.^3-2*x2+1; y3=x3.^5 +5*x3.^3-2*x3+1; y4=x4.^5 +5*x4.^3-2*x4+1; subplot(2,2,1),plot(x1,y1) ,title('图(1)') ,grid on, subplot(2,2,2),plot(x2,y2) ,title('图(2)'),grid on, subplot(2,2,3),plot(x3,y3) ,title('图(3)'),grid on, subplot(2,2,4),plot(x4,y4) ,title('图(4)') ,grid on,

实验习题-Matlab

1、 计算9.248.26107sin 369.12÷?? ? ??π+的值 2、 产生一个5阶魔术方阵，并执行如下操作： （1） 将矩阵的第2行3列元素赋值给变量c （2） 将由矩阵第2，3，4行第3，5列构成的子矩阵赋值给变量d 3、给出区间[0，1]上的6个等分点数据。 4、建立如下矩阵 （1）10102007000200700 02007??? ? ? ? ?? ? （2）1010010101001010100??? ? ? ? ??? 5、设有分块矩阵?? ????=????2232233 3S O R E A ，其中E,R,O,S 分别为单位阵、随机阵、零阵和对角阵，试通过数值计算比较A*A 与A.*A 的区别。 实验二、Matlab 程序设计 1、 已知函数?? ???≤≤<≤<≤-+=2110, 101,1)(2x x x x x x f 计算)5.1(),5.0(),1(f f f -. 2、 用for-end 循环语句求：100！和∑=100 1i i 。 3、 用while-end 循环语句求不超过1000的偶数之和与奇数之和。 4、 建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立方和等于该数本身。例如，153是一个水仙花数，因为153=13+53+33。 5、建立如下矩阵 （1）10101000120011 100??? ? ? ? ??? （2）20072000002007300000200740000020075000002007?? ? ? ? ? ? ???

1、 在同一坐标系下面画出)5.0sin(2.01.0x e y x +=和)5.0cos(2.01.0x e y x +=在区间]2,0[π上的曲线图。 2、 绘制三维螺旋线：].,0[,5.0,sin 2,cos 2π∈===t t z t y t x 3、 画出曲面)sin(xy z =的网线图。 4、 画出曲面)(22y x xe z +-=的图形。 5、作出下列曲面的3维图形， 1） )sin(22y x z +=π； 2）环面：?????=+=+=,sin ,sin )cos 1(,cos )cos 1(u z v u y v u x )2,0()2,0(ππ∈∈v u 。 实验四、使用Matlab 解决微积分问题 1、 求下列函数的极限： （1） 4 202 cos lim x e x x x -→- （2） x x 1lim 0+→ 2、 按要求实现下面的求导运算： （1） 已知)tan()1ln(22x x e y x -+=，求)3(,y y '； （2） 已知xy y x e y x z 22)(22++=，求y x z x z x z ???????222,,。 3、 已知函数]3,2[,2sin )(2π∈=x x e x f x 。使用Matlab 软件，完成下面的实验任务： （1） 求出函数)(x f 的一阶导数，二阶导数，并画出它们相应的曲线。 （2） 观察函数的单调区间，凹凸区间，以及极值点和拐点。 4、 使用Matlab 软件，完成下列积分运算： （1） 求不定积分??+-1,232x x dx dx e x x ；

matlab综合实验4

实验4、输入输出的应用 目的与要求： 掌握matlab 图形的绘制 掌握自定义函数的使用 掌握matlab 基础运用 要求：请将程序和运行结果的截屏附在题目后面。 一、编程题 1、产生一个1x10的随机矩阵，大小范围为[-5 5]，并且按照从大到小的顺序排列好！ 2、设 y=cos[0.5+((3sinx)/(1+x^2))]把x=[0,2π]间分为101点，画出以x 为横坐标，y 为纵坐标的曲线；线条颜色为红色，线型为虚线。 3、在[-10，10；-10，10]范围内画出函数z 的三维网格图形。 4、将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。 ①]2 2[)cos(π π，，-∈x x （曲线图）； ②4)y 2,-4x (-242),(2222≤≤≤≤+=；y x y x f （曲面图）。 5、请修改下面的程序，让他们没有for 循环语句！修改后的程序与源程序得到相当的结果。 A=[1 2 3; 4 5 6; 7 8 9]; [r c]=size(A); for i=1:1:r for j=1:1:c if (A(i,j)>8 | A(i,j)<2) A(i,j)=0; end end

end 6、请自己定义函数，求出递归调用（即1*2*3*…*n） （1）求前n项积，n为自然数。（2）给出调用命令，如n=100时的结果。 二、填空题 1、已知a=[4 5 6],b=[7 8],[x,y]=meshgrid(a,b),则x= y= 。 2、定义x=[0 1 2 3 4]的用冒号表达式表示为。 3、在算术、关系、逻辑三种运算符中，优先级最高。 4、设x=[4 8 1;9 7 6],[a,b]=size(x)得到的结果是a= b= 。 5、x=find(A)命令的含义是。 6、元素取整的命令有、。 7、a=[1 2 2]，求a的转置的matlab命令。 8、b=[1 4 5 2 3 7 8 11 9 13 15 6]，求所有元素和的matlab 命令。 9、syms t k的含义。 10、在线形图型格式的设置中，字符y表示什么颜色，,栅格的命令是。 11、一个字符串由多个字符组成，用符号来界定。

matlab练习题和答案

matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233，，1443678，i，，,,2357,,,,2335542，i,,,, A,1357B,,,2675342， i,,3239,,,,189543，，,,1894，， 再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];

MATLAB数学实验报告

Matlab 数学实验报告

一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150

plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1];

for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图

2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏，长满了草。他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？ 2.2.2问题分析 如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积，2a÷π×πx2=2aπ2，再求AB的面积，用扇形ABE的面积减去三角形ABE的面积即可。

matlab综合大作业(附详细答案)

m a t l a b综合大作业(附详细 答案) 标准化文件发布号：（9312-EUATWW-MWUB-WUNN-INNUL-DQQTY-

《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）： 1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5) 实验结果： A = 2)将矩阵A按列拉长得到矩阵B； 实验程序：B=A(:) 实验结果： B =

3)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的 矩阵C； 实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)] 实验结果： C = 4)寻找矩阵A中大于0的元素；] 实验程序：G=A(find(A>0)) 实验结果： G = 5)求矩阵A的转置矩阵D； 实验程序：D=A' 实验结果： D = 6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E； 实验程序：E=flipud(fliplr(A)) 实验结果： E =

7)删除矩阵A的第2列和第4列得到矩阵F； 实验程序：F=A; F(:,[2,4])=[] 实验结果： F = 8)求矩阵A的特征值和特征向量； 实验程序：[Av,Ad]=eig(A) 实验结果： 特征向量Av = + - + - - + + - 特征值Ad = 0 0 0 0 0 + 0 0 0 0 0 - 0 0 0 0 0 0 0 0 0 0 9)求矩阵A的每一列的和值； 实验程序：lieSUM=sum(A) 实验结果： lieSUM = 10)求矩阵A的每一列的平均值； 实验程序：average=mean(A) 实验结果： average = 2．符号计算（10分，每小题5分）： 1)求方程组20,0 uy vz w y z w ++=++=关于,y z的解； 实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z'); y= S. y, z=S. z

Matlab数学实验一2015(标准答案版)

Matｌab数学实验一——matlab初体验 一、实验目的及意义 ［1] 熟悉MAＴＬＡB软件的用户环境； ［2] 了解MAＴLAB软件的一般目的命令； ［３] 掌握MATLＡB数组操作与运算函数； 通过该实验的学习,使学生能熟悉matｌab的基础应用，初步应用ＭATLAB软件解决一些简单问题。 二、实验内容 １.认识ｍatｌab的界面和基本操作 2.了解mａtlab的数据输出方式（ｆormat) 3. MATLAＢ软件的数组(矩阵）操作及运算练习; 三、实验任务 根据实验内容和步骤,完成以下具体实验，要求写出实验报告（实验目的→问题→原理→算法与编程→计算结果或图形→心得体会) 完成如下题目,并按照实验报告格式和要求填写实验报告 1．在ｃommanｄwinｄow中分别输入如下值,看它们的值等于多少，并用maｔｌａb的ｈelp中查询这些缺省预定义变量的含义,用中文写出它们的意义。 iｊeps inf nａn pi realｍaｘrealｍin 2.分别输入一个分数、整数、小数等，(如:a=1/9)，观察显示结果，并使用forｍaｔ函数控制数据的显示格式，如：分别输入ｆｏｒmat sｈorｔ、fｏｒmat loｎg、format short e、forｍａt lｏng g、foｒmａt bank、forｍａt heｘ等，然后再在命令窗口中输入ａ,显示a的值的不同形式，并理解这些格式的含义。 3．测试函数clｅａr、clc的含义及所带参数的含义（利用ｍatlａb的ｈelp功能）。 4. 写出在命令窗口中的计算步骤和运行结果。 （1)计算 1.22 10 (ln log) 81 e ππ +- ； >＞(log(pi)+lｏｇ(pi)/lｏg(10)-ｅxp(1.２)）^２/81 >>ａns = ０．０3４8 (2) >> x=2;y＝４； >> z=x^２+exｐ(x+y)-y*lｏg(x)-3 z ＝ 40１.65６2 （3)输入变量 13 5.3, 25 a b ?? ==?? ?? ，在工作空间中使用who，whos，并用ｓave命令将变量存入”D：\exe0 １.ｍat”文件。测试clear命令，然后用load命令将保存的”D：\ｅｘe０１.ｍat”文件载入＞> ａ＝5.3 ａ＝

matlab数学实验练习题

Matlab 数学实验 实验一 插值与拟合 实验内容： 预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。 （1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,cos 10≤≤-=x x y （4）、22),ex p(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(τt e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面 一组t ，V 数据确定0V 和τ。 实验二 常微分方程数值解试验 实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。 实验内容： 实验三 地图问题 1. 下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量：以由西向东方向为x 轴，由南到北方向为y 轴，选择方便的原点，并将从最西边界点到最东边界点在x 轴上的区间适当地划分为若干段，在每个分点的y 方向测出南边界点和北边界点的y 坐标y1和y2，这样就得到了表中的数据（单位mm ）。

根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土 2 实验四狼追兔问题 狼猎兔问题是欧洲文艺复兴时代的著名人物达.芬奇提出的一个数学问题。当一个兔子正在它的洞穴南面60码处觅食时，一只恶狼出现在兔子正东的100码处。当两只动物同时发现对方以后，兔子奔向自己的洞穴，狼以快于兔子一倍的速度紧追兔子不放。狼在追赶过程中所形成的轨迹就是追击曲线。狼是否会在兔子跑回洞穴之前追赶上兔子？ 为了研究狼是否能够追上兔子，可以先考虑求出狼追兔子形成的追击曲线，然后根据曲线来确定狼是否能够追上兔子。 试验五：开放式基金的投资问题 某开放式基金现有总额为15亿元的资金可用于投资，目前共有8个项目可供投资者选择。每个项目可以重复投资，根据专家经验，对每个项目投资总额不能太高，且有个上限。这些项目所需要的投资额已经知道，在一般情况下，投资一年后各项目所得利润也可估计出来（见表一）， 表一: 投资项目所需资金及预计一年后所得利润(单位：万元)

MATLAB数学实验A

clear; clc; a=1;b=1; ezplot(sprintf('x^2/%f-y^2/%f',a^2,b^2)); hold on; ezplot(sprintf('x^2/%f-y^2/%f-1',a^2,b^2)); ezsurf('sin(a)*cos(b)','sin(a)*sin(b)','cos(a)',[0,pi,0,2*pi],60); hold on; ezsurf('x^2+y^2',[-1,1,-1,1],60);

clear all; x=-8:0.1:8; y=-8:0.1:8; [X,Y]=meshgrid(x,y); Z=sin(sqrt(X.^2+Y.^2))./sqrt(X.^2+Y.^2+2); [X,Y,Z]=peaks(50); surf(X,Y,Z)

syms x y; y=2*x^3-6*x^2-18*x+7; solve(diff(y,x),x) x=-1;eval(y) x=3;eval(y)

syms x y; z='x*y'; dblquad(z,1,4,-1,2) 结果 ans = 11.2500 求函数1+x -exp(2*x)+5的原函数clear all syms x C; f=int(1+x -exp(2*x)+5,'x')+C syms x y; >> x=0:0.01:1; >> y=sin(sin(x)); >> trapz(x,y)

x=0:0.05:1; y=[1.97687 2.17002 2.34158 2.46389 2.71512 3.06045 3.27829 3.51992 3.8215 4.2435 4.55188 4.88753 5.15594 5.698 6.04606 6.42701 7.00342 7.50192 7.89178 8.49315 9.0938] cftool 解常微分方程y’=-0.9y/(1+2x)的数值解y(0)=1 从0到0. 1的数值解，取步长0.02 clear all x1=0; x2=0.1; h=0.02; y(1,1)=1;

MATLAB验练习题(计算机) 南邮 MATLAB 数学实验大作业答案

“MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6)

1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408

3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）3 0sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans =