下载文档原格式
历年考试题目（附答案）一、填空题1、MATLAB常用操作界面包括、工作空间窗口、、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。
2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。
3、MATLAB中逗号主要用作；用作输入量与输入量之间的分隔符；。
4、工作空间浏览器主要用于内存变量的、和。
5、MATLAB实现将全下标转换为单下标的指令为、据单下标换算出全下标的指令为。
6、二维数组的标识有、、“逻辑1”标识。
二、简述下列命令在MATLAB中的作用1.clf %清除图对象clear %清除工作空间内的所有变量clc %清除当前屏幕上显示的所有内容，但不清除工作空间中的数据2.ceil %沿+∞方向取整factor %符号计算的因式分解3.sph2cart %球坐标变为直角坐标pow2 %2的幂4.digits(d) %设置当前的可变算术精度的位数为整数 d 位vpa(pi,100) %设置当前pi的小数位数为1005.box on %打开框状坐标轴开grid off %网格关6.logspace %对数分度向量cat %串接成高维数组7.sym2poly %符号多项式转变为双精度多项式系数向量poly2sym %双精度多项式系数转变为向量符号多项式8.plot3 %三维线图poly2str %以习惯方式显示多项式9.bar %二维直方图pie %二维饼图10.zoom on %打开图形缩放模式 edit %M文件编辑11.whos %对当前工作空间变量的信息进行列表figure %生成图形窗口12.cart2sph %直角坐标变为球坐标pol2cart % 极或柱坐标变为直角坐标13.diff %数值差分、符号微分dsolve %符号计算解微分方程14.ezplot3 %画三维曲线的简捷指令fix %向零取整factor %符号计算的因式分解三简答题1.叙述MATLAB的主要窗口包括哪些？其作用是什么？命令窗口（Commend Window）——供输入命令，实现计算或绘图功能。
综合练习matlab《MATLAB程序设计与⼯程应⽤》综合练习答题要求：1.直接将答案写在本⽂档中，可通过执⾏matlab指令来验证结果；2.答案⼀律⽤蓝⾊字体标明；3.发还给我的⽂件以“学号+姓名”的⽅式命名，否则拒收。
⼀、选择题(每题2分，共40分)1.最初的MATLAB核⼼程序是采⽤（）语⾔编写的。
A. FORTRANB. CC. BASICD. PASCAL2.已于2014年3⽉发布的MATLAB新版本的编号为（）。
A.MATLAB 2014RaB.MATLAB R2014aC.MATLAB 2014RbD.MATLAB R2014b3.当在命令窗⼝执⾏命令时，如果不想⽴即在命令窗⼝中输出结果，可以在命令后加上（）。
A. 冒号B. 逗号C. 分号D. 空格4.如果要重新执⾏以前输⼊的命令，可以使⽤（）键。
A. 左箭头←B. 右箭头→C. 下箭头↓D. 上箭头↑5.下列变量名中（）是合法的。
A. end,lbcxB. char_1,i,jC. x\y,a1234D. xy,a.16.已知数组123456789101112a=??，则要以单下标⽅式访问元素8时，应该使⽤的命令是（）。
D. a(8)7.已知数组123456789a=??，则a(:,end)是指（）。
A. 第⼀⾏元素B. 所有元素C. 第三⾏元素D. 第三列元素8.已知数组123456789a=??，则运⾏mean(a)命令是（）。
A. 计算a数组每列的平均值B. 计算a数组每⾏的平均值C. a数组增加⼀⾏平均值D. a数组增加⼀列平均值9.有以下程序代码：a=eye(5);for n=a(2:end,:)则for循环的循环次数是（）。
A. 3B. 4C. 5D. 110.有以下程序代码：x=reshape(1:12,3,4);for k=1:4if x(:,k)<=6commands1;elsecommands2;则代码commands1执⾏的次数是（）。
Matlab考试总结（100分）之前拷课件的时候，老师给了一个cpt00的文件,里面有9道题，比较简单，感觉像是考题，我就把他全部弄出了，并添加2012年的部分有参考价值的题目和练习题，把他们整理在一起，希望对大家有帮助。
由于时间比较紧,可能会有问题,大家可以跟我联系QQ:181816322。
一、矩阵基本运算（约10分）（1）求矩阵a 的逆阵、对角线及其行列式的值。
a =2 1 -3 -13 1 0 7-1 2 4 -21 0 -1 5a=[2 1 -3 -1;3 1 0 7;-1 2 4 -2;1 0 -1 5]；inv(a)det(a)（2）.产生一个5阶的随机矩阵,大小位于(-5,5),并求出方阵的逆、行列式的值、特征值与特征矢量及正交矩阵。
（8分，2012年考试题）情况一:产生的矩阵元素全为整数A= randint(5,5,[-5,5])情况二:产生的矩阵元素为小数A=-5+10*rand(5)矩阵的逆：B=inv(A)或者B=A^(-1)行列式的值：D=det(A)特征值和特征矢量：[m n]=eig(A)正交矩阵：E=orth(A)扩展：矩阵的伪逆：pinv(A)矩阵的秩rank(A)条件数：cond(A)LU分解：[l,u]=lu(A)QR分解：[q, r]=qr(A)1、2、无穷范数：norm(A,1) norm(A,2) norm(A,inf)产生5阶的魔方（魔术）矩阵：magic(5)2.求垂直于向量A=（1，2，3）和B=（3，4，5）的向量，并计算三个向量的混合积（平行六面体的体积）。
A=[1 2 3]；B=[3，4，5]；C=cross(A,B)S=dot(A,cross(B,C))二、求导，积分，微分，极限15分1．求下列函数的极限(作业题)（1）lim(x 2/sin 2(x/3)) x=0（2）lim((tanx-sinx)/sin 3x) x=0程序代码如下：syms x;a= limit(x^2/(sin(x/3))^2,0)b= limit((tan(x)-sin(x))/(sin(x))^3,0)（3）32lim(1)x x t x→∞+（3分） （2012年考题） syms t xlimit((1+2*t/x)^(3*x),x,inf)2．求下列函数的定积分（1） (x+sinx)/(1+cosx) [0,pi/2]（2） cos 5xsinx [0,pi/2]syms x;int((x+sin(x))/(1+cos(x)),x,0,pi/2)int(cos(x)^5*sin(x), x,0,pi/2)3．求导数(2012年考题)2222,,.(4)cos ln a t df d f d f f dx dt dxdt t x x ⎡⎤=⎢⎥⎣⎦。
Matalab课后作业学院：电气信息工程及其自动化班级：学号：姓名：完成日期： 2012年12月23日1、 matlab 软件主要功能是什么？电气工程及其自动化专业本科生主要用到哪些工具箱，各有什么功能？答：（1）主要功能：工业研究与开发； 数学教学，特别是线性代数；数值分析和科学计算方面的教学与研究；电子学、控制理论和物理学等工程和科学学科方面的教学与研究； 经济学、化学和生物学等计算问题的所有其他领域中的教学与研究；符号计算功能；优化工具；数据分析和可视化功能；“活”笔记本功能；工具箱；非线性动态系统建模和仿真功能。
（2）常用工具箱：（a ） MATLAB 主工具箱：扩充matlab 的数值计算、符号运算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能。
（b ）符号数学工具箱：符号表达式、符号矩阵的创建；符号可变精度求解；因式分解、展开和简化；符号代数方程求解；符号微积分；符号微分方程。
（c ） SIMULINK 仿真工具箱： Simulink 是用于动态系统和嵌入式系统的多领域仿真和基于模型的设计工具。
对各种时变系统，包括通讯、控制、信号处理、视频处理和图像处理系统，Simulink 提供了交互式图形化环境和可定制模块库来对其进行设计、仿真、执行和测试。
（d ）信号处理工具箱：数字和模拟滤波器设计、应用及仿真；谱分析和估计；FFT 、DCT 等变换；参数化模型。
（e ）控制系统工具箱：连续系统设计和离散系统设计；状态空间和传递函数以及模型转换；时域响应（脉冲响应、阶跃响应、斜坡响应）；频域响应（Bode 图、Nyquist 图）；根轨迹、极点配置。
2、设y=23e t 4-sin(43t+3),要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
程序如下：clcclearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3);y1=diff(y);subplot(2,1,1)plot(x,y)subplot(2,1,2)plot(x(1:150),y1)曲线如下图所示：3、设A是一个维数为的矩阵。
上机题A在Editor中按如下格式编程：%%%%%%%%%第1题%%%%%%%%%%%%%%%%%clc;clear;close all;%编写第1题程序%%%%%%%%%第2题%%%%%%%%%%%%%%%%%clc;clear;close all;%编写第2题程序所有题编完后用设置断点的方式逐题演示程序结果，时间30分钟。
1.(30分)在Matlab的work目录下新建文件夹MYPROG(2分)，并将其添加到搜索路径上（4分）,将当前工作目录切换到MYPROG下(2分)。
输入A为3×3的单位阵（4分），B为3×3的全1阵（4分），C为A加B的和（4分），D为A 点乘B（4分）。
存储C、D到数据文件data1.dat并放到MYPROG目录下（6分）。
2. (30分)t∈[0,1]，步长0.01，求t序列的维数(4分)，y=sin(2πt) (4分)，以t为横坐标，y为纵坐标，用黑色虚线画出图形(10分)，并标注x轴为t(4分)，标注y 轴为y(4分)，标题为正弦信号(4分)。
3.(40分)产生服从高斯分布方差为1的3×4×2多维阵列E(5分)，F=sum(E) (5分)，G=’I have learned how to use matlab.’ (5分)，H是一结构阵列(10分)，其形式如下Patient(1).Name.Bill ‘John’127.Name .Bill ‘Ann’109Patient(2)I是2×2单元阵列，其四个元素分别为E，F，G，H，并用cellplot显示单元阵列内容(15分)。
机题A 答案1. 用path 命令查看是否存在MYPROG 文件夹。
2. 1 10100.20.40.60.81-1-0.50.51ty正弦信号3.C =2 1 1 1 2 1 1 1 2 D =1 0 0 0 1 0 0 0 1A 1 0 0 0 1 0 0 0 1B = 1 1 1 1 1 1 1 1 1。
承诺书我们仔细阅读了全国大学生数学建模的竞赛规（）。
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们愿意承担由此引起的一切后果。
我们的参赛（报名）队号为： 6我们选择的题号为（A或B）：B参赛组别（研究生或本科）：本科参赛队员 (先打印,后签名,并留联系电话) ：可持续利用森林资源的策略摘要森林资源是人类赖以生存的重要资源，保护森林资源并对其进行可持续性利用是当代科学的重要课题，本文研究了森林的可持续利用问题，并对于森林砍伐、种植策略等问题进行了较深入的分析与讨论。
为了解决问题，我们首先建立了描述一棵树生长过程的含材体积模型。
通过查找资料确定树木生长速度的二次曲线模型bt at dtdv+-=2，一棵树木的含材体积⎰+-==232131bt at dt dt dv v (假设t =0时，含材体积为0)通过观察树木含材体积的曲线确立成材年限m 与成长停滞年限max 的关系。
a b m 43=，ab=max 接着，我们建立了单维离散动态模型，先考虑稳定状态(即时间充分长)下森林中只有一种树的各年龄段树数量变化情况，引入砍伐强度变量，认为木材的需求按年计算，即树木每年按需求砍伐一次。
我们认为种植策略是将砍伐过后的空地上种上幼苗，并及时将死掉的幼苗清除并重栽幼苗。
为了保持系统稳定，这里认为树木自我繁殖率与死亡率大致相等。
随后给出了该树种各年龄段的动态差分方程组并构造Leslie 矩阵：)()1(t X L t X *=+，计算矩阵最大特征根为1，表明该树种各年龄段的在确定的砍伐强度下分布情况),...,,...,,(max 21*X X X X X m =最终趋于稳定（即每年砍伐前*X 稳定）以3.0,2.01==k k 为例（种植面积100万公顷） 树木年龄（年）12345678稳定值（810⨯棵）2.460 2.463 2.462 2.458 2.454 2.452 1.473 1.773 这样，对于该树种的利用做到了可持续发展。
m a t l a b综合大作业(附详细答案)-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）：1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5)实验结果：A =0.1349 3.3818 0.6266 1.2279 1.5888-2.3312 3.3783 2.4516 3.1335 -1.67241.2507 0.9247 -0.1766 1.11862.42861.5754 1.6546 5.3664 0.8087 4.2471-1.2929 1.3493 0.7272 -0.6647 -0.38362)将矩阵A按列拉长得到矩阵B；实验程序：B=A(:)实验结果：B =0.1349-2.33121.25071.5754-1.29293.38183.37830.92471.65461.34930.62662.4516-0.17665.36640.72721.22793.13351.11860.8087-0.66471.5888-1.67242.42864.2471-0.38363)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C；实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)]实验结果：C =3.3783 3.13350.9247 1.11864)寻找矩阵A中大于0的元素；]实验程序：G=A(find(A>0))实验结果：G =0.13491.25071.57543.38183.37830.92471.65461.34930.62662.45165.36640.72721.22793.13351.11860.80871.58882.42864.24715)求矩阵A的转置矩阵D；实验程序：D=A'实验结果：D =0.1349 -2.3312 1.2507 1.5754 -1.29293.3818 3.3783 0.9247 1.6546 1.34930.6266 2.4516 -0.1766 5.3664 0.72721.2279 3.1335 1.1186 0.8087 -0.66471.5888 -1.67242.4286 4.2471 -0.38366)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；实验程序：E=flipud(fliplr(A))实验结果：E =-0.3836 -0.6647 0.7272 1.3493 -1.29294.2471 0.80875.3664 1.6546 1.57542.4286 1.1186 -0.1766 0.9247 1.2507-1.6724 3.1335 2.4516 3.3783 -2.33121.5888 1.2279 0.6266 3.3818 0.13497)删除矩阵A的第2列和第4列得到矩阵F；实验程序：F=A;F(:,[2,4])=[]实验结果：F =0.1349 0.6266 1.5888-2.3312 2.4516 -1.67241.2507 -0.17662.42861.5754 5.3664 4.2471-1.2929 0.7272 -0.38368)求矩阵A的特征值和特征向量；实验程序：[Av,Ad]=eig(A)实验结果：特征向量Av =-0.4777 0.1090 + 0.3829i 0.1090 - 0.3829i -0.7900 -0.2579 -0.5651 -0.5944 -0.5944 -0.3439 -0.1272-0.2862 0.2779 + 0.0196i 0.2779 - 0.0196i -0.0612 -0.5682 -0.6087 0.5042 - 0.2283i 0.5042 + 0.2283i 0.0343 0.6786 0.0080 -0.1028 + 0.3059i -0.1028 - 0.3059i 0.5026 0.3660 特征值Ad =6.0481 0 0 0 00 -0.2877 + 3.4850i 0 0 00 0 -0.2877 - 3.4850i 0 00 0 0 0.5915 00 0 0 0 -2.30249)求矩阵A的每一列的和值；实验程序：lieSUM=sum(A)实验结果：lieSUM =-0.6632 10.6888 8.9951 5.6240 6.208710)求矩阵A的每一列的平均值；实验程序：average=mean(A)实验结果：average =-0.1326 2.1378 1.7990 1.1248 1.24172．符号计算（10分，每小题5分）：1)求方程组20,0++=++=关于,y z的解；uy vz w y z w实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z');y= S. y, z=S. z实验结果：y =[ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w] [ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w] z =[ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))] [ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))]2)利用dsolve 求解偏微分方程,dx dyy x dt dt==-的解； 实验程序：[x,y]=dsolve('Dx=y','Dy=-x')实验结果：x =-C1*cos(t)+C2*sin(t)y = C1*sin(t)+C2*cos(t)3．数据和函数的可视化（20分，每小题5分）：1)二维图形绘制：绘制方程2222125x y a a +=-表示的一组椭圆，其中0.5:0.5:4.5a =；实验程序：t=0:0.01*pi:2*pi; for a=0.5:0.5:4.5; x=a*cos(t); y=sqrt(25-a^2)*sin(t); plot(x,y) hold on end实验结果：2) 利用plotyy 指令在同一张图上绘制sin y x =和10x y =在[0,4]x ∈上的曲线；实验程序：x=0:0.1:4; y1=sin(x); y2=10.^x;[ax,h1,h2]=plotyy(x,y1,x,y2); set(h1,'LineStyle','.','color','r'); set(h2,'LineStyle','-','color','g'); legend([h1,h2],{'y=sinx';'y=10^x'});实验结果：3)用曲面图表示函数22z x y =+；实验程序：x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; surf(X,Y,Z)实验结果：4)用stem 函数绘制对函数cos 4y t π=的采样序列；实验程序：t=-8:0.1:8;y=cos(pi.*t/4); stem(y)实验结果：4. 设采样频率为Fs = 1000 Hz ，已知原始信号为)150π2sin(2)80π2sin(t t x ⨯+⨯=，由于某一原因，原始信号被白噪声污染，实际获得的信号为))((ˆt size randn x x+=，要求设计出一个FIR 滤波器恢复出原始信号。
DSP题库详解1已知3阶椭圆IIR数字低通滤波器的性能指标为：通带截止频率0.4π，通带波纹为0.6dB，最小阻带衰减为32dB。
设计一个6阶全通滤波器对其通带的群延时进行均衡。
绘制低通滤波器和级联滤波器的群延时。
%Progranm 1% Group-delay equalization of an IIR filter.%[n,d] = ellip(3,0.6,32,0.4);[GdH,w] = grpdelay(n,d,512);plot(w/pi,GdH); gridxlabel('\omega/\pi'); ylabel('Group delay, samples');title('Original Filter');F = 0:0.001:0.4;g = grpdelay(n,d,F,2); % Equalize the passbandGd = max(g)-g;% Design the allpass delay equalizer[num,den,tau] = iirgrpdelay(6, F, [0 0.4], Gd); %设计六阶的全通滤波器[b,a]=iirgrpdelay(6,F,[0 0.4],Gd);He1=dfilt.df2(b,a);He=dfilt.df2(n,d);He_all=dfilt.cascade(He,He1);grpdelay(He_all) %DFILT:Digital Filter Implementation.%[GdA,w] = grpdelay(num,den,512);%figure(2);%plot(w/pi,GdH+GdA); grid%xlabel('\omega/\pi');ylabel('Group delay, samples');%title('Group Delay Equalized Filter');2设计巴特沃兹模拟低通滤波器，其滤波器的阶数和3-dB截止频率由键盘输入，程序能根据输入的参数，绘制滤波器的增益响应。
matlab秋季上机作业秋季是一个美丽而多彩的季节，天空湛蓝，阳光明媚，正是适合户外活动和学习的好时光。
而在大学里，秋季也意味着新学期开始，课程开始渐渐进入正轨。
对于学习MATLAB这门课程的同学们来说，秋季上机作业是一个重要的环节。
MATLAB是一种高级技术计算语言和环境，广泛应用于各个学科的科学计算、数据分析和可视化等领域。
秋季上机作业就是要求同学们运用MATLAB进行一系列的计算和实验，以增加对该语言的熟悉和应用能力。
在MATLAB秋季上机作业中，同学们通常会遇到各种各样的问题和挑战。
比如，如何使用MATLAB进行矩阵运算、如何画出漂亮的图形、如何编写函数等等。
这些问题都需要同学们充分理解课程内容，灵活运用MATLAB的各种功能和工具。
在完成秋季上机作业时，同学们可以按照以下步骤进行操作。
首先，要仔细阅读作业要求，明确任务的要求和目标。
然后，根据要求编写MATLAB代码，实现相应的功能。
在编写代码的过程中，同学们要注意代码的规范和可读性，命名要具有描述性，注释要清晰明了。
接着，运行代码并进行调试，确保程序的正确性。
最后，根据作业要求整理实验报告，详细记录实验过程和结果，并对实验结果进行分析和总结。
在MATLAB秋季上机作业中，同学们可以学到很多有用的知识和技能。
首先，通过编写代码和运行程序，同学们可以加深对MATLAB语言和环境的理解和掌握。
其次，通过解决问题和完成作业任务，同学们可以提高自己的逻辑思维和问题解决能力。
此外，通过实验报告的撰写和总结，同学们还可以提高自己的表达和写作能力。
当然，在完成秋季上机作业时，同学们也会遇到一些困难和挫折。
比如，代码运行出错、结果不符合预期等等。
这时，同学们要保持耐心和冷静，仔细检查代码和参数设置，找出问题所在，并进行适当的修改和调试。
如果遇到无法解决的问题，可以向同学、助教或老师寻求帮助，共同探讨和解决。
总的来说，MATLAB秋季上机作业是一个很好的学习机会，通过实践和实验，同学们可以更好地理解和掌握MATLAB语言和工具，提高自己的计算和分析能力。
电子信息科学与技术09级Matlab 期末考试试卷(2012年5月)院系： 物理与电子 年级： 电信09级 班级： 3班 学号： 1404090311 姓名： 徐云霄说明：请保留题目，在每个题目解答部分的空白处依次作答，并写清楚每个小题的题号。
作答要给出程序代码、仿真结果。
1、(30分) 如下图所示的弹簧振子系统在光滑的水平面运动，已知弹簧振子的质量kg m 2.0=，弹簧的倔强系数m N K /2.0=。
在0=t 时弹簧振子位于m x 1)0(=，速度为s m v /2)0(-=。
假设该系统在0≥t 时受到沿正方向的策动力t e F 2.01.0-=的作用，并且弹簧所受空气等阻力为cv f -=阻，其中m Ns c /1.0=，v 为弹簧振子的运动速度。
因此，该运动方程可以表示为：t e kx x c x m 2.01.0-+--=。
(1) 利用Dsolve 函数，编写程序代码求解弹簧振子运动方程的通解，并仿真s t ]20,0[∈弹簧振子的位置图像。
(5分)(2) 利用常微分方程的数值解法编写程序代码，利用子图绘制函数subplot(3,2,n)在同一窗口分别绘制s t ]20,0[∈弹簧振子：(10分)(a) 位置图像（要求标注：线的颜色（红色）、线型（-.点划线）、线宽（2）、数据点标注形状（* ））；(b) 速度图像；(c) 加速度图像；(d) 弹簧振子动能221mv E k =图像；(e) 弹簧振子势能221kx E p =图像；(f) 弹簧振子能量图像。
以上绘图分别要求标注对应的标题（如：位置图像、速度图像……）。
(3) 构建simulink 仿真模块，在同一个窗口利用三个坐标系按顺序分别仿真s t ]20,0[∈弹簧振子的加速度图像、速度图像、位置图像。
(5分)(4) 利用s-function ，构建simulink 仿真模块，在同一个窗口利用三个坐标系按顺序分别仿真s t ]20,0[∈弹簧振子的位置图像、速度图像、加速度图像。
数字信号处理matlab实验学院：电子与信息专业班级： 5学号：201130302068姓名：梁敏健实验名称：离散时间信号与系统分析、滤波器实验日期：2013.10实验1 离散时间信号与系统分析（2课时）题目一：以下程序中分别使用conv和filter函数计算h[n] = [3 2 1 -2 1 0 -4 0 3]和x[n] = [1 -2 3 -4 3 2 1]的卷积y和y1，并分析y和y1是否有差别代码：h = [3 2 1 -2 1 0 -4 0 3];x = [1 -2 3 -4 3 2 1];y = conv(h,x);n = 0:14;subplot(2,1,1);stem(n,y);xlabel('Time index n'); ylabel('Amplitude');title('Output Obtained by Convolution'); grid;x1 = [x zeros(1,8)];y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('Time index n'); ylabel('Amplitude');title('Output Generated by Filtering'); grid;可见，y与y1没有差别。
题目二：编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。
要求分别用filter、conv、impz三种函数完成。
,，源代码：对,，单位冲激响应a1=[1,0.75,0.125];b1=[1,-1];n=0:20;x1=[1 zeros(1,20)];y1filter=filter(b1,a1,x1);subplot(2,2,1);stem(n,y1filter);title('y1filter');xlabel('x');ylabel('y');title('impulse response generated by filter')x1_1=[1 zeros(1,10)];[h]=impz(b1,a1,10);y1conv=conv(h,x1_1);n=0:19;subplot(2,2,2);stem(y1conv) ;title('impulse response generated by conv')h=impz(b1,a1,21);subplot(2,2,3);stem([0:20], h);title('impulse response generated by impz') 结果：单位阶跃：a1=[1,0.75,0.125];b1=[1,-1];n=0:20;x2=ones(1,21);y1filter=filter(b1,a1,x2);subplot(2,2,3);stem(n,y1filter);title('y1filter_step');xlabel('x');ylabel('y');title('impulse response generated by filter')x2=ones(1,21);[h]=impz(b1,a1,20);y1=conv(h,x2);y1conv=y1(1:21);n1=0:20;subplot(2,2,2);stem(n1,y1conv,'filled');title('y1conv');xlabel('n');ylabel('y1[n]');title('impulse response generated by conv')a=[1,0.75,0.125];b=1;h=impz(b,a,21);subplot(2,2,1);stem([0:20], h);title('impulse response generated by impz') 结果：对单位冲激响应：a2=1;b2=[0 0.25*ones(1,4)];n=0:9;x1=[1 zeros(1,9)];y2filter=filter(b2,a2,x1);subplot(2,2,1);stem(n,y2filter);title('y2filter');xlabel('x');ylabel('y')title('impulse response generated by filter')x2=[1 zeros(1,5)];[h]=impz(b2,a2,5);y2conv=conv(h,x2);subplot(2,2,2);stem(n,y2conv,'filled')title('impulse response generated by conv')h=impz(b2,a2,10);subplot(2,2,3);stem([0:9], h);title('impulse response generated by impz') 结果：单位阶跃a1=1;b1=[0 0.25*ones(1,4)];n=0:20;x1=ones(1,21);y2filter=filter(b1,a1,x1);subplot(2,2,1);stem(n,y2filter);title('y2filter_step');xlabel('x');ylabel('y');title('impulse response generated by filter')y2=conv(b1,x1);y2conv=y2(1:21);subplot(2,2,2);stem(n,y2conv,'filled');title('y2conv');xlabel('n');ylabel('y[n]');title('impulse response generated by conv')n=0:20;b=[0,0.25,0.5,0.75,ones(1,17)];a=1;h=impz(b,a,21);subplot(2,2,3);stem([0:20], h);title('impulse response generated by impz')题目三2N点实数序列N=64。
MatLab考试题题库（必做题）（带答案）一，1.请登陆美国 MathWorks 公司的网站 ()，查看看现在大概有多少本 MATLAB-based books （以 MATLAB 为基本软件，来说明各个专业领域的教科书或工具书）。
哪一个领域的 MATLAB-based books 最多？中文书共有几本？答：1612本，数学方面的最多，中文书共有37本。
2.请在 MATLAB中直接输入下列常数，看它们的值是多少：a.ib.jc.epsd.infe.nanf.pig.realmaxh.realmin依次解为：ans = 0 + 1.0000i ans = 0 + 1.0000i ans =2.2204e-016 ans =Inf ans = NaN ans =3.1416 ans =1.7977e+308 ans =2.2251e-3083.试写一函数 regPolygon(n)，其功能为画出一个圆心在 (0, 0)、半径为 1 的圆，并在圆内画出一个内接正 n 边形，其中一顶点位于 (0, 1)。
例如 regPolygon(8) 可以画出如下之正八边型：解：新建regPolygon.m文件如下：function y=regPolyfon(n)n=8;%要画的n边形R=1; %圆的半径t=0:0.01:2*pi;x=R*cos(t);y=R*sin(t);m=linspace(pi/2,5/2*pi,n+1);xz=R*cos(m);yz=R*sin(m);hold onplot(x,y,xz,yz);axis 'equal';4.一条参数式的曲线可由下列方程式表示：x = sin(t), y = 1 - cos(t) + t/10当 t 由 0 变化到 4*pi 时，请写一个 MATLAB 的脚本 plotParam.m，画出此曲线在 XY 平面的轨迹。
解：新建plotParam.m ：t = linspace(0, 4*pi);x = sin(t);y = 1-cos(t)+t/10;plot(x, y, '-o'); -1-0.8-0.6-0.4-0.200.20.40.60.8100.511.522.535. 当一个小圆轮沿着一条曲线行进时，轮缘任一点的轨迹就会产生变化丰富的摆线。
（温馨提示：实验课结束后，请将所有作业（题目、代码、结果）利用word 整理成一个完整的实验报告，加上封面，打印，纸质档于18周周一交）第一次上机作业目的：1. 掌握MATLAB 各种表达式的书写规则2. 运行课堂上讲过的例子，熟悉矩阵、表达式的基本操作和运算。
作业：1. 熟悉matlab 集成环境界面。
回答以下问题，并操作相关的指令：（1） 分别写出清除命令窗口和清除变量的指令。
答： clc 和clear（2）在命令行输入命令后，matlab 的搜索过程是怎样的？答： （1）检查该命令对象是不是一个变量。
（2）检查该命令对象是不是一个内部函数。
（3）检查该命令对象是否为当前目录下的程序文件。
（4）检查该命令对象是否为MATLAB 搜索路径中其他目录下的M 文件。
（3）什么是matlab 的当前工作目录？写出两种设置当前工作目录的方法？ 答： 就是matlab 当前文件读取和存储的默认路径（1）在当前目录窗口中更改（2）在MATLAB 桌面工具栏中更改（3）使用cd 命令:cd c:\mydir---将c ：\mydir 设置为当前目录（4）什么是matlab 的搜索路径？写出两种设置搜索路径的方法？答： 指Matlab 运行文件时进行搜索的目录。
（1）用path 命令设置：（2）用Set Path 对话框设置（5）help 命令和doc 命令有什么作用，它们有什么区别？答： help 命令：最基本的帮助命令，查询信息直接显示在命令窗口。
doc 命令：在帮助窗口中显示HTML 帮助文档，显示函数的详细用法及 例子，比help 命令更详细。
2. 在matlab 中输入下列表达式，并求各表达式的值，显示MATLAB 工作空间的使用情况并用两种方式保存全部变量，变量保存的文件名必须包含自己的学号后四位数：（1）)1034245.01(26-⨯+⨯=w w=sqrt(2)*(1+0.34245*10^-6)w =1.4142（2）,)tan(22ac b e abc c b a x ++-+++=ππ 其中a=3.5,b=5,c=9.8。
