Matlab实验
MATLAB实验报告
学校：湖北文理学院
学院：物理与电子工程学院
专业：电子信息工程
学号： 2013128182 姓名：张冲
指导教师：宋立新
实验一 MATLAB环境的熟悉与基本运算
一、实验目的：
1．熟悉MATLAB开发环境
2．掌握矩阵、变量、表达式的各种基本运算
二、实验内容
1、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，
学习使用指令eye（其它不会用的指令，依照此方法类推）
2、学习使用clc、clear，观察command window、command history和workspace
等窗口的变化结果。
3、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、
exerc2、exerc3……），学习使用MATLAB的基本运算符。
三、练习
1)help rand，然后随机生成一个2×6的数组，观察command window、
command history和workspace等窗口的变化结果。
2)学习使用clc、clear，了解其功能和作用。
3)用逻辑表达式求下列分段函数的值
4)求[100,999]之间能被21整除的数的个数。（提示：rem,sum的用法）
四、实验结果
1）
2）clc:清除命令窗口所有内容，数值不变；clear:初始化变量的值。3）
4）
实验二 MATLAB数值运算
一、实验目的
1、掌握矩阵的基本运算
2、掌握矩阵的数组运算
二、实验内容
1)输入C=1:2:20，则C（i）表示什么？其中i=1,2,3, (10)
2)输入A=[7 1 5；2 5 6；3 1 5]，B=[1 1 1; 2 2 2; 3 3 3]，在命令窗
口中执行下列表达式，掌握其含义：
A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2)
A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A
3)二维数组的创建和寻访，创建一个二维数组（4×8）A，查询数组A第2
行、第3列的元素，查询数组A第2行的所有元素，查询数组A第6列的所有
元素。
4)两种运算指令形式和实质内涵的比较。设有3个二维数组A
2×4，B
2×4
，C
2×2
，
写出所有由2个数组参与的合法的数组运算和矩阵指令。
5)学习使用表4列的常用函数（通过help方法）
6)学习使用表5数组操作函数。
7)生成一个3行3列的随机矩阵，并逆时针旋转90°，左右翻转，上下翻转。
8)已知a=[1 2 3]，b=[4 5 6]，求a.\b和a./ b
9)用reshape指令生成下列矩阵，并取出方框内的数组元素。
三、实验结果
1）C(i)表示C中的第i个的数值；
2）
3）
4）
5）略6）略7）8）
9）
10）数组和矩阵有何不同？
数组中的元素可以是字符等，矩阵中的只能是数，这是二者最直观的区别。从外观形状和数据结构上看，二维数组和数学中的矩阵没有区别。但是矩阵作为一种变换或映射算子的体现，矩阵运算有着明确而严格的数学规则。而数组运算是Matlab软件所定义的规则，其目的是为了数据管理方便、操作简单、指令形式自然和执行计算的有效。虽然数组运算尚缺乏严谨的数学推理，而且数组运算仍在完善和成熟中，但是它的作用和影响正随着matlab的发展而扩大。
实验三符号工具箱的应用
一、实验目的：
1、掌握创建符号表达式和矩阵的方法。
2、掌握符号函数极限及导数求解方法。
3、掌握符号表达式的微分和积分运算。
4、掌握数值插值与曲线拟合的方法及其应用
二、实验内容
1、求阶跃函数的傅里叶、拉普拉斯及Z 变换。
2、求下列变上限积分对变量x 的导数：
?
+2
x x
dx
x a 。
3、求解方程 1）??
?=+=+
133x 28
22121x x x
2）???
??=+=+=+3
43 23x 21221
2121x x x x x
3）??
?=++=++
243x 21
32321321x x x x x
4、用两种方法创建符号矩阵，A =[ a, 2*b]
[3*a, 0]
5、计算二重不定积分
6、对符号方程f = ax 2+bx+c 求解。1）对x 求解，2）对a 求解。
7、求解微分方程0)0(,1)0(,02222===++dx dy
y y dx dy dx
y d 。
8、绘制线性系统025.0=+x x &&的相平面图。（分别给出x(0)=1,x ’(0)=0; x(0)=2,x ’(0)=0; x(0)=3,x ’(0)=0;）
dxdy
xe xy ??-
三、实验结果
1）2）3）
4）5）6）7）8）
实验三 MATLAB 语言的程序设计
一、实验目的：
1、熟悉MATLAB 程序编辑与设计环境
2、掌握各种编程语句语法规则及程序设计方法
3、函数文件的编写和设计
4、了解和熟悉跨空间变量传递和赋值
二、实验内容：
1、熟悉MATLAB 程序编辑与设计环境。
2、用for 循环语句实现求1～100的和。
3、用for 循环语句实现编写一个求n 阶乘的函数文件。
4、一行100列的Fibonacci 数组a ，元素满足Fibonacci 规则：a(1)=a(2)=1， a(i)=a(i-1)+a(i-2)(i=3,4,…)，求该数组中第一个大于10000的元素。 练习B 思考题
4、已知一维数组A=[2,4,5,8,10],B=[4,6,9,3,4]，用for 循环语句实现∑=+-n
i i n i B A 11。
求和函数可用sum()。.
5、利用for 循环求1！+2！+3！+ … +5！的值。.
6、编写M 函数表示曲线y2=e -t/3sin(3t)以及它的包络线y1=e -t/3，并从命令窗口输入命令语句绘制曲线，t 的取值范围是[0,4π]。
7、设01
.0)3(1
1.0)2(1)(4
2+-++-=
x x x f ，编写一个M 函数文件，使得调用f(x)时，x 可用矩阵代入，得出的f(x) 为同阶矩阵
三、实验结果：
1）2）
3）4）
5）6）
实验五MATLAB图形绘制
一、实验目的：
1、学习MATLAB图形绘制的基本方法；
2、熟悉和了解MATLAB图形绘制程序编辑的基本指令；
3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面，并添加图形
的各种标注；
4、掌握plot、subplot的指令格式和语法。
二、实验内容：
1.二维曲线绘图基本指令演示。本例运作后，再试验plot(t), plot(Y),
plot(Y,t) ，以观察产生图形的不同。
2.用图形表示连续调制波形Y=sin(t)sin(9t)及其包络线。
3.在一个图形窗口绘制正弦和余弦曲线，要求给图形加标题“正弦和余弦
曲线”，X轴Y轴分别标注为“时间t”和“正弦、余弦”，在图形的某个
位置标注“sin(t)”“cos(t)”，并加图例，显示网格，坐标为正方形
坐标系。
三、实验结果：
1）
2）
3）
实验六 MATLAB 数值及符号运算
一、实验目的：
1、掌握矩阵的基本运算；
2、掌握矩阵的数组运算；
3、掌握多项式的基本运算； 5、掌握求解代数方程的方法。
二、实验内容：
1.已知a=[1 2 3;4 5 6;7 8 0]，求其特征多项式并求其根。
2.用MATLAB 求解下面的的方程组。
?????
???????-=?????????????????????????----01741323151122231592127
4321x x x x 3.已知多项式a(x)=x 2+2x+3，b(x)=4x 2+5x+6，求a ，b 的积。
4.求解微分方程0)0(,1)0(,02222===++dt
dy
y y dt dy dt y d ，并画出相平面图。
5.已知系统冲激响应为h(n)=[0.05,0.24,0.40,0.24,0.15,-0.1,0.1] ，系 统 输入u(n)由指令randn('state',1);u=2*(randn(1,100)>0.5)-1产生，该 输入信号的起始作用时刻为0。试用直杆图（提示：用stem 指令）画出分别
显示该系统输入、输出信号的两张子图。
三、实验结果：
1）2）3）
4）
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数学实验练习题
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．掌握二维图形的绘制。 2．掌握图形的标注 3．了解三维曲线和曲面图形的绘制。 二、实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 设计提示 1．Matlab 允许在一个图形中画多条曲线：plot(x1，y1，x2，y2，……) 指令绘制y 1 = f 1(x 1), y 2 = f 2 (x 2 )等多条曲线。Matlab 自动给这些曲线以不同颜色。标注可用text 函数。 2．绘图时可以考虑极坐标和直角坐标的转换。 3．三维曲线绘图函数为plot3，注意参考帮助中的示例。 三、实验内容 1．生成1×10 维的随机数向量a ，分别用红、黄、蓝、绿色绘出其连线图、 杆图、阶梯图和条形图，并分别标出标题“连线图”、“杆图”、“阶梯图”、“条形图”。 2、绘制函数曲线，要求写出程序代码。 (1) 在区间[0:2π]均匀的取50个点，构成向量t (2) 在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求y1曲 线为红色点划线，标记点为圆圈；y2为蓝色虚线，标记点为星号。 (3) 分别在靠近相应的曲线处标注其函数表达式。 3.将图形窗口分成两个绘图区域，分别绘制出函数： ???+-=+=1 352221x x y x y 在[0,3]区间上的曲线，并利用axis 调整轴刻度纵坐标刻度，使1y 在[0,12] 区间上，2y 在[-2,1.5]区间上。 4．用mesh 或surf 函数，绘制下面方程所表示的三维空间曲面，x 和y 的
取值范围设为[-3，3]。 10102 2y x z +-= 思考题： 1. 编写一个mcircle(r)函数，调用该函数时，根据给定的半径r ，以原点 为圆心,画一个如图所示的红色空心圆。（图例半径r=5）；左图参考 polar 函数的用法，右图绘制圆形的参数方程为x=sin （t ），y=cos （t ）。其中，t 的区间为0~2*pi ，步长为0.1。 2．（1）绘一个圆柱螺旋线（形似弹簧）图。圆柱截面直径为10，高度为5， 每圈上升高度为1。如左图所示。 （2）利用（1）的结果，对程序做少许修改，得到如右图所示图形。
Matlab实验三答案
实验三 MATLAB语言的程序设计 一、实验目的： 1、熟悉MATLAB程序编辑与设计环境 2、掌握各种编程语句语法规则及程序设计方法 3、函数文件的编写和设计 4、了解和熟悉跨空间变量传递和赋值 二、实验基本知识： 1、程序流程控制语句 for循环结构 语法：for i＝初值：增量：终值 语句1 …… 语句n end 说明：1．i＝初值：终值，则增量为1。 2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。 while 循环结构 语法：while 逻辑表达式 循环体语句 end 说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达式的值为假，则程序执行end之后的语句。 2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位置加break语句、 以便程序能正常执行。（执行循环体的次数不确定；每一次执行循环体后，一定会改变while 后面所跟关系式的值。） 3、while循环也可以嵌套、其结构如下： while逻辑表达式1 循环体语句1 while逻辑表达式2 循环体语句2 end 循环体语句3 end if-else-end分支结构 if 表达式1 语句1 else if 表达式2(可选) 语句2
else(可选) 语句3 end end 说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。 2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。 3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。 switch-case结构 语法：switch表达式 case常量表达式1 语句组1 case常量表达式2 语句组2 …… otherwise 语句组n end 说明：1．switch后面的表达式可以是任何类型，如数字、字符串等。 2．当表达式的值与case后面常量表达式的值相等时，就执行这个case后面的语句组如果所有的常量表达式的值都与这个表达式的值不相等时，则执行otherwise后的执行语句。程序流程控制指令 break、return、pause 2、函数文件的结构和编写方法 函数只能在函数体内对变量进行操作，也就是只能访问函数本身工作空间中的变量。 M函数文件的结构 （1）函数定义行（function） （2）H1行（函数帮助文本的第一行） （3）函数帮助文本 （4）函数体 （5）注释 函数文件编写后，保存时，其文件名必须与函数名相同。注意：函数名不要与MATLAB 自身的函数命令相同。 三、实验内容： 2、用for循环语句实现求1～100的和 编程： sum=0; for i=1:100 sum=sum+i; end sum
MATLAB实验题答案
result5 = ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么如合法， 结果是多少 >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 58 11 >> result4=b*d result4 = 31 22 22 40 49 13 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 415 30 >> result7=a./b result7 = >> result8= Attempt to reference field of non-structure array. >> result9=a.\b result9 = >> result10=a92 result10 = 1 4 9 16 25 36 >> resultl 1=29a result11 = 2 4 8 16 32 64 >> result5=[b;c']*d 3、用MATLAB求解下面的的方程组。 1、求以下变量的值，并在MATLAB^验证。
1 2 x1 3 2 x2 11 5 x3 2 1 3 x4 >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 7 2 1 2 9 15 3 2 2 2 11 5 1 3 2 13 (1)求矩阵A的秩(rank) (2)求矩阵 A 的行列式(determinant) (3)求矩阵 A 的逆(inverse) (4)求矩阵 A 的特征值及特征向量 (eigenvalue and eigenvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3) >> a=det(A3) >> [V,D]=eig(A3) 10 n 10 查看y 的值) m1=0; for m=-10:10 m仁m1+2^m; end m1 m1 = 6、求分段函数的值。 用if 语句实现，算出下列表中x 对应的y 值。 x=input('enter x='); if x<0 y=x A2+x-6; elseif x>=0&&x<5 y=xA2-5*x+6; else y=xA2-x-1; end y 7、分别用if 和switch 语句实现，将百分 制成绩转换为成绩等级A、B、C、D、E。 其中90~1 00分为A，80~89 分为B，70~79 分为C，60~69 分为D，60 分以下为E。 对超出百分制范围的成绩，给出错误提示 信息。 if 结构程序： x=input('please enter score='); if x>=90&&x<=100 9 2 10 disp('A') 7 2 9 15 (1) 2 2 1 3 4 7 1 0 A 4、已知 2n 2 10 29
matlab数学实验复习题(有标准答案)
复习题 1、写出3 2、i ｎv(A)表示Ａ的逆矩阵； 3、在命令窗口健入 clc,４、在命令窗口健入clea 5、在命令窗口健入６、x=-1:０.2:17、det(A)表示计算Ａ的行列式的值；8、三种插值方法：拉格朗日多项式插值,分段线性插值，三次样条插值。 9、若A=123456789?? ????????,则fliplr （Ａ）＝321654987?????????? Ａ-３=210123456--??????????A ．＾２=149162536496481?????????? tril(A)=100450789?????????? tri ｕ（Ａ，－1)=123456089??????????ｄiag(A ）＝100050009?????????? Ａ(:,2),＝2 58A(3,:)=369 10、nor ｍcd ｆ（1，1，2)=0.5%正态分布mu=1,s ｉgm ａ＝２,x ＝1处的概率 ｅ45(@ｆ,［ａ,b ］,x0),中参数的涵义是@fun 是求解方程的函数M 文 件,[a,b ］是输入向量即自变量的范围a 为初值,ｘ0为函数的初值，t 为输出指定的[ａ,b],x 为函数值 １5、写出下列命令的功能：te ｘt （1,2,‘ｙ=s ｉｎ(x)’
hold on 16fun ｃｔion 开头; 17 ,4) ３，４） 21、设x 是一向量，则）的功能是作出将Ｘ十等分的直方图 ２2、interp １(［1,2,3]，[3,４,5]，2.5) Ans=4.５ 23、建立一阶微分方程组? ??+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来） 二、写出运行结果: １、>＞ｅy ｅ（3，4)=1000 01000010 2、>>s ｉze([1,2,3])=1;3 3、设b=ro ｕnd （unifrnd(-５，5,１，4)),则=3 ５ 2 -5 >>[x,m]=ｍin(b)；x ＝－5；m=4 ，［x,n ］=sort(b ） －５ 2 3 ５ 4 3 １ 2 ｍea ｎ(ｂ）=１.2５，m ｅdian(b)＝2.5,range(ｂ）=10 ４、向量ｂ如上题,则 ＞>ａn ｙ(b),all(ｂ<2),ａｌl(b<6） Ans ＝１ 0 １ ５、＞>[5 6;7 8]>[7 8；5 6]=00 11 6、若1234B ??=???? ,则 7、>>ｄｉag(d ｉag （B ）)=10 04 8、＞>[4:-2:１].＊[-1,６]=-4 12 9、>>acos(0.5),a ｔan(1） ans ＝ 1.６5９8 ans=
matlab实验(西北农林科技大学)
课程实验报告 学年学期2015-2016学年第2学期 课程名称MATLAB与科学计算 实验名称实验一 MATLAB环境与数值计算实验室测量测绘实验中心计算机室专业年级土木13级 学生姓名 学生学号 20130 提交时间 2016-06-02 成绩 任课教师王增红 水利与建筑工程学院
实验一 MATLAB环境与数值计算 1、目的和要求 （1）熟练掌握MATLAB的启动和退出、MATLAB的命令窗口；常用选单和工具栏，通过实例初步认识对MATLAB的命令和程序运行等。 （2）熟练掌握MATLAB变量的使用； （3）熟练掌握矩阵的创建、掌握MATLAB的矩阵和数组的运算； （4）熟悉MATLAB多项式的运用，掌握MATLAB的拟合和插值方法、积分运算和统计运算。2、内容和步骤 参见教材实验一、二。 3、实验报告提交要求 （1）用MALTAB命令和M程序两种方式求解三元一次方程组。 答：以求解三元一次方程组{2a?b+3c=5 3a+b?5c=5 4a?b+c=9 为例： （2）用MALTAB命令和M程序两种方式求解3×3复数矩阵的模和相角。答：
（3）创建一个4╳4矩阵，求逆矩阵、产生对角矩阵、求矩阵的秩等。 （5）创建一个一元四次多项式，进行多项式求值，求根计算。 432
（6）创建一个一元三次多项式，取十个不同自变量值，计算多项式的值，并分别进行一阶、二阶拟合，图示原始数据和拟合结果 （7）自建一个一元三次多项式，取十个不同自变量值，计算多项式的值，并用线性、相邻点和三次插值，图示其插值结果。 解： （8）给出多组温度的测定值（内含平均温度、最低、最高温度），并求其平均温度、最低、最高温度的平均值、标准差等。 解： >> a=fix(40*rand(6,3)+1)
浅析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实验题答案
1、求以下变量的值，并在MATLAB 中验证。 ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法，结果 是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 >> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB 求解下面的的方程组。 (1)????????????-=?????????????????????????----01741323151122231592127 4321x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2)???????=-++=--=-++=++56533332821 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知????????????----=1323151122231592127A (1)求矩阵A 的秩(rank)
matlab实验十--周期函数
实验十周期函数 【实验目的】 1.了解周期函数的基本概念。 2.了解周期函数经过四则运算、复合运算、求导运算、积分运算后的周期性。 3.学习、掌握MATLAB软件有关命令。 【实验内容】 从图形上观察六个三角函数的周期性。 【实验准备】 1.周期函数的基本概念 函数() f x是以T为周期的周期函数是指对任何x，有 += ()() f x T f x 使得上式成立的最小正数T称为函数的最小正周期。 2.周期函数的四则运算 3.周期函数的最小正周期 【实验重点】 1、周期函数的四则运算与复合 2、周期函数的求导与积分运算 【实验难点】 1、最小正周期的确定 【实验方法与步骤】 练习 1 图形上观察六个三角函数
sinx,cosx,tanx,cotx,secx,cscx的周期性。绘制正弦函数y=sinx 在区间[-6π，6π]的图形，相应的MATLAB代码为 >>x=-6*pi:2*pi/30:6*pi; >>y=sin(x); >>plot(x,y); >>xlabel('x');ylabel('y'); 运行结果见图10.1。 从图形中可以看出y=sinx为周期函数，最小正周期T≈6。实际上，最小正周期T=2π=6.28…。同样，可以画出余弦函数y=cosx的图形，见图10.2，其最小正周期也为T=2π。
画正切函数y=tanx 的图形时，要注意函数在,0,1,2,...2x k k π π=+=±±是不连续，所以我们只能分别绘出函数在区间(,),0,1,2,...22k k k ππ ππ-+=±±的图形。相应的MATLAB 代码为 >>x=-1.5:0.01:1.5; >>x1=x-pi;x2=x+pi; >>y=tan(x);y1=tan(x1);y2=tan(x2); >>plot(x,y,x1,y1,x2,y2); >>xlabel('x');ylabel('y'); 运行结果见图10.3。
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章MATLAB基础 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 答：在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 1.5 如何启动M文件编辑/调试器？ 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选
定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB 中有几种获得帮助的途径？ 答：（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。 （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 2.1 在MATLAB 中如何建立矩阵?? ? ? ??194375，并将其赋予变量a ？ 答：在Command Window 窗口输入操作:
Matlab实验第一次实验答案
实验一Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找sqrt（开方）函数的使用方法； 解：sqrt Square root Syntax B = sqrt(X) Description B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results. Remarks See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i
1.0000 1.4142 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B 解：A=[1 2;3 4 ]; B=[5 5;7 8 ]; A^2*B （2）矩阵除法 已知A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B 解：A=[1 2 3;4 5 6;7 8 9 ]; B=[1 0 0;0 2 0;0 0 3 ]; A\B,A/B （3）矩阵的转置及共轭转置
MATLAB实验三-选择结构程序设计方案
MATLAB实验报告 学院：光电学院 班级：073-1 姓名：刘颖 学号：2
实验三 选择结构程序设计 1. 求分段函数的值。 y=?????--≠≠<≤+--≠<-+其他 及且且,13250,653x 0,6222x x x x x x x x x x 用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的值。 程序设计： x=[-5.0,-3.0,1.0,2.0,2.5,3.0,5.0] if x<0&x~=-3 y=x.*x+x-6; elseif x>=0&x<5&x~=2&x~=3 y=x.*x-5*x+6; else y=x.*x-x-1 end y 运行结果： x = -5.0000 -3.0000 1.0000 2.0000 2.5000 3.0000 5.0000 y = 29.0000 11.0000 -1.0000 1.0000 2.7500 5.0000 19.0000 y = 29.0000 11.0000 -1.0000 1.0000 2.7500 5.0000 19.0000 2.输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。其中90分~100分为A ，80分~89分为B ，70分~79分为C ，60分~69分为D ，60分一下为E 。 要求： （1）分别用if 语句和switch 语句实现。 用switch 语句实现 程序设计： score=input('请输入百分制成绩:'); switch fix(score/10) case {9,10} grade='A'; case {8} grade='B'; case {7} grade='C'; case {6} grade='D'; otherwise grade='E'; end grade
MATLAB实验题答案
1、求以下变量的值，并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 >> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A
matlab实验报告
MATLAB 数学实验报告 指导老师： 班级： 小组成员： 时间：201_/_/_
Matlab 第二次实验报告 小组成员： 1 题目：实验四，MATLAB 选择结构与应用实验 目的：掌握if 选择结构与程序流程控制，重点掌握break，return , pause语句的应用。 问题：问题1：验证“哥德巴赫猜想” ，即：任何一个正偶数（n>=6）均可表示为两个质数的和。要求编制一个函数程序，输入一个正偶数，返回两个质数的和。 问题分析：由用户输入一个大于6 的偶数，由input 语句实现。由if 判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数，再向屏幕输出两个质数即可。 编程：function [z1,z2]=gede(n); n=input('please input n')
if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n- i) break end
end end 结果分析 如上图，用户输入了大于6的偶数返回两个质数5和31,通过 不断试验，即可验证哥德巴赫猜想。 纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰，更快的解决问题。 2题目：实验四，MATLAB选择结构与应用实验 目的：用matlab联系生活实际，解决一些生活中常见的实际问 题。
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练习题和答案
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)课后实验答案[1]
实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin851z e =+ (2) 221ln(1)2 z x x =++，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 2242011 122123t t z t t t t t ?≤=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)
4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果： m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 (2). 建立一个字符串向量例如： ch='ABC123d4e56Fg9';则要求结果是： ch='ABC123d4e56Fg9'; k=find(ch>='A'&ch<='Z'); ch(k)=[] ch =