MATLAB综合练习1 - 360文档中心

Matlab与数学建模综合练习(1)

《Matlab与数学建模》综合练习1．按顺序进行如下的操作：（1）产生一个5阶魔术方阵A；并计算A'与A-1(即inv(A))；>> A=magic(5)A =17 24 1 8 1523 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9 >> B=A'B =17 23 4 10 1124 5 6 12 181 7 13 19 258 14 20 21 215 16 22 3 9 >> inv(A)ans =-77/15600 133/2600 -23/650 3/2600 53/15600 89/2063 -97/2600 -3/650 33/2600 23/15600 -59/1950 1/325 1/325 1/325 71/1950 73/15600 -17/2600 7/650 113/2600 -577/15600 43/15600 1/200 27/650 -9/200 98/8837 （2）求A的特征值；>> eig(A)ans =65-2383/112-3846/2932383/1123846/293（3）计算A 的各列的总和与平均值；（4）计算A 的各行的总和与平均值；（5）若b=[1 2 3 4 5] '，求方程组 Ax=b 的解；（6）验证你的结论的正确性．2．产生行向量S =[1.0, 1.2, 1.4, …, 20]，并计算S * S' 与 S' * S ，你有何“发现”？3．设A=⎪⎪⎭⎫ ⎝⎛4321；B=⎪⎪⎭⎫⎝⎛5055；求C=A * B – B * A ，你有何“发现”？4．若设矩阵A=⎪⎪⎭⎫⎝⎛4321；B=⎪⎪⎭⎫⎝⎛5005；求C=A * B – B * A ，你又有何“发现”？ 5．如何建立如下的矩阵（命令方式和程序方式）？（１）1010200400020040002004⨯⎪⎪⎪⎪⎪⎭⎫⎝⎛；（２）1010010101001010100⨯⎛⎫⎪ ⎪⎪ ⎪ ⎪⎝⎭；>> a=ones(10,1)a =1 1 1 1 1 1 1 1 1 1>> A=0*a A =0 0 0 0 0 0 0>> b=diag(A)b =Columns 1 through 70 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 0Columns 8 through 100 0 00 0 00 0 00 0 00 0 00 0 00 0 00 0 00 0 00 0 0>> c=ones(10)c =Columns 1 through 71 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 1Columns 8 through 101 1 11 1 11 1 11 1 11 1 11 1 11 1 11 1 11 1 11 1 1>> C=10*cC =Columns 1 through 710 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 10Columns 8 through 1010 10 1010 10 1010 10 1010 10 1010 10 1010 10 1010 10 1010 10 1010 10 1010 10 10>> m=triu(C,1)m =Columns 1 through 70 10 10 10 10 10 100 0 10 10 10 10 100 0 0 10 10 10 100 0 0 0 10 10 100 0 0 0 0 10 100 0 0 0 0 0 100 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 0Columns 8 through 1010 10 1010 10 1010 10 1010 10 1010 10 1010 10 1010 10 100 10 100 0 100 0 0>> n=tril(C,-1)n =Columns 1 through 70 0 0 0 0 0 010 0 0 0 0 0 010 10 0 0 0 0 010 10 10 0 0 0 00 010 10 10 10 10 0 010 10 10 10 10 10 010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 10Columns 8 through 100 0 00 0 00 0 00 0 00 0 00 0 00 0 00 0 010 0 010 10 0>> K=m+n+bK =Columns 1 through 70 10 10 10 10 10 1010 0 10 10 10 10 1010 10 0 10 10 10 1010 10 10 0 10 10 1010 10 10 10 0 10 1010 10 10 10 10 0 1010 010 10 10 10 10 10 1010 10 10 10 10 10 1010 10 10 10 10 10 10Columns 8 through 1010 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 0 10 10 10 0 10 10 10 0（３）1010200411120041112004⨯⎪⎪⎪⎪⎪⎭⎫⎝⎛；>> c=ones(10,1)c =1 1 1 1 1 1 1 1 1 1>> C=1024*c C =1024102410241024102410241024102410241024>> F=diag(C)F =Columns 1 through 61024 0 0 0 0 00 1024 0 0 0 00 0 1024 0 0 00 0 0 1024 0 00 0 0 0 1024 00 0 0 0 0 10240 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 0Columns 7 through 100 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 01024 0 0 00 1024 0 00 0 1024 00 0 0 1024>> a=ones(10)a =Columns 1 through 61 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 1Columns 7 through 101 1 1 11 1 1 11 1 1 11 1 1 11 1 1 11 1 1 11 1 1 11 1 1 11 1 1 11 1 1 1>> U=triu(a,1)U =Columns 1 through 60 1 1 1 1 10 0 1 1 1 10 0 0 1 1 10 0 0 0 1 10 0 0 0 0 10 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 00 0 0 0 0 0Columns 7 through 101 1 1 11 1 1 11 1 1 11 1 1 11 1 1 11 1 1 10 1 1 10 0 1 10 0 0 10 0 0 0>> D=tril(A,-1)D =Columns 1 through 60 0 0 0 0 01 0 0 0 0 01 1 0 0 0 01 1 1 0 0 01 1 1 1 0 01 1 1 1 1 01 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 1Columns 7 through 100 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 00 0 0 01 0 0 01 1 0 01 1 1 0>> Q=D+U+FQ =Columns 1 through 61024 1 1 1 1 11 1024 1 1 1 11 1 1024 1 1 11 1 1 1024 1 11 1 1 1 1024 11 1 1 1 1 10241 1 1 1 1 11 1 1 1 1 11 1 1 1 1 11 1 1 1 1 1Columns 7 through 101 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1024 1 1 1 1 1024 1 1 1 1 1024 1 1 1 1 1024（４）1010101010101010101010⨯--⎛⎫⎪-- ⎪⎪ ⎪ ⎪--⎝⎭（５）1000120011100⎛⎫⎪⎪⎪⎪⎪⎝⎭a=10:10:100a =Columns 1 through 710 20 30 40 50 60 70Columns 8 through 1080 90 100>> A=diag(a)A =Columns 1 through 710 0 0 0 0 0 00 20 0 0 0 0 00 0 30 0 0 0 00 0 0 40 0 0 00 0 0 0 50 0 00 0 0 0 0 60 00 0 0 0 0 0 700 0 0 0 0 0 00 0 0 0 0 0 00 0 0 0 0 0 0Columns 8 through 100 0 00 0 00 0 00 0 00 0 00 0 00 0 080 0 00 90 00 0 100>> a=ones(10)a =Columns 1 through 71 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 1Columns 8 through 101 1 11 1 11 1 11 1 11 1 11 1 11 1 11 1 11 1 1>> D=tril(a,-1)D =Columns 1 through 70 0 0 0 00 01 0 0 0 00 01 1 0 0 00 01 1 1 0 00 01 1 1 1 00 01 1 1 1 10 01 1 1 1 1 1 01 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 1Columns 8 through 100 0 00 0 00 0 00 0 00 0 00 0 00 0 00 0 01 1 0>> W=A+DW =Columns 1 through 710 0 0 0 00 01 20 0 0 00 01 1 30 0 00 01 1 1 40 00 01 1 1 1 500 01 1 1 1 1 60 01 1 1 1 1 1 701 1 1 1 1 1 11 1 1 1 1 1 11 1 1 1 1 1 1Columns 8 through 100 0 00 0 00 0 00 0 00 0 00 0 00 0 080 0 01 90 01 1 100 （6）11123111113412111111220⎛⎫ ⎪⎪ ⎪ ⎪ ⎪⎪ ⎪⎪⎝⎭（７）20042000002004300000200440000020045000002004⎛⎫ ⎪ ⎪⎪ ⎪ ⎪ ⎪⎝⎭6．绘制下列曲线的图形（散点图与折线图）：]2,1[,1323-∈+--=x x x x y7．绘制下列曲面的图形：2222yx z +=（提示：曲面由两部分构成）8．在同一个图形上作下列两个函数的图象：（1）]2,0[),(sin 2π∈=x x y ；（2）]2,0[),(cos 2π∈=x x y9．假如你有一组实测数据，例如：x=[53 56 60 67.5 75 90 110]; y=[109 120.5 130 141.1 157.5 180 185]; 求其回归直线，画回归直线图形并计算最小误差平方和． 10．假如你有一组实测数据，例如：x=[75 86 95 108 112 116 135 151 155 160 163 167 171 178 185]; y=[10 12 15 17 20 22 35 41 48 50 51 54 59 66 75]; 求其回归直线，画回归直线图形并计算最小误差平方和．11．随机产生500个0到100的整数FS 作为学生的考试分数．（1）画出FS 的简单直方图；（2）画出每个分数段（0～10、10～20、…，90～100）的统计频数直方图； 12．求下列各结果：（1）用Matlab 因式分解：1001x-．（2）用Matlab 求极限：xx L xx e-+=→10)1(lim．（3）用Matlab 求积分：⎰+∞2sin x x d ．（4）用Matlab 求幂级数：∑∞=+---1212114)1(n n n n x 的和函数（化简结果）．13．非线性回归尝试说明：用线性回归方法将得到：x y 5961.00089.11+-=，但当18=x 时，2794.0-=y ，这是非常荒唐的结果！显然，一个基本要求是当0=x 时0=y ．试尝试使用非线性回归模型: bax y =.请尝试以下的命令：ezplot3('sin(t)', ' cos(t)', 't', [0,6*pi]) ezcontour('x*exp(-x^2 - y^2)') ezcontourf('x*exp(-x^2 - y^2)')ezmesh('(s-sin(s))*cos(t)','(1-cos(s))*sin(t)','s',[-2*pi,2*pi]) ezmeshc('(s-sin(s))*cos(t)','(1-cos(s))*sin(t)','s',[-2*pi,2*pi]) ezsurf('x*exp(-x^2 - y^2)') ezsurfc('x*exp(-x^2 - y^2)')。

数学实验matlab练习1

《数学实验》第一讲实验习题1.执行下列指令，观察其运算结果，理解其意义：(1)[1 2;3 4]+10-2ians =9 1011 12(2)[1 2;3 4].*[0.1 0.2;0.3 0.4]ans =0.1000 0.40000.9000 1.6000(3)[1 2;3 4].\[20 10;9 2]ans =20.0000 5.00003.0000 0.5000(4)[1 2;3 4].^2ans =1 49 16(5)exp([1 2;3 4])ans =2.7183 7.389120.0855 54.5982(6)log([1 10 100])ans =0 2.3026 4.6052(7)prod([1 2;3 4])ans =3 8(8)[a,b]=min([10 20;30 40])a =10 20b =1 1(9)abs([1 2;3 4]-pi)ans =2.1416 1.14160.1416 0.8584(10) [1 2;3 4]>=[4 3;2 1]ans =0 01 1(11) find([10 20;30 40]>=[40 30;20 10])ans =24(12) [a,b]=find([10 20;30 40]>=[40 30;20 10])a =22b =12(13) all([1 2;3 4]>1)ans =0 1(14) linspace(3,4,5)ans =3.0000 3.2500 3.5000 3.7500(15) A=[1 2;3 4];A(:,2)ans =242. 执行下列指令，观察其运算结果、变量类型和字节数，理解其意义：(1) clear;a=1,b=num2str(a),c=a>0,a= =b,a= =c,b= =cans =(2) clear;fun=’abs(x)’,x=-2,eval(fun),double(fun)ans =2ans =97 98 115 40 120 413. 本金K 以每年n 次，每次%p 的增值率(np 为每年增值额的百分比)增加，当增加到rK 时所花费的时间为(单位：年) 1001ln ln(.)T r n p += ln()ln(10.01)r T n p =+ 用MATLAB 表达式写出该公式并用下列数据计算：20512,.,r p n ===。

MATLAB复习题（答案仅供参考）1

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试题及答案

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)

1.输入20个数，求其中最大数和最小数。

要求分别用循环结构和调用MA TLAB的max函数、min函数来实现。

a=rand(1,20);max1=max(a) %用max函数求最大值min1=min(a) %用max函数求最大值%以下用循环方法求解：max2=a(1);min2=a(1);for i=2:20if a(i)>max2max2=a(i);endif a(i)<min2min2=a(i);endendmax2min2调用MATLAB的max函数、min函数来实现:>> a=rand(1,20)zuidazhi=max(a)zuixiaozhi=min(a)a =Columns 1 through 80.7513 0.2551 0.5060 0.6991 0.8909 0.9593 0.5472 0.1386 Columns 9 through 160.1493 0.2575 0.8407 0.2543 0.8143 0.2435 0.9293 0.3500 Columns 17 through 200.1966 0.2511 0.6160 0.4733zuidazhi =0.9593zuixiaozhi =0.13862. 求Fibonacci数列(1)大于4000的最小项。

(2)5000之内的项数。

(1)function t=fibonacc1f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendt=f(n);return(2)function n=fibonacc2f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendn=n-1;return3. 写出下列程序的输出结果:s=0;a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=afor j=1:4if rem(k(j),2)~=0s=s+k(j);endendendSAns=1081、解方程组Ax ＝b ，分别用求逆解法与直接解法求其解。

matlab练习题

matlab练习题Matlab是一种高级的技术计算软件，被广泛应用于科学与工程领域。

通过编写Matlab代码，我们可以解决各种数学问题、进行数据分析和可视化、进行仿真等。

本文将提供一些针对初学者的Matlab练习题，帮助读者熟悉Matlab的基本语法和常用操作。

练习一：矩阵运算在Matlab中，矩阵运算是一项重要的操作。

请编写Matlab代码完成以下任务：1. 创建两个3×3的矩阵A和B，元素的值自行选择。

2. 计算矩阵A和B的和，并将结果存储在变量C中。

3. 计算矩阵A和B的乘积，并将结果存储在变量D中。

4. 计算矩阵A的转置，并将结果存储在变量E中。

提示：使用Matlab的矩阵运算符和内置函数可以简化代码编写过程。

练习二：函数编写与图形绘制Matlab是一个功能强大的数学软件，可以用来编写自定义的函数并进行图形绘制。

1. 编写一个Matlab函数，该函数可以计算任意输入整数n的阶乘，并将结果输出。

2. 调用该函数，计算并输出5的阶乘。

3. 使用Matlab绘制函数y=2x+1的图像，其中x的取值范围为-10到10。

提示：使用Matlab的plot函数可以实现图形的绘制，并可使用grid函数添加刻度线。

练习三：数据分析Matlab提供了丰富的数据分析工具和函数，可以用来处理和分析实验数据。

1. 编写一个Matlab脚本，实现对给定实验数据进行统计分析的功能。

2. 实验数据保存在文本文件data.txt中，每行包含一个数据点。

读取数据并将其存储在一个向量中。

3. 计算数据的均值、方差和标准差，并将结果输出。

4. 使用Matlab的直方图函数绘制数据的频率分布直方图。

提示：使用Matlab的readmatrix函数可以读取文本文件中的数据。

调用mean、var和std函数可以计算均值、方差和标准差。

利用histogram函数可以绘制直方图。

练习四：数值解与绘图Matlab具有强大的数值计算和绘图功能，可以用来求解方程和绘制函数图像。

Matlab试题库1

一、填空1、在MATLAB命令窗口中的“>>”标志为MATLAB的_______提示符，“│”标志为_______提示符。

2、MATLAB的工作空间中只有三个变量v1, v2, v3，写出把它们保存到文件my_data.mat中的指令_______；3、设x是一维数组，x的倒数第3个元素表示为；设y为二维数组，要删除y的第34行和48列，可使用命令; ；4、fix(-1.5)= , round(-1.5)= .5、x为0～4pi，步长为0.1pi的向量，使用命令_______创建。

6、A=[1,2,3;4,5,6]; A(4)=__________, A(3,2)=__________________7、输入矩阵A=[1 3 2;3 -5 7;5 6 9]，使用全下标方式用_______取出元素“-5”，使用单下标方式用_______取出元素“-5”。

8、在Matlab中执行语句C=rem(25,4)的结果为。

9、Matlab的运算符分为算术运算符、关系运算符和。

10、在Matlab中圆周率π用来表示，非数值用来表示。

11、在Matlab中对数值2.3进行向∞方向取整的语句是。

12、在Matlab中命令可以在命令窗口中显示MATLAB函数或者命令的帮助信息。

13、在Matlab中__ 用于括住字符串。

14、Matlab通过数据类型把一组不同类型但同时又是在逻辑上相关的数据组成一个有机的整体，以便于管理和引用。

15、A=[1,2;3,1]；B=[1,0;0,1]；A~=B= 。

16、是Matlab的主要交互窗口，用于输入命令并显示（除图形以外）的执行结果。

17、在Matlab中引入矩阵除法的概念，有左除右除两种除法，若AX=B，则X= ，若XA=B，则X= 。

18、在Matlab语言中变量的命名应遵循如下规则：变量名必须以开头，大小写，变量名长度不超过位。

19、Matlab中Inf或inf表示、eps表示、NaN表示。

matlab习题及答案

matlab习题及答案《Matlab习题及答案：提升编程技能，解决实际问题》Matlab是一种强大的数值计算和数据分析工具，它被广泛应用于工程、科学和其他领域的计算和模拟。

为了帮助大家提升编程技能，解决实际问题，我们为大家准备了一些Matlab习题及答案，希望能够帮助大家更好地掌握Matlab的使用。

1. 习题一：编写一个Matlab程序，计算斐波那契数列的前20个数字，并将结果打印出来。

答案：下面是一个简单的Matlab程序，用于计算斐波那契数列的前20个数字。

```matlaba = 0;b = 1;fib = zeros(1, 20);fib(1) = a;fib(2) = b;for i = 3:20fib(i) = fib(i-1) + fib(i-2);enddisp(fib);```2. 习题二：编写一个Matlab程序，求解一个二次方程ax^2 + bx + c = 0的根。

答案：下面是一个简单的Matlab程序，用于求解二次方程的根。

```matlaba = 1;b = -3;c = 2;delta = b^2 - 4*a*c;if delta > 0x1 = (-b + sqrt(delta))/(2*a);x2 = (-b - sqrt(delta))/(2*a);disp(['The roots are ', num2str(x1), ' and ', num2str(x2)]);elseif delta == 0x = -b/(2*a);disp(['The root is ', num2str(x)]);elsedisp('The equation has no real roots');end```通过以上两个习题及答案的示例，我们可以看到Matlab的强大功能和灵活性。

通过练习这些习题，我们可以更好地掌握Matlab的基本语法和常用函数，从而在实际问题中更快更准确地解决数值计算和数据分析的挑战。

matlab试题1

例3-5 输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。
c=input('请输入一个字符','s'); if c>='A' & c<='Z'
disp(setstr(abs(c)+abs('a')-abs('A'))); elseif c>='a'& c<='z'
素数是大于1，且除了1和它本身以外，不能被其他任何整数所整除的整数。用筛选法求素数的基本思想是：要找出2~m之间的全部素数，首先在 2~m中划去2的倍数(不包括2)，然后划去3的倍数 (不包括3)，由于4已被划去，再找5的倍数 (不包括5)，…，直到再划去不超过的数的倍数，剩下的数都是素数。
例2-3 分别建立3×3、3×2和与矩阵A同样大小的零矩阵。
(1) 建立一个3×3零矩阵。
zeros(3) (2) 建立一个3×2零矩阵。
zeros(3,2) (3) 设A为2×3矩阵，则可以用zeros(size(A))建立一个与矩阵A同样大小零矩阵。 A=[1 2 3;4 5 6]; %产生一个2×3阶矩阵A zeros(size(A)) %产生一个与矩阵A同样大小的零矩阵
3.4 程序举例
例3-19 猜数游戏。首先由计算机产生 [1,100]之间的随机整数，然后由用户猜测所产生的随机数。根据用户猜测
的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You won”，同时退出游戏。用户最多可以猜7次。

MATLAB-考试试题-(1)汇总

MATLAB 考试试题 (1)产生一个1x10的随机矩阵，大小位于（-5 5），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏）答案：a=10*rand(1,10)-5;b=sort(a,'descend')1.请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5]2. 已知变量：A=’ilovematlab’；B=’matlab’, 请找出：（A） B在A中的位置。

（B）把B放在A后面，形成C=‘ilovematlabmatlab’3. 请修改下面的程序，让他们没有for循环语句！A=[1 2 3; 4 5 6; 7 8 9];[r c]=size(A);for i=1:1:rfor j=1:1:cif (A(i,j)>8 | A(i,j)<2)A(i,j)=0;endendend4. 请把变量A=[1 2 3; 4 5 6; 7 8 9]写到文件里(output.xls)，写完后文件看起来是这样的1 2 3 4 5 6 7 8 95.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。

6.编写M文件，从Yahoo网站批量读取60000.SH至600005.SH在2008年9月份的每日收盘价（提示：使用字符串函数）。

7. 将金牛股份（000937）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。

8.已知资产每日回报率为0.0025，标准差为0.0208，资产现在价值为0.8亿，求5％水平下资产的10天在险价值(Var)。

9.a=[1 2 3 4 5],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB 中最简单的方法计算b，注意最简单哦。

1、求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+ 2w=-3x +8z- 5w=9-6x+5y-2z+10w=-8求系数矩阵的秩；求出方程组的解。

matlab习题与答案

matlab习题与答案
MATLAB习题与答案
MATLAB是一种强大的数学软件，广泛应用于工程、科学和金融等领域。

通过MATLAB，用户可以进行数据分析、图像处理、模拟建模等多种操作。

为了帮
助大家更好地掌握MATLAB的应用，我们为大家准备了一些习题与答案，希望
能够帮助大家更好地理解和掌握MATLAB的使用。

习题一：编写一个MATLAB程序，实现对给定矩阵的转置操作。

解答：可以使用MATLAB中的transpose函数来实现矩阵的转置操作。

例如，
对于一个3x3的矩阵A，可以使用以下代码实现转置操作：
A = [1 2 3; 4 5 6; 7 8 9];
B = transpose(A);
习题二：编写一个MATLAB程序，实现对给定矩阵的逆矩阵计算。

解答：可以使用MATLAB中的inv函数来实现对矩阵的逆矩阵计算。

例如，对
于一个3x3的矩阵A，可以使用以下代码实现逆矩阵计算：
A = [1 2 3; 4 5 6; 7 8 9];
B = inv(A);
习题三：编写一个MATLAB程序，实现对给定矩阵的特征值和特征向量计算。

解答：可以使用MATLAB中的eig函数来实现对矩阵的特征值和特征向量计算。

例如，对于一个3x3的矩阵A，可以使用以下代码实现特征值和特征向量计算：A = [1 2 3; 4 5 6; 7 8 9];
[V, D] = eig(A);
通过以上习题与答案的学习，相信大家对MATLAB的应用已经有了更深入的了
解。

希望大家能够多加练习，不断提升自己的MATLAB技能，为今后的工作和学习打下坚实的基础。

matlab练习一、二答案

练习一答案1、一个三位整数各位数字的立方和等于该数本身，则称该数为水仙花数。

试编写程序找到,并输出全部水仙花数。

for n=100:999a=floor(n/100);b=floor(rem(n,100)/10);c=rem(n,10);if n==a^3+b^3+c^3s=nendend2、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6 是完数。

求[1 500]之间的全部完数。

for i=2:500s=1;for j=2:(i-1)if (rem(i,j)==0)s=s+j;endendif s==isendend3、写一个函数pifun 计算下列级数：F(n)=4*(1-1/3+1/5-1/7+...)并使用tic 和toc 指令来测量pifun(1000000)的计算时间。

n=input('请输入n值：');tic;f=0;for i=1:nf=f+4*((-1)^(i+1))/(2*i-1) ;endftoc4、输入一个百分制成绩，要求输出成绩等级为A~E，即90~100 为A，80~89 为B，70~79为C，60~69 为D，60 以下为E。

m=input('请输入一个百分制成绩：');m=floor(m/10);switch mcase {10,9}disp('A');case 8disp('B');case 7disp('C');case 6disp('D');disp('E');otherwisedisp('E');end5、输入x,y 的值，并将它们的值互换后输出。

x = input('x=');y = input('y=');disp ([x,y]);t=x;x=y;y=t;disp ([x,y]);6、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price 来表示)：price<200 没有折扣；200≤price<500 3%折扣；500≤price<1000 5%折扣；1000≤price<2500 8%折扣；2500≤price<5000 10%折扣；5000≤price 14%折扣输入所售商品的价格，求其实际销售价格。

MATLAB练习题和答案

MATLAB练习题和答案MATLAB练习题和答案控制系统仿真实验Matlab部分实验结果⽬录实验⼀MATLAB基本操作 (1)实验⼆Matlab编程 (8)实验三Matlab底层图形控制 (10)实验四控制系统古典分析 (23)实验五控制系统现代分析 (28)实验六PID控制器的设计 (35)实验七系统状态空间设计 (40)实验九直流双闭环调速系统仿真 (44)实验⼀ MATLAB 基本操作1 ⽤MATLAB 可以识别的格式输⼊下⾯两个矩阵12332357135732391894A ??=??144367823355422675342189543ii B i +??+?=+?再求出它们的乘积矩阵C ，并将C 矩阵的右下⾓2×3⼦矩阵赋给D 矩阵。

赋值完成后，调⽤相应的命令查看MATLAB ⼯作空间的占⽤情况。

A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];B=[1+4i,4,3,6,7,8;2,3,3,5,5,4+2i;2,6+7i,5,3,4,2;1,8,9,5,4,3]; C=A*B;D=C(4:5,4:6); whosName Size Bytes Class AttributesA 5x4 160 doubleB 4x6 384 doublecomplexC 5x6 480 double complexD 2x3 96 double complex2 选择合适的步距绘制出下⾯的图形(,)t∈-sin(/)t，其中111t=[-1:0.1:1];y=sin(1./t);plot(t,y)3 对下⾯给出的各个矩阵求取矩阵的⾏列式、秩、特征多项式、范数、特征根、特征向量和逆矩阵。

75350 083341009103150037193......A =??-????，5765710876810957910B ??=12345678910111213141516C=，33245518118575131D --??-?=----A=[7.5,3.5,0,0;8,33,4.1,0;0,9,103,-1.5;0,0,3.7,19.3];B=[5,7,6,5;7,10,8,7;6,8,10,9;5,7,9,10];C=[1:4;5:8;9:12;13:1rtf6];D=[3,-3,-2,4;5,-5,1,8;11,8,5,-7;5,-1,-3,-1];det(A);det(B);det(C);det(D); rank(A); rank(B); rank(C); rank(D); a=poly(A); b=poly(B); c=poly(C); d=poly(D);norm(A);norm(B);norm(C);norm(D);[v,d]=eig(A,'nobalance');[v,d]=eig(B,'nobalance');[v,d]=eig(C,'nobalance');[v,d]=eig(D,'nobalance');m=inv(A);n=inv(B);p=inv(C);q=inv(D);4 求解下⾯的线性代数⽅程，并验证得出的解真正满⾜原⽅程。

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考试试题及答案一、选择题（每题5分，共25分）1. 在MATLAB中，下列哪个命令用于创建一个行向量？A. v = [1; 2; 3]B. v = [1 2 3]C. v = [1, 2, 3]D. v = (1, 2, 3)答案：B2. 在MATLAB中，下列哪个命令用于计算矩阵A的行列式？A. det(A)B. det(A')C. det(inv(A))D. det(A^2)答案：A3. 在MATLAB中，下列哪个命令用于计算矩阵A的逆？A. inv(A)B. A^(-1)C. pinv(A)D. A\B答案：A4. 在MATLAB中，下列哪个命令用于求解线性方程组Ax= b？A. A\bB. A/BC. B/AD. A^-1b答案：A5. 在MATLAB中，下列哪个命令用于绘制二维图形？A. plot(x, y)B. scatter(x, y)C. bar(x, y)D. pie(x, y)答案：A二、填空题（每题5分，共25分）6. 在MATLAB中，可以使用______命令创建一个等差数列。

答案：linspace7. 在MATLAB中，可以使用______命令创建一个等比数列。

答案：logspace8. 在MATLAB中，可以使用______命令计算矩阵A的特征值。

答案：eig(A)9. 在MATLAB中，可以使用______命令计算矩阵A的特征向量。

答案：eigenvector(A)10. 在MATLAB中，可以使用______命令计算矩阵A的奇异值。

答案：svd(A)三、解答题（每题25分，共75分）11. 编写MATLAB程序，求解以下线性方程组：2x + 3y - z = 1x - y + 2z = 03x + 2y - 4z = -3答案：```A = [2 3 -1; 1 -1 2; 3 2 -4];b = [1; 0; -3];x = A\b;disp('解为：');disp(x);```12. 编写MATLAB程序，绘制以下函数的图形：y = sin(x) + cos(x)，x ∈ [0, 2π]答案：```x = linspace(0, 2pi, 100);y = sin(x) + cos(x);plot(x, y);title('y = sin(x) + cos(x)');xlabel('x');ylabel('y');grid on;```13. 编写MATLAB程序，计算以下矩阵的特征值和特征向量：A = [1 2 3; 4 5 6; 7 8 9]答案：```A = [1 2 3; 4 5 6; 7 8 9];[V, D] = eig(A);disp('特征值：');disp(diag(D));disp('特征向量：');disp(V);```14. 编写MATLAB程序，使用牛顿迭代法求解方程f(x) = x^3 - 4x + 2 = 0在x = 1附近的根。

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练习题

综合matlab练习题1. 计算矩阵897473535与638976242之和。

>> a=[5 3 5;3 7 4;7 9 8]; >> b=[2 4 2;6 7 9;8 3 6]; >> a+b ans =7 7 7 9 14 13 15 12 14 2. 求??+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。

>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i]; >> x’ ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i 3.计算=572396a 与??=864142b 的数组乘积。

>> a=[6 9 3;2 7 5]; >> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 404. 对于B AX =，如果=753467294A ，=282637B ，求解X 。

>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’; >> X=A\B X = -0.5118 4.0427 1.33185.已知：=987654321a ，分别计算a 的数组平⽅和矩阵平⽅，并观察其结果。

>> a=[1 2 3;4 5 6;7 8 9]; >> a.^2 ans =1 4 9 16 25 36 49 64 81 >> a^2 ans =30 36 42 66 81 96 102 126 150 6.-=463521a ，??-=263478b ，观察a 与b 之间的六种关系运算的结果。

MATLAB综合练习题汇总

一、填空题1、在MATLAB中，清除所有变量的命令为____________。

2、在MATLAB中，保存和载入数据分别使用命令____________、____________。

3、矩阵A的条件数记为___________，定义为___________。

4、在MATLAB中,对图形窗口进行灵活分割的命令为____________，对函数自适应采样的绘图函数为___________，绘制三维曲面的函数为____________、____________。

5、在MATLAB中，求矩阵特征值及对应的特征向量的命令为___________。

6、[2 1 3],[1 2 1]+A B==，A B ＝，A.*B=。

7、在MATLAB中，zeros(1,3)= ，ones(1,2)= 。

8、给定1n+个点可以唯一确定一个n阶多项式，在MATLAB 中，利用命令_________ 可以确定多项式的系数，对多项式求导的命令为____________。

9、在MATLAB中，使用命令____________可做线性插值。

10、在MATLAB中，输入命令：sum([2 1 5；9 8 5])，输出结果：ans= 。

11、在MATLAB 中,画二维图形的常用函数为____________，画三维图形常用函数为__________、___________、_______________。

12、在MATLAB 中，求矩阵A 的逆矩阵可用命令__________。

13、[7 1 5;2 5 6;3 1 5]A ，(2, 3)A ＝，A.^2= 。

14、在MATLAB 中，可以用命令计算多项式的值，用命令求多项式的根。

15、将区间[,]a b 三等分，步长为h ，各个节点的函数值为1234,,,f f f f ，则3/8辛普森法的公式为。

1/3辛普森法的公式为。

16、在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综合大作业(附详细答案)

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