matlab实验内容答案

1、实验报告说明：matlab 课程实验需撰写 8 个实验报告， 每个实验报告内容写每次 实验内容中标号呈黑体大号字显示的题目。第一次实验内容：实验一 MATLAB 运算基础一、实验目的1. 熟悉启动和退出 MATLAB 的方法。2. 熟悉 MATLAB 命令窗口的组成。3. 掌握建立矩阵的方法。4. 掌握 MATLAB 各种表达式的书写规则以及常用函数的使用。二、实验内容1.先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。1) z1 2s1ine8252)z21 ln( x 1 x2 ) ，其中212i20.4550.3a0.3 a3)z3eegsin( a0.3),

2、a3.0,2.9,2.8,L ,2.8, 2.9,3.02提示：利用冒号表达式生成 a 向量，求各点的函数值时用点乘运算。t0t14) z4 t 211t2 ，其中 t 0t22t 12t3提示：用逻辑表达式求分段函数值。2.已知12 34 41 3 1A 34 7 87 ， B2033 65 73 2 7求下列表达式的值：（1）A+6=B和 A-B+I（其中 I 为单位矩阵）。（2）A*B 和 A.*B 。（3）A3 和 A.3 。（4）A/B 和 BA 。（5）A， B和 A(1,3,;);B2 。3.设有矩阵 A 和 B1234530166789101769A 1112131415 ，

3、B023416171819209702122232425413111）求它们的乘积 C 。2）将矩阵 C的右下角 3×2 子矩阵赋给 D3）查看 MATLAB工作空间使用情况。4.完成下列操作：（1）求100，999之间能被 21 整除的数的个数。 提示：先利用冒号表达式，再利用 find 和 length 函数。（2）建立一个字符串向量，删除其中的大写字母。提示：利用 find 函数和空矩阵。第二次实验内容：实验三 选择结构程序设计一、实验目的1. 掌握建立和执行 M 文件的方法。2. 掌握利用 if 语句实现选择结构的方法。3. 掌握利用 switch 语句实现多分支选择结构的方

4、法。4. 掌握 try 语句的使用。二 、实验内容1. 求下列分段函数的值。2xx 6,x 0且 x 32y x25x 6,0 x 10,x 2且x 32 xx1其他要求：(1) 用 if 语句实现，分别输出 x= , , , , , , 时的 y值。 提示： x 的值从键盘输入，可以是向量。( 2) 仿照实验一第 1 题第 4 小题，用逻辑表达式实现，从而体会MATLAB逻辑表达式的一种应用。clearx=input('input the value of x :');if x<0 & x=3y=x*x+x-6;elseif x>=0 & x<

5、;10 & x=2 &x=3 y=x*x-5*x+6;elsey=x*x-x-1;endy2输入一个百分制成绩，要求输出成绩等级A,B,C,D,E 。其中 90100分为A，8089分为 B，7079 分为 C， 60 69 分为 D，60 分以下为 E。要求：( 1) 用 switch 语句实现。( 2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出错误信息。if 语句实现：clearx=input('input score:');if x<0disp('Sorry,The fail score!')elseif x<

6、60disp('E')elseif x<70disp('D')elseif x<80disp('C')elseif x<90disp('B')elseif x<=100disp('A')elsedisp('Sorry,The fail score!')end switch 语句实现：clear score=input('input the score:');switch score case num2cell(0:59) disp('E')ca

7、se num2cell(60:69)disp('D')case num2cell(70:79)disp('C')case num2cell(80:89)disp('B')case num2cell(90:100) disp('A') otherwise disp('Sorry!It''s the fail score,input the score again') end 3. 建立 5×6 矩阵，要求输出矩阵第 n 行元素。当 n 值超过矩阵的行数时，自动转为输出矩阵最 后一行元素，并给

8、出出错信息。clear A=1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;n=input(' 请输入输出元素行数 n=');tryB=A(n,:);catchB=A(5,:);lasterrendB实验四 循环结构程序设计一、实验目的1. 掌握利用 for 语句实现循环结构的方法。2. 掌握利用 while 语句实现循环结构的方法。3. 熟悉利用向量运算来代替循环操作的方法。、实验内容1.根据6122232 L12，求 的近似值。当 n分别取 100、1 000、10 000 时，结果是 n多少要求：分别

9、用循环结构和向量运算（使用 sum 函数）来实现。应用 sum 的方法： clear n=100;i=1:n; y=1./i.2; y1=sum(y); pai=sqrt(6*y1)第二种方法： clear y=0;n=input(' 输入数据 n： ') for i=1:n y=y+1/(i*i);end pai=sqrt(6*y)1 1 12.根据 . y 1 L ，求：3 5 2n 1(1) y<3时的最大 n 值。(2) 于( 1)的 n 值对应的 y值。clear y=0;n=1000;for i=1:1000 y=y+1/(2*i-1);if y>=3b

10、reak;endi=i+1;endi3. 已知f1 1f2 0,n 3f3 1fn fn 1 2fn 2 fn 3求 f1 f100 中：（ 1） 最大值、最小值、各数之和。（ 2） 正数、零、负数的个数。f=1,0,1;for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3); end k=length(find(f<0) k=length(find(f>0) k=length(find(f=0) 4若两个连续自然数的乘积减 1 是素数， 则称这两个连续自然数是亲密数对， 该素数是亲密素数。 求2，50区间内：(1) 亲密数对的对数；(2) 与上述亲密数对对应

11、的所有亲密素数之和。实验五 函数与文件一、实验目的1. 掌握定义和调用 ,MATLAB 函数的方法。2. 掌握 MATLAB文件的基本操作。二、实验内容1. 定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。function a,b,c,d=susb(x) a=exp(x);b=sin(x);c=log10(x);d=cos(x);调用：x=1+2ia,b,c,d=susb(x)2. 一个自然数是素数，且它的各位数字位置经过任意对换后仍为素数，则称是绝对素数。 例如 13 是绝对素数。试求所有两位的绝对素数。要求：定义一个判断素数的函数文件。functionq

12、=apnum(p); % 以下为函数文件，需单独保存for i=2:9n=find(rem(p,i)=0);p(n)=;% 找出能被 i 整除的 p 的位置% 消去能被 i 整除的 pend q=10*rem(p,10)+fix(p/10);for j=2:9% p 的个位和十位调换，调换后为 qm=find(rem(q,j)=0);q(m)=;% 找出能被 j 整除的 q 的位置% 消去能被 j 整除的 qend调用： clear all p=10:99;q=apnum(p)3. 已知 yf(40)f (30) f (20)1)当 f (n) n 10ln(n2 5)时，求 y 的值。2)当

13、 f (n) 1 2 2 3 3 4 n (n 1) 时，求 y的值。第三次实验内容：实验六 高层绘图操作一、实验目的1. 掌握绘制二维图形的常用函数。2. 掌握绘制三维图形的常用函数。二、实验内容21. 已知 y1 x ,y2 cos(2x), y3 y1 y2 完成下列操作：( 1)在同一坐标系下用不同的颜色和线型绘制三条曲线。(2) 以子图形式绘制 3 条曲线。(3) 分别用条图形、阶梯图、杆图和填充图绘制3 条曲线。12. 分别用 plot 和 fplot 函数绘制函数 y sin 的曲线，分析两曲线的差别。3. 绘制函数的曲面图和等高线。 已知 z (x2 2x)e x y xy 提

14、示：绘制三维曲面图，首先要选定一平面区域并在该区域产生网络坐标矩阵。在做本题前，先分 析并上机验证下列命令的执行结果，从中体会产生网络坐标矩阵的方法。x,y=meshgrid(-1:2,1:5)3. 已知 6<x<30， 15<y<36，求不定方程 2x+5y=126 的整数解。x=7:29; y=16:35;x,y=meshgrid(x,y); %在7,29 ×16,35区域生成网格坐标z=2*x+5*y; k=find(z=126); x(k),y(k)%找出解的位置%输出对应位置的x,y 即方程的解5. 裁掉表现函数 z=sin(y)cos(x)三维曲面

15、图中 z>部分。 x=0:2*pi;x,y=meshgrid(x);z=sin(y).*cos(x);I,J=find(z>for ii=1:length(I)z(I(ii),J(ii)=NaN;endsurf(x,y,z);实验七低层绘图操作一、实验目的1. 掌握图形对象属性的基本操作。2. 掌握利用图形对象进行绘图操作的方法。二、实验内容1. 建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在按下鼠标器 的左键之后显示出 Left Button Pressed 字样。hf=figure('Color',1,0,0,'WindowBu

16、ttonDownFcn','disp(''Left Button Pressed.'')');2 2x2. 先利用缺省属性绘制曲线 y x2e2x ，然后通过图形句柄操作来改变曲线的颜色、线型和线宽， 并利用文字对象给曲线添加文字标注 y x2e2x 。实验十二 菜单设计一、实验目的1. 了解图形用户界面的特点。2. 掌握菜单设计的方法。3. 掌握对话框设计的方法。二、实验内容1. 设计菜单。菜单条仅有 File菜单项， File下有 New 、Plot和Exit等 3个选项。选择 New时利用 Edit命令建立一 个新的 M 文件。选择

17、 Plot 将显示下一级菜单，其中有 Sine Wave和 Cosine Wave两个子菜单项，且 若选择了其中的 Sine Wave子菜单项，则将打开一个新的图形窗口并显示出正弦曲线。若选择了其 中的 Cosine Wave 子菜单项，则将打开一个新的图形窗口并显示出余弦曲线。 如果选择 Exit 菜单项， 则将关闭窗口并推出用户系统回到MATLAB命令窗口。2. 为图形窗口建立快捷菜单，用以控制窗口的背景颜色和大小。4. 设计一个对话框，其中有一个编辑框和按钮，当单击按钮时，使编辑框的内容加5。第四次实验内容：实验二 矩阵分析一、实验目的1. 掌握生成特殊矩阵的方法。2. 掌握矩阵分析的方

18、法。二、实验内容1. 产生 5阶希尔伯特矩阵 H和5阶帕斯卡矩阵 P，且求其行列式的值 Hh和 Hp以及它们的条件数 Th和 Tp，判断哪个矩阵性能更好，为什么H=hilb(5)%产生 55 阶希尔波特矩阵P=pascal(5)%产生 55 阶帕斯卡矩阵Hh=det(H);%计算行列式的值%计算条件数Ph=det(P);Th=cond(H);Tp=cond(P);HhPhThTp2. 建立一个 5×5 矩阵，求它的行列式值、迹、秩和范数。A=1 2 3 4 5;2 3 4 5 6;3 4 5 6 7;4 5 6 7 8;5 6 7 8 9;det(A)trace(A)rank(A)

19、norm(A)3. 已知29 6 18A 20 5 128 8 5求 A 的特征值及特征向量，并分析其数学意义。实验八 数据处理和多项式计算、实验目的2.3.掌握数据统计和分析的方法。 掌握数值插值与曲线拟合的方法及其应用。 掌握多项式的常用运算。、实验内容1. 将 100 个学生 5 门功课的成绩存入矩阵 P 中，进行如下处理：(1)分别求每门课的最高分、最低分及相应学生序号。 (2)分别求每门课的平均分合标准方差。(3) 5 门课总分的最高分、最低分及相应学生序号。(4) 将 5 门课总分按从大到小顺序存入 zcj 中，相应学生序号存入 xsxh。 提示：上机调试时，为避免输入学生成绩的麻

20、烦，可用取值范围在 45,95 之间的随机矩阵来表示学 生成绩。A=rand(5,100);P=100*A;% 求每门功课的最高分、序号% 求每门功课的最低分、序号% 求每门功课的平均值% 求每门功课的标准方差% 求五门功课的总成绩% 求总成绩的最高分% 求总成绩的最低分Y,U1=max(P,2);X,U2=min(P,2);YU1XU2Z1=mean(P,2)Z2=std(P,1,2)B1=cumsum(P);B=B1(5,:);YZ,UZ1=max(B,2) XZ,UZ2=min(B,2) zcj=-sort(-B,2)2. 已知 lg(x)在1,101区间 11个整数采样点的函数值如表

21、2所示。表 2 lg(x)在 10 个采样点的函数值x1112131415161718191101Lg(x)0试求 lg(x)的 5次拟合多项式 p(x),并分别绘制出 lg(x)和 p(x)在1,101区间的函数曲线。 x=1:10:101;y=0,;P,S=polyfit(x,y,5)plot(x,y,'k*',x,polyval(P ,x),'k-')4 3 223. 有 3个多项式 p1(x)x42x34x25,p2(x)x 2, p3(x)x22x 3,试进行下列操作：求 p(x) p1(x) p2(x) p3(x) 。1) 求 p(x) 的根。2)

22、当 x 取矩阵 A 的每一元素时，求 P( x)的值。其中1 1.2 1.4 A 0.75 23.5052.5(3) 当以矩阵 A 为自变量时，求 p1=1,2,4,0,5 p22=1,2;p2=0,0,0,p22; p33=1,2,3;p3=0,0,p33; p=p1+p2.*p3 x=roots(p) A=-1,;,2,;0,5,;G=poly(A)第五次实验内容：P（x）的值。其中 A 的值与（ 3）相同。% 求多项式的根% X=A 时 P（X）的值实验九 符号计算基础与符号微积分一、实验目的1. 掌握定义符号对象的方法。2. 掌握符号表达式的运算法则以及符号矩阵运算。3. 掌握求符号函

23、数极限及导数的方法。4. 掌握求符号函数定积分和不定积分的方法。、实验内容1. 已知 x=6,y=5,利用符号表达式求 zx13 x y提示：定义符号常数 x=sym( 6' ),y=sym(。 5') x=sym('6');y=sym('5');z=(x+1)/(sqrt(3+x)-sqrt(y)2. 用符号方法求下列极限或导数。1) lxim0sinx x(etanx1) 2(e 1)3sin xsyms x;f=(x*(exp(sin(x)+1)-2*(exp(tan(x)-1)/(sin(x)3; limit(f,x,0)2）已知 Atc

24、osxt3lnxdA，分别求 dA dxd2A,d2A2dt 2 dxdt%求 X 一阶导数%求 X 二阶导数syms a t x;A=ax,t3;t*cos(x),log(x); diff(A,x) diff(A,x,2)3. 用符号方法求下列积分。1)dx1 x4 x8x=sym('x'); f=1/(1+x4+x8); int(f)x=sym('x'); f=(x2+1)/(x4+1); int(f,0,inf)实验十 级数与方程符号求解一、实验目的1. 掌握级数求和的方法。2. 掌握将函数展开为泰勒级数的方法。3. 掌握微分方程符号求解的方法。4. 掌握

25、代数方程符号求解的方法。、实验内容1. 级数符号求和。10计算 Sn 1 2n 1n=sym('n');S=symsum(1/(2*n-1),n,1,10)21) 求级数 n2xn 1 之和函数，并求nn 之和。nn 1n 1%求 N的从 1 到无穷的和syms n x;S2=symsum(n2)*x(n-1),n,1,inf)函数S3=symsum(n2)/5n)2. 将 ln(x) 在 x=1 处按 5 次多项式展开为泰勒级数。x=sym('x'); f=log(x); taylor(f,x,6,1) 2 求微分方程的符号解。d2ydx2k2 y 0y(0)

26、 a y'(0) ba,b,k为任意常数syms x y a b k; m=dsolve('D2y+k2*y','y(0)=a','Dy(0)=b','x')4. 求下列方程和方程组的符号解。1) 3xex 5sin x 78.5 0x=solve('3*x*exp(x)+5*sin(x)','x')2)x2 y2 100 03x 5y 8 0x y=solve('sqrt(x2+y2)-100','3*x+5*y-8','x,y')53) l

27、n(1 x) 21 sinxx=solve('log(1+x)-5/(1+sin(x)=2','x')4) x2 9 x 1 1 0x=solve('x2+9*sqrt(x+1)-1','x')第六次实验内容： 实验十三 Simulink 的应用 一、实验目的1. 熟悉 Simulink 的操作环境并掌握构建系统模型的方法。2. 掌握 Simulink 中子系统模块的建立与封装技术。3. 对简单系统所给出的数学模型能转化为系统仿真模型并进行仿真分析。二、实验内容1. 假设从实际应用领域(力学、电学、生态或社会)中，抽象出有初始状态

28、为 0 的二阶微分方程x 0.2x 0.4 0.2u(t) ，u(t)是单位阶跃函数。 用积分器直接构造求解微分方程的模型并仿真。2. 利用传递函数模块建模。仍以下面的二阶微分方程为例 x 0.2x 0.4 0.2u(t)x 0.2x 0.4 0.2u(t)3. 利用状态方程模块建模。仍以下面的二阶微分方程为例第七次实验内容：实验十四 综合实验（ 1）MATLAB解决控制系统仿真分析与设计的方法。一、实验目的 综合运用所学知识，掌握利用二、实验内容1用 Matlab 表示传递函数为2s 3324s 3s 2s 1 的系统。2用 Matlab 表示传递函数为2(s 3)s(s 1)(s 2)的系统。3. 用 Matlab 表示状态空间表达式为 x3010x10001x20u123x31 的系统。x1y 1 0 0 x2x3A=0 1 0;0 0 1;-1 -2 -3;B=0;0;1;C=1 0 0;D=0;Gf=ss(A,B,C,D)4. 已知某控制系统的传递函数为G(s)12 s2 3s 2，求Matlab 描述的传递函数模型及零极点增益模型。num=1; den=1 3 2; Gtf=tf(num,den) Gzpk=zpk(Gtf