(完整word版)Matlab实验报告
MATLAB实验报告

MATLAB实验报告一、实验目的本次 MATLAB 实验旨在深入了解和掌握 MATLAB 软件的基本操作和应用,通过实际编程和数据处理,提高解决问题的能力,培养编程思维和逻辑分析能力。
二、实验环境本次实验使用的是 MATLAB R2020a 版本,运行在 Windows 10 操作系统上。
计算机配置为英特尔酷睿 i5 处理器,8GB 内存。
三、实验内容(一)矩阵运算1、矩阵的创建使用直接输入、函数生成和从外部文件导入等方式创建矩阵。
例如,通过`1 2 3; 4 5 6; 7 8 9` 直接输入创建一个 3 行 3 列的矩阵;使用`ones(3,3)`函数创建一个 3 行 3 列元素全为 1 的矩阵。
2、矩阵的基本运算包括矩阵的加减乘除、求逆、转置等。
例如,对于两个相同维度的矩阵`A` 和`B` ,可以进行加法运算`C = A + B` 。
3、矩阵的特征值和特征向量计算通过`eig` 函数计算矩阵的特征值和特征向量,加深对线性代数知识的理解和应用。
(二)函数编写1、自定义函数使用`function` 关键字定义自己的函数,例如编写一个计算两个数之和的函数`function s = add(a,b) s = a + b; end` 。
2、函数的调用在主程序中调用自定义函数,并传递参数进行计算。
3、函数的参数传递了解值传递和引用传递的区别,以及如何根据实际需求选择合适的参数传递方式。
(三)绘图功能1、二维图形绘制使用`plot` 函数绘制简单的折线图、曲线等,如`x = 0:01:2pi; y = sin(x); plot(x,y)`绘制正弦曲线。
2、图形的修饰通过设置坐标轴范围、标题、标签、线条颜色和样式等属性,使图形更加清晰和美观。
3、三维图形绘制尝试使用`mesh` 、`surf` 等函数绘制三维图形,如绘制一个球面`x,y,z = sphere(50); surf(x,y,z)`。
(四)数据处理与分析1、数据的读取和写入使用`load` 和`save` 函数从外部文件读取数据和将数据保存到文件中。
MATLAB实验报告(1-4)

信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂,并在MATLAB软件中练习例题,最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算:算数运算向量运算:1.向量元素要用”[ ]”括起来,元素之间可用空格、逗号分隔生成行向量,用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值,step表示步长或者增量,xn 为结束值。
矩阵运算:1.矩阵”[ ]”括起来;矩阵每一行的各个元素必须用”,”或者空格分开;矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例:计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算:定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语:>>plot(x,y)2).输出多个图像表顺序:例如m和n表示在一个窗口中显示m行n列个图像,p表示第p个区域,表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围:axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母:xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前:title(‘……’)6).输出:grid on举例1:举例2:3.matlab程序流程控制1).for循环:for循环变量=初值:增量:终值循环体End2).while循环结构:while 逻辑表达式循环体End3).If分支:(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号:y=k*exp(a*t)举例:2).正弦信号:y=k*sin(w*t+phi)3).复指数信号:举例:4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号:y=tripuls(t,width,skew)(skew的取值在-1~+1之间,若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算:语句格式:d iff(function,’variable’,n)Function:需要进行求导运算的函数,variable:求导运算的独立变量,n:求导阶数2).连续时间信号的积分运算:语句格式:int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件,并且能够输入简单的语句并输出相应的结果和波形图,也在一定程度上巩固了c语言的一些语法。
MATLAB实验报告(四份!!!)

0380160815880001449451490168779983287916534749293454515731] >> double(x) ans = -33.5325 60.6078 -54.7098 27.3192 实验三 MATLAB 绘图
一 实验环境 计算机 MATLAB 软件 二 实验目的 1. 掌握 MATLAB 的基本绘命令。 2. 掌握运用 MATLAB 绘制一维、二维、三维的图形方法。 3. 给图形加以修饰。 三 预备知识 特殊的二维图形函数有:直方图、柄图、阶梯图、饼图、频数累计柱状 图、极坐标图 四 实验内容和步骤
3.创建一个 3x3 矩阵,并求其转置,逆矩阵。 >> A=magic(3) A= 8 3 4 >> A' ans = 8 1 6 >> inv(A) 3 5 7 4 9 2 1 5 9 6 7 2
ans = 0.1472 -0.0611 -0.0194 -0.1444 0.0222 0.1889 0.0639 0.1056 -0.1028
3. 画出 y=的曲线(x(-5,5)) 。在这曲线上加入相同区间里的 y=的曲线,并且采用绿色折 线标识。 >> x=linspace(-5,5,100); >> y=x.^2; >> plot(x,y) >> hold on >> z=x.^(1/3); >> plot(x,z,'g--')
4. 在同一窗口不同坐标系里分别绘出 y1=sinx,y2=cosx,y3=cinh(x),y4=cosh(x)4 个图形。 >> x=linspace(0,2*pi,30); >> subplot(2,2,1);plot(x,sin(x)); >> subplot(2,2,2);plot(x,cos(x)); >> subplot(2,2,3);plot(x,sinh(x)); >> subplot(2,2,4);plot(x,cosh(x));
MATLAB实验报告

实验一、Matlab 软件的安装与基本使用一、实验目的和要求:学会MA TLAB 软件的安装方法,熟悉MA TLAB 的工作环境及其基本功能的设置,以及基本指令的使用,了解MA TLAB 数值数组及其矩阵操作的使用、复数的输入方法。
二、实验内容和步骤:1、完成MA TLAB 软件的安装。
2、安装完成后启动MA TLAB 软件,熟悉MA TLAB 的工作环境。
打开MA TLAB 7,默认打开的窗口包括: (1)命令窗口(Command Window); (2)命令历史窗口(Command History); (3)工作间管理窗口(Workspace);(4)当前路径窗口(Current Directory)此外,还有捷径键、编译窗口、图形窗口和帮助窗口等其他种类的窗口,结合课本实例,熟悉每个窗口的功能和基本使用。
3、(1)求23)]47(212[÷-⨯+的运算结果,注意命令后有无分号的作用。
命令及结果如下:t=[12+2*(7-4)]/3^2t = 2(2)求20012x x v t at =++的值,其中t=5s 、010x m =、015/v m s =、29.81/sec a m =-。
命令及结果如下:syms t x0 v0 a; t=5,x0=10,v0=15,a=-9.18; x=x0+v0*t+1/2*a*t^2 t = 5 x0 = 10 v0 = 15 x = -29.7500(3) 设3x =, 4v =,用MA TLAB 命令求下列各式的值:1)22log()x v + 命令及结果如下:syms x v x=3,v=4; log(x^2+v^2)x = 3 ans =3.21892)2(2)x v - 命令及结果如下:syms x vx=3,v=4;f=(sqrt(x-3))/(x-2*v)^2 x = 3 f = 0 3) |sin 2|vx e命令及结果如下: syms x v x=3,v=4;t=abs(2*x)*exp(v) x = 3t = 327.58894、用三种不同的方法输入矩阵1234246836912a ⎛⎫ ⎪= ⎪ ⎪⎝⎭,运行并观察实验结果,矩阵的三种基本输入方法如下:a=[1,2,3,4;3,4,6,8;3,6,9,12] a =1 2 3 4 3 4 6 8 3 6 9 12>> a=[1,2,3,43,4,6,8 3,6,9,12]a =1 2 3 4 3 4 6 8 3 6 9 12>> >>利用M 文件创建 1 2 3 4 3 4 6 8 3 6 9125、输入复数矩阵⎥⎦⎤⎢⎣⎡-+++=i ii i24356223B ,并求其实部、虚部、模和幅角(分别用弧度和角度表示),了解MA TLAB 向量化运算的特点。
MATLAB实验报告(1-4)

信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂,并在MATLAB软件中练习例题,最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算:算数运算向量运算:1.向量元素要用”[ ]”括起来,元素之间可用空格、逗号分隔生成行向量,用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值,step表示步长或者增量,xn为结束值。
矩阵运算:1.矩阵”[ ]”括起来;矩阵每一行的各个元素必须用”,”或者空格分开;矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例:计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算:定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语:>>plot(x,y)2).输出多个图像表顺序:例如m和n表示在一个窗口中显示m行n列个图像,p表示第p个区域,表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围:axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母:xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前:title(‘……’)6).输出:grid on举例1:举例2:3.matlab程序流程控制1).for循环:for循环变量=初值:增量:终值循环体End2).while循环结构:while 逻辑表达式循环体End3).If分支:(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号:y=k*exp(a*t)举例:2).正弦信号:y=k*sin(w*t+phi)3).复指数信号:举例:4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号:y=tripuls(t,width,skew)(skew的取值在-1~+1之间,若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算:语句格式:d iff(function,’variable’,n)Function:需要进行求导运算的函数,variable:求导运算的独立变量,n:求导阶数2).连续时间信号的积分运算:语句格式:int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件,并且能够输入简单的语句并输出相应的结果和波形图,也在一定程度上巩固了c语言的一些语法。
(完整word版)Matlab数学实验报告

Matlab 数学实验报告一、实验目的通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令,通过建立数学模型解决理论或实际问题。
了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。
二、实验内容2.1实验题目一2.1.1实验问题Feigenbaum曾对超越函数y=λsin(πx)(λ为非负实数)进行了分岔与混沌的研究,试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图2.1.2程序设计clear;clf;axis([0,4,0,4]);hold onfor r=0:0.3:3.9x=[0.1];for i=2:150x(i)=r*sin(3.14*x(i-1));endpause(0.5)for i=101:150plot(r,x(i),'k.');endtext(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end加密迭代后clear;clf;axis([0,4,0,4]);hold onfor r=0:0.005:3.9x=[0.1];for i=2:150x(i)=r*sin(3.14*x(i-1));endpause(0.1)for i=101:150plot(r,x(i),'k.');endend运行后得到Feigenbaum图2.2实验题目二2.2.1实验问题某农夫有一个半径10米的圆形牛栏,长满了草。
他要将一头牛拴在牛栏边界的桩栏上,但只让牛吃到一半草,问拴牛鼻子的绳子应为多长?2.2.2问题分析如图所示,E为圆ABD的圆心,AB为拴牛的绳子,圆ABD为草场,区域ABCD为牛能到达的区域。
问题要求区域ABCD等于圆ABC的一半,可以设BC等于x,只要求出∠a和∠b就能求出所求面积。
Matlab数学实验报告模板

数学实验报告
日期:20年月日
班级姓名学号
实验(一)猜数字游戏
任务描述:
首先由计算机随机产生一个[1,100]之间的一个整数,然后由用户猜测所产生的这个数。
根据用户猜测的情况给出不同的提示,如果猜测的数大于产生的数,则显示“High”,小于则显示“Low”,等于则显示“You won!”,同时退出游戏。
用户最多有7次机会。
实验原理与数学模型:
实验所用软件及版本:
主要内容(要点):
1
实验过程记录(含基本步骤、主要程序清单及异常情况记录等):实验结果报告与实验总结:
思考与深入:
教师评语:
实验(二)
问题背景描述:
实验目的:
实验原理与数学模型:
实验所用软件及版本:
主要内容(要点):
实验过程记录(含基本步骤、主要程序清单及异常情况记录等):
实验结果报告与实验总结:思考与深入:
教师评语:。
(完整word版)武汉大学MATLAB报告

Matlab实训任务书学生姓名:好人专业班级:…………班指导教师:好老师工作单位:信息工程学院题目: 离散时间信号和系统时域分析及MATLAB实现课题要求:深入研究离散时间信号和系统时域分析的理论知识。
利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能,实现离散时间信号和系统时域分析的仿真波形。
课题内容:一、用MATLAB绘制常用信号的时域波形(通过改变参数分析其时域特性)1、单位序列,2、单位阶跃序列,3、正弦序列,4、离散时间实指数序列,5、离散时间虚指数序列,6、离散时间复指数序列。
二、用MATLAB实现信号的时域运算1、相加,2、相乘,3、数乘。
三、用MATLAB实现信号的时域变换(参数变化,分析波形的变化)1、反转,2、时移(超时,延时),3、展缩,4、倒相。
四、用MATLAB实现离散时间系统卷积和仿真波形给出几个典型例子,对每个例子要求画出e(k),h(k),e(i),h(i),h(-i),Rzs(k)波形。
五、用MATLAB实现离散时间系统的单位响应,阶跃响应的仿真波形给出几个典型例子,四中调用格式。
六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形给出几个典型例子,要求可以改变激励的参数,分析波形的变化。
指导教师签名:2013年 1 月24 日系主任(或责任教师)签名:2012年 1 月24 日摘要 (I)Abstract (II)1 绪论 (1)2 MATLAB简介 (1)3 MATLAB实验内容 (2)3.1用MATLAB绘制常用信号的时域波形............................................ 错误!未定义书签。
3.1.1单位序列 (2)3.1.2单位阶跃序列标 (4)3.1.3正弦序列 (4)3.1.4 离散时间实指数序列 (5)3.1.5离散时间虚指数序列 (6)3.1.6离散时间复指数序列 (7)3.2用MATLAB实现信号的时域运算 (8)3.2.1 相加 (8)3.2.2 相乘 (8)3.2.3数乘 (9)3.3用MATLAB实现信号的时域变换 (11)3.4用MATLAB实现离散时间系统卷积和仿真波形 (12)3.5 实现离散时间系统的单位响应,阶跃响应的仿真波形 (13)3.6实现离散时间系统对实指数序列信号的零状态响应的仿真波形 (16)4实验小结 (17)5参考文献 (18)随着计算机软硬件技术的飞速发展,集成电路被广泛应用,电路越来越复杂,集成电路越来越高,加之元器件层出不穷,使得越来越多的工作已经无法依靠手工来完成。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
图形加注
Title画题头
Xlabel x轴标注
Ylabel y轴标注
Text任意定位的标注
Gtext鼠标定位标注
Dgrid网格
关于坐标轴尺寸的选择和图形处理等控制命令:
图形控制命令
Axis人工选择坐标轴尺寸
Clg清除图形窗口
figure(1)
subplot(311)
plot(t,s_am);holdon;
plot(t,A + mt,'r--');
title('AM调制信号及其包络');
xlabel('t');
rt = s_am.*cos(2 * pi * fc * t);
rt = rt - mean(rt);
[f,rf] = T2F(t,rt来自;plot(f,psf);
axis([ -2*fc 2*fc 0 max(psf)]);
title(AM信号功率谱');
xlabel('f');
实验内容
<设计性实验>
1、用Matlab产生一个频率为2Hz、功率为1的正弦信源m(t),设载波频率为10Hz,A分别等于0.5、1、1.5,试画出:
(1)AM调制信号;
实验原理
实验1程序:
dt = 0.001; fm = 2; fc = 10; T = 5; t = 0:dt:T;
mt = sqrt(2) * cos(2 * pi * fm * t);
A = 0.5;
s_am = (A + mt).* sin(2 * pi * fc * t);
B = 2 * fm;
(2)调制信号的功率谱密度;
(3)相干解调后的信号波形(选做)
(4)研究调制深度与失真的关系。(选做)
<设计性实验>(选做)
2、消息 是周期为2s的周期信号,它在区间 的定义如下:
该信号用DSB方案调制50Hz的载波。
(1)绘出调制信号。
(2)绘出调制信号的频谱。
(3)比较调制信号和未调制信号的频谱。
disp(str1)
disp(str2)
实验4结果为:
实验5程序:
n=input('n');
sum=1;
if(n==0)
break;
else
form=1:n
sum=sum*m;
end
end
sum
仿真结果:
n5
sum= 120
实验名称
实验二:Matlab在通信原理中的应用
实验目的
1、掌握通信原理课程中基本知识的Matlab编程、仿真方法
实验环境
MATLAB7.0
实验名称
实验一:熟悉MATLAB环境及Matlab基本操作
实验目的
1、熟悉MATLAB环境及各种基本Matlab操作
实验原理
MATLAB是以复杂矩阵作为基本编程单元的一种程序设计语言。它提供了各种矩阵的运算与操作,并有较强的绘图功能。
利用FILE菜单可方便对文件或窗口进行管理。其中FILE|NEW的各子菜单,M-FILE(M文件)、FIGURE(图形窗口)、或MODEL(SIMULINK编辑界面)分别可创建对应文件或模块。EDIT菜单允许用户和WINDOWS的剪切板交互信息。
end
disp(a)
end
实验结果及分析
实验3仿真结果:
实验4程序如下:
clearall
n=0;
sum=0;
while(sum<1000)
n=n+1;
sum=sum+n;
end
str1 = ['最后一个自然数为:',num2str(n)];
str2 = ['计算结果为:',num2str(sum)];
<设计性实验>
3、用MATLAB程序建立九九乘法表。
4、求连续自然数的和,当和大于等于1000时,求最后一个自然数以及自然数的和。
5、MATLAB函数是可以递归调用的,即在函数的内部调用函数自身,试用递归方式编写n!函数。
实验结果及分析
实验2:程序如下:
clearall
sum=0;
forn=1:2:97
在MATLAB中把数据绘成图形可有多种命令一供选择。下面列出了这些命令:
绘图命令
Plot线性X-Y坐标图
Loglog双对数坐标图
Semilogx X轴对数半对数坐标图
Semilogy Y轴对数半对数坐标图
Polar极坐标图
Mesh三维消隐图
Contour等高线图
Bar条形图
Stairs阶梯图
实验原理
MATLAB语言最基本的赋值语句结构为:变量名列表=表达。表达式由操作符或其它字符,函数和变量名组成,表达式的结果为一个矩阵,显示在屏幕上,同时输送到一个变量中并存放于工作空间中以备调用。如果变量名和“=”省略,则ANS变量将自动建立,例如键入:1900/81,得到输出结果:ans =23.4568。
sum=sum+n;
end
sum=sum+99
disp('程序运行结束!')
实验结果如下。
实验2结果分析:
这道题有多种方法,例如如下程序:
clearall
sum=0;
fori=1:2:100
sum=sum+i;
end
sum
这道题需要注意的一点是语句sum=sum+i后要加分号。这样就不会出现i每赋值一次窗口就显示一次sum值。
Ginput利用鼠标的十字准线输入
Hold保持图形
Shg显示图形窗口
Subplot将图形窗口分成N块子窗
实验内容
<验证性实验>
1、将Matlab当前目录变更到D:\自己姓名拼音文件夹下,并建立新文件画出y=2 * x的曲线,并将曲线图形保存在当前目录下。
2、使用for语句及while语句求解1至100整数中奇数的和。
实验结果及分析
1、A=0.5
2、A=1
实验结果及分析
A=1.5
实验名称
实验三Matlab在信号与系统中的应用
实验目的
1、掌握信号与系统课程中基本知识的Matlab编程、仿真方法
实验3:程序如下:
clearall
forx=1:9
a=[];
fory=1:x
term=([num2str(y),'*',num2str(x),'=',num2str(x*y,'%2d')]);
ifx*y<10
term=[term,' '];
else
term=[term,' '];
end
a=[a,term];
[t,rt] = lpf(f,rf,2 * fm);
subplot(312)
plot(t,rt);holdon;
plot(t,mt/2,'r--');
title('相干解调后的信号波形与输入信号的比较');
xlabel('t')
subplot(313)
[f,sf]=T2F(t,s_am);
psf = (abs(sf).^2)/T;