matlab.docx
《matlab.docx》由会员分享,可在线阅读,更多相关《matlab.docx(23页珍藏版)》请在冰豆网上搜索。
matlab
专业基础实践任务书
学生姓名:
专业班级:
指导教师:
工作单位:
信息工程学院
题目:
专业基础实践第3套综合题
初始条件:
(1)提供实验室机房及其Matlab6.5以上版本软件;
(2)《MATLAB教程》学习。
要求完成的主要任务:
(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求):
(1)选择一本《MATLAB教程》,认真学习该教程的全部内容,包括基本使用方法、数组运算、矩阵运算、数学运算、程序设计、符号计算、图形绘制、GUI设计等内容;
(2)对该套综合题的10道题,进行理论分析,针对具体设计部分的原理分析、建模、必要的推导和可行性分析,画出程序设计框图,编写程序代码(含注释),上机调试运行程序,记录实验结果(含计算结果和图表)。
(3)对实验结果进行分析和总结;
(4)要求阅读相关参考文献不少于5篇;
(5)根据课程设计有关规范,按时、独立完成专业基础实践说明书。
时间安排:
(1)布置课程设计任务,查阅资料,学习《MATLAB教程》2.5天;
(2)进行编程设计2天;
(3)完成专业基础实践报告书1天;
(4)答辩0.5天;
指导教师签名:
年月日
系主任(或责任教师)签名:
年月日
题目:
Matlab应用课程设计第3套题
1.算y1=
和y2=
;
2.绘制函数
在
时的曲线。
3.用图形表示离散函数
。
4.分析下面每条指令的功能并运行,观察执行结果。
(1)X=0:
0.5:
6;
Y=X.*exp(-X);
plot(X,Y),xlabel(‘x’),ylabel(‘y’),title(‘y=x*exp(-x)’);
(2)A=zeros(3,5)
A(:
)=-4:
10
L=abs(A)>4
islogical(L)
X=A(L)
(3)A=[1:
5;6:
10]
B=rand(2,5)
A*B’
(4)A=zeros(3,2)
A(:
)=-1:
4
A=A*(1+i)
A1=A.’;
B1=A’;
(5)A=ones(2,3)
B=eye(3)
C=magic(4)
D=diag(C)
E=repmat(C,1,3)
5.计算
在x=0.1与10处的值。
6.求函数
n=100的值。
7.求1500个元素的随机数向量A中大于0.5的元素个数。
8.用图形表示连续调制波形
运行下面的程序,观察结果,并简单说明。
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)')
9.计算下列矩阵A的标准正交基
10.有如下数据:
x
1
1.1
1.2
1.3
1.4
y
1.00000
1.23368
1.55271
1.99372
2.61170
利用三次样条插值方法对其进行插值,得到每隔0.05的结果。
目录
摘要1
1MATLAB概述1
1.1MATLAB简介1
1.2MATLAB的功能1
1.3MATLAB应用2
2.Matlab应用课程设计第3套题3
2.13
2.23
2.34
2.45
2.510
2.611
2.712
2.813
2.914
2.1015
3.体会心得16
4.参考文献17
摘要
MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其它编程语言的程序等。
MATLAB主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
本次课程设计对该套综合题的10道题,进行了理论分析,针对具体设计部分作了原理分析、建模、必要的推导和可行性分析,编写了程序代码(含注释),上机调试了运行程序,并且记录下了实验结果(含计算结果和图表)。
最后对实验结果进行分析和总结。
关键词:
MATLAB;课程设计;理论分析;上机调试;总结
1MATLAB概述
1.1MATLAB简介
MATLAB是集数值计算、符号运算及图形处理等强大功能于一体的科学计算语言,它是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。
MATLAB具有友好的工作平台和编程环境、简易易用的程序语言、强大地科学计算及数据处理能力、出色的图形处理功能、应用广泛的模块集和工具箱、实用的程序接口和发布平台以及模块化的设计和系统级的仿真。
在美国及其他发达国家的理工科院校里,MATLAB已经成为了一门必修的课程,在科研院所、大型公司或企业的工程计算部门,MATLAB也是最为普遍的计算工具之一。
1.2MATLAB的功能
(1)MATLAB以矩阵作为数据操作的基本单位,但无需预先指定矩阵维数
(2)按照IEEE的数值计算标准进行计算
(3)提供十分丰富的数值计算函数,方便计算,提高效率
(4)MATLAB命令与数学中的符号,公式非常接近,可读性强,容易掌握
(5)MATLAB提供丰富的绘图命令,很方便实现数据的可视化
(6)MATLAB具有程序结构控制,函数调用,数据结构,输入输出,面向对象等程序语言特征,而且简单易学,编程效率高
1.3MATLAB应用
MATLAB的应用范围非常广,包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。
附加的工具箱(单独提供的专用MATLAB函数集)扩展了MATLAB环境,以解决这些应用领域内特定类型的问题。
MATLAB产品族可以用来进行以下各种工作:
(1)数值分析
(2)数值和符号计算
(3)工程与科学绘图
(4)控制系统的设计与仿真
(5)数字图像处理技术
(6)数字信号处理技术
2.Matlab应用课程设计第3套题
2.1计算
和
理论分析:
利用MATLAB的基本函数cos()和sqrt()就可以处理这两个算式。
程序代码如下:
>>20*cos(0.5*pi)/(2+sqrt(3))
ans=
3.2814e-016
>>5*cos(0.4*pi)/(4+sqrt(3))
ans=
0.2696
小结:
要注意乘法要用‘*’而不能用‘·’或者不写。
分母要用括号括起来,而‘π’直接写'pi'就可以了。
这与平时的一些书写习惯有所差别。
Matlab在处理类似的数学算式方面很有优势,以后要多加运用。
2.2绘制函数
在
时的曲线。
理论分析:
利用matlab的绘图功能绘制曲线图,我用的是最基本的plot()函数.
程序代码:
x=0:
0.01:
1;%在0~1中以0.01位步进取点
y=x.*exp(-x);%函数
plot(y);%显示图形
title('电信1304班张云波');%标题
xlable('X');ylable('Y');%X,Y轴标题
运行结果如图1所示:
小结:
在二维绘图函数中,最重要、最基本的函数就是plot函数。
其可以在二维坐标系中绘制线性图形。
其调用格式为
Plot(x)
Plot(x,y1,x,y2......),并且可以使用title()函数为图像加上标题,使用xlable()和ylable()函数,为X轴和Y轴设置标题。
2.3用图形表示离散函数
。
理论分析:
在这道题目中,运用到绝对值函数abs(),以及离散点绘图函数stem()
程序代码:
n=-10:
22;%产生一组自变量数据
y=(abs(n-6)).^(-1);%产生相应的函数值
stem(y);%产生离散图
title('电信1304班张云波');%显示标题
运行结果如图2所示:
2.4分析下面每条指令的功能并运行,观察执行结果。
(1)题目
X=0:
0.5:
6;
Y=X.*exp(-X);
plot(X,Y),xlabel(‘x’),ylabel(‘y’),title(‘y=x*exp(-x)’);
理论分析:
X=0:
0.5:
6;表示自变量X的变化范围是0到6,每次自增0.5;Y=X.*exp(-X)表示应变量与自变量之间的函数关系;plot(X,Y)表示用plot函数在二维坐标系里绘制线性函数;xlabel(‘x’),ylabel(‘y’)表示在当前坐标轴上加注x轴,y轴的表示依次为x,y。
而title(‘y=x*exp(-x)’)表示将"y=x*exp(-x)"标注为图形的标题。
或者注释如下:
X=0:
0.5:
6;%长度为13的x轴采样点
Y=X.*exp(-X);%函数表达式
plot(X,Y);%绘制一条曲线,如下所示
xlabel(‘x’);%将水平轴标注为x轴
ylabel(‘y’);%将竖直轴标注为y轴
title(‘y=x.*exp(-x)’);%将函数图形取名为y=x.*exp(-x)
程序代码:
X=0:
0.5:
6;
Y=X.*exp(-X);
plot(X,Y);
xlabel('x');
ylabel('y');
title('y=x.*exp(-x)');
运行结果如图3所示:
(2)题目
A=zeros(3,5)
A(:
)=-4:
10
L=abs(A)>4
islogical(L)
X=A(L)
理论分析:
A=zeros(3,5)%取0数组A,3行5列
A(:
)=-4:
10%将-4至10的15个数值按列依次给A数组赋值
L=abs(A)>4%判断A中元素的绝对值与4的大小,将结果赋值给L
islogical(L)%将L变为L的逻辑数组
X=A(L)%将逻辑值为1的元素列写出来
程序代码:
A=zeros(3,5)
A(:
)=-4:
10
L=abs(A)>4
islogical(L)
X=A(L)
运行结果如下
A=
00000
00000
00000
A=
-4-1258
-30369
-214710
L=
00011
00011
00011
ans=
1
X=
5
6
7
8
9
10
(3)题目:
A=[1:
5;6:
10]
B=rand(2,5)
A*B’
理论分析:
A=[1:
5;6:
10]%输入数组A
B=rand(2,5)%生成一个2行,5列随机数组B
A*B’%矩阵A和B的转置矩阵相乘
程序代码:
A=[1:
5;6:
10]
B=rand(2,5)
A*B'
运行结果如下:
A=
12345
678910
B=
0.19340.30280.15090.37840.8537
0.68220.54170.69790.86000.5936
ans=
7.033310.2671
16.428827.1440
(4)题目
A=zeros(3,2)
A(:
)=-1:
4
A=A*(1+i)
A1=A.’;
B1=A’;
理论分析:
A=zeros(3,2)%取0数组A,3行2列
A(:
)=-1:
4%将-1至4的6个数值按列依次给A数组赋值
A=A*(1+i)%将矩阵A转化为复数矩阵
A1=A.’%将矩阵A转置
B1=A’%将矩阵A共轭转置
程序代码:
A=zeros(3,2)
A(:
)=-1:
4
A=A*(1+i)
A1=A.'
B1=A'
运行结果:
A=
00
00
00
A=
-12
03
14
A=
-1.0000-1.0000i2.0000+2.0000i
03.0000+3.0000i
1.0000+1.0000i4.0000+4.0000i
A1=
-1.0000-1.0000i01.0000+1.0000i
2.0000+2.0000i3.0000+3.0000i4.0000+4.0000i
B1=
-1.0000+1.0000i01.0000-1.0000i
2.0000-2.0000i3.0000-3.0000i4.0000-4.0000i
(5)题目
A=ones(2,3)
B=eye(3)
C=magic(4)
D=diag(C)
E=repmat(C,1,3)
理论分析:
A=ones(2,3)%2行3列全1矩阵A
B=eye(3)%3阶单位方阵B
C=magic(4)%生成一个4阶幻方矩阵,并赋值给C
D=diag(C)%提取对角矩阵C
E=repmat(C,1,3)%生成一个包含3个C的复制的大矩阵E
程序代码:
A=ones(2,3)
B=eye(3)
C=magic(4)
D=diag(C)
E=repmat(C,1,3)
运行结果如下:
A=
111
111
B=
100
010
001
C=
162313
511108
97612
414151
D=
16
11
6
1
E=
162313162313162313
511108511108511108
976129761297612
414151414151414151
小结:
利用矩阵函数可以方便快捷的生成我们所想要的矩阵,并且可以进行矩阵运算,很方便,很实用。
Ones()函数生成全1矩阵,zeros()生成全0矩阵,eye()函数是生成一个单位矩阵,magic()函数生成幻方矩阵,其对角线上元素的和相等,并且幻方矩阵的元素值从1到n^2,repmat()函数按指定维上的数目复制矩阵。
2.5计算
在x=0.1与10处的值。
理论分析1:
可以直接输入方程式,然后就可以计算出结果,如下所示。
程序代码:
x=0.1
y=2*x^5+0.1*x^3-10
x=10
y=2*x^5+0.1*x^3-10
运行结果:
x=
0.1000
y=
-9.9999
x=
10
y=
200090
理论分析2:
这道题目是运用MATLAB计算多项式。
这里运用到的是polyval()函数,通过polyval()函数可以很快的计算出在相应值下的多项式的值。
程序代码:
p=[200.100-10];
polyval(p,0.1)
polyval(p,10)
运行结果如下:
ans=
-9.9988
ans=
20000
总结:
对于多项式的计算,matlab提供了专用的函数polyval(),通过polyval()函数可以很快的计算出在相应值的下的多项式结果。
2.6求函数
n=100的值。
理论分析1:
在这道题目中我用了两种方法,这是第一种即编写程序算出结果。
代码如下:
y=0;
n=100;
fork=1:
n
y=y+1/(k^3);
end
y
结果:
y=1.2020
理论分析2:
在matlab中有专门的计算级数和的函数即symsum()函数。
其格式为:
symsum(x,y,a,b)计算级数x的和,其中y为自变量,x为符号表达式,[a,b]为参数y的取值范围。
代码如下:
symsn;
s=1/n^3;
symsum(s,n,1,100)
运行结果
ans=
81473483330743503583074181861672511931518122336172216406894149391331289704097519580221863303145356050828007873151451209887/67781182783092495848656345091844021571734190630914590229332161379950257170828098031102950264769178652556660142954086400000
而
>>81473483330743503583074181861672511931518122336172216406894149391331289704097519580221863303145356050828007873151451209887/67781182783092495848656345091844021571734190630914590229332161379950257170828098031102950264769178652556660142954086400000
ans=
1.2020
小结:
对于级数求和函数symsum()要注意先定义变量,而计算结果一开始是一大串数数字,我以为错了,最后才发现中间还有一个'/'号。
我把结果运算一下,和上面所得结果一样。
2.7求1500个元素的随机数向量A中大于0.5的元素个数。
理论分析:
对于列写随机数需要用到函数rand().得到矩阵后需要用程序控制,有循环语句,选择语句和嵌套循环语句。
程序代码:
A=rand(50,30);
n=0;
fori=1:
50;
forj=1:
30;
ifA(i,j)>0.5,n=n+1;
end,end,end
n
运行结果如下:
n=
783
小结:
在写代码的过程中需要注意,i循环的次数数组A的行数50,j的循环次数为数组A的列数30。
为而最后三个end与for,for,if三者依次对应,不能缺少,否则得不到最终结果。
2.8用图形表示连续调制波形
运行下面的程序,观察结果,并简单说明。
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)')
解答:
t1=(0:
11)/11*pi;%定义一组自变量(0:
11),步长为pi/11
y1=sin(t1).*sin(9*t1);%数组的乘法
t2=(0:
100)/100*pi;%定义一组自变量(0:
100),步长为pi/100,较y1点密y2=sin(t2).*sin(9*t2);
subplot(2,2,1),%将窗口分为4个,当前为第一个
plot(t1,y1,'r.'),%将点描绘在坐标轴上
axis([0,pi,-1,1]),%控制x,y轴的数值范围,axis[xmin,xmax,ymin,ymax]
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.')%点少时显示的是折线,r表示显示的点是红色
axis([0,pi,-1,1]),title('子图(3)')
subplot(2,2,4),plot(t2,y2)%点多时才显示曲线
axis([0,pi,-1,1]),title('子图(4)')
结果如图4所示:
总结:
plot()函数和subplot()函数是matlab中非常有用的绘图工具,可以很方便很直观的展示复杂函数的图形图像,plot()函数里面可以修改显示图像的颜色,而且,我通过axis()函数来调整X轴和Y轴的输出范围,在以后的工作学习中,要多加运用和学习。
2.9计算下列矩阵A的标准正交基
理论分析:
运用orth()函数处理即可
程序代码:
A=[114;225;336]
[Q,R]=qr(A)
运行结果如下:
A=
114
225
336
Q=%Q即为所求标准正交基
-0.26730.95620.1195
-0.5345-0.0439-0.8440
-0.8018-0.28950.5228
R=
-3.7417-3.7417-8.5524
0-0.00001.8685
00-0.6050
总结:
orth函数的功能是求矩阵正交基。
B=orth(A),返回矩阵A的正交基,B的列与A的列具有相同的空间,B的列向量是正交向量,满足B'*B=eye(rank(A)),B的列数是A的秩。
但是正交基和标准正交基不是同一个概念。
求一个矩阵的标准正交基需要用[Q,R]=qr(A),结果中Q即为所求标准正交基。
2.10有如下数据:
x
1
1.1
1.2
1.3
1.4
y
1.00000
1.23368
1.55271
1.99372
2.61170
利用三次样条插值方法对其进行插值,得到每隔0.05的结果。
理论分析:
利用三次样条插值函数即可,或者利用一维插值函数也可以。
(1)程序代码:
x=