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

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

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

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

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

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

（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 a e ez a a --=+=---提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。

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

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 ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦， 3161769023497041311B ⎡⎤⎢⎥-⎢⎥⎢⎥=-⎢⎥⎢⎥⎢⎥⎣⎦ （1） 求它们的乘积C 。

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

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 xx x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且其他要求： （1） 用if 语句实现，分别输出x= -5.0, -3.0, 1.0, 2.0, 2.5, 3.0, 5.0 时的y 值。

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

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

clearx=input('input the value of x :');if x<0 & x~=3y=x*x+x-6;elseif x>=0 & x<10 & x~=2 &x~=3y=x*x-5*x+6;elsey=x*x-x-1;endy2．输入一个百分制成绩，要求输出成绩等级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<=100disp('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,:);catchB=A(5,:);lasterrendB实验四 循环结构程序设计一、实验目的1. 掌握利用for 语句实现循环结构的方法。

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

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

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

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

应用sum 的方法：clearn=100;i=1:n;y=1./i.^2;y1=sum(y);pai=sqrt(6*y1)第二种方法：cleary=0;n=input('输入数据n ：')for i=1:ny=y+1/(i*i);endpai=sqrt(6*y)2.根据. 11113521y n =++++- ，求：（1）y<3时的最大n 值。

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

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

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

-------------------------------------------------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.掌握定义和调用,MA TLAB函数的方法。

2.掌握MA TLAB文件的基本操作。

二、实验内容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 整除的qend调用：clear allp=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 函数绘制函数1siny x =的曲线，分析两曲线的差别。

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

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

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

[x,y]=meshgrid(-1:0.5: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>0.25部分。