数学实验--期末复习提要VIP免费

《数学实验》复习数据处理命令与常用函数数学实验课典型问题实验报告及思考题max求数组最大值min求数组最小值mean求数组的平均数median求数组的中间数sort数组按单增排序sum数组元素求和prod数组元素求积hist统计频率fix向零方向取整zeros全零矩阵.ones全壹矩阵.eye单位矩阵rand均匀分布随机数length数组的长度det求行列式eig特征值和特征向量roots求多项式的零点.polyval计算多项式的值.数据处理命令与常用函数abs(x)绝对值sqrt(x)开平方exp(x)自然指数log(x)e为底的对数rat(x)分数表示sin正弦函数asin反正弦函数cos余弦函数acos反余弦函数tan正切函数atan反正切函数ezplot函数简单绘图fplot函数绘图.quad数值积分.dblquad数值重积分inline构造临时函数fzero求函数零点.fmin求一元函数最小值fmins求多元函数最小值diff符号求导数int符号求积分mod求余数polyfit数据多项式拟合.数据处理命令与常用函数v0=515;alpha=45*pi/180;T=2*v0*sin(alpha)/9.8;t=T*(0:16)/16;x=v0*t*cos(alpha);y=v0*t*sin(alpha)-0.5*9.8*t.^2;plot(x,y,x,y,’r*’)Xmax=x(17)例1.发射角为450的抛射曲线绘制%计算飞行时间%飞行时刻%计算航点00.511.522.53x10402000400060008000Xmax=2.7064e+004例5生肖问题——今年是鼠年还是狗年？n=input('inputn:=');S='鼠牛虎兔龙蛇马羊猴鸡狗猪';k=rem(n-4,12)+1;%取余，符号与n-4相同s=S(k);s=strcat(int2str(n),‘年是’,s,'年')inputn:=2006s=2006年是狗年Inputn:=2008s=2008年是鼠年例6天干地支问题？year=input('inputyear:=');S1='甲乙丙丁戊已庚辛壬癸';S2='子丑寅卯辰巳午未申酉戍亥';k1=mod(year-4,10)+1;%取余，符号与10相同s1=S1(k1);k2=mod(year-4,12)+1;s2=S2(k2);calendar=strcat(int2str(year),'年是',s1,s2,'年')例7程序功能如下：把边长为2以原点为中心的正方形旋转pi/24，将其压缩(r=0.89),重复24次并绘图xy=[-1-1;1-1;11;-11;-1-1];A=[cos(pi/24)-sin(pi/24);%创建正交矩阵sin(pi/24)cos(pi/24)];x=xy(:,1);y=xy(:,2);%提取坐标数据line(x,y),pause(1)%画线并暂停一秒fork=1:24%for循环，重复24次xy=.89*xy*A’;%旋转并压缩x=xy(:,1);y=xy(:,2);line(x,y),pause(1)end-1-0.500.51-1-0.500.51year=input('inputyear:=');n1=year/4;n2=year/100;n3=year/400;ifn1==fix(n1)&n2~=fix(n2)disp('是润年')elseifn1==fix(n1)&n3==fix(n3)disp('是润年')elsedisp('不是润年')end润年条件有二①能被4整除，但不能被100整除；②能被4整除，又能被400整除。例8判润年程序%fix（x）舍去小数取整；‘~=’不等于号，&逻辑与。while(B-A)>0.2iff==1tk=(B-A)/(vb+vc);elsetk=(B-A)/(vc+va);endA=A+va*tk;B=B-vb*tk;plot(A,0,'R.',B,0,'g.'),pause(1)f=-f;k=k+1;endk,tk,A,BA=0;B=100;va=10;vb=8;vc=60;f=1;k=0;plot(A,0,'ro',B,0,'go'),holdon例9摩托车问题求解的函数文件例10.计算曲线段f(x)=exp(ax)sin(bx),绕X轴旋转的旋转曲面体积symsabx%定义多个符号变量f=exp(a*x)*sin(b*x);f1=subs(f,a,-0.2);f2=subs(f1,b,0.5);V=pi*int(f2*f2,x,0,2*pi)double(V)20xV=pi*(-125/116*exp(-4/5*pi)+125/116)ans=3.1111202)]([dxxfV%subs为f，f1中参数a，b的代换；int（）为函数积分。f=inline(‘exp(-0.2*x).*sin(0.5*x)’);%构造临时函数t=(0:20)*pi/10;theta=t;r=f(t);x=t'*ones(size(t));y=r'*cos(theta);z=r'*sin(theta);mesh(x,y,z)colormap([000])axisoffview(-17,54)例11.绕X轴旋转的旋转曲面图形绘制例12.牟合方盖问题h=2*pi/100;t=0:h:pi;r=0:0.05:1;x=r‘*cos(t);y=r’*sin(t);%极坐标zz=sqrt(1-x.^2);meshz(x,y,zz)colormap([001])axisoffN=2000;P=rand(N,3);%N*3随机矩阵x=P(:,1);y=P(:,2);z=P(:,3);II=find(x.^2+y.^2<=1&x.^2+z.^2<=1);%寻找非零元素坐标m=length(II);%求向量维数V=8*m/N例13.取地球半径6400km,输入经差绘三维地球。functionearthface(Dtheta)ifnargin==0,Dtheta=15;endR=6400;theta=(-180:Dtheta:180)*pi/180;fai=(-90:Dtheta:90)*pi/180;X=R*cos(theta)'*cos(fai);Y=R*cos(theta)'*sin(fai);Z=R*sin(theta)'*ones(size...