资源描述：
《《数学软》实验指导书.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、《数学软件》实验指导书实验一matlab运算基础实验目的：熟悉matlab 编程环境，掌握建立矩阵的方法。实验内容：1.设有矩阵A和B，（1）求他们的乘积C。（2）将矩阵C的右下角3*2子矩阵赋给DclcclearA=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25];B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3)2.完成下列操作：（1）求[100,999]之间能被21整除的数的个数。提示：先利用冒号表达式，再
2、利用find和length函数。（2）建立一个字符串向量，删除其中的大写字母。提示：利用find函数和空矩阵。%P286T4clcclearx=100:999;y1=find(rem(x,21)==0);length(y1)s='ABCabcD12f';l=s>'A'&s<'Z';s(l)=[]实验二选择结构程序设计实验目的：掌握建立和执行M文件的方法，实现选择程序设计的方法。实验内容：求下列分段的值：要求：（1）用if语句实现，分别输出x=-5.0，-3.0，1.0，2.0，2.5，3.0，5.0时的y值。提示：x的值可以从键盘输入，可以是向量。（2）用逻辑表
3、达式实现。x=input('pleaseenterxvalue:');L=(x<0&x~=-3);Y(L)=x(L).^2+x(L)-6;L=(x>=0&x~=2&x~=3&x<10);Y(L)=x(L).^2-5*x(L)+6;y=x.^2-x-1;实验三循环结构程序设计实验目的：掌握实现循环程序设计的方法，利用向量运算来代替循环操作的方法。实验内容：根据，求的近似值。当n分别取100、1000、10000时，结果是多少？要求：分别用循环结构和向量运算（使用sum函数）来实现。%P288T1clcclearn=input('Pleaseentern:');y1
4、=0;fori=1:n;y1=y1+1/i/i;y=sqrt(6*(y1));end%P288T1n=input('Pleaseenterthenumbern:');i=1:n;f=1./i.^2;y1=sum(f);y=sqrt(6*(y1))实验四函数与文件实验目的：掌握定义和调用matlab函数的方法，掌握matlab文件的基本操作。实验内容：一个自然数是素数，且它的各位数字位置经过任意对换后仍为素数，则称是绝对素数。例如13是绝对素数。试求所有两位的绝对素数。要求：定义一个判断函数的函数文件。fori=10:99a=floor(i/10);b=mod(i
5、,10);ii=b*10+a;ifsushu(i)&sushu(ii)iendendfunctionflag=sushu(i)flag=1;forj=2:i/2ifmod(i,j)==0flag=0;endend实验五高层绘图操作实验目的：掌握绘制二维、三维图形的常用函数。实验内容：1.已知，完成下列操作：i.在同一坐标系下用不同的颜色和线型绘制3条曲线。ii.以子图形式绘制3条曲线。iii.分别用条形图、阶梯图、杆图和填充图绘制3条曲线。2.绘制函数的曲面图和等高线。i.ii.%-----------------------------------------
6、--------------------------------%P290_1x=(0:pi/10:2*pi);y1=x.^2;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'r-.',x,y2,'b-',x,y3,'c:');figure;subplot(2,2,1);plot(x,y1,'r-.');subplot(2,2,2);plot(x,y2,'b-');subplot(2,2,3);plot(x,y3,'c:');figure;subplot(2,2,1);bar(x,y1,'r');title('bar');subplot(2,2
7、,2);stairs(x,y1,'g');title('stairs');subplot(2,2,3);stem(x,y1,'b');title('stem');subplot(2,2,4);fill(x,y1,'y');title('fill');%-------------------------------------------------------------------------%P290_4figure;%x=-10:0.1:10;x=-3:0.1:3;[x,y]=meshgrid(x);z=(x.^2-2*x).*exp(-x.^2-y.^2
8、-x.*y);mesh(
