matlab实验内容答案

实验报告说明：matlab 课程实验需撰写8个实验报告，每个实验报告内容写每次实验内容中标号呈黑体大号字显示的题目。

第一次实验内容：实验一 MATLAB 运算基础一、实验目的1.熟悉启动和退出MATLAB 的方法。

2.熟悉MATLAB 命令窗口的组成。

3.掌握建立矩阵的方法。

4.掌握MATLAB 各种表达式的书写规则以及常用函数的使用。

二、实验内容1.先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。

（1）22sin 8511z e ︒=+（2）12ln(2z x =，其中2120.455i +⎡⎤=⎢⎥-⎣⎦（3）0.30.33sin(0.3), 3.0, 2.9, 2.8,,2.8,2.9,3.02a ae e z a a --=+=---g L 提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。

（4）2220141122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪--≤<⎩，其中t ＝0：：提示：用逻辑表达式求分段函数值。

2.已知12344347873657A -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦，131203327B -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦求下列表达式的值：（1） A+6=B 和A-B+I(其中I 为单位矩阵)。

（2） A*B 和A.*B 。

（3） A^3和A^.3 。

（4） A/B 和B\A 。

（5）[A ，B]和[A([1,3],;);B^2] 。

3.设有矩阵A 和B12345678910111213141516171819202122232425A ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦， 30161769023497041311B ⎡⎤⎢⎥-⎢⎥⎢⎥=-⎢⎥⎢⎥⎢⎥⎣⎦（1） 求它们的乘积C 。

（2） 将矩阵C 的右下角3×2子矩阵赋给D （3） 查看MATLAB 工作空间使用情况。

4.完成下列操作：（1）求[100，999]之间能被21整除的数的个数。

提示：先利用冒号表达式，再利用find 和length 函数。

（2）建立一个字符串向量，删除其中的大写字母。

提示：利用find 函数和空矩阵。

第二次实验内容：实验三 选择结构程序设计一、实验目的1. 掌握建立和执行M 文件的方法。

2. 掌握利用if 语句实现选择结构的方法。

3. 掌握利用switch 语句实现多分支选择结构的方法。

4. 掌握try 语句的使用。

二 、实验内容1.求下列分段函数的值。

2226,0356,010,231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且其他要求：（1） 用if 语句实现，分别输出x= , , , , , , 时的y 值。

提示：x 的值从键盘输入，可以是向量。

（2） 仿照实验一第1题第4小题，用逻辑表达式实现，从而体会MATLAB 逻辑表达式的一种应用。

clearx=input('input the value of x :'); if x<0 & x~=3 y=x*x+x-6;elseif x>=0 & x<10 & x~=2 &x~=3 y=x*x-5*x+6; elsey=x*x-x-1; end y2．输入一个百分制成绩，要求输出成绩等级A,B,C,D,E 。

其中90～100分为A ，80～89分为B ，70～79分为C ，60～69分为D ，60分以下为E 。

要求：（1） 用switch 语句实现。

（2） 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出错误信息。

if 语句实现： clearx=input('input score:'); if x<0disp('Sorry,The fail score!') elseif x<60disp('E') elseif x<70disp('D') elseif x<80disp('C') elseif x<90disp('B') elseif x<=100 disp('A') elsedisp('Sorry,The fail score!') endswitch 语句实现： clearscore=input('input the score:'); switch scorecase num2cell(0:59) disp('E')case num2cell(60:69) disp('D')case num2cell(70:79) disp('C')case num2cell(80:89) disp('B')case num2cell(90:100) disp('A') otherwisedisp('Sorry!It''s the fail score,input the score again') end3. 建立5×6矩阵，要求输出矩阵第n 行元素。

当n 值超过矩阵的行数时，自动转为输出矩阵最 后一行元素，并给出出错信息。

clearA=[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,:); catch B=A(5,:); lasterr end B实验四 循环结构程序设计一、实验目的1. 掌握利用for 语句实现循环结构的方法。

2. 掌握利用while 语句实现循环结构的方法。

3. 熟悉利用向量运算来代替循环操作的方法。

二、实验内容1.根据2222211116123n π=++++L ，求π的近似值。

当n 分别取100、1 000、10 000时，结果是多少要求：分别用循环结构和向量运算（使用sum 函数）来实现。

应用sum 的方法： clearn=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); endpai=sqrt(6*y)2.根据. 11113521y n =++++-L ，求： （1）y<3时的最大n 值。

（2）于（1）的n 值对应的y 值。

cleary=0;n=1000; for i=1:1000 y=y+1/(2*i-1); if y>=3 break; end i=i+1; end i3. 已知1231231012,3n n n n f f f f f f f n ---=⎧⎪=⎪⎨=⎪⎪=-+>⎩求1100f f -中：（1） 最大值、最小值、各数之和。

（2） 正数、零、负数的个数。

-------------------------------------------------f=[1,0,1];for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3);endk=length(find(f<0))k=length(find(f>0))k=length(find(f==0))-------------------------------------------------4．若两个连续自然数的乘积减1是素数，则称这两个连续自然数是亲密数对，该素数是亲密素数。

求[2，50]区间内：（1）亲密数对的对数；（2）与上述亲密数对对应的所有亲密素数之和。

实验五函数与文件一、实验目的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+2i[a,b,c,d]=susb(x)2. 一个自然数是素数，且它的各位数字位置经过任意对换后仍为素数，则称是绝对素数。

例如13是绝对素数。

试求所有两位的绝对素数。

要求：定义一个判断素数的函数文件。

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

（2）当)1(433221)(+⨯++⨯+⨯+⨯=n n n f Λ时，求y 的值。

第三次实验内容：实验六 高层绘图操作一、实验目的1. 掌握绘制二维图形的常用函数。

2. 掌握绘制三维图形的常用函数。

二、实验内容1. 已知213),2cos(2,12y y y x y x y ⨯===完成下列操作： （1）在同一坐标系下用不同的颜色和线型绘制三条曲线。

（2）以子图形式绘制3条曲线。

（3）分别用条图形、阶梯图、杆图和填充图绘制3条曲线。

2. 分别用plot 和fplot 函数绘制函数1sin y x=的曲线，分析两曲线的差别。

3. 绘制函数的曲面图和等高线。

已知 222(2)xy xyz x x e ---=-提示：绘制三维曲面图，首先要选定一平面区域并在该区域产生网络坐标矩阵。

在做本题前，先分析并上机验证下列命令的执行结果，从中体会产生网络坐标矩阵的方法。

[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)三维曲面图中z>部分。