MATLAB上机模拟试题
Matlab上机模拟试题
Matlab 上机模拟试题1、 计算当1,2,5x =-时，21()21x f x x e x -=+--的值。
2、计算0tan 60+3、 求解矩阵方程：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡315432343122321X 4、 求解矩阵方程：⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡212101343122321X 5、 计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。
6、 用MATLAB 计算差分方程()0.7(1)0.45(2)0.6(3)0.8()0.44(1)0.36(2)0.02(3)y n y n y n y n x n x n x n x x +-----=--+-+-当输入序列为()()x n n δ=时的输出结果()y n ,n 0,1,...40 =。
7、用MA TLAB 计算差分方程()0.7(1)0.45(2)0.6(3)0.8()0.44(1)0.36(2)0.02(3)y n y n y n y n x n x n x n x x +-----=--+-+-所对应的系统函数的频率响应。
8、用FFT 计算序列(){13112331}x n =-的频谱。
9、求下列直接型系统函数的零、极点。
1234123410.10.30.30.2()10.10.20.20.5zz z z H Z z z z z ------------=++++10、画出序列()()()()()() x n n 2n 13n 24n 35n 4δδδδδ=+-+-+-+- 的图形并补零，序列长度为20.11、产生一个1024点的正弦加白噪声序列，并绘其频谱。
12、已知信号()2(cos 2)cos100f x x x = ，求其FFT 变换的幅值和相位值，并绘出结果。
13、给定()()()()()()x n n 2n 13n 24n 35n 4δδδδδ=+-+-+-+-， ()()()()()y n n 2n 1n 22n 3δδδδ=+-+-+-。
matlab试题及答案
matlab试题及答案Matlab 上机考试试题考试要求:1、从10道题目中随机抽取3道独立完成，时间1小时。
(输入randperm(10)，取前三个数)2、每个题目兴建一个.m的文件，命名方式ks+N.m(N为题号)，然后将所选三个题目放入一个文件夹，文件名为学号+姓名。
考试完成后将文件夹通过FTP提交。
3、考试完成后要写一份报告，内容包括以下:(建一个.Doc的文档,文件名为学号+姓名)(1) 题号，题目;(2) 运行结果及其分析;(3) 图也要粘贴在文档中。
4、查阅资料写一篇2000字左右的关于matlab在电子信息中的应用的小论文或综述，也可以具体的写matlab在电子信息中某一个方面或某一个点的应用。
(打印或手写都可，打印版要交电子文档)5、所有要交的东西在1月3号之前必须交齐。
(由学习委员统一收齐交给我，电子文档也拷到学习委员处，统一拷给我)。
所交项目包括:考试报告打印版，小论文打印版(两个装订在一起，考试报告在上，小论文在下，最好做一个统一的封皮)，考试报告doc文档，小论文doc文档。
Matlab 上机考试试题1.求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+ 2w=-3X +8z- 5w=9-6x+5y-2z+10w=-8(1)求系数矩阵的秩; (2)求出方程组的解。
22sinx，yz,2.在[-10，10;-10，10]范围内画出函数的三维图形。
22x，y1Hs(),3.试画出系统的零极点分布图，判断系统是否稳定，同时求其单位32sss，，，221冲激响应和频率响应(幅频特性和相频特性)。
4. 将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。
,,cos(x)，x,[,，]22(1)(曲线图);22xy(2)(曲面图)。
f(x,y),，;(-2,x,2,-4,y,4)2224,11，z5.系统传递函数为，按照以下要求求解: Hz(),,,1210.81,，zz(1)求其极零点图，判断系统的稳定性，画出系统的频谱特性;050,,n(2)当系统输入信号为:，时，画出xnnn()[5cos(0.2)2sin(0.7)],，，,, 系统的输出。
MATLAB上机模拟考试题
M A T L A B上机模拟考试题(一)(总8页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--(1)在MATLAB的命令窗口中执行_____命令，将命令窗口的显示内容清空。
（）off(2)在MATLAB的命令窗口中执行_____命令，使数据输出显示为十六进制表示。
（） long rat hex short e(3)下列变量名中_____是合法的。
（） *y,a,1 \y,a1234 ,1 bcx ,i,j(4)已知x=0:5,则x有_____个元素。
（）(5)一下运算符中哪个的优先级最高_____。
（） A./ B.^ C.~=D.&(6)使用检测函数isnumeric(10)的结果是_____。
（）(7)三维图形中默认视角是_____度。
（） A.方位角=0 俯仰角=90B.方位角=90 俯仰角=0C.方位角= 仰俯角=30D.方位角=0 仰俯角=180(8)将符号表达式化简为因式分解因式分解因式分解因式分解形式，使用_____函数。
（）(9)运行以下命令，则_____描述是正确的。
（）>>syms a b cd >>A=[a b;c d] 占用的内存小于100B B.创建了5个符号变量占用的内存是a b c d的总和 D.不存在（10）已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。
（11）运行命令bitor（8,7）的结果是_____。
（12）运行以下命令：>>x=0:10; >>y1=sin(x); >>y2=5*sin(x);>>y3=[10*sin(x);20*sin(x)]; >>plot(x,y1,x,y2,x,y3) 则在一个图形窗口中，可以看到_____条曲线。
（13）符号表达式“g=sym(sin(a*z)+cos(w*v))”中的自由符号变量是_____。
MATLAB上机实验练习题及答案
MATLAB上机实验练习题及答案09级MATLAB上机实验练习题1、给出一个系数矩阵A[2 3 4;5 4 1;1 3 2]，U=[1 2 3]，求出线性方程组的一个精确解。
2、给出两组数据x=[0 0.3 0.8 1.1 1.6 2.3]’y=[0.82 0.72 0.63 0.60 0.55 0.50]’，我们可以简单的认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。
3、解线性方程4、通过测量得到一组数据：5、已知一组测量值6、从某一个过程中通过测量得到：分别采用多项式和指数函数进行曲线拟合。
7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近插值法，双线性插值法，二重三次方插值法）8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。
分别使用不同的颜色，线形和标识符。
9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额占总销售额的比例，并分离第三种产品的切片。
X= 19.3 22.1 51.634.2 70.3 82.4 61.4 82.9 90.8 50.5 54.9 59.1 29.4 36.3 47.010、对应时间矢量t ，测得一组矢量y采用一个带有线性参数的指数函数进行拟合，y=a 0+a 1e -t +a 2te -t ，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。
11、请创建如图所示的结构数组（9分）12、创建如图所示的元胞数组。
（9分）13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。
要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。
14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。
MATLAB上机练习题(一)
MATLAB 上机练习题（一）
1. 计算下式
7
2009/42008/41949222?+-⨯-⨯=
2. 试计算d
3461.1,30,5b d a b c
a b c +=-+===其中
要求只输出d 的值，其他变量的值不输出。
**************************************************************
3. 计算下式：
181(2)sin()cos()6,3
y
x x y π-+-==（，其中 0.1104
2log 805(3)2009
6(4)1,5;
2e
e x x x -⨯+-=其中
4. 完成以上题目后清除工作空间中的所有变量。
**************************************************************
5. 求下列方程的根：
（1）325x x -=
（2）32340x x x -+=
6. 根据课本任意找出4个矩阵在MA TLAB 中输入和输出。
附1：增加的一些常用函数pi：常数
round(x)：四舍五入至最近整数
fix(x) ：无论正负，舍去小数至最近整数min(x): 向量x的元素的最小值
max(x): 向量x的元素的最大值
diff(x): 向量x的相邻元素的差
length(x): 向量x的元素个数
sin(x)：正弦函数
cos(x)：余弦函数
tan(x)：正切函数
asin(x)：反正弦函数
acos(x)：反余弦函数
atan(x)：反正切函数。
matlab上机练习（附答案）
1.以下两种说法对吗？（1）MATLAB进行数值的表达精度与其指令窗中的数据显示精度相同。
（2）MATLAB指令窗中显示的数据有效位数不超过七位。
2.历史指令窗所记录的内容与diary指令所产生的“日志”内容有什么不同？DIARY filename causes a copy of all subsequent command window inputand most of the resulting command window output to be appended to thenamed file. If no file is specified, the file 'diary' is used.DIARY OFF suspends it.DIARY ON turns it back on.DIARY, by itself, toggles the diary state.Use the functional form of DIARY, such as DIARY('file'),when the file name is stored in a string.3.如何把用户自己的“工作目录”永久地设置在MATLAB的搜索路径上？“位于搜索路径上的目录”与“当前目录”在MATLAB中的功用相同吗？4.如何向MATLAB工作空间输入一个含有100个左右元素的一维或二维数值数组？用直接键入法？用数组编辑器？用M文件编辑器？5.运用数组算术运算符去掉下面程序里的for/end循环：x=11:15for k=1:length(x)z(k)=x(k)^2+2.3*x(k)^0.5;endx=11:15 z1=x.^2+2.3*x.^0.56.不使用数组算术运算符，重写下面的程序代码：x=[2 1 4]z=1./(1+x.^2)x=2;k=1;while i<=4,z2(k)=1/(1+i^2);i=i+1;x=x+1;end7．某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
MATLAB期末上机试题带答案
MATLAB 期末上机考试试题带答案版姓名：学号：成绩：1.请实现下图：50100150200250x y x=linspace(0,8*pi,250);y=sin(x);plot(x,y)area(y,-1)xlabel('x')ylabel('y')title('y=sin(x)')2.请实现下图：x=linspace(0,2*pi,100);y1=sin(x);subplot(2,2,1)plot(x,y1,'k--')grid onxlabel('x')ylabel('y')title('sin(x)')legend('y=sin(x)')y2=cos(x);subplot(2,2,2)plot(x,y2,'r--')grid onxlabel('x')ylabel('y')title('cos(x)')legend('y=cos(x)')y3=tan(x);subplot(2,2,3)plot(x,y3,'k-')grid onxlabel('x')ylabel('y')title('tan(x)')legend('y=tan(x)')y4=cot(x);subplot(2,2,4)plot(x,y4)grid onxlabel('x')ylabel('y')title('cot(x)')legend('y=cot(x)')3.解方程组：a=[321;1-13;24-4]；b=[7;6;-2]；x=a\b4.请实现下图：x y x=linspace(0,4*pi,1000);y1=sin(x);y2=sin(2*x);plot(x,y1,'--',x,y2,'b*')grid onxlabel('x');ylabel('y');title('耿蒙蒙')legend('sin(x)','sin(2*x)')5.请在x ，y 在（-2，2）内的z=xexp (-x 2-y 2)绘制网格图[x,y]=meshgrid(-2:0.1:2);z=x.*exp (-x.^2-y.^2);mesh(x,y,z)6.请实现peaks 函数：-33x Peaksy [x,y]=meshgrid(-3:1/8:3);z=peaks(x,y);mesh(x,y,z)surf(x,y,z)shading flataxis([-33-33-88])xlabel('x');ylabel('y');title('Peaks')7.请在x=[0,2],y=[-0.5*pi,7.5*pi],绘制光栅的振幅为0.4的三维正弦光栅。
Matlab语言与计算仿真上机考试试卷
2008级自动化Matlab 语言与计算仿真上机考试试卷要求：1、 只准带书查阅，不许带其它辅导资料，每个人独立作答；2、 请将试题解答程序和仿真结果拷贝到word 文件中并加以整理，写清题号；3、 Word 文件以自己“学号+姓名”命名，注意答题时时刻保存，答题结束后等待老师将word 文件拷贝下来方可离开。
1（10分）、计算矩阵a=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡987654321的特征值分解和Lu 分解。
2（10分）、求解下述方程组的解：⎪⎪⎩⎪⎪⎨⎧=-++=--=-++=++56533332921w z y x w y x w z y x z y x 3（10分）、已知 A =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--231134101，求：方阵A 的行列式值、秩、逆、特征值及对应的特征向量。
4（15分）、已知某系统的输入电压为x 时，输出电流为y ，且y=F(x)。
现测得数据如下：根据测得数据，求出四次拟合多项式(x)p 4。
5（15分）、计算积分⎰--51x dxe2。
6（10分）、求二阶微分方程的解：cos622=+xdx yd，初始条件为：y(0)=5，y'(0)=0。
7（15分）、有一组测量数据满足-at e=y，t的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线，并添加标题-at e=y和图例框。
8（15分）、有一组学生的考试成绩（见表），成绩在100分时为满分，成绩在90-99之间为优秀，成绩在80-89分之间为良好，成绩在60-79分之间为几个，成绩在60分以下时为不及格。
编制一个根据成绩划分等级的程序。
matlab上机练习及答案
第二、三次上机练习：目的：运行课本第四章及课堂上讲过的例子，掌握Matlab 的流程控制语句、函数及脚本文件的编程、调试方法。
作业：1、 完成下列操作：1） 求[100，999]之间能被21整除的个数。
2） 建立一个字符串向量（要求字符串向量中必须包含自己的姓名首字母，大小写均可），删除其中的大写字母2. 编写脚本文件，实现用magic(6)产生一矩阵，用for 循环指令求解其所有元素的和。
3. 定义一个函数文件，求∑=ni m i 1，要求在函数文件中包含能够通过help 查询到的说明；然后调用该函数文件求∑∑∑===++101501210011k k k k k k 的值。
4. 已知)7.1cos(12ln )7.1sin(++++=x xx y π，当x 取-3.0，-2.9，-2.8，…，2.8，2.9，3.0时， 1） 求各点的函数值；2） 求这些数据的平均值；5、求分段函数的值。
222603565231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩， 且， 0且及， 其它用if 语句实现，分别输出x=-5.0，-3.0，1.0，2.0，2.5时的值。
6、输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。
其中90分~100分为A ，80分~89分为B ，70分~79分为C ，60分~69分为D ，60分以下为E 。
要求：分别用if 语句和switch 语句实现。
7、根据222221111...,6123n ππ=++++求的近似值。
当n 分别取100、1000、10000时，结果是多少？（要求：分别用循环结构和向量运算来实现）8、已知n=1时，f 1=1；n=2时，f 2=0；n=3时，f 3=1；n>3时，f n =f n-1 -2f n-2+ f n —3； 求f 1~ f 100中，最大值、最小值以及各数之和。
Matlab上机练习题及答案
Matlab 上机练习题及答案---------------------------------------------------------------------1、 矩阵Y= ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3472123100451150425，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。
取出子矩阵⎥⎦⎤⎢⎣⎡21301，并求该矩阵的维数。
解：命令为：Y=[5,2,4;0,15,1;45,100,23;21,47,3] Y(2,3) Y(10)sub2ind([4 3],2,3)[i,j]=ind2sub([4 3],10)find(Y==100) sub2ind([4 3],3,2)B=Y(2:2:4,3:-2:1) 或 B=Y([2 4],[3 1]) [m n]=size(Y)---------------------------------------------------------------------2、已知矩阵A=[1 0 -1 ；2 4 1； -2 0 5]，B=[0 -1 0；2 1 3；1 1 2] 求2A+B 、A 2-3B 、A*B 、B*A 、A .*B ，A/B 、A\B解：命令为：A=[1 0 -1 ；2 4 1； -2 0 5] B=[0 -1 0；2 1 3；1 1 2] E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B---------------------------------------------------------------------3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵，并计算两者的乘积。
解：命令为： A=eye(3,4) B=8*ones(4)C=A*B---------------------------------------------------------------------4、创建矩阵a=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------7023021.5003.120498601,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b 构成矩阵d ，计算a*b 、c<d ,c&d, c|d ,~c|~d解：命令为：a=[-1,0,-6,8;-9,4,0,;0,0,,-2;0,-23,0,-7] b=a(:,[1 2]) c=a([1 2],:) d=b ’ e=a*b f=c<d g=c&d h=c|d i=~c|~d---------------------------------------------------------------------5、求!201∑=n n解：命令文件为 sum=0; s=1;for n=1:20 s=n*s; sum=sum+s; end sum---------------------------------------------------------------------6、求a aa aaa aa a S n ++++=得值，其中a 是一个数字，由键盘输入，表达式中位数最多项a 的个数，也由键盘输入。
MATLAB上机模拟考试题(一)
(1)在MATLAB的命令窗口中执行_____命令，将命令窗口的显示内容清空。
（） A.clearB.clcC.echo offD.cd(2)在MATLAB的命令窗口中执行_____命令，使数据输出显示为十六进制表示。
（）A.format longB.format ratC.format hexD.format short e(3)下列变量名中_____是合法的。
（） A.x*y,a,1 B.x\y,a1234 C.end,1 bcx D.char_1,i,j(4)已知x=0:5,则x有_____个元素。
（） A.5 B.6 C.7 D.8(5)一下运算符中哪个的优先级最高_____。
（） A./ B.^ C.~= D.&(6)使用检测函数isnumeric(10)的结果是_____。
（） A.1 B.0 C.false D.true(7)三维图形中默认视角是_____度。
（） A.方位角=0 俯仰角=90 B.方位角=90 俯仰角=0 C.方位角=37.5 仰俯角=30 D.方位角=0 仰俯角=180(8)将符号表达式化简为因式分解因式分解因式分解因式分解形式，使用_____函数。
（）A.collectB.expandC.hornerD.factor(9)运行以下命令，则_____描述是正确的。
（）>>syms a b c d >>A=[a b;c d]A.A占用的内存小于100BB.创建了5个符号变量C.A占用的内存是a b c d的总和D.不存在（10）已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。
（11）运行命令bitor（8,7）的结果是_____。
（12）运行以下命令：>>x=0:10; >>y1=sin(x); >>y2=5*sin(x); >>y3=[10*sin(x );20*sin(x)]; >>plot(x,y1,x,y2,x,y3)则在一个图形窗口中，可以看到_____条曲线。
Matlab上机题库及详细答案
解：c=input('请输入一个字符','s');
if c>='A' & c<='Z'
disp(setstr(abs(c)+abs('a')-abs('A')));
elseif c>='a'& c<='z'
disp(setstr(abs(c)- abs('a')+abs('A')));
(2)均值为0.6方差为0.1的5阶正态分布随机矩阵
解：>>x=20+(50-20)*rand(5)；
>>y=0.6+sqrt(0.1)*randn(5)
例13：将101~125等25个数填入一个5行5列的表格中，使其每行每列及对角线的和均为565。
解：M=100+magic(5)
M =117 124 101 108 115
-0.4606
0.3848
例2、用简短命令计算并绘制在0x6范围内的sin(2x)、sinx2、sin2x。
解：x=linspace(0,6)
y1=sin(2*x),y2=sin(x.^2),y3=(sin(x)).^2;
plot(x,y1,x, y2,x, y3)
例3：画出指数衰减曲线y1=exp(-t/3)*sin(3*t)和它的包络y2=exp(-t/3),t的取值范围是（0，4pi）。
MATLAB上机模拟考试题(一)
(1)在MATLAB的命令窗口中执行_____命令，将命令窗口的显示内容清空。
（） A.clearB.clcC.echo offD.cd(2)在MATLAB的命令窗口中执行_____命令，使数据输出显示为十六进制表示。
（）A.format longB.format ratC.format hexD.format short e(3)下列变量名中_____是合法的。
（） A.x*y,a,1 B.x\y,a1234 C.end,1 bcx D.char_1,i,j(4)已知x=0:5,则x有_____个元素。
（） A.5 B.6 C.7 D.8(5)一下运算符中哪个的优先级最高_____。
（） A./ B.^ C.~= D.&(6)使用检测函数isnumeric(10)的结果是_____。
（） A.1 B.0 C.false D.true(7)三维图形中默认视角是_____度。
（） A.方位角=0 俯仰角=90 B.方位角=90 俯仰角=0 C.方位角=37.5 仰俯角=30 D.方位角=0 仰俯角=180(8)将符号表达式化简为因式分解因式分解因式分解因式分解形式，使用_____函数。
（）A.collectB.expandC.hornerD.factor(9)运行以下命令，则_____描述是正确的。
（）>>syms a b c d >>A=[a b;c d]A.A占用的内存小于100BB.创建了5个符号变量C.A占用的内存是a b c d的总和D.不存在（10）已知数组a=[1 2 3;4 5 6;7 8 9],则a(:,end)是指_____元素。
（11）运行命令bitor（8,7）的结果是_____。
（12）运行以下命令：>>x=0:10; >>y1=sin(x); >>y2=5*sin(x); >>y3=[10*sin(x );20*sin(x)]; >>plot(x,y1,x,y2,x,y3)则在一个图形窗口中，可以看到_____条曲线。
matlab上机题汇总(必做)
上机题汇总1设置matlab 的工作环境，将工作目录设置为d:\work ，添加搜索目录d:\example2在matlab 的命令窗口里完成如下计算，其中t 的值分别取-1,0,1，表达式如下：4/3)2t y e π-=3自行产生一个5行5列的数组，得到最中间的三行三行矩阵。
4用magic 产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置5求方程组的根 x 1+4x 2-3x 3=22x 1+5x 2-x 3=11x 1+6x 2+x 3=126已知：一个多项式的系数向量是p=[1 -6-72 -27]，求这个多项式的根。
7已经两个多项式的系数分别是：[1 2 3 4]和[1 4 9 16]，请求这两个多项式的乘积，及商和余数。
8给定一个多项式的根是[-5 -3+4i -3-4i],求原来的多项式9 A=[2 3 4；1 5 7；6 2 5]用什么函数，保证第一列排序的时候，其他列跟着变化。
最后的结果是：[1 5 7；2 3 4；6 2 5]10补充题：电路分析电路分析常常涉及到对方程组的求解，这些方程常常是利用描述电流进入和离开节点的电流方程，或者描述电路中网络回路上的电压的电压方程得出的。
下图描述了3个回路电压的方程方程式如下：11121222132423432532()0()()0()0V R I R I I R I I R I R I I R I I R I V -++-=⎧⎪-++-=⎨⎪-++=⎩假设5个电阻值为已知，2个电压值也为已知，求3个电流值。
11自行产生一个5行5列的数组，用两种方法得到最中间的三行三行矩阵。
12根据a=reshape(-4:4,3,3)做一个矩阵，然后（1）取出所有大于0的元素构成一个向量（可推广到从一个矩阵里筛选出符合条件的元素组成一个向量）（2）将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示(可推广到将原矩阵中不符合条件的全用0来表示，符合条件的值不变)。
Matlab上机题库及详细答案
case num2cell(25:49) %价格大于等于2500但小于5000
Байду номын сангаасrate=10/100;
otherwise %价格大于等于5000
rate=14/100;
end
price=price*(1-rate) %输出商品实际销售价格
例19已知，当n=100时，求的值。
解：程序如下：
y=0;n=100;for i=1:n；y=y+1/(2*i-1);End
例20：一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数
解：for m=100:999
m1=fix(m/100);m2=rem(fix(m/10),10);m3=rem(m,10);
if m==m1*m1*m1+m2*m2*m2+m3*m3*m3；disp(m)；end
D=diag(1:5);D*A %用D左乘A，对A的每行乘以一个指定常数
例15：输入x,y的值，并将它们的值互换后输出。
程序如下：
x=input('Input x please.')；y=input('Input y please.');z=x;x=y;y=z;
disp(x);disp(y);
例16：:求一元二次方程ax2 +bx+c=0的根。
switch fix(price/100)
case {0,1} %价格小于200
rate=0;
case {2,3,4} %价格大于等于200但小于500
rate=3/100;
case num2cell(5:9) %价格大于等于500但小于1000
MATLAB机考样题(带答案)
MATLAB 机考样题： （1）Generateandplotsequence121[]2cos() [][8x n n and x nx n π==-n1=-20:20; x1=2*cos(pi/8*n1) n2=n1-4;x2=2*cos(pi/8*n2) subplot(211) plot(n1,x1); subplot(212) plot(n2,x2)？（2）Write a MATLAB program to compute and plot the impulse response of a causal finite-dimensional discrete-time system characterized by a difference equation of the following form:3]-x [n 86.0 2]-1.32x [n 1]-0.34x [n x [n]8.13]-0.72y[n -2]-0.5y[n 1]-0.3y[n y[n]--+=++N=input('请输入你要求的点数N='); num=[1.8 0.34 -1.32 -0.86]; den=[1 0.3 0.5 -0.72]; x=[1 zeros(1,N-1)]; y=filter(num,den,x) plot(0:N-1,y)（3）Write a MATLAB program to compute and display the poles and zeros ,generate the pole-zero plot of a polynomials in z -13213211768.018.052.115.1082.2393.61.8)(-------+++-+=z z z z z z z Hnum=[8.1 6.93 -23.82 10.5]; den=[1 1.52 0.18 -0.1768]; sos=tf2sos(num,den) zplane(num,den)（4）Try to give a program to evaluate the following DTFT in the rangeπω≤≤0 :ωωωωωωωω4324321245535952)(j j j j j j j j e e e e e e e e z G --------++++++++=%由于用freqz 计算频点至少是2个，所以至少输入两个频点w1=input('请输入你要计算的频点w1='); w2=input('请输入你要计算的频点w2='); w=[w1 w2]; num=[2 5 9 5 3]; den=[5 45 2 1 1];h=freqz(num,den,w)（6）Write a MATLAB program to compute and plot the magnitude response of a causal LTI discrete-time system with a transfer function given by2127.05.01)1(15.0)(---+--=z z z z Hnum=0.15*[1 0 -1]; den=[1 -0.5 0.7];[h,w]=freqz(num,den,512); plot(w/pi,abs(h))（7）Consider the following FIR transfer function:123456()10.60.490.480.140.120.09H z z z z z z z ------=++---+Using MATLAB to determine its zero locations and plot its magnitude and phase response.h=[1 0.6 .49 -0.48 -0.14 -0.12 0.09];figure(1) zplane(h,1);[H,w]=freqz(h,1,512); figure(2)plot(w/pi,abs(H)) figure(3)plot(w/pi,angle(H))（8）Givena signal()4cos0.1x t t t π=+, when using a samplingfrequency f T = 20KHz, plot the magnitude and phase spectrum of the sampled sequence(given length-64).fs=2e4; n=(0:63)/fs;x=4*n+cos(0.1*pi*n); h=fft(x,512); figure(1)plot(0:2/511:2,abs(h))figure(2)plot(0:2/511:2,angle(h))（9）design an IIR butterworth digital lowpass filter with the following specifications: sampling rate of 40kHz, passband edge frequency of 4kHz, stopband edge frequency of 8kHz, passband ripple of 0.5dB, and a minimum stopband attenuation of 40dB,plot frequency-magnitude and check if your design fits the specification.fs=40;wp=4*2/fs;ws=8*2/fs;ap=0.5;as=40;[num,den]=butter(n,wn);[h,w]=freqz(num,den,512);figure(1)plot(w/pi,20*log10(abs(h)))axis([0 1 -50 0])figure(2)subplot(211)plot(w/pi,20*log10(abs(h)))axis([0 wp -0.5 0]);title('通带纹波')subplot(212)plot(w/pi,20*log10(abs(h)))axis([ws 1 -50 -30])title('阻带纹波')（10）Design a Hanning FIR lowpass filter meeting the followingspecifications: passband edge frequency=2kHz, stopband edgefrequency=2.5kHz, passband ripple δp=0.005, stopband rippleδs=0.005,and sampling rate of 10kHz.Plot its gain and phase responses andcheck if it meets the specifications?ft=10;fp=2;wp=2*pi*fp/ft;ws=2*pi*fs/ft;ds=0.005;ap=20*log10(1-ds);as=20*log10(ds)wc=(wp+ws)/2;dw=ws-wp;M=ceil(3.11*pi/dw);N=2*M;b=fir1(N,wc/pi,hann(N+1));[h,w]=freqz(b,1,512);figure(1)plot(w/pi,20*log10(abs(h)));axis([0 1 -50 0])title('magnitude response');figure(2)plot(w/pi,unwrap(angle(h)));title('phase response');figure(3) %局部放大，观察通带与阻带subplot(211)plot(w/pi,20*log10(abs(h)));axis([0 wp/pi ap 0])title('通带纹波')subplot(212)plot(w/pi,20*log10(abs(h)));axis([ws/pi 1 as 0])title('阻带纹波')%从图中可以看出，通带和阻带中纹波都不满足要求，所以不满足指标%as= -46.0206<43.9 所以不能用hanning窗设计%应当用hamming或blackman窗设计。
完整版本matlab上机练习试题含答案
1.计算 a6 9 3 2 4 1 的数组乘积1 71618 82 7 与 b6 8 5 计算多项式除法 (3x 3+13 x 2+6x+8)/( x+4)54>> a=[6 9 3;2 7 5]; >> d=deconv([3 13 6 8],[1 4])>> b=[2 4 1;4 6 8]; d =>> a.*b 31 2ans =6 求欠定方程组2 4 7 4x 8 的最小范数解12 3639 35 6 584240>> a=[2 4 7 4;9 3 5 6];4 9 2 37>> b=[8 5]';>> x=pinv(a)*b2.对于 AXB ，假如 A7 6 4 ，B 26，求解 X 。
3 5 728x =-0.2151>> A=[4 9 2;7 6 4;3 5 7]; 0.4459 >> B=[37 26 28] ’;0.7949 >> X=A\B 0.2707X =7 用符号函数法求解方程 at 2+b* t+c=0-0.5118 >> r=solve('a*t^2+b*t+c=0','t') 4.0427 r =1.3318[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] 1 2 5 ， b8 7 4 [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]3. a64 362 ，察看 a 与 b 之间的3a 11 a 12六种关系运算的结果8求矩阵 A的队列式值、逆和特点根a 21 a 22 >> a=[1 2 3;4 5 6];>> b=[8 –7 4;3 6 2]; >> syms a11 a12 a21 a22; >> a>b >> A=[a11,a12;a21,a22]ans =>> AD=det(A) % 队列式 0 1 0 >> AI=inv(A) % 逆 11>> AE=eig(A) % 特点值>> a>=b A =ans =[ a11, a12] 0 1 0 [ a21, a22] 11AD =>> a<b a11*a22-a12*a21 ans =AI =1 0 1 [ -a22/(-a11*a22+a12*a21), a12/(-a11*a22+a12*a21)]1[ a21/(-a11*a22+a12*a21), -a11/(-a11*a22+a12*a21)] >> a<=b AE =ans =1 0 1 [ 1/2*a11+1/2*a22+1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)]1>> a==b [ 1/2*a11+1/2*a22-1/2*(a11^2-2*a11*a22+a22^2+4*a12*a21)^(1/2)]ans =9 因式分解： x 45x 3 5 x 2 5x 60 0 0 >> syms x;>> f=x^4-5*x^3+5*x^2+5*x-6; >> a~=b >> factor(f) ans =ans =1 1 1 (x-1)*(x-2)*(x-3)*(x+1) 111ax 214 计算多项式乘法 (x 2+2x+2)(x 2+5x+4)，用符号微分求 df/dx 。
MATLAB机考样题（带答案）
MATLAB机考样题（带答案）MATLAB 机考样题：（1）Generate and plot sequence121[]2cos() [][4]8x n n and x n x n π==-, with 20n 20-≤≤.n1=-20:20;x1=2*cos(pi/8*n1);n2=n1-4;x2=2*cos(pi/8*n2);subplot(2,1,1);plot(n1,x1);subplot(212);plot(n2,x2); （2）Write a MATLAB program to compute and plot the impulse response of a causal finite-dimensional discrete-time system characterized by a difference equation of the following form:3]-x[n 86.0 2]-1.32x[n 1]-0.34x[n x[n]8.13]-0.72y[n -2]-0.5y[n 1]-0.3y[n y[n]--+=++N=input('请输⼊你要求的点数N=');num=[1.8 0.34 -1.32 -0.86];den=[1 0.3 0.5 -0.72]; x=[1 zeros(1,N-1)];（单位冲击）y=filter(num,den,x);plot(0:N-1,y);（3）Write a MATLAB program to compute and display the poles and zeros , to compute and display the second-order factored form , and to generate the pole-zero plot of a z-transform that is a ratio of two polynomials in z -1. Using this program to analyze the following G(z):3213211768.018.052.115.1082.2393.61.8)(-------+++-+=z z z z z z z Hnum=[8.1 6.93 -23.82 10.5];den=[1 1.52 0.18 -0.1768];sos=tf2sos(num,den) %tf2sos 表⽰为1/z 的升幂zplane(num,den)（4）Try to give a program to evaluate the following DTFT in the range πω≤≤0 :ωωωωωωωω4324321245535952)(j j j j j j j j e e e e e e e e z G --------++++++++=%由于⽤freqz 计算频点⾄少是2个，所以⾄少输⼊两个频点w1=input('请输⼊你要计算的频点w1=');w2=input('请输⼊你要计算的频点w2=');w=[w1 w2];num=[2 5 9 5 3];den=[5 45 2 1 1];h=freqz(num,den,w)（6）Write a MATLAB program to compute and plot the magnitude response of a causal LTI discrete-time system with a transfer function given by2127.05.01)1(15.0)(---+--=z z z z Hnum=0.15*[1 0 -1];den=[1 -0.5 0.7];[h,w]=freqz(num,den,512);plot(w/pi,abs(h)); （7）Consider the following FIR transfer function:123456()10.60.490.480.140.120.09H z z z z z z z ------=++---+Using MATLAB to determine its zero locations and plot its magnitude and phase response.h=[1 0.6 .49 -0.48 -0.14 -0.12 0.09];figure(1)zplane(h,1);[H,w]=freqz(h,1,512);figure(2)plot(w/pi,abs(H));figure(3)plot(w/pi,angle(H));（8）Given a signal ()4cos0.1x t t t π=+, when using a sampling frequency f T = 20KHz, plot the magnitude and phase spectrum of the sampled sequence(given length-64).fs=2e4;n=0:63;x=4*n/fs+cos(0.1*pi*n/fs);h=fft(x,1024);figure(1);plot(0:2/1023:2,abs(h));figure(2);plot(0:2/1023:2,angle(h));（9）design an IIR butterworth digital lowpass filter with the following specifications: sampling rate of 40kHz, passband edge frequency of 4kHz, stopband edge frequency of 8kHz, passband ripple of 0.5dB, and a minimum stopbandattenuation of 40dB,plot frequency-magnitude and check if your design fits the specification.fs=40;wp=4*2/fs; %wp<1,没有乘piws=8*2/fs; %ws<1,没有乘piap=0.5;as=40;[n,wn]=buttord(wp,ws,ap,as);[num,den]=butter(n,wn);[h,w]=freqz(num,den,512);figure(1);plot(w/pi,20*log10(abs(h)));axis([0 1 -50 0]);figure(2);subplot(2,1,1);plot(w/pi,20*log10(abs(h)))axis([0 wp -0.5 0]);title('通带纹波');subplot(2,1,2);plot(w/pi,20*log10(abs(h)));axis([ws 1 -50 -30]);title('阻带纹波');（10）Design a Hanning FIR lowpass filter meeting the following specifications: passband edge frequency=2kHz, stopband edge frequency=2.5kHz, passband ripple δp=0.005, stopband rippleδs=0.005, and sampling rate of 10kHz.Plot its gain and phase responses and check if it meets the specifications?ft=10;fp=2;fs=2.5;wp=2*pi*fp/ft;ws=2*pi*fs/ft;ds=0.005;ap=20*log10(1-ds)as=20*log10(ds)wc=(wp+ws)/2;dw=ws-wp;M=ceil(3.11*pi/dw);N=2*M;b=fir1(N,wc/pi,hann(N+1));[h,w]=freqz(b,1,512);figure(1);plot(w/pi,20*log10(abs(h)));axis([0 1 -50 0]);title('magitude response');figure(2);plot(w/pi,unwrap(angle(h)));title('phase response');figure(3);subplot(211);plot(w/pi,20*log10(abs(h)));axis([0 wp/pi ap 0]);title('通带纹波')subplot(212);plot(w/pi,20*log10(abs(h)));axis([ws/pi 1 as 0]);title('阻带纹波');%从图中可以看出，通带和阻带中纹波都不满⾜要求，所以不满⾜指标%as= -46.0206<43.9 所以不能⽤hanning 窗设计%应当⽤hamming 或blackman 窗设计(11)Writing a MATLAB program to compute 128-point DFT of the followingsequence, you must firstly use DFT definition (directly computing DFT) to compute and use MATLAB function to test the result. Plot the two results in one figure.n=0:31;k=1:128;x=sin(pi*n/4);Xk1=zeros(1,128);for t=1:128 for m=1:32Xk1(t)=Xk1(t)+x(m)*exp(-1i*2*pi*(t-1).*(m-1)/128);endendsubplot(2,1,1);plot(k,Xk1);Xk2=fft(x,128);subplot(2,1,2);plot(k,Xk2);(12)Using the function fir1 and window of Kaiser, design a linear-phase FIR lowpass filter meeting the followingspecifications: passband edge frequency=2kHz, stopband edge frequency=2.5kHz, passband ripple δp=0.005, stopband ripple δs=0.005, and sampling rate of 10kHz.Plot its gain and phase responses and check if it meets the specifications?Design a Type 1 Chebyshev IIR lowpass filter meeting the specifications as below: sampling rate of 12kHz, passband edge frequency of 2.1kHz, stopband edge frequency of 2.7kHz, passband []sin(4),031g n n n π=≤≤ripple of 0.6dB, and a minimum stopband attenuation of 45dB. Write down the exact expression for the transfer function generated. Does your design meet the specifications?Fp=2100;Fs=2700;Ft=12000;Rp=0.6;Rs=45;Wp=2*Fp/Ft;Ws=2*Fs/Ft;[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs)[B,A]=cheby1(N,Rp,Wn)[h,w]=freqz(B,A,512);figure(1);plot(w/pi,20*log10(abs(h)));axis([0 1 -50 0]);figure(2);subplot(2,1,1);plot(w/pi,20*log10(abs(h)))axis([0 Wp -0.6 0]);title('通带纹波');subplot(2,1,2);plot(w/pi,20*log10(abs(h)));axis([Ws 1 -50 -30]);title('阻带纹波');(13)Using the function fir1 and window of Kaiser, design a linear-phase FIR lowpass filter meeting the following specifications: passband edge frequency=2kHz, stopband edge frequency=2.5kHz, passband ripple δp=0.005, stopband rippleδs=0.005, and sampling rate of 10kHz.Plot its gain and phase responses and check if it meets the specifications?[n,wn,beta,typ]=kaiserord([2000 2500],[1 0],[0.005 0.005],10000);b=fir1(n,wn,kaiser(n+1,beta),'noscale');[h,omega]=freqz(b,1,256);subplot(2,1,1);plot(omega/pi,20*log10(abs(h)));xlabel('\omega/\pi'); ylabel('Gain, dB');subplot(2,1,2)plot(omega/pi,angle(h));gridtitle('Phase Spectrum');xlabel('\omega/\pi'); ylabel('Phase, radians');(14)Given a signal when using a sampling frequency fT= 20KHz, plot the magnitude and phase spectrum of the sampled sequence(givenlength-64). ft=20000;n=1:64;x=sin(0.1*pi*n/ft)+2*cos(0.3*pi*5*n/ft)+3*sin(0.5*pi*n/ft); subplot(2,1,1);[h,w]=freqz(x,1,256);plot(w/pi,abs(h));title('Magnitude spectrum of the sampled samples');xlabel('\omega/\pi'); ylabel('Magnitude');subplot(2,1,2);plot(w/pi,angle(h));title('Phase spectrum of the sampled samples');xlabel('\omega/\pi'); ylabel('Phase');(15)Write a MATLAB program to compute the first L samples of the inverse of rational z-transforms where the value of L is provided by the user through the command input. Using this program to compute and plot the first 50 samples of the inverse of following G(z). Use the command stem for plotting the sequence generated by the inverse transform L=input('input the L=');r=[10/4 -8/2];p=[-1/4 -1/2];k=-2;[B,A]=residuez(r,p,k)[h t]=impz(B,A,L)stem(t,h);Writing a MATLAB program to compute the circular convolution of two length-N sequences via the DFT-based approach. Using this program to determine the following pair of sequences:g[n]={7, 4, -9, 0, 2, -5}, h[n]={1, -1, 2, 0, 10, 5} orAnd plot the result sequencex1=[7 4 -9 0 2 -5];x2=[1 -1 2 0 10 5];L=length(x1);y=zeros(1,L);x2tr = [x2(1) x2(L:-1:2)];for k = 1:L,sh = circshift(x2tr', k-1)';h = x1.*sh;y(k) = sum(h);enddisp(y);()sin(0.1)2cos(0.3)3sin(0.5)x t t t t πππ=++5.0||,284102)(113>+-++-=--z z z z Gn=0:length(x1)-1;stem(n,y);n=0:50;x=sin(5*pi*n/16);stem(n,x); %或者如下syms n1; %x1=x’x1=sin(5*pi*n1/16).^2; %y=x*x1 y=symsum(x1,n1,0,50)y1=double(y)。
(完整版)matlab上机题目一
上机一1、（1（2） 查询所有人的成绩。
2、已知：A=[27,24,0,-1;23,35,1,14;e 4,-8,13,20;10,12,19,21;11,18,ln2,0]， 求：（1）求A 中第4、5、2行与第2、4、1列组成的C 阵，并将C 中第2行元素全改成8，得到D 阵；（2） 创建与A 同维的G 阵，使C 阵左右翻转后位于G 中第2行至第4行、第2列至第4列位置，且G 中第1行第1列元素为1，其余元素为0。
（3） 求A 阵除以3的余数阵H 。
3、已知：A=[5,17,-2,3;-1,0,21,11;14,9,-10,8]，B=[2,5,8,-1;12,-2,7,8;11,13,16,4] 求：（1）将A 除以B 阵的对应元素，并将结果朝负无穷大方向取整。
（2）求A./B ， B.\A ，A.*B 和A*B ’的值。
（3）求2*ones(3,4)-eye(3,4)|A-B>B&B-8的值.4、要求在闭区间[0，2π ]上产生50个等距采样的一维数组。
试用两种不同的指令实现。
5、对矩阵，⎥⎦⎤⎢⎣⎡=4321A 问如下MATLAB 指令：)(),(),5.0(^),5.0(.^A sqrtm A sqrt A A 所得的结果相同吗？他们中的哪个结果为复数矩阵？ 6、说出MATLAB 指令A （3，1，2，： ）=1：4所产生数组的维数、大小和长度；然后对A 进行降维处理；最后指出所有非零元素的“全下标”位置。
7、现有五个数组： (2×3)的实数数组A ； (2×3)的复数数组B ； (2×3)的字符串数组C ； (2×3)的异质（元胞）数组D ； (2×3)的结构（构架）数组E 。
问，它们各占多少字节？8、写出使以下这段文字成为字符串的MATLAB 程序。
注意保持这段文字的格式。
在英式用法中，引号通常是单引号，如‘Fire!’。
Matlab上机题代码及结果4题
例1下图描述了六个城市之间的航空航线图，其中1、2、......、6表示六个城市，带箭头线段表示两个城市之间的航线。
用MATLAB软件完成以下操作：（1）构造该图的邻接矩阵A；（2）若某人连续乘坐五次航班，那么他从哪一个城市出发到达哪一个城市的方法最多？（3）若某人可以乘坐一次、二次、三次或四次航班，那么他从哪一个城市出发总是不能达到哪一个城市？航空航线图（六城市）解：（1）构造邻接矩阵；（2）计算矩阵可达矩阵，找出该矩阵的最大元素，并确定它所在的位置；（3）计算可达矩阵，找出该矩阵中零元素的位置。
在MATLAB软件的M编辑器中编写m文件：% 图与矩阵clearA=[0,1,0,0,0,1;0,0,1,1,0,0;0,0,0,1,1,0;0,1,0,0,0,0;1,0,1,0,0,0;0,1,0,0,1,0]; % 构造邻接矩阵B=A^5;C=A+A^2+A^3+A^4;disp('邻接矩阵A为：');disp(A);disp('矩阵A^5为：');disp(B);m=max(max(B)); % 计算矩阵B的最大值[m_i,m_j]=find(B==m); % 寻找矩阵B中元素等于m的位置fprintf('矩阵A^5最大值%d的位置在：\n',m);disp([m_i,m_j]);disp('矩阵A＋A^2＋A^3＋A^4为：');disp(C);[z_i,z_j]=find(C==0); % 寻找矩阵C中零元素的位置disp('矩阵A＋A^2＋A^3＋A^4零元素的位置在：');disp([z_i,z_j]);在MATLAB命令窗口中输入m文件名称,计算结果为：邻接矩阵A为：0 1 0 0 0 10 0 1 1 0 00 0 0 1 1 00 1 0 0 0 01 0 1 0 0 00 1 0 0 1 0矩阵A^5为：2 5 5 53 12 4 43 2 02 3 5 5 2 10 2 1 3 2 12 6 4 5 4 11 4 4 7 4 2矩阵A^5最大值7的位置在：6 4矩阵A＋A^2＋A^3＋A^4为：2 6 5 6 4 21 4 4 6 3 12 5 4 5 4 11 3 3 3 1 03 5 6 64 23 6 6 54 1矩阵A＋A^2＋A^3＋A^4零元素的位置在：4 6从计算结果中可以看出，矩阵A^5最大值出现在矩阵的第六行第四列，说明：这个人如果从城市6出发连续乘坐五次航班后到达城市4，他可以选择的乘机路线最多，共有7种不同的方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
页脚内容1Matlab 上机模拟试题1、 计算当1,2,5x =-时，21()21x f x x e x -=+--的值。
2、 计算300tan 605cos30+-3、 求解矩阵方程：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡3542343122321X 4、求解矩阵方程：⎥⎦⎤⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡212101343122321X页脚内容25、 计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。
6、 用MATLAB 计算差分方程 P138 157()0.7(1)0.45(2)0.6(3)0.8()0.44(1)0.36(2)0.02(3)y n y n y n y n x n x n x n x x +-----=--+-+-当输入序列为()()x n n δ=时的输出结果()y n ,n 0,1,...40?=。
1231230.80.440.360.02()10.70.450.6z z z H z z z z -------++=+--页脚内容3N=41;a=[0.8 -0.44 0.36 0.22];b=[1 0.7 -0.45 -0.6];x=[1 zeros(1,N-1)];k=0:1:N-1;y=filter(a,b,x);stem(k,y)xlabel('n');ylabel('幅度')7、用MATLAB 计算差分方程 P156()0.7(1)0.45(2)0.6(3)0.8()0.44(1)0.36(2)0.02(3)y n y n y n y n x n x n x n x x +-----=--+-+-所对应的系统函数的频率响应。
B=[0.8 -0.44 0.36 0.02];A=[1 0.7 -0.45 -0.6];freqz(B,A);23230.80.440.360.02()10.70.450.6j j j j j j j e e e H ee e e ωωωωωωω--------++=+--8、用FFT 计算序列(){13112331}x n =-的频谱。
页脚内容4xn=[1 3 -1 1 2 3 3 1];Xn=fft(xn);figure(1);plot(abs(Xn));9、求下列直接型系统函数的零、极点。
1234123410.10.30.30.2()10.10.20.20.5z z z z H Z z z z z ------------=++++ b=[1,-0.1,-0.3,-0.3,-0.2];a=[1,0.1,0.2,0.2,0.5];disp('H(z)');printsys(b,a,'z');[z,p,k]=tf2zp(b,a)10、画出序列()()()()()() x n n 2n 13n 24n 35n 4δδδδδ=+-+-+-+-的图形并补零，序列长度为20.xn=[1 2 3 4 5];b=zeros(1,15);页脚内容5xn=[xn,b];stem(xn)11、产生一个1024点的正弦加白噪声序列，并绘其频谱。
x=[0:1023];y=sin(x);y1=awgn(y,10);freqz(y1)12、已知信号()2(cos 2)cos100f x x x =，求其FFT 变换的幅值和相位值，并绘出结果。
y=2.*(cos(2.*x)).*(cos(100.*x));y1=fft(y);freqz(y1)13、给定()()()()()()x n n 2n 13n 24n 35n 4δδδδδ=+-+-+-+-，()()()()()y n n 2n 1n 22n 3δδδδ=+-+-+-。
分别求两序列相加、相卷积。
用4个子图显示。
n1=0:4;x1=[1 2 3 4 5];subplot(4,1,1),stem(n1,x1);ylabel('x1');n2=0:3;x2=[ 1 2 1 2];subplot(4,1,2),stem(n2,x2);ylabel('x2');n=0:4;x1=[x1 ,zeros(1,5-length(n1))];x2=[ x2,zeros(1,5-length(n2))];x=x1+x2;y=conv(x1,x2);subplot(4,1,3),stem(n,x);ylabel('x');subplot(4,1,4),stem(y);ylabel('y')14、生成一个长度为50的汉宁窗，并观察期频率特性。
n=51;window=hanning(n);[h,w]=freqz(window,1);subplot(1,2,1)页脚内容6页脚内容7stem(window);subplot(1,2,2);plot(w/pi,20*log(abs(h)/abs(h(1))));15、已知有限长序列x(n)=[7,6,5,4,3,2],求DFT 和IDFT ，要求：画出序列傅立叶变换对应的幅度谱和相位谱；画出原信号与傅立叶逆变换IDFT[X(k)]的图形进行比较。
xn=[7,6,5,4,3,2];N=length(xn);n=0:(N-1);k=0:(N-1);Xk=xn*exp(-j*2*pi/N).^(n'*k);x=(Xk*exp(j*2*pi/N).^(n'*k))/N;subplot(2,2,1),stem(n,xn);title('x(n)');subplot(2,2,2),stem(n,abs(x));title('IDFT|X(k)|');subplot(2,2,3),stem(k,abs(Xk));16、已知周期序列的主值x(n)=[7,6,5,4,3,2],求x(n)周期重复次数为3次时的DFS 和IDFS 。
要求：画出原信号序列的主值和周期序列的图形；画出离散傅立叶变换对应的幅度谱和相位谱. xn=[7,6,5,4,3,2];title('|X(k)|'); subplot(2,2,4),stem(k,angle(Xk)); title('arg|X(k)|');N=length(xn);n=0:3*N-1;k=0:3*N-1;xn1=xn(mod(n,N)+1);Xk=xn1*exp(-j*2*pi/N).^(n'*k);subplot(2,2,1),stem(xn);title('原主值信号x(n)');subplot(2,2,2),stem(n,xn1);title('周期序列信号');subplot(2,2,3),stem(k,abs(Xk));title('|X(k)|');subplot(2,2,4),stem(k,angle(Xk));title('arg|X(k)|');17、求x(n)=[7,6,5,4,3,2], 0=<n<=5的DTFT，将（-2*pi，2*pi）区间分成500份。
要求：画出原信号；画出离散傅立叶变换对应的幅度谱和相位谱；求有限长序列x(n)=[7,6,5,4,3,2]在N=100时的DFT，并与DTFT 进行对比。
xn=[7,6,5,4,3,2];N=length(xn);n=0:N-1;页脚内容8w=linspace(-2*pi,2*pi,500);X=xn*exp(-j*n'*w);subplot(3,1,1),stem(n,xn,'k');ylabel('x(n)');subplot(3,1,2),plot(w,abs(X),'k');axis([-2*pi,2*pi,1.1*min(abs(X)),1.1*max(abs(X))]);ylabel('幅度谱');subplot(3,1,3),plot(w,angle(X),'k');axis([-2*pi,2*pi,1.1*min(angle(X)),1.1*max(angle(X))]);ylabel('相位谱');（2）N=100;xn=[7,6,5,4,3,2,zeros(1,N-6)];n=0:(N-1);k=0:(N-1);Xk=xn*exp(-j*2*pi/N).^(n'*k);x=(Xk*exp(j*2*pi/N).^(n'*k))/N;subplot(2,1,1),stem(k,abs(Xk));title('|X(k)|');页脚内容9页脚内容10 subplot(2,1,2),stem(k,angle(Xk));title('arg|X(k)|');18、绘出下列直接型系统函数的频率响应。
1234123410.10.30.30.2()10.10.20.20.5z z z z H Z z z z z ------------=++++ B=[1 -0.1 -0.3 -0.3 -0.2];A=[1 0.1 0.2 0.2 0.5];freqz(B,A)19、设采样周期T=250μs （采样频率f s =4kHz ），用脉冲响应不变法和双线性变换法设计一个三阶巴特沃兹滤波器，并绘出频率响应，其3dB 边界频率为f c =1kHz 。
[B,A]=butter(3,2*pi*1000,'s');[num1,den1]=impinvar(B,A,4000);[h1,w]=freqz(num1,den1);[B,A]=butter(3,2/0.00025,'s');[num2,den2]=bilinear(B,A,4000);[h2,w]=freqz(num2,den2);f=w/pi*2000;plot(f,abs(h1),'-.',f,abs(h2),'-');grid;xlabel('频率/Hz ')ylabel('幅值/dB')程序中第一个butter的边界频率2π×1000，为脉冲响应不变法原型低通滤波器的边界频率；第二个butter的边界频率2/T=2/0.00025，为双线性变换法原型低通滤波器的边界频率.图3.1给出了这两种设计方法所得到的频响，虚线为脉冲响应不变法的结果；实线为双线性变换法的结果。
脉冲响应不变法由于混叠效应，使得过渡带和阻带的衰减特性变差，并且不存在传输零点。