科研实践报告格式l副本.docx

科研实践报告格式l副本.docx

《科研实践报告格式l副本.docx》由会员分享，可在线阅读，更多相关《科研实践报告格式l副本.docx（23页珍藏版）》请在冰豆网上搜索。

科研实践报告格式l副本

论文题目:

单边带幅度调制系统仿真

学院：

计算机与信息学院

专业年级：

电子信息工程2008级

学号：

081151007

姓名：

钟季水

任课教师、职称：

谢秀娟讲师

2010年12月30日

目录

1、引言………………………………………………………………………………3

2、matlab基础知识…………………………………………………………………3

3、matlab数值分析…………………………………………………………………8

4、matlab的符号运算………………………………………………………………13

5、matlab的绘图功能………………………………………………………………17

6、simulink基础……………………………………………………………………21

7、单边带幅度调制系统仿真………………………………………………………23

8、结束语……………………………………………………………………………24

题目（黑体小二居中。

中文标题必须严格限制在20字内，如超出可采用副标题形式）

电子信息工程专业08年级钟季水

1引言

2一级标题（宋体小四靠左）

2.1MATLAB基础知识

2．实现下列运算，并观察它们的结果，理解运算规则。

（1）矩阵加减：

已知矩阵A=[111;222;333];B=A;求A+B=？

、A-B=？

、A+2=？

、A-1=？

>>A=[111;222;333;];

>>B=A;

>>A+B

ans=

222

444

666

>>A-B

ans=

000

000

000

>>A+2

ans=

333

444

555

>>A-1

ans=

000

111

222

（2）矩阵乘：

已知矩阵C=ones（1,3）;D=ones（3,1）;求C*D=?

、D*C=?

>>C=ones（1,3）;D=ones（3,1）;

>>C*D

ans=

3

>>D*C

ans=

111

111

111

（3）矩阵乘和数组乘：

已知矩阵E=[11111;22222;33333];F=E;求矩阵乘E*F=？

数组乘E.*F=?

、F.*E=?

矩阵乘和数组乘有何区别？

>>E=[11111;22222;33333];

>>F=E;

>>E.*F

ans=

11111

44444

99999

>>F.*E

ans=

11111

44444

99999

（4）矩阵除：

已知矩阵A=[1111;2222];B=[1111];求A/B=？

、A/2=？

数组除：

已知矩阵A=[223344;112233;445566];B=[113322;111111;225533];求A./B=？

、B./A=？

>>A=[1111;2222];

>>B=[1111];

>>A/B

ans=

1.0000

2.0000

>>A/2

ans=

0.50000.50000.50000.5000

1.00001.00001.00001.0000

>>A=[223344;112233;445566];

>>B=[113322;111111;225533];

>>A./B

ans=

221122

112233

221122

>>B./A

ans=

0.50000.50001.00001.00000.50000.5000

1.00001.00000.50000.50000.33330.3333

0.50000.50001.00001.00000.50000.5000

3、已知A=[223344;112233;445566];

（1）用size（）函数计算矩阵A的维数，调用格式：

d=size（A）%将矩阵A的行数和列数赋给变量d

[m,n]=size（A）%将矩阵A的行数赋给变量m、列数赋给变量n

>>A=[223344;112233;445566];

>>d=size（A）

d=

36

>>[m,n]=size（A）

m=

3

n=

6

（2）用rand（）函数产生随机矩阵，，调用格式：

rand（n）%产生值在0~1之间随机分布的n*n的随机方阵

rand（m,n）%%生值在0~1之间随机分布的n*m的随机矩阵

>>rand（n）

ans=

0.95010.45650.92180.41030.13890.0153

0.23110.01850.73820.89360.20280.7468

0.60680.82140.17630.05790.19870.4451

0.48600.44470.40570.35290.60380.9318

0.89130.61540.93550.81320.27220.4660

0.76210.79190.91690.00990.19880.4186

>>rand（m,n）

ans=

0.84620.67210.68130.50280.30460.6822

0.52520.83810.37950.70950.18970.3028

0.20260.01960.83180.42890.19340.5417

（3）计算矩阵长度（列数）的函数length（），调用格式：

a=length（B）%将矩阵B的列数赋值给变量a

>>a=length（B）

a=

6

（4）矩阵元素的求积运算函数prod（），调用格式：

prod（A）%若A为向量，将计算矩阵A所有元素之积；若A为矩阵，将产生一行向量，其元素分别为矩阵A的各列元素之积。

prod（A,k）%将对矩阵A按k定义的方向进行示积运算，若k=1则按列的方向求积，若k=2则按行的方向求积。

>>prod（A）

ans=

8830307272

>>prod（A,2）

ans=

576

36

14400

（5）矩阵元素的求和运算函数sum（），调用格式同prod（）函数。

>>sum（A）

ans=

7710101313

>>sum（A,2）

ans=

18

12

30

三级标题（宋体小四靠左）

……

2.2matlab的数值分析

1、实现下列连续时间信号

（1）

（画t=-1~3）

程序如下:

clearall;

t=-1:

0.00001:

3;

t0=0;

u=stepfun（t,t0）;

plot（t,u）;

axis（[-21-0.21.2]）

title（'\itx（t）\rm=u（t）'）

title（'\itx（t）\rm=u（t）','FontSize',12）

（2）

程序如下

clearall;

t=-1:

0.0001:

3;

t1=0;

u1=stepfun（t,t1）;

t2=2;

u2=stepfun（t,t2）;

u=u1-u2;

plot（t,u）;

axis（[-13-11.1]）;

title（'\itx（t）\rm=u（t）-u（t-2）'）

title（'\itx（t）\rm=u（t）-u（t-2）','FontSize',12）

（3）

（画t=-1~3）

程序如下：

clearall

t=-1:

0.00001:

3;

t1=0;

u1=stepfun（t,t1）;

u=t.*u1;

plot（t,u）;

axis（[-13-15]）;

title（'\itx（t）\rm=tu（t）'）

title（'\itx（t）\rm=tu（t）','FontSize',12）

（4）

（画t=-1~3）

clearall;

t=-1:

0.0001:

3;

b1=10;

a1=-1;

b2=-5;

a2=-2;

u1=a1*exp（a1*t）;

u2=a2*exp（a2*t）;

u=u1+u2;

plot（t,u）;

axis（[-1.23.1-182]）

title（'\ity\rm=10e^{-\itt}-5e^{-\it2t}'）

title（'\ity\rm=10e^{-\itt}-5e^{-\it2t}','FontSize',12）

提示：

画阶跃信号可先创建如下函数，并保存为Heaviside.m，再调用即可。

functionf=Heaviside（t）

f=（t>0）;

2、已知x（t）的波形如下图，

（1）编写实现该信号波形的函数x（t）的MATLAB函数；

functionyt=x（t）

yt=1*t.*（t>=0&t<=2）+2*（t>2&t<=3）+（-1）*（t>3&t<=5）;

（2）画出x（t），x（0.5t）和x（2-0.5t）。

2.3matlab的符号运算

1、用符号表达式表示信号

，并绘制波形。

符号表达式表示如下：

functionyt=x（t）

yt='sin（pi/4*t）';

绘制波形程序如下：

closeall;

t=linspace（0,2*pi,100）;%100个点的x坐标

y=sin（pi/4*t）;%对应的y坐标

plot（x,y）;

2、实现下列连续时间信号x（t）的微分和积分。

程序如下：

functionyt=x（t）

yt=1*t.*（t>=0&t<=2）+2*（t>2&t<=3）+（-1）*（t>3&t<=5）;

t=0:

0.001:

5;

subplot（2,1,1）

plot（t,x（t））

title（'x（t）'）

axis（[0,5,-1,2]）

quad（'x2_1（t）',0,5）

ans=

2.0000

>>h=0.0001;

>>t=0:

h:

5;

>>y=diff（1*t.*（t>=0&t<=2）+2*（t>2&t<=3）+（-1）*（t>3&t<=5））/h;

>>plot（t（1:

length（t）-1）,y）

3、信号

，用MATLAB符号运算的相关命令求：

、

、

、

、

，并绘出其时域波形。

程序如下：

symst%定义符号变量t

f=sym（’（t/2+1）*（heaviside（t+2）-heaviside（t-2））’）%创建符号表达式f，heaviside函数见第二讲

subplot（2,3,1）,ezplot（f,[-3,3]）%subplot把图形窗口分成2行3列的6个小窗口，当前波形画在第1个小窗口中，下同

y1=subs（f,t,t+2）%subs（）函数指将连续信号f中的时间变量t用t+2替换，下同

subplot（2,3,2）,ezplot（y1,[-5,1]）

y2=subs（f,t,t-2）

subplot（2,3,3）,ezplot（y2,[-1,5]）

y3=subs（f,t,-t）

subplot（2,3,4）,ezplot（y3,[-3,3]）

y4=subs（f,t,2*t）

subplot（2,3,5）,ezplot（y4,[-2,2]）

y5=-f

subplot（2,3,6）,ezplot（y5,[-3,3]）

请注意：

上述命令执行后，f、y1、y2、y3、y4、y5的符号表达式（在Command窗口中显示）。

]

2.4MATLAB的绘图功能

1）绘制曲线y=x3+x+1，x的取值范围为[-5,5]。

closeall;

symsx

x=-5:

0.0001:

5;

y='x^3+x+1';

ezplot（y）

2）有一组测量数据满足y=e-at，t的变化范围为0~10，用不同的线型和标记画出a=0.1,a=0.2,a=0.5三种情况下的曲线；在结果图中添加标题“y=e-at”，并用箭头线标识出各曲线a的取值；在结果图中添加图例框。

提示：

（1）添加标题语句：

title（'\ity\rm=e^{-\itat}'）

title（'\ity\rm=e^{-\itat}','FontSize',12）

（2）箭头线标识语句：

text（t（6）,y1（6）,'\leftarrow\ita\rm=0.1','FontSize',11）

（3）图例框语句：

legend（'a=0.1','a=0.2','a=0.5'）

closeall;

t=linspace（0,10,100）;

y1=1*exp（0.1*t）;

y2=1*exp（0.2*t）;

y3=1*exp（0.5*t）;

plot（t,y1,'b',t,y2,'g',t,y3,'m'）

closeall;t=linspace（0,10,100）;

y1=1*exp（-0.1*t）;

y2=1*exp（-0.2*t）;

y3=1*exp（-0.5*t）;

plot（t,y1,'b',t,y2,'g',t,y3,'m'）

title（'\ity\rm=e^{-\itat}'）

title（'\ity\rm=e^{-\itat}','FontSize',12）

text（t（6）,y1（6）,'\leftarrow\ita\rm=0.1','FontSize',11）

text（t（6）,y2（6）,'\leftarrow\ita\rm=0.2','FontSize',11）

text（t（6）,y3（6）,'\leftarrow\ita\rm=0.5','FontSize',11）

legend（'a=0.1','a=0.2','a=0.5'）

3）用subplot语句在一个图形窗口上开多个大小不等的子窗口进行绘图并添加注释，见图。

subplot（'position',[0.1,0.15,0.3,0.65]）

hist（randn（1,1000）,20）;

xlabel（'直方图'）

subplot（'position',[0.45,0.52,0.25,0.28]）

[xp,yp,zp]=peaks;

contour（xp,yp,zp,15,'k'）

holdon

pcolor（xp,yp,zp）

shadinginterp

holdoff

axisoff

text（-1.2,-4,'伪彩色图'）

subplot（'position',[0.72,0.5,0.25,0.3]）

sphere（25）;

axisequal,axis（[-0.75,0.75,-0.75,0.75,-0.75,0.75]）

light（'Position',[132]）;

light（'Position',[-3-13]）;

materialshiny

axisoff

text（-0.8,-0.7,-1,'三维图'）

subplot（'position',[0.45,0.15,0.5,0.25]）

t=0:

pi/15:

pi;

y=sin（4*t）.*sin（t）/2;

plot（t,y,'-bs','LineWidth',2,...%设置线型

'MarkerEdgeColor','k',...%设置标记点边缘颜色

'MarkerFaceColor','y',...%设置标记点填充颜色

'MarkerSize',5）

axis（[0,3.14,-0.5,0.5]）

xlabel（'带标记点的线图'）

subplot（'position',[0.1,0.9,0.8,0.1]）

text（0.25,0.2,'多窗口绘图示例',...

'fontsize',25,'fontname','隶书','color','b'）

axisoff

2.5Simlink基础

仿真一个连续三阶系统（自定义），当输入为阶跃信号时，查看示波器（Scope）的输出；当输为为其它信号时，查看示波器的输出。

3结束语

正文排版注意事项：

1.符号：

正文、图表中的变量都要用斜体，英文缩写、计量单位、函数名称、运算符号、括号等都要用正体；相同参数与变量用统一的形式进行表示（注意大小写与上下标）；向量、矢量与矩阵应为粗斜体；文中的计量单位一律使用《中华人民共和国法定计量单位》。

2.公式：

公式须用字符格式录入。

（正文中公式务必用公式编辑器编辑）

3.插图：

要有图序、图名，插图要精绘，线条清晰，粗细一致，图中文字书写清楚，最好是可活动字体（即图和字可以单独分开修改）；图的整体效果要紧凑，大小合理。

4.表格：

要有表序、表名，用三线表（表格的左、右端不封），表中数字采用小数点对齐式。

5.论文采用WORD文档，A4页面，正文采用宋体、小四号字，行距1.25倍。

字数应以3000~5000字为宜。

参考文献（5号宋体）

[1]作者.题名[M].译者.版本.出版地:

出版社,出版年:

页次.（图书）

[2]作者.题名[J].期刊名,出版年,卷号（期号）:

页次.（期刊）

[3]作者.题名[C]//编者.文集名.会议名,会议地址,会议时间.出版地:

出版社,出版年,卷号（期号）:

页次.（会议录）

[4]作者.题名[D].学位授予单位所在城市:

学位授予单位,年.（学位论文）

[5]作者.标准号题名[S].出版地:

出版社,出版年.（标准）

[6]作者.题名:

其他题名,编号[R].出版地:

出版社,出版年.（科学技术报告）

[7]专利所有者.题名,专利号[P].公告日期.（专利）

[8]作者.题名[文献类型标志/文献载体标志].出版地:

出版社,出版年（更新或修改日期）[引用日期].访问路径.（电子文献）

注：

作者采用姓在前名在后的形式，欧美作者的名可用缩写字母。

作者不超过3个人时，需全部照录，超过3个人时，只需著录前3个作者，其后加“，等”或“，etal”。

上交材料：

1、报告电子版；

2、报告纸质材料。

以上材料由各班学委于12月30日之前收齐，再交给任课教师，其中电子版发至任课教师邮箱，纸质材料于12月31日上午10:

00交到田C402。