MATLAB上机题汇总(较难题,抽1题)合肥学院
MATLAB 上机 习题及答案
15、今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5，要求先求得P(x)=P1(x)+P2(x)，然后计算xi=0.2*i各点上的P(xi)（i=0,1,2,…,5）值。
p1=[1.0 0.0 0.0 -2.0 1.0];>> p2=[0.0 0.0 1.0 4.0 -0.5];>> p1x=poly2sym(p1);p2x=poly2sym(p2);>> p=p1x+p2xp =x^4+2*x+1/2+x^2>> x=0:5;>> x.^4+2*x+1/2+x.^2ans =0.5000 4.5000 24.5000 96.5000 280.5000 660.50001、试个MATLAB的工作空间中建立以下2个矩阵：A＝[1 2]1234B⎡⎤=⎢⎥⎣⎦，求出矩阵A和B的乘积，并将结果赋给变量C。
>> A=[1 2]A =1 2>> B=[1 23 4]B =1 23 4>> C=A*BC =7 102、利用MATLAB提供的帮助信息，了解inv命令的调用格式，并作简要说明。
help invINV Matrix inverse.INV(X) is the inverse of the square matrix X.A warning message is printed if X is badly scaled ornearly singular.See also SLASH, PINV, COND, CONDEST, LSQNONNEG, LSCOV. Overloaded methodshelp gf/inv.mhelp zpk/inv.mhelp tf/inv.mhelp ss/inv.mhelp lti/inv.mhelp frd/inv.mhelp sym/inv.mhelp idmodel/inv.m3、使用help命令查询函数plot的功能以及调用方法，然后利用plot命令绘制函数y=sin(x)的图形，其中0xπ≤≤。
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([43],2,3)[i,j]=ind2sub([43],10)find(Y==100)sub2ind([43],3,2)B=Y(2:2:4,3:-2:1)或B=Y([24],[31])[m n]=size(Y)---------------------------------------------------------------------2、已知矩阵A=[10-1；241；-205]，B=[0-10；213；112]求2A+B 、A 2-3B 、A*B 、B*A 、A.*B ，A/B 、A\B 解：命令为：A=[10-1；241；-205]B=[0-10；213；112]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,12.3;0,0,5.1,-2;0,-23,0,-7]b=a(:,[12])c=a([12],:)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:20s=n*s;sum=sum+s;end sum---------------------------------------------------------------------6、求a aa aaa aa a S n ++++=得值，其中a 是一个数字，由键盘输入，表达式中位数最多项a 的个数，也由键盘输入。
Matlab 上机题目及答案
1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
for m=100:999m1=fix(m/100); %求m的百位数字m2=rem(fix(m/10),10); %求m的十位数字m3=rem(m,10); %求m的个位数字if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend2.从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。
sum=0;n=0;val=input('Enter a number (end in 0):');while (val~=0)sum=sum+val;n=n+1;val=input('Enter a number (end in 0):');endif (n > 0)summean=sum/nend3. 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。
求[1,500]之间的全部完数。
for m=1:500s=0;for k=1:m/2if rem(m,k)==0s=s+k;endendif m==sdisp(m);endend4. 从键盘上输入数字星期，在屏幕上显示对应英文星期的单词。
function weekn=input('input the number:');if isempty(n)errror('please input !!')endif n>7|n<1error('n between 1 and 7')endswitch ncase 1disp('Monday')case 2disp('Tuesday')case 3disp('Wednesday')case 4disp('Thursday')case 5disp('Friday')case 6disp('Saturday')case 7disp('Sunday')end5. 某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
Matlab试题和答案
M a t l a b试题和答案(总7页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--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)求出方程组的解。
2.在[-10，10；-10，10]范围内画出函数2222 sinyx yx z++=的三维图形。
3.试画出系统321()221H s s s s =+++的零极点分布图，判断系统是否稳定，同时求其单位冲激响应和频率响应(幅频特性和相频特性)。
4. 将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。
(1)]22[)cos(ππ，，-∈x x （曲线图）； (2)4)y 2,-4x (-242),(2222≤≤≤≤+=；y x y x f （曲面图）。
matlab上机考试题及答案
matlab上机考试题及答案1. 题目：编写一个MATLAB函数，计算并返回一个向量中所有元素的平方和。
答案：函数定义如下：```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目：使用MATLAB的内置函数，找出一个矩阵中的最大元素及其位置。
答案：可以使用`max`函数来找出矩阵中的最大元素，同时使用`find`函数来获取其位置。
示例代码如下：```matlabA = [1, 2, 3; 4, 5, 6; 7, 8, 9];[maxValue, linearIndex] = max(A(:));[row, col] = ind2sub(size(A), linearIndex);```3. 题目：给定一个向量，使用MATLAB编写代码，实现向量元素的逆序排列。
答案：可以使用`flip`函数来实现向量的逆序排列。
示例代码如下：```matlabvector = [1, 2, 3, 4, 5];reversedVector = flip(vector);```4. 题目：编写一个MATLAB脚本，计算并绘制一个正弦波的图像。
答案：可以使用`sin`函数生成正弦波数据，并使用`plot`函数绘制图像。
示例代码如下：```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);xlabel('x');ylabel('sin(x)');title('Sine Wave');```5. 题目：给定一个3x3的矩阵，使用MATLAB编写代码，计算其行列式。
答案：可以使用`det`函数来计算矩阵的行列式。
示例代码如下：```matlabmatrix = [1, 2, 3; 4, 5, 6; 7, 8, 9];determinant = det(matrix);```结束语：以上是MATLAB上机考试的题目及答案，希望能够帮助大家更好地掌握MATLAB的编程技巧和函数使用。
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上机题目4
上机题目4一、单元数组和结构体1. 分别使用直接输入法和cell 函数创建 2×2 单元数组A ，第 1个元素为字符串array example ，第 2个元素为数组[1 2 3；4 5 6；7 8 9]，第三个为5+6i ，第四个元素为一维行向量[0:2:10]，并在workspace 中查看该单元数组。
2. 取出单元数组A 中每个单元的全部内容。
3. 把单元数组A 中第三个单元的内容换为字符串replacing 。
4. 分别使用直接输入法和struct 指令两种方法创建一个结构数组，名字为student 。
该数组中包括学生的姓名（name ）、学号（number ）、MATLAB 成绩（score ）。
要求最少使用两个学生的信息。
在workspace 中查看该结构体数组。
二、绘图1. 编写程序，该程序在同一窗口中绘制在 []0,2π之间的正弦曲线和余弦曲线，步长为/10π，正弦曲线设置为蓝色实线，余弦曲线颜色设置为红色虚线，两条曲线交点处，用红色星号标记。
并在生成的图形中添加如下注释：坐标轴横轴（xlabel ）标注：角度；纵轴（ylabel ）标注：曲线；图例（legend ）两条曲线分别命名为：正弦曲线、余弦曲线；标题（title ）标注：正弦曲线和余弦曲线。
2. 绘制三维曲线：226621z x xy y x y =++++-，1010x -<<，1010y -<<3. 绘制双曲抛物面：22164x y z =-，1616x -<<，44y -<<， 并在完成画图后，熟悉图形窗口工具栏（包括放大、缩小、三维旋转、光标等）和Insert 菜单栏。
比较shading flat 和shading interp 两种图形颜色修饰方法。
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上机题目2
上机题目二（1）数据类型1.建立一个类型为int8的整形变量a，并将其赋值为100；并通过class函数查看变量a的实际类型2.建立一个类型为int8的整形变量h，并将其赋值为200，察看结果，理解数据溢出3.建立一个类型为single的浮点数变量b, 并将其赋值为12.342, 并通过class函数查看变量b的实际类型4.建立一个类型为double的浮点数变量c，并将其赋值为13.342，并通过class函数查看变量c的实际类型5.实验将a与b和c分别相加，查看计算是否能够进行，并将计算结果保存到变量d中，并查看d的实际类型6.分别用笔和matlab计算表达式uint8(2^8-5)/uint8(int8(2^9))的值，理解整数除法和数据溢出的概念。
（2）基本逻辑、关系运算a=2，b=5，c=6，开展如下关系和逻辑运算，输入程序前可以先口算，然后和程序结果对比，理解关系运算和逻辑运算。
1.c>b||a<b2.c>b&&a<b3.c>b&&a>b4.c>b||a>b5.c==a&&c~=b6.c==a||c~=b（3）数组1.输入数组。
2.取出A的第二行和第三行组成一个新的数组。
3.取出A的第1-2行和第3-4列交点处的数据组成一个新的数组4.C=A([1 2], 2:3)，写出C的结果。
5.取出数组A的第四行第五列处的值。
6.把数组A变为4行9列的新的数组，并存储给变量B。
7.进行关系运算A==8，察看返回结果。
8.查找数组A中等于-1的元素的位置，分别返回这些元素的单下标索引和双下标索引。
9.查找数组A中大于4的元素的位置，并把这些元素取出来组成一个新的数组。
10.分别使用直接输入方法、冒号操作符、和linspace函数，创建一个行向量d，其元素为1, 3, 5, 7, 9,11。
并将向量d添加到数组A中，使其成为A的第一行。
Matlab考试题库及答案(教师出卷参考专用)
Matlab考试题库及答案(教师出卷参考专用)一、选择题1.以下哪个函数用于在Matlab中创建一个图形窗口？A. figureB. plotC. graphD. window答案：A2.在Matlab中，以下哪个选项可以用来定义一个矩阵？A. A = [1 2 3; 4 5 6]B. A = (1, 2, 3, 4, 5, 6)C. A = {1, 2, 3, 4, 5, 6}D. A = 1 2 3; 4 5 6答案：A3.以下哪个函数用于求解线性方程组Ax=b？A. solveB. linsolveC. solve(A, b)D. linsolve(A, b)答案：D4.在Matlab中，如何计算矩阵A和矩阵B的乘积？A. A BB. A \ BC. A . BD. A .\ B答案：A5.以下哪个函数用于在Matlab中绘制三维散点图？A. scatterB. scatter3C. plot3D. bar3答案：B二、填空题1.在Matlab中，要创建一个名为"myfig"的图形窗口，可以使用______函数。
答案：figure('Name', 'myfig')2.在Matlab中，要计算矩阵A的行列式，可以使用______函数。
答案：det(A)3.在Matlab中，若要计算变量x的平方，可以使用______运算符。
答案：.^24.在Matlab中，若要计算矩阵A的逆矩阵，可以使用______函数。
答案：inv(A)5.在Matlab中，要绘制一个正弦波形，可以使用______函数。
答案：plot(sin(x))三、判断题1.在Matlab中，矩阵的索引从1开始计数。
（）答案：正确2.在Matlab中，可以使用逻辑运算符"&&"和"||"。
（）答案：错误3.在Matlab中，矩阵乘法满足交换律。
MATLAB期末上机试题带答案
MATLAB 期末上机考试试题带答案版姓名： 学号： 成绩：1.请实现下图：50100150200250-1-0.8-0.6-0.4-0.200.20.40.60.81xyy=sin(x)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 on xlabel('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=[3 2 1;1 -1 3;2 4 -4]；b=[7;6;-2] ；x=a\b4.请实现下图：2468101214-1-0.8-0.6-0.4-0.200.20.40.60.81xy'自己的姓名'sin(x)sin(2*x)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 函数：-3-2-10123-22-55xPeaksy[x,y]=meshgrid(-3:1/8:3); z=peaks(x,y); mesh(x,y,z)surf(x,y,z)shading flataxis([-3 3 -3 3 -8 8])xlabel('x');ylabel('y');title('Peaks')7.请在x=[0,2],y=[-0.5*pi,7.5*pi],绘制光栅的振幅为0.4的三维正弦光栅。
Matlab试题和答案
Matlab试题和答案以下是为大家整理的matlab试题和答案的相关范文，本文关键词为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)求出方程组的解。
2.在[-10，10；-10，10]范围内画出函数z?sinx2?y2x?y22的三维图形。
3.试画出系统h(s)?1的零极点分布图，判断系统是否稳定，同时求其单s3?2s2?2s?1位冲激响应和频率响应(幅频特性和相频特性)。
4.将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。
完整版本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.求解方程组：x x2x21233x x x6123x3x4x4123 2.用magic产生一个5*5的矩阵，将这个矩阵的第2列与第5列互换位置。
3.编写一个M函数，求和S= 1+2+3+…+1000 4.定义一个矩阵，放入相应的数据，然后保存到MAT文件中，清除内存后，然后再调入内存。
查看数据是否改变。
5.根据下面的方程绘制图形，y=cos(x)+x，x 取值-pi-pi，然后用多项式拟合函数进行拟合。
6.创建一个5行5列的magic矩阵，删除其中的第二、三行与第三、四列。
使其变成3行3列的矩阵。
7.建一个2行2列的元胞数组，a={zeros(2,2,2),’hello’;1758,int16(1:10)},请在此基础上创建一个新的8行8列元胞数组。
8.要求编写一个函数文件，将一个给定图片的所有颜色都变深（变深多少自己斟酌）,其中黑色为0，白色为255 9.根据下面的方程绘制图形，y=cos(2x)+sin(x),然后用5次多项式进行拟合。
10.要求编写一个M函数，完成下列函数的求解：y＝sin（x）；y＝sin（x＋pi）；y＝sin（x＋pi／3）；y＝sin（x＋pi/2）自变量x的变化范围从－pi到pi，要求四个图像画在同一个图形窗体中。
11.矩阵的使用多项式在[-1 ,3 ]之间求值并画出相应的曲线。
32x4x7x10 12.元胞和结构的定义创建一个元胞数组A{1,1}=’James Bond’ ;A{1,2}=[1 2 ;3 4 ;5 6] ;A{2,1}=pi ;A{2,2}=zeros(5) ;请读取A元胞数组中的数字5。
13.M文件的编写计算分段函数y2cos(x1)x1x10xx xx10 14.函数文件的编写编写函数文件，求半径为r的圆的周长和面积。
15.图形绘制和分析2y x绘制在[－3，3]上以0.3为步长各数据点的条形图。
e 16.多项式拟合用一个三次多项式在区间[0,20]内逼近函数sin(x)。
(完整版)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基础考试题及答案
matlab基础考试题及答案1. MATLAB中，如何创建一个名为"myMatrix"的3x3矩阵，其元素为1到9的连续整数？答案：可以使用以下命令创建矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用`magic`函数：```matlabmyMatrix = magic(3);```2. 在MATLAB中，如何计算一个向量的元素平方？答案：可以使用`.^`运算符来计算向量元素的平方：```matlabv = [1 2 3 4];v_squared = v.^2;```3. MATLAB中，如何将一个向量逆序？答案：可以使用`flip`函数或者`end`索引来逆序一个向量：```matlabv = [1 2 3 4];v_reversed = flip(v);```或者```matlabv_reversed = v(end:-1:1);```4. 在MATLAB中，如何计算矩阵的行列式？答案：可以使用`det`函数来计算矩阵的行列式：```matlabA = [1 2; 3 4];detA = det(A);```5. MATLAB中，如何使用循环结构来计算1到100的自然数之和？答案：可以使用`for`循环来计算总和：```matlabsum = 0;for i = 1:100sum = sum + i;end```或者使用`sum`函数直接计算：```matlabsum = sum(1:100);```6. 在MATLAB中，如何创建一个复数？答案：可以通过指定实部和虚部来创建一个复数：```matlabcomplexNumber = 3 + 4i;```7. MATLAB中，如何绘制一个正弦波？答案：可以使用`plot`函数和`sin`函数来绘制正弦波：```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);```8. 如何在MATLAB中计算两个矩阵的点积？答案：可以使用`dot`函数或`*`运算符来计算两个矩阵的点积：```matlabA = [1 2; 3 4];B = [5 6; 7 8];dotProduct = dot(A, B);```或者```matlabdotProduct = A * B;```9. MATLAB中，如何使用条件语句来检查一个数是否为偶数？答案：可以使用`if`语句和模运算符来检查一个数是否为偶数：```matlabnum = 10;if mod(num, 2) == 0disp('The number is even.');elsedisp('The number is odd.');end```10. 在MATLAB中，如何将一个字符串转换为小写？答案：可以使用`lower`函数来将字符串转换为小写：```matlabstr = 'Hello World';lowerStr = lower(str);```结束语：以上是MATLAB基础考试题及答案，希望这些题目能帮助你更好地理解和掌握MATLAB的基础知识。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1.根据下面的方程绘制图形，y=cos(x)+x ，x 取值-pi-pi ，然后用多项式拟合函数进行拟合。
2.要求编写一个函数文件，将一个给定图片的所有颜色都变深（变深多少自己斟酌）,其中黑色为0，白色为255
3.根据下面的方程绘制图形，y=cos(2x)+sin(2x),然后用5次多项式进行拟合。
4.要求编写一个M 函数，完成下列函数的求解：y ＝sin （x ）；y ＝sin （x ＋pi ）；y ＝sin （x ＋pi ／3）；y ＝sin （x ＋pi/2）自变量x 的变化范围从－pi 到pi ，要求四个图像画在同一个图形窗体中。
5.M 文件的编写
计算分段函数 ⎪⎩⎪⎨⎧≠+=+++=10101)1cos(2x x x x x x x y
6.函数文件的编写
编写函数文件，求半径为r 的圆的周长和面积。
7.题目见教材p85，其中的数据可以给出文本文件或者由下面的式子表达
a=[10*rand(24,1)+15,1*rand(24,1)+30,2*rand(24,1)+45]，这是24行3列的一个矩阵，第一列表示温度，第二列表示湿度，第三列表示气压。
要求如下：
（1） 将温度数据取出
（2） 找出一天中最高温度
（3） 按温度进行从小到大（默认）排序
（4） 求出最高温度的时刻
（5） 按温度进行从大到小排序
8.已知Fibonacci （斐波那）数列由式子⋯⋯=+=--4,3,21k a a a k k k 可以生成，其中初值为121==a a ，试编写出生成某项Fibonacci 数值的MATLAB 函数，并求20a 的值。
9.编写M 函数，判断给定的三边能够构成一个三角形？并求其周长
10.给定一个声音文件，如1.wav ，引入这段声音，将音量改为2倍之后，再保存成2.wav
11.所表示的三维曲面，x,y的取值范围是[-8,8]
12.绘制y=sin(2x)的曲线（x从0到2π），采用指令方式在图形中：
线条颜色为红色；
线型为虚线
添加网格
添加标题“正弦曲线”
添加图例“y=sin(2x)”
添加坐标轴标签“x或y”
13.某一电力系统信号包含四种类型的信号分量，分别是：
基频正弦分量y1=10 sin(2πf1t)，其中f1=50 Hz；
3次谐波分量y2=2 sin(2πf2t+j)，其中f2=150 Hz，j=0.25π；
直流分量y3=2；
随机扰动分量。
试利用SIMULINK 建模仿真，并在同一示波器中观察四种类型信号分量及它们叠加信号的波形。
14.作出函数])6,0[(53434∈++-=x x x x y 的图形，用红十字标出其在
[0，6]之间的最小值点，并在最小值点附近标出该最小值点的坐标值。
15.编写函数，利用plot 等语句绘制y=at 2图像，其中a=[1 2 5 10],t 范围[-2，5]。