物理与电子信息工程学院实验报告册姓名陈楠学号10110003306班级10通信本课程名称通信系统仿真与实践指导教师李理敏第2012/2013学年第2学期上机作业一：1、写出可以完成下列任务的Matlab命令(函数):（1）Matlab中的圆周率约等于多少？evalpi：圆周率，其值为imag(log(-1)).约3.1415926535897（2）Matlab中的浮点运算相对精度是多少？eps：浮点运算相对精度（3）浮点数表示范围为：10-308～10308，但具体是多少？（4）Matlab中的最大整数、最小整数分别是多少？（5）如何知道命令realmax的具体用法？（6）如何知道函数max的作用？并举例加以说明。C=max(A)C=max(A,B)C=max(A,[],dim)[C,I]=max(...)2、回答以下问题：（1）Matlab是什么含义？Matlab:MatrixLaboratory矩阵实验室（2）Matlab命令who与whos有什么区别？Who：显示工作空间中的所有变量Whos：查看工作空间中变量的详细属性（3）如何使用Matlab命令补全功能？命令补全功能：Tab键（4）Matlab命令clear、clc与home有什么区别？clear：清除工作空间的所有变量；clc：清除命令窗口的内容，对工作环境中的全部变量无任何影响；home：将光标移到左上角，腾出命令窗口，以前的代码不会删除，只会上移；第二讲：5阶的单位阵8阶均匀分布的随机矩阵及其下三角矩阵生成列向量x=[1,3,5,7,9,…,99]生成以x的前5个元素为对角线的矩阵A生成一个与A同阶的正态分布的随机矩阵B计算A的转置与B的kronecker乘积矩阵C生成由B与其同阶魔方矩阵点乘得到的矩阵D生成由D的第2、4、5行和第4、1列组成的子矩阵E求出矩阵E中绝对值最大的元素的绝对值第三讲：第三章习题3、设有一系统，其输入/输出关系由差分方程确定，即y(n)=0.95y(n-1)+x(n)-0.5x(n-1)(1)求该系统的单位取样响应(2)若x(n)=,n,求该系统的输出并画图。代码如下：clearalla=[10.95];b=[1-0.5];n=0:30;impz(b,a,30),gridontitle('系统单位取样响应h(n)')clearalla=[10.95];b=[1-0.5];n=0:30;x=0.3.^n;y=filter(b,a,x);stem(n,y,'fill'),gridonxlabel('n'),title('输出响应)')结果如下；5.，求该信号的Fourier变换。代码如下：clearallsymst;f=exp(-t*t);F=fourier(f);subplot(1,1,1);ezplot(F);结果如下：6、一个信号的Fourier变换是F(W)=1，试求该信号的时域表达式。代码如下：clearallsymstw;F=1;f=ifourier(F,t)subplot(1,1,1);ezplot(f);9、一个带通信号的定义如下：X(t)=6(cos(20πt)-3sin(30πt))cos(400πt+3/8πt)画出x(t)的普函数。画出该信号解析函数的普函数。载波频率为f=200Hz，求该信号的低通等效信号。clearallts=0.001;fs=1/ts;t=0.1:ts:0.2;df=fs/length(t);f=-500:df:500-df;x=6*(cos(20*pi*t)-3*sin(30*pi*t)).*cos(400*pi*t+3/8*pi*t);X=fft(x)/fs;xa=hilbert(x);Xa=fft(xa)/fs;subplot(2,1,1);plot(t,x);title('信号x');xlabel('时间t')subplot(2,1,2);plot(f,fftshift(abs(X)));title('信号x幅度谱');xlabel('频率f')figuresubplot(2,1,1);plot(t,abs(xa));title('信号xa包络');xlabel('时间t')subplot(2,1,2);plot(f,fftshift(abs(Xa)));title('信号xa幅度谱');xlabel('频率f')clearallts=0.001;fs=1/ts;t=0