实验1 MATLAB运算基础.doc
实验1 MATLAB运算基础1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。(1) (2) ，其中 (3) ， (4) 其中2.已知,求下列表达式的值：（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,:);B23设有矩阵A和B，（1） 求他们的乘积C。（2） 求矩阵C的右下角3x2子矩阵赋给D。（3） 查看MATLAB 工作空间的使用情况。4完成下列操作：（1） 求100，999之间能被21整除的数的个数。(2)建立一个字符串向量，删除其中的大写字母。实验2 选择结构程序设计2输入一个百分制成绩，要求输出成绩等级A，B，C，D，E。其中90100分为A，8089分为B，7079分为C，6069分为D，60分以下为E。要求：（1） 分别用if语句和switch语句实现。（2） 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出错误成绩。程序：3假定某地区电话收费标准为：通话时间3分钟以下，收费0.5元；3分钟以上，则每超过分钟加收0.15元；通话时间在7：0022：00之间者，按上述标准全价收费，在其他通话时间者，按上述标准半价收费。计算某人在t1到t2通话时间，应缴多少电话费程序：price=input('请输入通话时间'); 参3章27 switch (time) case 0,1,2,3 price=0.5; case 2,3,4 price =3/100; case num2cell(5:9) rate=5/100; case num2cell(10:24) rate=8/100; case num2cell(25:49) rate=10/100; otherwise rate=14/100; end price=price*(1-rate)4建立一个5×6矩阵，要求输出矩阵第n行元素。当n值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出错误信息。程序：实验三 循环结构程序设计1.根据，求的近似值。当n分别取100、1000、1000时，结果是多少？2根据,求：（1）y<3时的最大值。（2）与（1）的n值对应的y 值3一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。试输出全部的水仙花数。要求（1）用循环结构实现。（2）用向量运算实现。4已知,求f1f100中：（1） 最大值、最小值、各数之和。（2） 正数、零、负数的个数。试验四 函数与文件1、一个自然数是素数，且它的各位数字位置经过任意对换后仍为素数，则称为绝对素数。试求所有的两位绝对素数。实验五 高层绘图操作1.画出空间曲线在范围内的图形，并画出相应的等高线。空间曲线2.11绘制参数曲线 的图形. 程序：t=0:0.1:10;x=3.*t.2;y=1./(t.2+2.*t);z=atan(t+2);plot3(x,y,z);grid on;2.13绘制函数的曲面图和等高线。（1） 程序：x=-10:0.1:10;x,y=meshgrid(x);z=(x.2-2.*x).*exp(-x.2-y.2-x.*y);mesh(x,y,z)x,y,z=peaks;contour3(x,y,z,12,'k');（2） 程序：x=-5:0.2:5;x,y=meshgrid(x);f=1./sqrt(x-1).2+y.2)-1./sqrt(x+1).2+y.2)+eps;mesh(x,y,f);contour3(x,y,f,12,'k');参4章68程序：x=-30:1:30;y=x;x,y=meshgrid(x,y)R=x.2+y.2;Z=(10*sin(sqrtR)/sqrt(1+R)mesh（x,y,z）contour3(x,y,f,12,'k');2.根据给定的参数方程，绘制下列曲面的图形。a) 椭球面作出椭球面的图形.解：令 () 程序： syms u v;u=0:0.2:2*pi;u,v=meshgrid(u);x=2.*sin(u).*cos(v);y=3.*sin(u).*sin(v);z=cos(u);mesh(x,y,z)程序：syms u v;u=0:0.2:2*pi;u,v=meshgrid(u);x=3.*cos(u).*sin(v);y=2.*cos(u).*cos(v);z=cos(u);mesh(x,y,z)b) 椭圆抛物面c) 单叶双曲面2.4作出单叶双曲面的图形解：令 ()程序：syms u v;u=-pi/2:0.2:pi/2;v=0:0.2:2*pi;u,v=meshgrid(u,v);x=sec(u).*sin(v);y=2.*sec(u).*cos(v);z=3*tan(u);mesh(x,y,z);axis(-10,10,-10,10,-10,10);view(-7,60);程序：syms u v;u=-pi/2:0.2:pi/2;v=0:0.2:2*pi;u,v=meshgrid(u,v);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u2;mesh(x,y,z);axis(-10,10,-10,10,-10,10);view(-7,60);d) 双曲抛物面e) 旋转面f) 圆锥面g)h)i) 待添加的隐藏文字内容3环面j) 正螺面3.已知，。完成下列操作：（1） 在同一坐标系下用不同颜色和线形绘制3条曲线。（2） 以子图形式绘制3条曲线（3） 分别用条形图、阶梯图、杆图和填充图绘制3条曲线。4绘制极坐标曲线，并分析参数对曲线形状的影响。5.绘制函数的曲面图和等高线。（1）（2）实验六 数据处理与多项式计算1、 将100个学生5门功课的成绩存入矩阵P中，进行如下处理：（1） 分别求5门课的最高分、最低分和相应学生的学号。（2） 分别求5门课的平均分和标准方差。（3） 5门课总分的最高分、最低分及相应学生学号。（4） 将5门课总分按从大到小的顺序存入zcj中，相应学生学号存入xsxh。2、某气象观测站测得某日6：0018：00之间每隔2小时室外温度如下表所示：时间h681012141618室外温度t118020022025.030.028.024.0室外温度t215.019.024.028.034.032.030.0试用三次样条差值分别求出该日室内外6：0017：30时之间每隔两小时各点的近似温度。3、已知lg（x）在1，101区间11个整数采样点的函数值如下表所示。x1112131415161718191101Lg（x）01.04141.32221.49141.61281.70761.78531.85131.90851.95902.0043试求lg（x）的5次拟合多项式p（x），并分别绘制出lg（x）和p（x）在1，101区间上的函数曲线。4、有三个多项式，试进行下列操作：（1）求。（2）求的根。（3）当x取矩阵A的每一元素时，求的值。其中（4）当以矩阵A为自变量时，求的值。实验七 符号计算基础和符号微积分1、已知，利用符号表达式求2、已知求：（1） B=P1*P2*A。（2） B的逆矩阵并验证结果。（3） 包括B矩阵主对角线元素的下三角阵。（4） B的行列式值。3、用符号法求下列极限和导数。（1）（2）已知，求4、用符号法求下列积分。（1）（2）实验八 级数和符号方程求解1、 级数符号求和。（1） 计算（2） 求级数之和函数，并求之和2、将ln（x）在x=1处按5次多项式展开为泰勒级数。3、求微分方程的符号解。4、求下列方程和方程组的符号解。（1）（2）
