2012秋季南理工matlab综合作业(课程考试)

合集下载

MATLAB历年考试题目(附答案)

MATLAB历年考试题目(附答案)

历年考试题目(附答案)一、填空题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

综合练习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总结+2012考题 (1)

Matlab总结+2012考题 (1)

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 ⎡⎤=⎢⎥⎣⎦。

Matlab上机作业部分参考答案

Matlab上机作业部分参考答案

-2
-1.7
-1.4
-1.1
-0.8
-0.5
-0.2
0.1
0.4
பைடு நூலகம்0.7
1
1.3
.10289 .11741 .13158 .14483 .15656 .16622 .17332 .1775
.17853 .17635 .17109 .16302
【求解】用下面的语句可以立即得出给定样本点数据的三次插值与样条 插值,得出的结果如,可见,用两种插值方法对此例得出的结果几乎一 致,效果均很理想。 >> x=[-2,-1.7,-1.4,-1.1,-0.8,-0.5,-0.2,0.1,0.4,0.7,1,1.3,... 1.6,1.9,2.2,2.5,2.8,3.1,3.4,3.7,4,4.3,4.6,4.9]; y=[0.10289,0.11741,0.13158,0.14483,0.15656,0.16622,0.17332,... 0.1775,0.17853,0.17635,0.17109,0.16302,0.15255,0.1402,... 0.12655,0.11219,0.09768,0.08353,0.07019,0.05786,0.04687,... 0.03729,0.02914,0.02236]; x0=-2:0.02:4.9; y1=interp1(x,y,x0,'cubic'); y2=interp1(x,y,x0,'spline'); plot(x0,y1,':',x0,y2,x,y,'o')
S 2i 1 2 4 8 262 263 用数值方法可以求出,
i 0
63
试不采用循环的形式求出和式的数值解。由于数值方法采用 double 形式进行计算的,难以保证有效位数字,所以结果不 一定精确。试采用符号运算的方法求该和式的精确值。 【求解】用符号运算的方式可以采用下面语句 >> sum(2.^[1:63]) ans = 1.844674407370955e+019 >> sum(sym(2).^[1:63]) ans = 18446744073709551614

matlab课程设计题目全

matlab课程设计题目全

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是一个维数为的矩阵。

MATLAB语言考试上机题A(2012)

MATLAB语言考试上机题A(2012)

上机题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。

2012年南京理工大学数学建模竞赛论文

2012年南京理工大学数学建模竞赛论文

承诺书我们仔细阅读了全国大学生数学建模的竞赛规()。

我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与本队以外的任何人(包括指导教师)研究、讨论与赛题有关的问题。

我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。

如有违反竞赛规则的行为,我们愿意承担由此引起的一切后果。

我们的参赛(报名)队号为: 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 这样,对于该树种的利用做到了可持续发展。

matlab综合大作业(附详细答案)

matlab综合大作业(附详细答案)

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 滤波器恢复出原始信号。

南理工数字信号处理matlab题库

南理工数字信号处理matlab题库

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截止频率由键盘输入,程序能根据输入的参数,绘制滤波器的增益响应。

南京理工大学课程考试试卷答案及评分标准

南京理工大学课程考试试卷答案及评分标准

_
u 注:该题为基本题,考核含理想二极管电路伏安特性曲线的绘制。 0 -2V
第 2

共 3

南京理工大学课程考试试卷答案及评分标准
课程名称:
试卷编号: 1002210466


学分:
4.5
教学大纲编号: 10022104 考试时间: 120 分钟 ┈┈ 2 分
∴ i (t ) =
30 + 6iL (t ) 54 36e 2.5t = = 6 4e 2.5t A (t ≥ 0+ ) 9 9
3Ω u1 4Ω
u1
i
_ u1 +
0.5u1 2Ω
+
0.5u1 u ┈┈ 4 分
2 U oc 3 × 3 = = 0.75W 4Req 4 × 3
i
iD
Req =
u 6 × (0.5u1 )6Ω = 1.5u1 = 3 i u1 + 0.5u1 0.5u1 _
+
当 u = uD = 0 时,D 导通, i =
南京理工大学课程考试试卷答案及评分标准
课程名称:
试卷编号: 1002210466 一、


学分:
4.5
教学大纲编号: 10022104 考试时间: 120 分钟
4. (共 7 分) S 闭合时: I A =
3∠ 30 o
380∠30 o 660∠0 o = = 6.6∠ 36.87 o A o 80 + j 60 100∠36.87
1.5u1
_
又: 3i1 + 6i2 = 3u1 + 3u1 = 6
┈┈ 4 分
U 2 = 2 I 1 + 4 I 2 + 2 I 1 + 2 I 2 = 4 I 1 + 6 I 2 (2) 1 3 由(2) : I 1 = U 2 + ( I 2 ) (3) 4 2 1 1 (3)代入(1) U 1 = U 2 + 3( I 2 ) + 2 I 2 = U 2 + ( I 2 ) : 2 2

南京理工大学课程考试试卷答案及评分标准 课程名称: 信号与系统 学分

南京理工大学课程考试试卷答案及评分标准 课程名称: 信号与系统 学分
二、(10分)已知某离散系统的单位样值响应为 ,激励信号为 ,
(1)求该系统的零状态响应 ;
(2)该系统是因果的,还是非因果的?请说明理由。
解:(1)
┈┈6分
或:
(2)该系统是因果的,因为 ┈┈4分
注:该题为基本题,考核离散系统的零状态响应的卷积求和法与系统的因果性等知识点。
三、(14分)系统如图2所示,已知x1(t)=Sa(t),x2(t)=Sa(2t),
2.系统微分方程式 ,则零输入响应分量为 ┈┈4分
注:该题为基本题,考核零输入响应等知识点。
3.已知f(t)的频带宽度为Δω,则f(2t-4)的频带宽度为2Δω
注:该题为基本题,考核傅氏变换的尺度变换性质等知识点。┈┈2分
4.已知 ,收敛域为 ,则 的z逆变换 =
注:该题为基本题,考核z变换的时移性质、逆z变换等知识点。┈┈4分
四、(14分)某因果离散时间系统由两个子系统级联而成,如图3所示,其中: ,
பைடு நூலகம்图3
(1)求该系统的系统函数H(z),并注明收敛域,说明系统是否稳定;
(2)画出该系统并联形式的框图或信号流图;
(3)求该系统的单位样值响应h[n];
解:(1) ,不稳定┈┈5分
(2)
┈┈5分
(3)
或: , ┈┈4分
注:该题为综合题,考核离散系统的系统函数及其收敛域的确定、稳定性、系统模拟、单位样值响应等知识点。
南京理工大学课程考试试卷答案及评分标准
课程名称:信号与系统学分:4.5教学大纲编号:04030803
试卷编号:001考试方式:闭卷、笔试满分分值:100分考试时间:120分钟
注意:所有答案均写在答题纸上,否则该题以零分记。
一、基本概念题(共40分)

matlab秋季上机作业

matlab秋季上机作业

matlab秋季上机作业秋季是一个美丽而多彩的季节,天空湛蓝,阳光明媚,正是适合户外活动和学习的好时光。

而在大学里,秋季也意味着新学期开始,课程开始渐渐进入正轨。

对于学习MATLAB这门课程的同学们来说,秋季上机作业是一个重要的环节。

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实验(1)

华南理工大学数字信号处理matlab实验(1)

数字信号处理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考试题题库(必做题)(带答案)

MatLab考试题题库(必做题)(带答案)

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. 当一个小圆轮沿着一条曲线行进时,轮缘任一点的轨迹就会产生变化丰富的摆线。

Matlab实验题目与答案

Matlab实验题目与答案

(温馨提示:实验课结束后,请将所有作业(题目、代码、结果)利用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。

南京理工大学课程考试试卷

南京理工大学课程考试试卷

L0 60.00 cm(不确定度可忽略); t1 10.0 0.5 oC、 t2 90.0 0.5 oC;
n1 1.00 0.05 cm, n2 3.00 0.05 cm; D 100.0 0.1cm; d 8.00 0.05cm。
已知线膨胀系数 l
n2 n1 d
2DL0 t2 t1
43.5001.0 22.00 1.000 15、已知 L L1 2L2 5L3 , L1 L1 32.4 0.1mm; L2 L2 4.580 0.004 mm;
L3 L3 5.02 0.02 mm,求 L L
20、将一块灵敏度 K H 1.50 10 3 V·A-1T-1 的霍耳片,通以工作电流 I H 2.00 mA,放在一个磁 场 B 中。当磁场方向和 I H 都为正时,测得U a 15.44 mV;当 B 和 I H 都为负时,测得 U b 3.41 mV;当 B 为正, I H 都为负时,测得U c 14.38 mV;当 B 为负, I H 为正时,测 得U d 3.52 mV。试求:
(1)、霍耳电压U H 值;
(2)、该处的磁感应强度 B。 21、 在灵敏电流计实验中
(1)、 R2 增大时,电磁阻尼是增大还是减小?
(2)、有人说短路开关 K 4 的作用是保护电流计,对吗?
(3)、测量灵敏电流计内阻 Rg 时,保持其他条件不变,当U 1.20V,R 8.0 KΩ时,电流计
偏转 100mm;当U 0.60V,若 R2 3.4 KΩ,电流计偏转 100mm,则 Rg 为多大?
12、已知 R U I ,U 3.00V, I 20 μA,则 R = (17) 。
(3)、 l l
18、用 6502 示波器观察其正弦电信号,屏幕上正好出现一个完整的波形。此时“电压灵敏度选择” 指向每格 0.5V,“电压灵敏度微调”处于“校准”位置,正弦波最高处到最低处相差 4 格。“时 基选择”指向每格 1ms,“时基调节”处于“校准”位置,试回答:

南理工matlab答案

南理工matlab答案

南理工matlab答案一、问答题(15分)1. 在安装MATLAB 软件程序时,哪些组件(模块)是必须选择的?否则,不能建立MATLAB 工作环境。

答:Toolboxes MA TLAB Complier2. 下面的5个变量名,哪些合法哪些非法,为什么?①Do_you_know_the_way_to_heaven_Jesus ②1dog ③ dog1 ④_help ⑤What’s_up 答:①非法,因为每个变量名最长只能包含19个字符;②非法,因为变量名的首字符必须是字母;③ 合法;④非法,因为变量名的首字符必须是字母;⑤非法,因为变量名中不能有标点符号。

3.分别写出下面数组的大小和内容. 注意后面的数组可能根据前面数组的定义.①. a=1:2:5;② . b=[a', a', a'];③ . c=b(1:2:3, 1:2:3);④ . d=a+b(2,:);⑤ . w=[zeros(1,3), ones(3,1)', 3:5'];⑥ . b([1, 3],2)=b([3, 1],2);答:① a =【 1 3 5】② b =【 1 1 13 3 35 5 5】③ c =【 1 15 5】④ d =【 4 6 8】⑤ w =【 0 0 0 1 1 1 3 4 5】⑥ b =【 1 5 13 3 35 1 5】4. 根据下列矩阵回答问题① A 的大小是多少?② A(4,1)的值是多少?③ A(:,1:2)的大小和值分别是多少?④ A([1 3],end)的大小和值分别是多少?1.10.02.13.5 6.00.00.6 6.6 2.84.52.1 1.20.3 1.4 1.31.45.19.0 1.30.8A --=??-??--??答:①size(A)=【4 5】②A(4,1)= -1.4③size(A(:,1:2))= 【4 2】A(:,1:2)=【 1.1000 00 0.60002.1000 1.2000-1.4000 5.1000】④size(A([1 3],end)=【2 1】A([1 3],end)=【 6.00001.3000】5. 变量a的赋值为64,分别用4个类型的数组表示为:①(1×1)的实数数组a;②(1×1)的字符串数组a;③ (1×1)的单元(cell)数组a;④ (1×1)的结构(structure)数组a.a。

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

各城市每月国民生产总值原始数据 四、 A 、 B 、 C 三个城市, 上半年每个月的国民 城市 1月 2月 3月 4月 5月 6月 生产总值(单位:亿 72 80 104 90 120 113 元)如右表。请编写 A M文件: B 60 90 115 102 110 95 (1)分别计算三个城市 C 51 40 70 85 96 110 这6个月的平均国民 生产总值。 (2)写出生成如右下图 的二维累加直方图的 M文件程序代码。 (3)计算这三个城市的 月国民生产总值最大 的城市及其最大值。
请编写实现上述运算的M程序代码,并给出每条语句的 注释说明。(10分)
三、球的位置和速度。 假设一静止小球在离地h0 的地方以初速度v0 做垂直运动,其等式为
1 2 gt + v0 t + h0 = h(t ) 2 ) gt + v0 v(t=
其中g 为重力加速度(−9.81m/s2),h(t)代表在t 时刻小球的高度,v(t)代 表在时刻t 小球的速度,设定初始高度和速度分别为h0=6m, v0=50m/s。 编写一个MATLAB程序,计算小球运动的速度和高度,并分别利用 二维曲线图和三维曲线图画出h、v 关于时间t 的函数(t从0秒至10秒, 时间间隔0.1s);并计算在t=10秒时小球的速度和高度值。 请编写上述运算的M程序代码,并给出每条语句的注释说明。(20分)
综合大作业
注意事项: • • • 共7道大题,请注明课程名称、学号和姓名 纸质作业统一使用学校信纸手写完成,标清题号,要求字迹清楚,尤其是程序代码中 的符号、变量,如出现错误、不规范或不清楚,均按0分计算。 除提交纸质作业外,必须同时通过Email提交对应程序代码的电子文档到下列邮箱地 址: Email address: matlab_zhang@; 该Email将于第4周末关闭。 电子文档为TXT格式,文件命名格式为:“matlab作业”+“姓名”+“学号”+“ .txt” ; Email邮件标题需包含姓名和学号信息,电子文档以附件形式发送,否则按0分计算。 我会对每位同学的正确提交进行回复。 必须保证纸质作业的程序代码和电子文档内容完全一致,代码可在matlab6.5版本下运 行,否则均按0分计算。 未严格按照上述要求交作业者,视为放弃该课程学分,作缺考处理。 作业提交时间:纸质作业和电子文档上交的截止日期在第4周日之前。 如有问题可第3周之前与我联系.
• • • •
一、问答题(15分) 1.写出MATLAB中常用的6个预定义变量名,并说明各自的含义。 2. 下面的5个变量名,哪些合法哪些非法,为什么? ①dog31 ②what’s_up ③ Do_you_know_the_way_to_heaven_Jesus ④5dog_HElp ⑤_sum2up 3.分别写出下面各变量的维数和内容. 注意后面的变量可能根据前面变 量的定义计算. ①. a=1:2:5; ② . b=[a', a', a']; ③ . c=a+b(2,:); ④ . b([1, 3],2)=b([3, 1],2); ⑤ . d=b(1:2:3, 1:2:3) ; ⑥ . x=[zeros(1,3), ones(3,1)', 3:5'];
二、求线性方程组Ax=B的解x,其中
3.0 4.0 −2.0 5.0 1.0 2.0 −1.0 −5.0 −2.0 6.0 −3.0 6.0 −4.0 −5.0 3.0 A = 7.0 3.0 −6.0 −5.0 −2.0 −3.0 6.0 4.0 2.0 −6.0 4.0 5.0 2.0 4.0 1.0 −1.0 4.0 −1.0 = , B −2.0 4.0 −4.0 0 1 −6 = , x 10 −5 −2 x1 x 2 x3 x4 x5 x6
六、George,Sam,Mary, Betty,Chung 和Suzie 是好 朋友,Chung要去英国留学,其他几个朋友商定去送 别他。 George,Sam,Mary, Betty和Suzie为买礼物分 别花了$16,$7,$15,$9 ,$12 。请创建一个二维饼 图,并标出每个人的名字及其花费占总数的百分比, 请问 5个人平均花费多少钱? 要求给出程序代码,以及每条语句的注释说明。(10 分)
4. 根据下列矩阵回答问题
1.1 2.1 0.0 −3.5 = A j 3.0 0.6 −6.6 2.8 −1.4 + j 5.3 5.1 −9.0 1.3
① ② ③ ④ A的维数是多少? A(3,1)的值是多少? A(:,1:2)的维数和取值分别是多少? A([1 3],end)的维数和取值分别是多少??
5. 假设变量a的赋值为108,分别用4种数据类型表示该变量为: ① (1×1)的字符 串数组a; ② (1×1)的实数数组a;③ (1×1)的单元(cell)数组a; ④ (1×1)的 结构(structure)数组x.a。 请问它们在内存空间中各占多少字节? 6. 求下列表达式的值 ①. −11/5+6 ②. (−11/5)+6 ③ . − 11/(5+6) ④. round(-11/5)+6 ⑤ . ceil(-11/5)+6 ⑥ . floor(-11/5)+6
7. 变量a,b,c ,d的定义分别如下,
1 1 1 a= 2 − 1 , b = , c = , d= eye(2) [ ] −1 2 −2
如果下面运算是合法 result=a+b; ②. result=a*d; ③. result=a.*d; ④. result=a*c; ⑤. result=a.*c; ⑥. result=a\b; ⑦. result=a.\b; ⑧. result=a .^ b; 8. M(脚本)文件和M函数文件的主要区别在哪些方面?
七、请谈谈你对这门课程的认识,要求不少于200字 (10分)。
May you all make success in your life!
要求给出每条语句的注 释说明。(15分)
五、某函数z(x, y)定义如下: 2+y2) -0.5(x z(x, y)=e 其中自变量x,y的取值区间均为[-4,4],数 据间隔均为0.2。请分别创建3个图形,画 出该函数的三维网格曲线图、三维曲面图 和等高线图,并标注每个图形的坐标轴名 称。要求给出程序代码,以及每条语句的 注释说明。 (20分)
相关文档
最新文档