《数学实验期末复习提要》由会员分享，可在线阅读，更多相关《数学实验期末复习提要（23页珍藏版）》请在金锄头文库上搜索。
1、数学实验复习数据处理命令与常用函数数学实验课典型问题实验报告及思考题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
2、) 分数表示 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 *c
3、os(alpha); y=v0 *t *sin(alpha) -0.5*9.8*t.2; plot(x,y,x,y,r*) Xmax=x(17)例1. 发射角为450的抛射曲线绘制%计算飞行时间 %飞行时刻 %计算航点Xmax=2.7064e+004例5 生肖问题今年是鼠年还是狗年？ n=input(input n: =); S=鼠牛虎兔龙蛇马羊猴鸡狗猪; k=rem(n-4,12)+1; %取余，符号与n-4相同 s=S(k); s=strcat(int2str(n),年是, s,年)input n:=2006s= 2006年是狗年 Input n:=2008s=2008年是鼠年例6 天干地
4、支问题？year=input(input year:=); 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;1 1;-1 1;-1 -1; A=cos(pi/24) -sin(pi/24); %创建正交矩阵si
5、n(pi/24) cos(pi/24); x=xy(:,1);y=xy(:,2); %提取坐标数据 line(x,y),pause(1) %画线并暂停一秒 for k=1:24 %for循环，重复24次 xy=.89*xy*A; %旋转并压缩x=xy(:,1);y=xy(:,2); line(x,y),pause(1) endyear=input(input year:=); n1=year/4; n2=year/100; n3=year/400; if n1=fix(n1)elsetk=(B-A)/(vc+va);endA=A+va*tk;B=B-vb*tk;plot(A,0,R.,B,0,
6、g.),pause(1)f=-f;k=k+1; end k,tk,A,BA=0; B=100; va=10;vb=8;vc=60; f=1;k=0; plot(A,0,ro,B,0,go),hold on例9 摩托车问题求解的函数文件例10. 计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积syms a b x %定义多个符号变量 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) V =pi*(-125/116*exp(-4
7、/5*pi)+125/116) ans = 3.1111%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(0 0 0) axis off view(-17,54)例11. 绕X轴旋转的旋转曲面图形绘制例12.牟合方盖问题 h=2*pi/100; t=0:h:pi; r=0:0.05:1; x=r*c
8、os(t); y=r*sin(t); %极坐标 zz=sqrt(1-x.2); meshz(x,y,zz) colormap(0 0 1) axis offN=2000;P=rand(N,3); %N*3随机矩阵 x=P(:,1);y=P(:,2);z=P(:,3); II=find(x.2+y.2 z = 0.05。即求逆累积函数在 x=0.95 处的值 mu=170;sigam=6; z=norminv(0.95,mu,sigam) data=mu+sigam*randn(1000,1); II=find(data=z); F=length(II)/10000 x=150:.1:190; y=normpdf(x,mu,sigam); plot(x,y,k) xx=z:.1:190; yy =normpdf(xx,mu,sigam); fill(z,xx,190,0,yy,0,r)z = 179.8691F = 0.0047
