电子科技大学成都学院MATLAB复习题

matlab期末复习题及答案1. MATLAB基础操作在MATLAB中，如何创建一个名为"myMatrix"的3x3矩阵，其元素分别为1, 2, 3, 4, 5, 6, 7, 8, 9？答案：在MATLAB中，可以通过直接输入矩阵元素来创建矩阵。

例如，要创建一个名为"myMatrix"的3x3矩阵，可以使用以下命令：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```这样，MATLAB就会创建一个3x3的矩阵，其元素按照行顺序排列。

2. 矩阵运算给定两个矩阵A和B，其中A = [1 2; 3 4]，B = [5 6; 7 8]，计算矩阵A和B的和。

答案：在MATLAB中，可以使用加号(+)来计算两个矩阵的和。

对于给定的矩阵A和B，可以使用以下命令来计算它们的和：```matlabA = [1 2; 3 4];B = [5 6; 7 8];C = A + B;```执行上述命令后，矩阵C的结果将是：```C = [6 8;10 12];```3. 条件语句编写一个MATLAB脚本，判断一个给定的数n是否为质数，并输出相应的信息。

答案：在MATLAB中，可以使用if-else语句来判断一个数是否为质数。

以下是一个简单的脚本示例：```matlabn = input('请输入一个数：');if mod(n, 2) == 0 && n > 2disp('该数不是质数');elseif n == 1disp('1不是质数');elseisPrime = true;for i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;break;endendif isPrimedisp('该数是质数');elsedisp('该数不是质数');endend```该脚本首先接收用户输入的数n，然后通过一个for循环检查n是否有除了1和它自身以外的因数，从而判断n是否为质数。

matlab考试题库及答案大学# MATLAB考试题库及答案大学一、选择题1. MATLAB中用于生成0到1之间均匀分布随机数的函数是： - A. rand- B. randi- C. randperm- D. randn答案：A2. 下列哪个命令可以用于绘制函数图像？- A. plot- B. text- C. title- D. xlabel答案：A3. MATLAB中，以下哪个不是矩阵的属性？- A. size- B. length- C. rank- D. transpose答案：D4. 以下哪个函数可以用于求解线性方程组？- A. solve- B. linsolve- C. linprog- D. fsolve答案：A5. MATLAB中，用于计算矩阵特征值的函数是：- A. eig- B. eign- C. eigen- D. eigenvalue答案：A二、填空题1. MATLAB的基本数据单位是________。

答案：矩阵2. 使用MATLAB进行数值计算时，可以利用________来存储数据。

答案：变量3. MATLAB中的向量可以是________或________。

- 答案：行向量；列向量4. 矩阵的转置可以通过________操作实现。

答案：.'5. MATLAB中，使用________函数可以计算矩阵的行列式。

答案：det三、简答题1. 简述MATLAB中使用循环结构的注意事项。

答案：在MATLAB中使用循环结构时，应注意以下几点：确保循环体内部逻辑正确，避免无限循环；使用for循环时，循环变量的初始化和步长设置要合理；使用while循环时，循环条件要明确，确保循环能够正常退出。

2. 描述MATLAB中如何实现数组的多维索引。

答案：在MATLAB中，多维数组的索引可以通过使用圆括号来实现，每个维度的索引用逗号分隔。

例如，对于一个三维数组A，可以使用A(i,j,k)来访问第i行、第j列、第k层的元素。

大学matlab考试题及答案大学MATLAB考试题及答案一、选择题1. MATLAB的全称是什么？A. Matrix LaboratoryB. Microprocessor Application Tool SetC. Microsoft Advanced Tool SetD. Microprocessor Application Test System答案：A2. 在MATLAB中，以下哪个命令用于绘制三维图形？A. plotB. scatterC. surfD. bar答案：C3. MATLAB中用于求解线性方程组的函数是？A. solveB. linsolveC. linprogD. fsolve答案：A二、简答题1. 简述MATLAB的基本数据类型有哪些，并给出至少两个每种类型的示例。

答案：MATLAB的基本数据类型包括数值数组、字符数组和单元数组。

数值数组可以是向量、矩阵或多维数组。

例如，向量 `v = [1 2 3]`，矩阵 `M = [1 2; 3 4]`。

字符数组是由单引号或双引号括起来的字符序列，如 `C = 'Hello'`。

单元数组可以包含不同类型的数据，如`{1, 'text', [1; 2; 3]}`。

2. 描述如何在MATLAB中实现矩阵的转置和翻转。

答案：矩阵的转置可以通过 `'T'` 来实现，例如 `A'` 表示矩阵A 的转置。

矩阵的翻转可以通过 `flip` 函数实现，例如 `flip(A)` 可以翻转矩阵A的所有行和列，`flipud(A)` 仅翻转矩阵A的行，而`fliplr(A)` 仅翻转矩阵A的列。

三、编程题1. 编写一个MATLAB函数，该函数接受一个向量作为输入，并返回向量中所有元素的和以及平均值。

```matlabfunction [sumVal, avgVal] = calculateSumAndAverage(V)sumVal = sum(V);avgVal = mean(V);end```2. 设计一个MATLAB脚本来解决以下问题：给定一个3x3的矩阵，找出其中最大的元素，并显示其位置。

MATLAB复习题（答案仅供参考）1⼀、填空题1. MATLAB命令窗⼝中可⽤clear 命令清除⼯作区中的变量；⽤clc 命令清除命令窗⼝中的内容。

2. MATLAB中的运算符包括数学运算符、矩阵运算符、数组运算符。

3 MATLAB的含义为矩阵运算。

4. MATLAB中的函数类型有嵌套函数、主函数、⼦函数、私有函数和匿名函数。

5. MATLAB中⽤系统预定义的特殊变量nargin 表⽰函数输⼊参数个数。

6. size命令的功能确定数组⼤⼩。

7. MATLAB中⽤linspace 符号定义单元数组。

8. MATLAB中⽤hold on 命令打开图形保持功能。

9 MATLAB中如果⼀个变量名为123，该变量名是否合法否（填是或否）。

10. MATLAB中⽤inv 函数完成矩阵的求逆运算，⽤det 函数求矩阵的⾏列式。

11.已知C={2+5i,’Hi’;1:4,[1 2]} 输⼊s=C(3) 即s=__'_Hi__'___.12.先输⼊P=struct(‘name’,’Jack‘,’num’,43) ，再输⼊s= 即s=___Jack___.13.当x=[ 1 2 3; 4 5 6] , y=[ 5 8; 5 8 ] 时，__不能__(能﹨不能)使⽤plot 命令绘图。

14．sym 或syms _能_______(能﹨不能)创建符号矩阵。

15.已知x=3:2:24,则x有__11_个元素。

结构的开始是“if”命令，结束是_end____命令。

17.已知x为⼀个向量，计算ln（x）的MATLAB命令是计算_log(x)___。

18.命令：>>subplot(3,3,x)中有__9_个⼦图。

19.已知s=‘显⽰“hey”’，则s的元素个数是___9__。

y=f()global WW=2y=f1(5)function y=f()global Wy=w+W创建以上函数⽂件，在命令窗⼝中运⾏“y=f”命令则显⽰__出错___。

matlab考试题目及答案1. 题目：编写一个MATLAB函数，实现计算并返回一个向量中所有元素的平方和。

答案：以下是一个简单的MATLAB函数，用于计算并返回一个向量中所有元素的平方和。

```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目：给定一个3x3的矩阵A，使用MATLAB编写代码，求出矩阵A 的转置。

答案：可以通过简单的转置操作来求得矩阵A的转置。

```matlabA = [1 2 3; 4 5 6; 7 8 9];A_transpose = A';```3. 题目：编写一个MATLAB脚本，实现对一个二维数组进行排序，并返回排序后的数组。

答案：以下是一个MATLAB脚本，用于对一个二维数组进行排序，并返回排序后的数组。

```matlabfunction sortedArray = sort2DArray(array)sortedArray = sort(array(:));end```4. 题目：给定一个向量x，使用MATLAB编写代码，计算并返回向量x的元素个数。

答案：可以通过内置函数`numel`来计算向量x的元素个数。

```matlabx = [1, 2, 3, 4, 5];numElements = numel(x);```5. 题目：编写一个MATLAB函数，实现计算并返回两个向量元素的点积。

答案：以下是一个简单的MATLAB函数，用于计算两个向量的点积。

```matlabfunction dotProduct = calculateDotProduct(vector1, vector2)dotProduct = dot(vector1, vector2);end```6. 题目：给定一个矩阵B，使用MATLAB编写代码，求出矩阵B的行列式。

大学matlab考试题及答案一、选择题（每题2分，共10分）1. MATLAB中用于创建向量的命令是？A. vectorB. arrayC. linspaceD. [ ]答案：D2. 以下哪个函数用于计算矩阵的行列式？A. detB. rankC. invD. eig答案：A3. 在MATLAB中，以下哪个命令用于绘制三维曲面？A. plot3B. surfC. meshD. contour答案：B4. MATLAB中，以下哪个函数用于求解非线性方程？A. fsolveB. fzeroC. linsolveD. solve答案：A5. 如何在MATLAB中创建一个全为1的3x3矩阵？A. ones(3,3)B. zeros(3,3)C. eye(3,3)D. rand(3,3)答案：A二、填空题（每题3分，共15分）1. 在MATLAB中，使用________函数可以计算矩阵的逆。

答案：inv2. 要将向量[1, 2, 3]转置成行向量，可以使用________运算符。

答案：.'3. MATLAB中，________函数用于计算两个矩阵的点乘。

答案：dot4. 若要在MATLAB中生成一个从0到1的100个等间隔的数，可以使用________函数。

答案：linspace5. 使用________函数可以在MATLAB中创建一个随机数矩阵。

答案：rand三、简答题（每题5分，共20分）1. 描述MATLAB中矩阵的基本操作有哪些？答案：矩阵的基本操作包括矩阵的创建、转置、加法、减法、乘法、点乘、求逆、行列式计算、特征值和特征向量的计算等。

2. MATLAB中如何实现循环结构？答案：MATLAB中可以使用for循环、while循环和switch-case结构来实现循环结构。

3. 简述MATLAB中函数定义的基本语法。

答案：MATLAB中函数定义的基本语法是：function [输出参数列表] = 函数名(输入参数列表)函数体end4. MATLAB中如何实现条件判断？答案：MATLAB中可以使用if-else语句、switch-case语句来实现条件判断。

matlab试题及答案# MATLAB试题及答案一、选择题1. MATLAB的基本数据单位是：A. 矩阵B. 向量C. 标量D. 数组答案：A2. 下列哪个命令可以用来绘制函数图形？A. `plot`B. `graph`C. `draw`D. `chart`答案：A3. MATLAB中，以下哪个是正确的矩阵转置操作？A. `transpose(A)`B. `A'`C. `A^T`D. `flip(A)`答案：B二、简答题1. 简述MATLAB中矩阵的基本操作。

答案：在MATLAB中，矩阵是最基本的数据结构，可以进行加、减、乘、除等基本运算。

矩阵的创建可以使用方括号`[]`，例如`A = [1 2;3 4]`。

矩阵的转置使用单引号`'`，例如`A'`。

矩阵的求逆使用`inv`函数，例如`inv(A)`。

2. MATLAB中如何实现循环结构？答案：MATLAB中实现循环结构主要有两种方式：`for`循环和`while`循环。

`for`循环用于已知迭代次数的情况，例如：```matlabfor i = 1:5disp(i);end````while`循环用于迭代次数未知的情况，例如：```matlabi = 1;while i <= 5disp(i);i = i + 1;end```三、计算题1. 给定矩阵A和B，请计算它们的乘积C，并求C的行列式。

A = [1 2; 3 4]B = [5 6; 7 8]答案：首先计算矩阵乘积C：```matlabC = A * B;```然后计算C的行列式：```matlabdetC = det(C);```结果为：```matlabC = [19 22; 43 50]detC = -16```2. 编写一个MATLAB函数，计算并返回一个向量的范数。

答案：```matlabfunction norm_value = vector_norm(v)norm_value = norm(v);end```四、编程题1. 编写一个MATLAB脚本，实现以下功能：- 随机生成一个3x3的矩阵。

matlab编程考试题及答案1. 编写一个MATLAB函数，该函数接受一个向量作为输入，并返回该向量中所有元素的平方和。

```matlabfunction sumOfSquares = vectorSquareSum(inputVector)sumOfSquares = sum(inputVector .^ 2);end```答案：该函数通过使用点乘运算符（`.^`）来计算向量中每个元素的平方，然后使用`sum`函数计算这些平方值的总和。

2. 给定一个3x3的矩阵A，编写MATLAB代码计算其行列式。

```matlabA = [1 2 3; 4 5 6; 7 8 9];determinantA = det(A);```答案：使用MATLAB内置函数`det`可以直接计算矩阵A的行列式。

3. 编写一个MATLAB脚本，实现对一个给定字符串进行反转。

```matlaboriginalString = 'HelloWorld';reversedString = fliplr(originalString);```答案：使用`fliplr`函数可以将字符串中的字符从左到右翻转，实现字符串的反转。

4. 给定两个向量x和y，编写MATLAB代码计算它们之间的欧几里得距离。

```matlabx = [1 2 3];y = [4 5 6];euclideanDistance = sqrt(sum((x - y) .^ 2));```答案：欧几里得距离可以通过计算两个向量对应元素差的平方和的平方根来得到。

5. 编写一个MATLAB函数，该函数接受两个参数，一个是矩阵，另一个是标量值，返回矩阵中所有元素与该标量值的差的绝对值。

```matlabfunction absDiffMatrix = matrixScalarDifference(matrix, scalar)absDiffMatrix = abs(matrix - scalar);end```答案：该函数通过从矩阵的每个元素中减去标量值，然后使用`abs`函数计算结果的绝对值。

matlab期末考试题目及答案1. 题目：编写一个MATLAB函数，实现矩阵的转置操作。

答案：可以使用`transpose`函数或`.'`操作符来实现矩阵的转置。

例如，对于一个矩阵`A`，其转置可以通过`A'`或`transpose(A)`来获得。

2. 题目：使用MATLAB求解线性方程组Ax=b，其中A是一个3x3的矩阵，b是一个3x1的向量。

答案：可以使用MATLAB内置的`\`操作符来求解线性方程组。

例如，如果`A`和`b`已经定义，求解方程组的代码为`x = A\b`。

3. 题目：编写MATLAB代码，计算并绘制函数f(x) = sin(x)在区间[0, 2π]上的图像。

答案：首先定义x的范围，然后计算对应的函数值，并使用`plot`函数绘制图像。

代码示例如下：```matlabx = linspace(0, 2*pi, 100); % 定义x的范围y = sin(x); % 计算函数值plot(x, y); % 绘制图像xlabel('x'); % x轴标签ylabel('sin(x)'); % y轴标签title('Plot of sin(x)'); % 图像标题```4. 题目：使用MATLAB编写一个脚本，实现对一个给定的二维数组进行排序，并输出排序后的结果。

答案：可以使用`sort`函数对数组进行排序。

如果需要对整个数组进行排序，可以使用`sort`函数的两个输出参数来获取排序后的索引和值。

代码示例如下：```matlabA = [3, 1, 4; 1, 5, 9; 2, 6, 5]; % 给定的二维数组[sortedValues, sortedIndices] = sort(A(:)); % 对数组进行排序sortedMatrix = reshape(sortedValues, size(A)); % 将排序后的值重新构造成矩阵disp(sortedMatrix); % 显示排序后的结果```5. 题目：编写MATLAB代码，实现对一个字符串进行加密，加密规则为将每个字符的ASCII码值增加3。

第二章1、x=[2,4];y=x八3+(x-0.98)八2./(x+1.35)八3-5*(x+1./x) 2、y=cos(pi/3)-(9-sqrt(2))A(1/3) 3、a=3;A=4; b=a.A2; B=b.A2-1; c=a+A-2*B;C=a+2*B+c4、x=[1 2 3;4 5 6;7 8 9]Desktop->Workspace, 双击变量x 5、clearx=magic(3) y=randn(3,3) xy=[x,y] yx=[x,y] z=xy(:,1:2) 6、clear x=eye(4,4)y=triu(x)7、clear x=rand(4,5) y=x([1,2],:)z=(y>=0.3).*y8、clear x=randn(5,5) y=inv(x)9、clear x=randn(5,5)z=x A510、clearA=[1,4,8;-3,6,-5;2,-7,-12];B=[5,4,3;6,-2,3;-1,3,-9];C=A*BD=A.*B11、clear x=linspace(0,2*pi,125);y=cos(x).*(0.5+3*sin(x)./(1+x.A2)); plot(x,y)12、z=-45:1:45;x=z.*sin(3*z);y=z.*cos(3*z);plot3(x,y,z);13、x=-2:0.1:2,y=x;[x,y]=meshgrid(x,y);z=x.A2.*exp(-x.A2-y.A2)surf(x,y,z);14、x=-2:0.1:2,y=x;[x,y]=meshgrid(x,y);z=x.A2.*exp(-x.A2-y.A2);surf(x,y,z);z1=0.05*x-0.05*y+0.1;hold on,mesh(x,y,z1);15、(1)n=2;alfa=0;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,1);plot(t,x,t,y);(2)n=2;alfa=0;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,1);plot(t,x,t,y);n=2;alfa=pi/3;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa); subplot(2,2,2);plot(t,x,t,y);n=2;alfa=pi/2;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa);subplot(2,2,3);plot(t,x,t,y);n=2;alfa=pi;t=0:0.1:10;x=cos(t);y=sin(n*t+alfa);subplot(2,2,4);plot(t,x,t,y);球体clearfor i=1:100for j=1:100x(i,j)=i/100*cos(j*2*pi/100);y(i,j)=i/100*sin(j*2*pi/100);z(i,j)=sqrt(1.001-x(i,jF2-y(i,j)八2);endendsurf(x,y,z);hold on;surf(x,y,-z);axisequal棱锥clearfor i=1:100for j=1:100x(i,j)=i/100*cos(j*2*pi/4);y(i,j)=i/100*sin(j*2*pi/4);z(i,j)=sqrt(x(i,j)A2+y(i,j)A2);endendsurf(x,y,-z);求最大值x=[67 87 56 99 43] max=0;for i=1:5 if max>=x(i) max=max; elsemax=x(i) endendmax求最小值x=[67 87 56 99 43] min=inf;for i=1:5if min<=x(i) min=min;else min=x(i) end end min求和x=[67 87 56 99 43]; sum=0; for i=1:5 sum=sum+x(i) end sum第三章1、h0=[446,714,950,1422,1634];t0=[7.04,4.28,3.40,2.54,2.13];t1=interp1(h0,t0,500,'linear') 2、x0=[1,0,-1];y0=[0,1,0]; p=polyfit(x0,y0,3); x=-1:0.1:1; y=polyval(p,x);plot(x,y,-x,-y),axis equal3、clear x0=[1.0,1.1,1.2,1.3,1.4];y0=[0.25,0.2268,0.2066,0.1890,0.1736];p=polyfit(x0,y0,3);p1=polyder(p); y=polyval(p1,[1.0,1.2])4、p=[3 4 7 2 9 12]; r=roots(p)5、r=[-3 -5 -8 -9]; p=poly(r)6、function ypie=fun1(x,y)ypie=x42./y-x.*cos(y);[x,y]=ode45( 'fun1' ,[0,5],1)plot(x,y)7、function y=fun2(x)y=x.A4-3*x A3+5*cos(x)+8;fplot( 'fun2' ,[1,5])hold on;x=0:5;y=0*x;plot(x,y)a=fzero(fun2,2)b=fzero(fun2,3)8、a=[2 4 9;4 2 4;9 4 18];[v,gama]=eig(a)9、function xpie=fun3(t,x)f=exp(-t)xpie=[0 1 0 0;1 0 0 -1;0 0 0 1;0 -1 -10]*x+[0;1;0;0]*f[t,x]=ode45( 'fun3' ,[0,5],[0 0 0 0]'); plot(t,x(:,1),'b' ,t,x(:,3), 'r' )10、function f=fun4(x,y)f=4*(x-y)-x八2-y八2;v=-2:0.2:2;[x,y]=meshgrid(v)z=fun4(x,y);surf(x,y,z)求极小值 (有误) function f=fun5(x,y)x=v(1);y=v(2)f=-(4*(x-y)-x.A2-y.A2);v=-50:2:50;[x,y]=meshgrid(v)z=fun4(x,y);surf(x,y,z)min=fminsearch('fun4',[0,0])11、function f=fun6(x,y)f=sin(y)+exp(x)-x.*y.A2ezplot( 'fun6' )12、function f=fun7(x,y)f=(x-y).A2.*(sin(x+y)).A2s=dblquad( 'fun7' ,pi,2*pi,0,pi)14、a=randn(4,4)[l,u]=lu(a)a=randn(4,4)[u,gama,v]=svd(a)15、b=[3 0.5 4];a=[3 4 5 6]; [r,p,k]=residue(b,a)b=[3 0.5 4];a=[3 4 5 6];[r,p,k]=residue(b,a)t=0:0.1:100;y=0*t;for i=1:3 y=y+r(i)*exp(p(i)*t) end plot(t,y) b=[3 0.5 4];a=[3 4 5 6]; impulse(b,a) b=[3 0.5 4];a=[3 4 5 6]; step(b,a)第四章1、n=1:11;x=cos(pi/6*n); subplot(2,1,1);stem(x); y=abs(fft(x)); subplot(2,1,2);stem(y);clear n=0:11;x=cos(pi/6*n); subplot(2,1,1);stem(x); y=abs(fft(x));subplot(2,1,2);stem(y);2、(有误)n=0:19;x=5*0.6.八n;subplot(3,1,1);stem(n,x);for i=-20:39xl(i+21)=x(mod(i+40,20)+1); end n1=-20:39; subplot(3,1,2);stem(n1,x1);x2=x1(10:29); subplot(3,1,3);stem(n,x2)新方法m=10; e=0:19;c=0.6;k=5;a=k*c.A e;a=a';b=circshift(a,m);L=length(a)-1;n=0:L;subplot(2,1,1);stem(n,a);axis([0,L,min(a),max(a)]);subplot(2,1,2);stem(n,b);axis([0,l,min(a),max(a)])3、x=[0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0 0 0 0 0 0 0 0];h=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]y=conv(x,h);stem(y)4、x=[0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0.8 0 0 0 0 0 0 0 0];h=[1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0]y=conv(x,h);subplot(2,1,1);stem(y)xx=fft(x);hh=fft(h);yy=ifft(xx.*hh);subplot(2,1,2);stem(yy)5、b=[2 3];a=[1 0.4 1];[z p k]=tf2zp(b,a) 6、b=[4 15.6 6 2.4 -6.4];a=[3 2.4 6.3 -11.4 6];[z p k]=tf2zp(b,a)7、b=[18];a=[18 3 -4 -1];[r p k]=residue(b,a)8、b=[0.2 0.3 1];a=[1 0.4 1]; freqs(b,a)9、Fs=2000;wc=[100 200]/(Fs/2);N=10;[b,a]=butter(N,wc, 'z' ); %freqz(b,a)dimpulse(b,a,100)10、wp=100;ws=200;rp=2;rs=15;Fs=500;wp=wp/(Fs/2);ws=ws/(Fs/2);[N,wc]=buttord(wp,ws,rp,rs); [b,a]=butter(N,wc, 'z' ) freqz(b,a)b=fir1(48,[0.35 0.65]);a=1; freqz(b,a)12、b=fir1(37,0.3);a=1; freqz(b,a)13、低通：F=[0:1/56:1];A=[ones(1,50),zeros(1,56-49)]; b=fir2(56,F,A);a=1; freqz(b,a)11、带通：F=[0:1/56:1];A=[ones(1,25),zeros(1,25),ones(1,7)];b=fir2(56,F,A);a=1;freqz(b,a)高通F=[0:1/56:1];A=[zeros(1,50),ones(1,56-49)]; b=fir2(56,F,A);a=1; freqz(b,a)%滤波器设计(椭圆滤波器)wp1=300;wp2=700;ws1=301;ws2=699; rp=0.1;rs=50;Fs=2000; wp=[wp1]/(Fs/2);ws=[ws1]/(Fs/2);[N,wc]=ellipord(wp,ws,rp,rs,'z')[b,a]=ellip(N,rp,rs,wc,'high')%滤波器特性分析freqz(b,a)%信号采样及时频变换Fs=2000;t=0:1/Fs:10;yt=sin(2*pi*200*t)+sin(2*pi*500*t)+sin(2*pi*800*t); figure(2);subplot(3,1,1);plot(t,yt)figure(2);subplot(3,1,2);plot(t,yt);axis([0,0.1,-5,5]) yf=abs(fft(yt));f=(1:length(yf))/length(yf)*Fs;figure(2);subplot(3,1,3);plot(f,yf)%信号滤波及时频分析yt1=filter(b,a,yt);figure(3);subplot(3,1,1);plot(t,yt1)figure(3);subplot(3,1,2);plot(t,yt1);axis([0,0.1,-5,5]) yf1=abs(fft(yt1));f=(1:length(yf1))/length(yf1)*Fs; figure(3);subplot(3,1,3);plot(f,yf1)。