MATLAB课设报告

课程设计任务书

学生姓名：董航专业班级：电信1006班

指导教师：阙大顺，李景松工作单位：信息工程学院

课程设计名称：Matlab应用课程设计

课程设计题目：Matlab运算与应用设计5

初始条件：

1.Matlab6.5以上版本软件；

2.课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应

用”、线性代数及相关书籍等；

3.先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。

要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）

1.课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成；

2.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，

针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括：

①初步了解Matlab、熟悉Matlab界面、进行简单操作；

②MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计；

③基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形

标注、简单颜色设定等；

④使用文本编辑器编辑m文件，函数调用；

⑤能进行简单的信号处理Matlab编程；

⑥按要求参加课程设计实验演示和答辩等。

3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括：

①目录；

②与设计题目相关的理论分析、归纳和总结；

③与设计内容相关的原理分析、建模、推导、可行性分析；

④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结；

⑤课程设计的心得体会（至少500字）；

⑥参考文献（不少于5篇）；

⑦其它必要内容等。

时间安排：1.5周（分散进行）

参考文献：

[1](美)穆尔,高会生,刘童娜,李聪聪．MA TLAB实用教程(第二版) . 电子工业出版社,2010.

[2]王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011.

[3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011.

[4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010.

指导教师签名：年月日

系主任（或责任教师）签名：年月日

目录

1 Matlab的简介................................... 错误!未定义书签。2课程设计内容..................................... 错误!未定义书签。

2.1

2.2

2.34

2.45

2.512

2.613

2.714

2.815

2.916

2.1017

3课程设计心得体会 (20)

4参考文献 (21)

1 Matlab 软件简介

1.1 MATLAB产生的历史背景

MATLAB 是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。

MATLAB是矩阵实验室（Matrix Laboratory）的简称，和Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用 MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB 成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。

MATLAB的应用：利用MATLAB可以进行以下各种工作

●数值计算

●数值和符号计算

●工程和科学绘图

●控制系统的设计与仿真

●数字图像处理

●数字信号处理

●通讯系统设计与仿真

●财务与金融工程

1．2 MATLAB的语言特点和开发环境

MATLAB作为一种科学计算的高级语言之所以受欢迎，就是因为它有丰富的函数资源和工具箱资源，编程人员可以根据自己的需要选择函数，而无需再去编写大量繁琐的程序代码，从而减轻了编程人员的工作负担。被称为第四代编程语言的MATLAB最大的特点就是简洁开放的程序代码和直观实用的开发环境。具体地说MATLAB主要有以下特点：

(1)库函数资源丰富

(2)语言精炼，代码灵活

(3)运算符多而灵活

(4)面向对象控制功能优良，使界面编程更方便、自由。

(5)程序设计自由

(6)图形功能强大

(7)程序的兼容性好

(8)源代码开放

(9)形形色色的工具箱

2．课程设计内容

2.1 已知t=linspace(0,2*pi,6)，求t的正弦和正弦的绝对值。

2.1.1理论分析

该函数linspace（）为赋值向量函数，还要用到简单的正弦sin（）和绝对值abs（）函数。

2.1.2程序代码

t = linspace(0,2*pi,6); %产生6维的从0到2*pi等分行向量

sin(t) %求t的正弦值

abs(sin(t)) %求sin（t）的的绝对值

2.1.3程序运行结果

ans =

0 0.9511 0.5878 -0.5878 -0.9511 -0.0000 ans =

0 0.9511 0.5878 0.5878 0.9511 0.0000 2.1.4总结

在matlab中乘号不能不写，也不能写成‘×’。只能写成‘*’。

2.2 制函数x

=在1

xe

y-

≤x时的曲线。

0≤

2.2.1理论分析：

这是一个简单的函数，但要了解该函数用程序代码书写的正确和规范，以及了解绘制二维曲线的函数。

2.2.2程序代码

x = 0:0.01:1; %x轴上以0.01的速度从0递增到1的数组y = x.*exp(-x); %y关于x的表达式

plot(x,y) %以x为横坐标y为纵坐标画出二维曲线

2.2.3程序运行结果

2.2.4总结：

在matlab中两数之间应该用‘*’符号，不能省略，而且注意如果是进行向量的运算得用‘.*’运算。

2.3 已知a=[1,3,0;2,0,1;4,6,5],a1=logical([1 0

1]),a2=logical([1 1 0]),b=a>1,求a(a1,a2)和a(b)。

2.3.1理论分析：

注意矩阵的赋值以及各种运算

2.3.2程序代码

a = [1,3,0;2,0,1;4,6,5]; %定义a的矩阵并赋值

a1 = logical([1 0 1]); %a1的逻辑值数组

a2 = logical([1 1 0]); %a2的逻辑值数组

b = a > 1; %a,b的逻辑表达式

a(a1,a2) %求矩阵a的1,3行和1,2列的子矩阵a(b) %求a中大于1的数

2.3.3程序运行结果

ans =

1 3

4 6

ans =

2

4

3

6

5

2.3.4总结

在matlab中矩阵的形式各种各样，还有各种赋值形式，我们要了解矩阵的各种形式以及运算

2.4分析下面每条指令的功能并运行，观察执行结果。

（1）X=0:0.25:3;

Y=X.*exp(-X);

plot(X,Y),xlabel(‘x’), ylabel(‘y’),title(‘y=x*exp(-x)’);

(2)A=zeros(3,5)

A(:)=-4:5

L=abs(A)>4

islogical(L)

X=A(L)

(3)A=[1:5;6:10]

pow2(A)

(4)A=zeros(3,2)

A(:)=1:6

A=A*(1+i)

A1=A.’;

B1=A’;

(5)A=ones(2,3)

B=ones(2)

C=eye(4)

D=diag(C)

E=repmat(C,1,3)

2.4.1理论分析：

该大题是让我们练习编写代码，观察执行结果，根据结果去了解各个函数的功能及用法，所以我们必须熟练了解一些常用函数的用法。

2.4.1.2程序代码

X=0:0.25:3; %生成数组

Y=X.*exp(-X); %y1的表达式

plot(X,Y),xlabel('x'), ylabel('y'),title('y=x*exp(-x)'); %标注x,y轴及标题

2.4.1.3程序运行结果

2.4.2.2程序代码

A = zeros(2,5) %生成一个2×5的全0矩阵

A(:) = -4:5 %将A矩阵赋值

L = abs(A) > 4 %将A中绝对值与4比较后的逻辑值矩阵islogical(L) %求L的逻辑值

X = A(L) %求出L中逻辑值为1的那个元素对应A的值

2.4.2.3程序运行结果

A =

-4 -2 0 2 4

-3 -1 1 3 5

L =

0 0 0 0 0

0 0 0 0 1

ans =

1

X =

2.4.

3.2程序代码

A = [1:5;6:10] %生成一个二维数组

pow2(A) %求该数组的每个元素的平方值

2.4.

3.3程序运行结果

A =

1 2 3 4 5

6 7 8 9 10

ans =

2 4 8 16 32

64 128 256 512 1024 2.4.4.2程序代码

A=zeros(3,2) %生成一个3行2列的全0 矩阵

A(:)=1:6 %将A从新赋值

A=A*(1+i) %将A中的每个元素进行运算

A1=A.' %求A的转置矩阵

B1=A' %求A的共轭转置矩阵

2.4.4.3程序运行结果

A =

0 0

0 0

0 0

A =

1 4

2 5

3 6

A =

1.0000 + 1.0000i 4.0000 + 4.0000i

2.0000 + 2.0000i 5.0000 + 5.0000i

3.0000 + 3.0000i 6.0000 + 6.0000i

A1 =

1.0000 + 1.0000i

2.0000 + 2.0000i

3.0000 + 3.0000i

4.0000 + 4.0000i

5.0000 + 5.0000i

6.0000 + 6.0000i

B1 =

1.0000 - 1.0000i

2.0000 - 2.0000i

3.0000 - 3.0000i

4.0000 - 4.0000i

5.0000 - 5.0000i

6.0000 - 6.0000i

2.4.5.2程序代码

A = ones(2,3) %生成一个2×3的全为1的矩阵

B = ones(2) %生成一个2×2的全为1的方阵

C = eye(4) %生成一个4阶单位矩阵

D = diag(C) %取C矩阵中的对角元素作为列向量

E = repmat(C,1,3) %将C矩阵复制为1×3块的矩阵

2.4.5.3程序运行结果

A =

1 1 1

1 1 1

B =

1 1

1 1

C =

1 0 0 0

0 1 0 0

0 0 1 0

0 0 0 1

D =

1

1

1

1

E =

1 0 0 0 1 0 0 0 1 0 0 0

0 1 0 0 0 1 0 0 0 1 0 0

0 0 1 0 0 0 1 0 0 0 1 0

0 0 0 1 0 0 0 1 0 0 0 1 2.4.3总结：

在该程序编程题中，有很多函数，要注意他们的用法以及他们之间的联系，比如A（：）= -4 : 5,是给矩阵数组赋值，在之前

zeros()是确定A 的矩阵形式，所以赋值是应与A 数组矩阵有相同多的元素个数，本来原题中是zeros(3,5),但运行时却有错误，是因为zeros （3,5）赋给A ，A 是3×5矩阵共有15个元素，而后来赋值时又只赋给 -4:5 只有10个元素所以就出错了；当改成zeros （2,5）就可以了。还有各种其他函数的功能和用法。

2.5计算101.0235-+=x x y 在x=0.1与10处的值。

2.5.1理论分析：

该题是简单的二元方程带值求解问题，只需用matlab 编写出该方程的表达式，带值即可求解。

2.5.2程序代码

x = 0.1; %将x 初始化为0.1

y = 2*x^5+0.1*x^3-10; %y 关于x 的表达式

y %求当x=0.1时y 的值 x = 10; %将x 赋值为10

y = 2*x^5+0.1*x^3-10; %y 关于x 的表达式

y %求当x=10时y 的值

2.5.3程序运行结果

y =

-9.9999

y =

200090

2.5.4总结：

在写出表达式之前一定要定义出表达式中的变量

2.6求函数333131211n

++++

Λ, n =100的值

2.6.1理论分析：

该题是一道用循环函数编写的数学题

2.6.2程序代码

s = 0; %为s变量赋初值为0

for i = 1:100 %循环函数，从1到100

s = s+1/i^3; %s的累加求和表达式

end

s %求s的值

2.6.3程序运行结果

s =

1.2020

2.6.4总结：

在这个题目中应该注意s = s+1/i^3的通式的表达，这是一个累加表达式，但是还要注意等比数列的写法。

2.7求1500个元素的随机数向量A中大于0.5的元素个数。

2.7.1理论分析：

该题看起来很简单，其实涉及到了很多方面，产生1500个元素得用的随机函数，依次比较1500个数得用到循环函数，还要看是否大于0.5，所以得用到if函数，再其次就是要用到累加表达式。

2.7.2程序代码

A = rand(1500); %随机函数产生1500个随机数向量赋给一维数组A s = 0; %给s赋初值0

for i = 1:1500 %循环语句，从1到1500循环1500次

if A(i)>0.5 %判断数组A中1500是否大于0.5

s = s+1; %如果大于0.5，则将s的值加1，累计个数

end

end

s %求出s最后的累计值

2.7.3程序运行结果

s =

760

2.7.4总结

要注意随机函数的用法，以及要分析好函数的嵌套，那个rand（）是随机函数所以每次产生的数都是不一样的，即每次运行的结果应该是不同的。应该巧妙运用到for函数与if函数的嵌套。

2.8用图形表示连续调制波形)9

t

y ,运行下面的程序，观察结

sin(t

)

sin(

果，并简单说明。

t1=(0:11)/11*pi;

y1=sin(t1).*sin(9*t1);

t2=(0:100)/100*pi;

y2=sin(t2).*sin(9*t2);

subplot(2,2,1),plot(t1,y1,'r.'),axis([0,pi,-1,1 ]),title('子图 (1)')

subplot(2,2,2),plot(t2,y2,'r.'),axis([0,pi,-1,1 ]),title('子图 (2)')

subplot(2,2,3),plot(t1,y1,t1,y1,'r.')

axis([0,pi,-1,1]),title('子图 (3)')

subplot(2,2,4),plot(t2,y2)

axis([0,pi,-1,1]),title('子图 (4)')

2.8.1理论分析：

这题是用图形表示连续调制函数y=sin（t）sin（9t）的波形，注意比较几个图得到不同。

2.8.2程序代码

t1 = (0:11)/11*pi; %生成一个0到1/π步距为1/（11*π）的数组

y1 = sin(t1).*sin(9*t1); %y1关于t1的数学表达式

t2 = (0:100)/100*pi; %生成一个0到1/π步距为1/（100*π）的数组

y2 = sin(t2).*sin(9*t2); %y2关于t2的数学表达式

subplot(2,2,1),plot(t1,y1,'r.'),axis([0,pi,-1,1]),title('子图（1）') %指定分成一个图到指定位置，以r的格式画出t1，y1点图形，x轴从0到π，y轴-1到1，标记为子图（1）

subplot(2,2,2),plot(t2,y2,'r.'),axis([0,pi,-1,1]),title('子图（2）') %指

定分成一个图到指定位置，以r的格式画出t2，y2点图形，x轴从0到π，y轴-1到1，标记为子图（2）

subplot(2,2,3),plot(t1,y1,t1,y1,'r.')

axis([0,pi,-1,1]),title('子图（3）') %指定分成一个图到指定位置，以r的格式画出t1，y1连线的折线图形，x轴从0到π，y轴-1到1，标记为子图（3）

subplot(2,2,4),plot(t2,y2)

axis([0,pi,-1,1]),title('子图（4）') %指定分成一个图到指定位置，以r的格式画出t2，y2默认的曲线图形，x轴从0到π，y轴-1到1，标记为子图（4）

2.8.4总结：

在这几个图中，子图（1）描写函数的能力最差，子图（3）也有较大的失真，因为数据太少。对于离散数据，Plot指令默认处理方法是：自动地把这些离散点用直线连接，使之成为连续曲线。Subplot是指定几号子图。Axis设置平面轴的范围、坐标分割线。

2.9已知方程组

123

123

123

3x x2x9

-5x x+3x5

x x4x8

-+=

?

?

+=

?

?-+=

?

，用矩阵除法来解线性方程组。

2.9.1理论分析：

本题是运用矩阵法来求解线性方程组，在matlab中有很简单实用的方法求解线性方程组

2.9.2程序代码

A = sym([3,-1,2;-5,1,3;1,-1,4]); %将方程组的各个元素的系数定义为符号数组

b = sym([9;5;8]); %将方程组等号后的值定义为符号列向量

X = A\b %解出方程组

2.9.4总结：

将各个方程组的各个元素的的系数赋给A，用到方程组的符号解，矩阵计算是求解线性方程组最简便有效的方法。用sym()符号函数，注意b要用sym（9；5；8），不能用sym（9,5,8）。

2.10已知-1

，其中t的范围是[0 10]，计算y的微分和积分并

y e cos(t)

给出曲线图。

2.10.1理论分析：

这是一个求积分和微分的题，以及也涉及到画图函数。有各种不同的画图函数，注意选择方便的。

2.10.2程序代码以及结果：

syms t,y = exp(-1)*cos(t); %定义符号变量

df = diff(y,'t') %将符号t视作变量求y对t的微分jf = int(y,t) %将符号t作为变量求y对t的积分subplot(2,1,1),ezplot(df,[0,10]) %在分成后的指定地方绘制微分在0到10之间的曲线图

subplot(2,1,2),ezplot(jf,[0,10]) %在分成后的指定地方绘制积分在0到10之间的曲线图

2.10.3程序运行结果

2.10.4总结

要了解积分和微分函数的运用，该题用到了符号积分和符号微分，对ezplot（）画图函数的运用；

3心得体会

通过这次MATLAB的学习，我对MATLAB有了一个基础的认识，matlab是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。它集图示和精确计算于一身，在应用数学、物理、化工、机电工程、医药、金融和其他需要进行复杂数值计算的领域得到了广泛应用。MATLAB是一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程的特点。用户可以在命令窗口中将输入语句与执行命令同步，也可以先编写好一个较大的复杂程序（M文件）后再一起运行。

在这短短的一周内从开始的一头雾水，到自己看书学习，到同学讨论，再进行整个题目的理论分析和计算，参考课程上的代码，写出自己的代码。

我们也明白了学无止尽的道理，在我们所查的很多参考书中，很多知识是我们从没有接触过的，我们对它的了解还仅限于皮毛，对它的很多功能以及函数还不是很了解，所以在这个学习的过程中我们穿越在知识的海洋中，一点一点吸取着它的知识。在MATLAB编程中需要很多的参考书，要尽量多的熟悉matlab自带的函数及其作用，因为matlab的自带函数特别多，基本上能够满足一般的数据和矩阵的计算，所以基本上不用你自己编函数。这一点对程序非常有帮助，可以使程序简单，运行效率高，可以节省很多时间。本次课设中用了很多MATLAB自带的函数，使程序变得很简单。

把基本的知识看过之后，就需要找一个实际的程序来动手编一下，不要等所有的知识都学好之后再去编程，你要在编程的过程中学习，程序需要什么知识再去补充，编程是一点一点积累的，所以你要需做一些随手笔记什么的。

在编写程序代码时，需要什么函数，需要什么模块就应该去着重看那个知识点，不要一步登天，一步一步学，如果太急于把所有东西都学到，也是不好的，更是实现不了的。所以那时一天一天积累的，慢慢地学通这个软件。

总之，通过这次学习，我了解了一下这个软件总体的功能，以及通过自己编写一些代码也学到了一些用法和知识。更了解到了，我们还有好多东西去学，学无止尽。

4参考文献

[1] （美）穆尔,高会生,刘童娜,李聪聪. MATLAB实用教程(第二版)[M] . 北京:

电子工业出版社,2010.

[2] 王正林,刘明. 精通MATLAB(升级版) [M]. 北京:电子工业出版社,2011.

[3] 陈杰. MATLAB宝典(第3版) [M]. 北京:电子工业出版社,2011.

[4] 刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版)[M] . 人民邮电

出版社,2010.

[5] 张德丰. MATLAB在电子信息工程中的应用[M].北京:电子工业出版社,2009.

本科生课程设计成绩评定表

指导教师签字：

年月日

MATLAB实验报告50059

实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1

例1-2 例1-3 例1-4

3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。

5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？

matlab实验报告

数学实验报告 班级: 学号： 姓名： 实验序号：1 日期：年 月 日 实验名称：特殊函数与图形 ◆ 问题背景描述：绘图是数学中的一种重要手段，借助图形，可以使抽象的对象得到 明白直观的体现，如函数的性质等。同时，借助直观的图形，使初学者更容易接受新知识，激发学习兴趣。 ◆ 实验目的：本实验通过绘制一些特殊函数的图形，一方面展示这些函数的特点属性， 另一方面，就 Matlab 强大的作图功能作一个简单介绍。 实验原理与数学模型： 1、 球2222x y z R ++= ，x=Rsin φcos θ, y= Rsin φsin θ, z= cos φ, 0≤θ≤2π , 0≤φ≤π 环面 222222222()4(),(cos )cos ,x y z a r a x y x a r φθ+++-=+=- (cos )sin ,sin ,02,02y a r z r φθφφπθπ=-=≤≤≤≤ 2、 平面摆线：2 22 31150,(sin ),(1cos ),0233 x y x a t t y a t t π+-==-=-≤≤ 3、 空间螺线：（圆柱螺线）x=acost , y=asint , z=bt ;（圆锥螺线）22 cos ,sin ,x t t y t t z t === 4、 椭球面sin cos ,sin sin ,cos ,02,0x a y b z c φθφθφθπφπ===≤<≤≤ 双叶双曲面3 tan cos ,tan sin ,sec ,02,22 x a y b z c π φθφθφθπφπ===≤<- << 双曲抛物面2 sec ,tan 2 u x au y bu z θθ=== 实验所用软件及版本：mathematica(3.0) 主要内容（要点）： 1、 作出下列三维图形（球、环面） 2、 作出下列的墨西哥帽子 3、 作出球面、椭球面、双叶双曲面，单叶双曲面的图形 4、 试画出田螺上的一根螺线 5、 作出如图的马鞍面

Matlab通信系统仿真实验报告

Matlab通信原理仿真 学号： 2142402 姓名：圣斌

实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);

Matlab数学实验报告一

数学软件课程设计 题目非线性方程求解 班级数学081 姓名曹曼伦

实验目的：用二分法与Newton迭代法求解非线性方程的根； 用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 编程实现二分法及Newton迭代法； 学会使用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。 实验内容： 比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方） (1)在区间（0，1)内用二分法； (2)用迭代法x=(2-exp(x))/10，取初值x=0 。 试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0，则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正，则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负，则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法； format long

syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果： （1）二分法： symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 （2）迭代法 由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次，而用迭代法求得同样精度的解仅用5次，但由于迭代法一般只具有局部收敛性，因此通常不用二分法来求得非线性方程的精确解，而只用它求得根的一个近似解，再用收敛速度较快的迭代法求得其精确解。

MATLAB仿真实验报告

MATLAB 仿真实验报告 课题名称：MATLAB 仿真——图像处理 学院：机电与信息工程学院 专业：电子信息科学与技术 年级班级：2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作，熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能，理解多维图像的相关性质 3、熟悉Help 命令的使用，掌握对相关函数的查找，了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用，能够处理多维图像。 二、实验条件

MATLAB调试环境以及相关图像处理的基本MATLAB语句，会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents ： ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating（联系起来）A1, A2 ... along the dimension DIM. Calls to CAT can be nested（嵌套）. (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts（整数下标）. (4)、Manipulating multi-dimensional arrays

数学实验“几种常见的求积分近似解的方法”实验报告(内含matlab程序)

西京学院数学软件实验任务书

实验二十一实验报告 一、实验名称：Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 二、实验目的：进一步熟悉Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 三、实验要求：运用Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。 四、实验原理： 1．Romberg 积分法： 龙贝格积分法是用里查森外推算法来加快复合梯形求积公式的收敛速度，它的算法如下，其中()i m T 是通过一系列逼近原定积分的龙贝格分值. 计算(0)1[()()]2 b a T f a f b -= + 对1,2,3,k n = ，计算下列各步： 21()(1)1 111 1(21)()[()]222k k k k k j b a j b a T T f a ---=---=++∑

对1,2,,m k = 和,1,2,,1i k k k =-- ，计算111 441 m i i i m m m m T T T --+-=- 随着计算的步骤的增加，()i m T 越来越逼近积分()b a f x dx ?。 2．Gauss 型积分法： 高斯积分公式的思想是用n 个不等距的节点123,,,n x x x x 对被积函数进行插值，然后对插值后的函数进行积分，其积分公式为： 1 1 1 ()()n k k k f x dx A f x -=≈∑? 如果积分区间不是[1,1]-，则需转换到此区间： 11()()222 b a b a b a b a f x dx f t dt ---+= +? ? 其中系数k A 、节点k x 与n 的关系如下表所示： 3．高斯-切比雪夫积分法： 第一类切比雪夫积分形式为： 1 1 ()()n k k k f x dx A f x -=≈∑? 其中k A n π= ，21cos 2k k x n π-= 4．高斯-拉盖尔积分法： 高斯-拉盖尔公式有两种形式： 1 ()()n x k k k e f x dx A f x +∞ -=≈∑?

自动控制原理MATLAB仿真实验报告

实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些？ 2、 如何判断系统稳定性？ 3、 系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间范围0---Tn 。 3、),(T sys step ；表示时间范围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：) ()()()(1)(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.

MATLAB Simulink系统建模与仿真 实验报告

MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************

实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：

1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：

matlab数学实验报告5

数学实验报告 制作成员班级学号 2011年6月12日

培养容器温度变化率模型 一、实验目的 利用matlab软件估测培养容器温度变化率 二、实验问题 现在大棚技术越来越好，能够将温度控制在一定温度范围内。为利用这种优势，实验室现在需要培植某种适于在8.16℃到10.74℃下能够快速长大的甜菜品种。为达到实验所需温度，又尽可能地节约成本，研究所决定使用如下方式控制培养容器的温度：1，每天加热一次或两次，每次约两小时； 2，当温度降至8.16℃时，加热装置开始工作；当温度达到10.74℃时，加热装置停止工作。 已知实验的时间是冬天，实验室为了其它实验的需要已经将实验室的温度大致稳定在0℃。下表记录的是该培养容器某一天的温度 时间(h)温度（℃）时间（h）温度（℃）09.68 1.849.31 0.929.45 2.959.13 3.878.981 4.989.65 4.988.811 5.909.41 5.908.691 6.839.18 7.008.5217.938.92 7.938.3919.048.66 8.978.2219.968.43 9.89加热装置工作20.848.22 10.93加热装置工作22.02加热装置工作10.9510.8222.96加热装置工作12.0310.5023.8810.59 12.9510.2124.9910.35 13.889.9425.9110.18 三、建立数学模型 1，分析：由物理学中的傅利叶传热定律知温度变化率只取决于温度

差，与温度本身无关。因为培养容器最低温度和最高温度分别是：8.16℃和10.74℃；即最低温度差和最高温度差分别是：8.16℃和10.74℃。而且，16.8/74.10≈1.1467，约为1，故可以忽略温度对温度变化率的影响2， 将温度变化率看成是时间的连续函数，为计算简单，不妨将温度变化率定义成单位时间温度变化的多少，即温度对时间连续变化的绝对值（温度是下降的），得到结果后再乘以一系数即可。 四、问题求解和程序设计流程1）温度变化率的估计方法 根据上表的数据，利用matlab 做出温度-时间散点图如下： 下面计算温度变化率与时间的关系。由图选择将数据分三段，然后对每一段数据做如下处理：设某段数据为{(0x ,0y ),(1x ,1y ),(2x , 2y ),…,(n x ,n y )},相邻数据中点的平均温度变化率采取公式： 温度变化率=（左端点的温度-右端点的温度）/区间长度算得即：v( 2 1i i x x ++)=(1+-i i y y )/(i i x x - +1). 每段首尾点的温度变化率采用下面的公式计算：v(0x )=(30y -41y +2y )/(2x -0x )v(n x )=(3n y -41+n y +2+n y )/(n x -2-n x )

matlab实验报告

实验一小球做自由落体运动内容：一小球竖直方向做自由落体，并无损做往返运动。程序： theta=0:0.01:2*pi x=cos(theta) y=sin(theta) l=1 v=1 while l<10 for t=1:10 y=y+(-1)^l*v*t plot(x,y,[-1,1],[-56,2],'.') axis equal pause(0.1) end l=l+1 end 结果：

-50 -40 -30 -20 -10 收获：通过运用小球自由落体规律，及（-1）^n 来实现无损往 返运动！ 实验二 旋转五角星 内容：一个五角星在圆内匀速旋转 程序：x=[2 2 2 2 2 2] y=[0 4/5*pi 8/5*pi 2/5*pi 6/5*pi 0] y1=2*sin(y) x1=2*cos(y) theta=0:4/5*pi:4*pi

x2=2*cos(theta) y2=2*sin(theta) plot(x,y,x1,y1,x2,y2) axis equal theta1=theta+pi/10 x2=2*cos(theta1) y2=2*sin(theta1) plot(x2,y2) axis equal theta=0:4/5*pi:4*pi for rot=pi/10:pi/10:2*pi x=2*cos(theta+rot) y=2*sin(theta+rot) plot(x,y) pause(0.1) end 结果：

-2 -1.5-1-0.500.51 1.52 -2-1.5-1-0.500.511.5 2 收获：通过theta1=theta+pi/10，我们可以实现五角星在圆内匀速 旋转！ 实验三 转动的自行车 内容：一辆自行车在圆内匀速转动 程序：x=-4:0.08:4; y=sqrt(16-x.^2); theta1=-pi/2:0.01*pi:3*pi/2; x3=0.5*cos(theta1); y3=0.5*sin(theta1); theta=-pi/2+0.02*pi for k=1:100

MATLAB实验报告

数字信号处理及MATLAB 实验报告 班级： 学号： 姓名：

4.7.2 例4,2 设x(n)是由两个正弦信号及白噪声的叠加，试用FFT文件对其作频谱分析。程序清单 %产生两个正弦加白噪声 N=256; f1=.1;f2=.2;fs=1; a1=5;a2=3; w=2*pi/fs; x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N); %应用FFT求频谱 subplot(2,2,1); plot(x(1:N/4)); title('原始信号'); f=-0.5:1/N:0.5-1/N; x=fft(x); y=ifft(x); subplot(2,2,2); plot(f,fftshift(abs(x))); title('频域信号'); subplot(2,2,3); plot(real(x(1:N/4))); title('时域信号');

例4.3 设x（n）为长度N=6的矩形序列，用MATLAB程序分析FFT取不同长度时x（n）频谱的变化。N=8,32,64,时x（n）的FFT MATLAB实现程序如下。 x=[1,1,1,1,1,1]; N=8; y1=fft(x,N); n=0:N-1; subplot(3,1,1);stem(n,abs(y1),'.k');axis([0,9,0,6]); N=32; y2=fft(x,N); n=0:N-1; subplot(3,1,2);stem(n,abs(y2),'.k');axis([0,40,0,6]); N=64; y3=fft(x,N); subplot(3,1,3);stem(n,abs(y3),'.k');axis([0,80,0,6]);

参考答案Matlab实验报告

实验一 Matlab基础知识 一、实验目的： 1.熟悉启动和退出Matlab的方法。 2.熟悉Matlab命令窗口的组成。 3.掌握建立矩阵的方法。 4.掌握Matlab各种表达式的书写规则以及常用函数的使 用。 二、实验内容： 1.求[100，999]之间能被21整除的数的个数。（rem） 2.建立一个字符串向量，删除其中的大写字母。（find） 3.输入矩阵，并找出其中大于或等于5的元素。（find） 4.不采用循环的形式求出和式 63 1 2i i= ∑ 的数值解。(sum) 三、实验步骤： ●求[100，199]之间能被21整除的数的个数。（rem） 1.开始→程序→Matlab 2.输入命令： ?m=100:999; ?p=rem(m,21); ?q=sum(p==0) ans=43 ●建立一个字符串向量，删除其中的大写字母。（find） 1.输入命令：

?k=input('’,’s’)； Eie48458DHUEI4778 ?f=find(k>=’A’&k<=’Z’)； f=9 10 11 12 13 ?k(f)=[ ] K=eie484584778 ●输入矩阵，并找出其中大于或等于5的元素。（find） 1.输入命令： ?h=[4 8 10;3 6 9; 5 7 3]； ?[i,j]=find(h>=5) i=3 j=1 1 2 2 2 3 2 1 3 2 3 ●不采用循环的形式求出和式的数值解。（sum） 1.输入命令： ?w=1:63； ?q=sum(2.^w) q=1.8447e+019

实验二 Matlab 基本程序 一、 实验目的： 1. 熟悉Matlab 的环境与工作空间。 2. 熟悉M 文件与M 函数的编写与应用。 3. 熟悉Matlab 的控制语句。 4. 掌握if,switch,for 等语句的使用。 二、 实验内容： 1. 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 2. 编程完成，对输入的函数的百分制成绩进行等绩转换，90~100为优，80~89为良，70~79为中，60~69为及格。 3. 编写M 函数文件表示函数 ，并分别求x=12和56时的函数值。 4. 编程求分段函数 2226;03 56;0532 1;x x x x y x x x x x x x +-<≠=-+≤<≠≠-+且且及其它，并求输入x=[-5.0,-3.0,1.0,2.0,2.5,3.0,3.5]时的输出y 。 三、 实验步骤： 根据y=1+1/3+1/5+……+1/（2n-1）,编程求：y<5时最大n 值以及对应的y 值。 1. 打开Matlab ，新建M 文件 2. 输入命令： 51022-+x

控制理论实验报告MATLAB仿真实验解析

实验报告 课程名称：控制理论（乙） 指导老师：林峰 成绩：__________________ 实验名称：MATLAB 仿真实验 实验类型：________________同组学生姓名：__________ 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 实验九 控制系统的时域分析 一、 实验目的： 1．用计算机辅助分析的办法，掌握系统的时域分析方法。 2．熟悉Simulink 仿真环境。 二、实验原理及方法： 系统仿真实质上就是对系统模型的求解，对控制系统来说，一般模型可转化成某个微分方程或差分方程表示，因此在仿真过程中，一般以某种数值算法从初态出发，逐步计算系统的响应，最后绘制出系统的响应曲线，进而可分析系统的性能。控制系统最常用的时域分析方法是，当输入信号为单位阶跃和单位冲激函数时，求出系统的输出响应，分别称为单位阶跃响应和单位冲激响应。在MATLAB 中，提供了求取连续系统的单位阶跃响应函数step ，单位冲激响应函数impulse ，零输入响应函数initial 等等。 二、实验内容： 二阶系统，其状态方程模型为 ? 1x -0.5572 -0.7814 1x 1 = + u ? 2x 0.7814 0 2x 0 1x y = [1.9691 6.4493] +[0] u 2x 四、实验要求： 1．编制MATLAB 程序，画出单位阶跃响应曲线、冲击响应曲线、系统的零输入响应、斜坡输入响应； （1）画出系统的单位阶跃响应曲线； A=[-0.5572 -0.7814;0.7814 0 ]; B=[1;0];

本科毕业设计__基于matlab的通信系统仿真报告

创新实践报告
报 告 题 目： 学 院 名 称： 姓 名：
基于 matlab 的通信系统仿真 信息工程学院 余盛泽 11042232 温 靖
班 级 学 号： 指 导 老 师：
二 O 一四年十月十五日

目录
一、引言 ....................................................................................................................... 3 二、仿真分析与测试 ................................................................................................... 4
2.1 随机信号的生成................................................................................................................ 4 2.2 信道编译码......................................................................................................................... 4 2.2.1 卷积码的原理 ......................................................................................................... 4 2.2.2 译码原理................................................................................................................. 5 2.3 调制与解调........................................................................................................................ 5 2.3.1 BPSK 的调制原理 ................................................................................................... 5 2.3.2 BPSK 解调原理 ....................................................................................................... 6 2.3.3 QPSK 调制与解调................................................................................................... 7 2.4 信道..................................................................................................................................... 8 2.4.1 加性高斯白噪声信道 ............................................................................................. 8 2.4.2 瑞利信道................................................................................................................. 8 2.5 多径合并............................................................................................................................. 8 2.5.1 MRC 方式 ................................................................................................................ 8 2.5.2 EGC 方式................................................................................................................. 9 2.6 采样判决............................................................................................................................. 9 2.7 理论值与仿真结果的对比 ................................................................................................. 9
三、系统仿真分析 ..................................................................................................... 11
3.1 有信道编码和无信道编码的的性能比较 ....................................................................... 11 3.1.1 信道编码的仿真 .................................................................................................... 11 3.1.2 有信道编码和无信道编码的比较 ........................................................................ 12 3.2 BPSK 与 QPSK 调制方式对通信系统性能的比较 ........................................................ 13 3.2.1 调制过程的仿真 .................................................................................................... 13 3.2.2 不同调制方式的误码率分析 ................................................................................ 14 3.3 高斯信道和瑞利衰落信道下的比较 ............................................................................... 15 3.3.1 信道加噪仿真 ........................................................................................................ 15 3.3.2 不同信道下的误码分析 ........................................................................................ 15 3.4 不同合并方式下的对比 ................................................................................................... 16 3.4.1 MRC 不同信噪比下的误码分析 .......................................................................... 16 3.4.2 EGC 不同信噪比下的误码分析 ........................................................................... 16 3.4.3 MRC、EGC 分别在 2 根、4 根天线下的对比 ................................................... 17 3.5 理论数据与仿真数据的区别 ........................................................................................... 17
四、设计小结 ............................................................................................................. 19 参考文献 ..................................................................................................................... 20

浅析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);