2012秋季南理工matlab综合作业(课程考试)
各城市每月国民生产总值原始数据 四、 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分)
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程序设计与⼯程应⽤》综合练习答题要求：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考试总结（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上机作业部分参考答案
-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课程设计题目全
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)
上机题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年南京理工大学数学建模竞赛论文
承诺书我们仔细阅读了全国大学生数学建模的竞赛规（）。
我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。
我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。
如有违反竞赛规则的行为，我们愿意承担由此引起的一切后果。
我们的参赛（报名）队号为： 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综合大作业(附详细答案)
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题库
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