Matlab与通信仿真实验报告册答案

合集下载

MATLAB通信建模实验仿真实验报告

实验1：上采样与内插一、实验目的1、了解上采样与内插的基本原理和方法。

2、掌握上采样与内插的matlab程序的设计方法。

二、实验原理上采样提高采样频率。

上采样使得周期降低M倍，即新采样周期Tu和原有采样周期Ts的关系是T u=T s/M,根据对应的连续信号x(t),上采样过程从原有采样值x(kT s)生成新采样值x(kT u)=x(kT s/M)。

操作的结果是在每两个采样值之间放入M-1个零值样点。

更实用的内插器是线性内插器，线性内插器的脉冲响应定义如下：上采样值x(kT u)=x(kT s/M)通过与线性内插器的脉冲响应的卷积来完成内插。

三、实验内容仿真正弦波采样和内插，通过基本采样x（k）,用M=6产生上采样x u(k)，由M=6线性内插得到样点序列x i（k）。

四、实验程序% File: c3_upsampex.mM = 6; % upsample factorh = c3_lininterp(M); % imp response of linear interpolatort = 0:10; % time vectortu = 0:60; % upsampled time vectorx = sin(2*pi*t/10); % original samplesxu = c3_upsamp(x,M); % upsampled sequencesubplot(3,1,1)stem(t,x,'k.')ylabel('x')subplot(3,1,2)stem(tu,xu,'k.')ylabel('xu')xi = conv(h,xu);subplot(3,1,3)stem(xi,'k.')ylabel('xi')% End of script file.% File: c3_upsample.mfunction out=c3_upsamp(in,M)L = length(in);out = zeros(1,(L-1)*M+1); for j=1:Lout(M*(j-1)+1)=in(j); end% End of function file.% File: c3_lininterp.m function h=c3_lininterp(M) h1 = zeros(1,(M-1)); for j=1:(M-1) h1(j) = j/M; endh = [0,h1,1,fliplr(h1),0]; % End of script file.四、实验结果012345678910x0102030405060x u01020304050607080-11x i实验二：QPSK、16QAM信号的散点图、正交、同相分量波形图一、实验目的1、了解QPSK、16QAM调制的基本原理。

中国地质大学基于MATLAB通信仿真实验报告

中国地质大学（武汉）通信原理matlab仿真报告姓名: 张彪_班号: 075112 _学号：20111002253院系：_机电学院专业：_通信工程指导教师：_赵娟老师联系方式：预留指导时间：1、选题为3.周期信号m(t)，周期为2，在区间[0,2]内定义为：m(t)=t , 1=<t=<2；m(t)=2-t , 1=<t<=2。

（1）以时间间隔为0.1 的步长对此信号进行8 电平的均匀量化，画出量化输出波形；（2）求各量化点的量化误差，画出量化误差波形；（3）通过计算误差信号的功率，求该系统的SNRq（以dB 计）。

2、解题原理与思路模拟信号的取值范围：a～b，量化电平数＝M则均匀量化时的量化间隔为：(b-a)/M且量化区间的端点mi=a+i v i=0,1,…,M 若量化输出电平qi取量化间隔的中点，则qi=[mi+m （i-1）]/2 i=1,2…,M .我的解题思路很简单，大致分为三部分，第一部分构造出题中所给的分段函数；第二部分自定义均匀量化函数，并一并写入计算信噪比等功能；第三部分将第一部分带入自定义函数并求出量化误差波形，最后输出波形，完成题目要求。

思路框图3、实验代码及注释：（1）分段函数的构造并代入自定义函数及求量化误差。

clear all;x1=0:0.1:1;y1=x1;x2=1:0.1:2y2=-x2+2;y3(1:11)=y1;y3(11:21)=y2; %构造出分段函数x=0:0.1:2;plot(x,y3); %输出分段函数w=juny(y3,1,8); %代入自定义函数求值plot(x,y3,x,w); %输出原波形与量化波形title('原波形与量化波形')figure(2)z=abs(y3-w); %求量化误差plot(x,z) %输出量化误差波形title('量化误差波形')（2）自定义函数：量化函数的实现function h=juny(f,v,L)%f是量化的信号序列%v是量化电平最大值%L是量化电平数n=length(f);t=2*v/L; %纵坐标取间隔p=zeros(1,L+1); %确定范围及初始化for i=1:L+1,p(i)=-v+(i-1)*t; %确定原信号抽样值endfor i=1:n %求量化值if f(i)>=v %在此范围内的量化取值（大于量化电平最大值时）h(i)=v;endif f(i)<=-v %在此范围内的量化取值（小于量化电平最小值时）h(i)=-v;endflag=0; %设立标志for j=2:L/2+1 %对量化电平数即纵坐标取半if flag==0if f(i)<p(j) %原信号与抽样信号循环比较h(i)=(p(j-1)+p(j))/2; %确定量化电平flag=1;endendendfor j=L/2+2:L+1 %另一部分的量化（同上）if flag==0if f(i)<p(j)h(i)=(p(j-1)+p(j))/2;flag=1;endendendendnq=v^2/(3*L^2); %计算量化噪声功率nq %显示量化噪声功率Am=2;snr=(3/2)*(Am/v)^2*L^2; %计算信噪比snr %显示信噪比snr=10*log(snr)/log(10); %将信噪比改为dB表示snr4、运行结果：运算结果：nq =0.0052 snr =384 snr =25.8433dB5、分析：在实验过程中出现许多问题，比如分段函数的构造写了一下午，本来是很简单，可以直接是一个三角函数，调用即可，不过如果是其他形式的分段函数就不能了，所以过程中用if 语句等等去做，不过都没做好，最后想起来用的这种方法构造的。

Matlab实验指导书(含答案)详解

实验一：Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境。

2．学习使用图形函数计算器命令funtool及其环境。

二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式;学会使用变量和矩阵的输入，并进行简单的计算;学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算:1．单函数运算操作。

➢求下列函数的符号导数(1) y=sin（x);（2）y=（1+x）^3＊（2-x）；➢求下列函数的符号积分(1) y=cos(x)；（2) y=1/（1+x^2）；（3）y=1/sqrt(1—x^2）；（4) y=（x-1）/(x+1)/（x+2)；➢求反函数(1) y=(x—1)/（2＊x+3)；(2) y=exp(x）;（3）y=log（x+sqrt（1+x^2）);➢代数式的化简（1) （x+1)*（x—1)＊（x—2）/(x—3)/（x-4)；（2) sin（x）^2+cos(x)^2;（3）x+sin(x)+2＊x—3＊cos（x)+4＊x*sin(x）;2．函数与参数的运算操作。

➢从y=x^2通过参数的选择去观察下列函数的图形变化（1）y1=（x+1)^2（2）y2=（x+2)^2（3) y3=2＊x^2（4）y4=x^2+2(5）y5=x^4（6）y6=x^2/23．两个函数之间的操作➢求和(1) sin（x）+cos(x)（2) 1+x+x^2+x^3+x^4+x^5➢乘积(1) exp（—x）＊sin(x）（2）sin(x)＊x➢商(1）sin（x）/cos(x）;（2）x/(1+x^2)；(3）1/(x-1）/(x-2）;➢求复合函数(1) y=exp（u） u=sin（x)（2）y=sqrt（u) u=1+exp（x^2)(3）y=sin（u）u=asin（x）（4) y=sinh(u）u=—x三、设计提示1．初次接触Matlab应该注意函数表达式的文本式描述。

MATLAB 课后实验答案

实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

(1) 0122sin851z e =+ (2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦(3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解： M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 运算结果：z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)z1 =0.2375z2 =0.7114 - 0.0253i 0.8968 + 0.3658i0.2139 + 0.9343i 1.1541 - 0.0044iz3 =Columns 1 through 40.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416iColumns 5 through 80.7822 + 3.1416i 0.7602 + 3.1416i 0.7254 + 3.1416i 0.6784 + 3.1416iColumns 9 through 120.6196 + 3.1416i 0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780 + 3.1416iColumns 13 through 160.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i -0.0771 + 3.1416iColumns 17 through 20-0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416iColumns 21 through 24-0.8536 + 3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416iColumns 25 through 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245 Columns 29 through 32-3.0017 -2.3085 -1.8971 -1.5978Columns 33 through 36-1.3575 -1.1531 -0.9723 -0.8083Columns 37 through 40-0.6567 -0.5151 -0.3819 -0.2561Columns 41 through 44-0.1374 -0.0255 0.0792 0.1766Columns 45 through 480.2663 0.3478 0.4206 0.4841Columns 49 through 520.5379 0.5815 0.6145 0.6366Columns 53 through 560.6474 0.6470 0.6351 0.6119Columns 57 through 600.5777 0.5327 0.47740.4126Column 610.3388z4 =0 0.2500 0 1.2500 1.0000 2.2500 2. 已知：求下列表达式的值：(1) A+6*B和A-B+I（其中I为单位矩阵）(2) A*B和A.*B(3) A^3和A.^3(4) A/B及B\A(5) [A,B]和[A([1,3],:);B^2]解：M 文件：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]运算结果：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*BA-B+eye(3)A*BA.*BA^3A.^3A/BB\A[A,B][A([1,3],:);B^2]ans =18 52 -1046 7 10521 53 49ans =12 31 -332 8 840 67 1ans =68 44 62309 -72 596154 -5 241ans =12 102 468 0 2619 -130 49ans =37226 233824 48604247370 149188 60076678688 454142 118820 ans =1728 39304 -6439304 343 65850327 274625 343 ans =16.4000 -13.6000 7.600035.8000 -76.2000 50.200067.0000 -134.0000 68.0000109.4000 -131.2000 322.8000-53.0000 85.0000 -171.0000-61.6000 89.8000 -186.2000ans =12 34 -4 1 3 -134 7 87 2 0 33 65 7 3 -2 7ans =12 34 -43 65 74 5 111 0 1920 -5 403. 设有矩阵A和B(1) 求它们的乘积C。

《MATLAB仿真技术》实验指导书2016附答案分析

实验项目及学时安排实验一 MATLAB环境的熟悉与基本运算 2学时实验二 MATLAB数值计算实验 2学时实验三 MATLAB数组应用实验 2学时实验四 MATLAB符号计算实验 2学时实验五 MATLAB的图形绘制实验 2学时实验六 MATLAB的程序设计实验 2学时实验七 MATLAB工具箱Simulink的应用实验 2学时实验八 MATLAB图形用户接口GUI的应用实验 2学时实验一 MATLAB环境的熟悉与基本运算一、实验目的1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识1.熟悉MATLAB环境:MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。

2.掌握MATLAB常用命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英文字母开头（3）长度不大于31个（4）区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。

MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符4.MATLAB的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式5.MATLAB的基本运算表7 两种运算指令形式和实质内涵的异同表6.MATLAB的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）2、学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。

3、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、exerc2、exerc3……），学习使用MA TLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。

《MATLAB与通信仿真》实验指导书-(通信2012级)

电子信息学院《MATLAB与通信仿真》实验指导书刘紫燕编写适用专业：通信工程贵州大学二O一四年二月前言本课程是通信工程专业的选修课程。

课程内容包含MATLAB的基本操作、MATLAB程序设计、函数文件、MATLAB的图形和数据处理、SIMULNK的基本应用及其在通信工程中的应用等。

通过本课程的学习，掌握MATLAB软件使用和编程方法，验证和深化书本知识，从而加强基础知识，掌握基本技能，提高MATLAB软件的编程能力，并用MATLAB软件解决通信工程专业中的实际问题，同时，针对通信工程专业的特点，要求学生掌握使用MATLAB来研究和开发与本专业相关的系统的方法。

本课程设置5个实验，均为设计性实验。

建议实验学时为12学时。

实验一是MATLAB软件的基本操作；实验二是MATLAB程序设计；实验三是MATLAB的图形绘制；实验四是MATLAB的数据处理；实验五是MATLAB/Simulink在电路中的仿真设计。

每个实验2个小时左右，为了使学生更好的掌握实验内容，学生务必要做到以下几点：(1)实验前认真预习实验。

明确实验目的，熟悉实验内容，理论分析实验结果，编写相应的程序代码，并撰写出预习报告。

(2)实验过程中积极思考，深入分析命令、程序的执行过程和执行结果，对比理论分析结果，分析评判实验结果，并把实验中出现的问题及解决方法记录下来。

(3)实验完成后，总结本次实验有哪些收获，还存在什么问题，撰写并提交最终的实验报告。

本指导书实验项目和要求明确，学生容易着手实验并得出实验结果。

本实验指导书适用于通信工程专业。

目录实验一MATLAB的基本操作 (1)实验二MATLAB程序设计 (9)实验三MATLAB的图形绘制 (12)实验四MATLAB的数据处理 (14)实验五MATLAB/Simulink在电路中的仿真设计 (16)实验报告的基本内容及要求 (23)贵州大学实验报告 (24)实验一 MATLAB 的基本操作实验学时：2实验类型：设计实验要求：必做一、实验目的1、熟悉MATLAB 的界面功能和基本操作；2、掌握MATLAB 的基本运算方法；3、掌握MATLAB 的矩阵运算方法；4、掌握MATLAB 中帮助命令的使用方法。

(完整版)MATLAB)课后实验答案[1]

实验一 MATLAB 运算基础1.先求下列表达式的值，然后显示 MATLAB 工作空间的使用情况并保存全部变量解:4.完成下列操作：(1) 求[100,999] 之间能被21整除的数的个数。

(1) z i 2sin 85° 1 e 2Z 2 卯(x L)，其中x2 0.45 2i 5Z 3 0.3a 0.3a e e 2sin (a 0.3)3.0, 2.9,L ,2.9, 3.0 Z 4t 2t 2 t 2 1 2t 其中 t=0:0.5:2.5(2) 建立一个字符串向量，删除其中的大写字母。

解: (1)结果：m=100:999;n=fin d(mod(m,21)==0);len gth( n)ans =43(2).建立一个字符串向量例如:ch二'ABC123d4e56Fg9：则要求结果是: ch二'ABC123d4e56Fg9：k=fi nd(ch>='A'&ch<='Z'); ch(k)=[] ch =123d4e56g9 实验二MATLAB 矩阵分析与处理1.设有分块矩阵A E 3 3 °2 3 R 3 2S 2 2 ，其中E 、R 、0、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证A 2 E R RS 0 S 22 3 4 1 1 1 3 4 5 1 1 1 4 5 6 x , 0.95 x 2 0.67 x 3 0.52(1) 求方程的解。

(2) 将方程右边向量元素 b 3改为0.53再求解，并比较 b 3的变化和解的相对变化。

(3)计算系数矩阵A 的条件数并分析结论。

解：M 文件如下：解：M 文件如下; 5.下面是一个线性方程组:实验三选择结构程序设计1. 求分段函数的值。

x2x 6 x 0且x 3y x2 5x 6 0 x 5且x 2 及x 3x2x 1 其他用if语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0 时的y值。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

1.在Command Window里面计算（1）(358)510++÷⨯； >>(3+5+8)/5*10 ans=32（2）sin(3)π >>sin(3*pi)/sqrt(9/5) ans=2.7384e-16（3）123456789A⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，789456123B⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，计算：,,\,/C A BD A B A C C B=⨯=+；>> a=[1 2 3;4 5 6;7 8 9];>> b=[7 8 9;4 5 6;1 2 3];>> c=a*bc =18 24 3054 69 8490 114 138>> d=a+bd =8 10 128 10 128 10 12>> e=a\cWarning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 1.541976e-018.e =34.0000 22.0000 62.0000-50.0000 -23.0000 -100.000028.0000 16.0000 56.0000>> f=c/bWarning: Matrix is singular to working precision.f =NaN NaN NaNNaN -Inf InfNaN NaN NaN（4）3 1.247.5 6.6 3.15.4 3.46.1A⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，求1,,A A A-'；(求矩阵的行列式)>> a=[3 1.2 4;7.5 6.6 3.1;5.4 3.4 6.1]; >> a'ans =3.0000 7.5000 5.40001.2000 6.6000 3.40004.0000 3.1000 6.1000>> inv（c）c =2.1555 0.4555 -1.6449-2.1040 -0.2393 1.5013-0.7354 -0.2698 0.7833>> det(a)ans =13.7880（5）12345678i iZi i++⎡⎤=⎢⎥++⎣⎦，输入复数矩阵；>> z=[1+2*i 3+4*i;5+6*i 7+8*i] or>>z=[1 3;5 7]+[2 4;6 8]i z =1.0000 +2.0000i3.0000 +4.0000i5.0000 +6.0000i7.0000 +8.0000i2.建立.m文件，用for循环语句生成5×5的矩阵A：125236569⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦,将A矩阵进行水平和垂直翻转得到矩阵B和C。

将A矩阵的前3行，3列变成0并赋值给D。

>> a=zeros(5);for m=1:5for n=1:5a(m,n)=m+n-1endenda =1 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 00 0 0 0 0a =1 2 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 01 2 3 4 52 3 4 5 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 0 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 0 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 0 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 0 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 7 0 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 0 0 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 0 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 0 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 8 0 0 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 7a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 85 6 0 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 85 6 7 0 0a =1 2 3 4 52 3 4 5 63 4 5 6 74 5 6 7 85 6 7 8 0a =1 2 3 4 52 3 4 5 65 6 7 8 9c=flipud(a)c =5 6 7 8 9 4 5 6 7 8 3 4 5 6 7 2 3 4 5 6 1 2 3 4 5 d=fliplr(a)d =5 4 3 2 16 5 4 3 27 6 5 4 38 7 6 5 49 8 7 6 5 >> d=a;>> d(1:3,1:3)=0d =0 0 0 4 5 0 0 0 5 6 0 0 0 6 7 4 5 6 7 8建立的m文件a=zeros(5)for m=1:5for n=1:5a(m,n)=m+n-1endendd=a;d(1:3,1:3)=0建立的m文件，存盘时，不能存成1.m。

不运行，运行就为1.我存为mydesign.m默认工作目录是work文件夹edit pathdef 改路径最后一句加上cd d:\matlabBFile---Set path3.建立.m文件，随机产生一个5×5的矩阵，元素值为从0到255，要求用0和255对该矩阵进行标记，元素值大于等于128的标记为255，元素值小于128的标记为0。

a=255*rand(5);a =193.2237 180.0418 209.9817 111.8798 124.8899189.4988 8.1174 177.1813 97.2974 113.6245100.0179 70.6154 80.8604 195.2068 164.8098167.1469 11.7737 242.3066 202.7760 180.888043.6526 24.7686 8.7838 47.6525 192.4451a=255*rand(5);for i=1:5for j=1:5if a(i,j)>128a(i,j)=255;elsea(i,j)=0 end end enda =*********** 255 0 255 0 0 *********** 255 255 255 0 255 ***********4.编写程序使用0.618搜索法（即：近似黄金分割法）求给定函数的极值。

搜索法求解min ()t t ϕ≥的基本过程如下：给出[a,b]，使得t 在[a,b]中，其中[a,b]称为搜索区间，通过迭代缩短[a,b]的长度。

当[a,b]的长度小于某个预设的值，或者导数的绝对值小于某个预设的正数，则迭代终止。

以函数3min ()21t t t t ϕ≥=-+作为处理对象，其中搜索区间定为[0,3]，精度ε定为0.1。

图1-1 程序流程图b=3;er=0.5;while(1)m=a+0.382*(b-a); n=a+0.618*(b-a); p=m^3-2*m+1;q=n^3-2*n+1;if p<=qif n-a<ermbreak; elseb=n;endelseif b-m<ernbreak; elsea=m;endendn =0.8753实验二绘图和确知信号分析实验1.假设N=12，对于M=4、5、7、10时，在0≤n≤2N -1区间使用plot 和stem 分别绘制信号)2sin(][NnM n x M π=，并添加适当标注。

clc;clear;M=[4,5,7,10]; N=12; n=0:2*N-1; for i=1:length(M) xn=sin(2*pi*M(i)*n/N); subplot(length(M),1,i) plot(n,xn); hold on ; stem(n,xn,'r'); s=int2str(M(i)) title(['M=' s]); endclc,clear; N=12; M=[4,5,7,10]; n=0:2*N-1;y1=sin(2*pi*M(1)*n/N); y2=sin(2*pi*M(2)*n/N); y3=sin(2*pi*M(3)*n/N); y4=sin(2*pi*M(4)*n/N); subplot(221)plot(n,y1);title('M=4');xlabel('n');ylabel('y1'); subplot(222)plot(n,y2);title('M=5');xlabel('n');ylabel('y2'); subplot(223)plot(n,y3);title('M=7');xlabel('n');ylabel('y3'); subplot(224)plot(n,y4);title('M=10');xlabel('n');ylabel('y4'); figure subplot(221)stem(n,y1);title('M=4');xlabel('n');ylabel('y1'); subplot(222)stem(n,y2);title('M=5');xlabel('n');ylabel('y2'); subplot(223)stem(n,y3);title('M=7');xlabel('n');ylabel('y3'); subplot(224)stem(n,y4);title('M=10');xlabel('n');ylabel('y4');2.设周期信号一个周期[0,T]的波形为()1,0/20/2t T s t T t T≤≤⎧=⎨<≤⎩，，其中T=1。