matlab数学实验复习题有答案

上传人:公**** 文档编号:457552929 上传时间:2023-06-01 格式:DOC 页数:12 大小:125.50KB
返回 下载 相关 举报
matlab数学实验复习题有答案_第1页
第1页 / 共12页
matlab数学实验复习题有答案_第2页
第2页 / 共12页
matlab数学实验复习题有答案_第3页
第3页 / 共12页
matlab数学实验复习题有答案_第4页
第4页 / 共12页
matlab数学实验复习题有答案_第5页
第5页 / 共12页
点击查看更多>>
资源描述

《matlab数学实验复习题有答案》由会员分享,可在线阅读,更多相关《matlab数学实验复习题有答案(12页珍藏版)》请在金锄头文库上搜索。

1、-复习题1、写出3个常用的绘图函数命令:plot、ezplot、fplot2、invA表示A的逆矩阵;3、在命令窗口健入clc,作用是去除工作间管理窗口的所有容4、在命令窗口健入clear,作用:去除存中所有变量5、在命令窗口健入figure,作用是翻开一个新的图形;6、*=-1:0.2:1表示在区间-1,1以0.2为步长等距取值7、detA表示计算A的行列式的值;8、三种插值方法:拉格朗日多项式插值,分段线性插值,三次样条插值。9、假设A=,则fliplrA=A-3=A2=trilA=triuA,-1=diagA=A(:,2),=A(3,:)=10、 normcdf1,1,2=0.5%正态分

2、布mu=1,sigma=2,*=1处的概率11、 unifpdf(5,7,2,6)=【0.25;0】11、 命令format short的作用保存小数点后四位而format long:保存小数点后14位12、 format rat的作用是最接近的有理数12、 interp1(*0,y0,*)的作用是求以*0,y0为节点数组,*为插值点数组的分段线性插值13、 13、a,b,c,d=fzerofun,*0中参数的涵义是a是变号点的近似值,b是对应,的函数值,c是停顿运行的原因(c=1即为找到该点,c=0就是没有找到)d是一个构造变量,fun是求解方程的函数M文件,*0是零点或变号点附近的值。14

3、、龙格-库塔方法可用如下MATLAB命令求解微分方程t,*=ode45(f,a,b,*0),中参数的涵义是fun是求解方程的函数M文件,a,b是输入向量即自变量的围a为初值,*0为函数的初值,t为输出指定的a,b,*为函数值15、写出以下命令的功能:a*is equal纵、横坐标轴采用等长刻度te*t1,2,y=sin*在*=1,y=2处加上字符串y=sin(*);hold on把新的plot产生的图形画在原来的图形上。titley=sin*在图形正上方加上字符串y=sin(*)16、Matlab中自定义函数M文件的第一行必须以function开头;17、二种数值积分的库函数名为:quad;q

4、uadl18、unifrnd1,2,3,4的功能是:随机生成3行4列均匀分布,每个元素服从(1,2)的矩阵19、binornd20,0.3,3,4的功能是随机生成3行4列服从(20,0.3)的二项分布的矩阵20、eigA的功能是矩阵A的特征值21、设*是一向量,则hist*的功能是作出将*十等分的直方图22、interp1(1,2,3,3,4,5,2.5)Ans=4.523、 建立一阶微分方程组的函数M文件。(做不出来)二、写出运行结果:1、eye(3,4)=2、size(1,2,3)=1;33、设b=roundunifrnd-5,5,1,4,则=3 5 2 -5*,m=min(b);*=-5

5、;m=4,*,n=sort(b)-5 2 3 54 3 1 2mean(b)=1.25,medianb=2.5,rangeb=104、向量b如上题,则any(b),all(b2),all(b5 6;7 87 8;5 6=6、假设,则7、diag(diag(B)=8、4:-2:1.*-1,6=-4 129、acos(0.5),atan(1)ans=1.6598ans=0.744810、 norm(1,2,3)Ans=3.394111、length1,3,-1=312、*=0:0.4:2;plot(*,2*,k*)13、zeros(3,1);ans=00014、ones(3)=,vander(2,

6、3,5)=16、floor1:0.3:3=1 1 1 1 2 2 218、subplot(2,2,1); fplot(sin,0,2*pi);subplot(2,2,2);plot(1,2,-1);*=linspace(0,6*pi);subplot(2,2,3);plot3(cos(*),sin(*),*);subplot(2,2,4);polar(*,5*sin(4*/3);19、 t=linespace0,2,110.0 0.20.4 0.60.8 1.01.2 1.41.6 20、 a,b=binostat(15,0.2)a=3 b=2.4y1=binopdf(5,10,0.7)=0.

7、1029,y2=binocdf(5,10,0.7)=0.150321、log10(1,10,100)=0 1 222、p=1;for k=2:3:9 p=p*k;end;p p=8023、s=0;for k=2:3:9 s=s+k;end;s s=1524、Ans=3.864725、a1=norminv(0.6,3,4)a1=4.013426、unifinv(0.4,1,5),unifpdf(0.4,1,5),unifpdf(2,1,5)Ans=2.6 0 0.2527、 A=0 1-1;2 1 0;1-1 1;0 1-11-1 1A(1,3,:)1-1 10 1-1A(3,1,:)=1 -1

8、 1 0 1 -1A(2,:)=2 1 0-2*A(1,:)= 0 -2 228、quad(sin(*),0,pi/2)=1.000029、trapz(3,4,6,1,2,3)=6.500030、 int(*-sin(*),0,1)Ans=cos(1) - 1/231、round(3:0.4:5),ceil(3:0.4:5);floor(3:0.4:5)3 3 4 4 5 53 3 3 4 4 5limit(1+1/(3*)*,inf)=1diff(sin(3*)+*3,2)=6*-9*sin(3*) taylor(e*p(3*),5,1):命令输入: y=taylor(e*p(3*),*,1

9、,Order,5)Ans=e*p(3) + 3*e*p(3)*(* - 1) + (9*e*p(3)*(* - 1)2)/2 + (9*e*p(3)*(* - 1)3)/2 + (27*e*p(3)*(* - 1)4)/8a1=mod(15,4),b1=rem(15,4)=3,3a2=mod(-15,-4),b2=rem(-15,-4)=-3,-3a3=mod(15,-4),b3=rem(15,-4)=-1,-3a4=mod(-15,4),b4=rem(-15,4)=1,-334、 *=binornd(20,0.4,2,4)8 7 10 810 7 9 12sign(*),1 1 1 11 1

10、 1 1y=-poissrnd(8,2,4)-16-108-7-7-8-6-9sign(y)-1-1-1-1-1-1-1-135、a1,b1=binostat(20,0.4)a1=8 b1=4.8a2,b2=poisstat(8)ans=8,8a3,b3=chi2stat(15)ans=15 3036、运行M文件:chi2fign=5;a=0.9;*a=chi2inv(a,n);*=0:0.1:15;y=chi2pdf(*,n);plot(*,y,b);hold on;*f=0:0.1:*a;yf=chi2pdf(*f,n);fill(*f,*a,yf,0,g);te*t(*a*1.01,0.

11、005,num2str(*a);te*t(2.5,0.05,alpha=0.9,fontsize,20);te*t(9,0.09,*chi2(4),fontsize,16);37、t=linspace(0,2*pi);polar(t,3*t,g*)38、 quadl(e*p(2*).*log(3*),1,3)ans = 398.635239、*0=0:2*pi/6:2*pi;y0=sin(*0).*cos(*0);*=linspace(0,2*pi,100);y=sin(*).*cos(*);y1=spline(*0,y0,*);*;y;y1plot(*,y,k,*,y1,b-)注:此处省略1

12、00组数据. z-40、A=round(unifrnd(0,100,3,3);L,U=lu(A)L = 0.9897 0.4699 1.0000 0.1649 1.0000 0 1.0000 0 0U = 97.0000 80.0000 92.0000 0 35.8041 26.8247 0 0 -89.656841、a=sparse(1 3 3,2 3 5,1 2 3,4,5);s=full(a)s = 0 1 0 0 0 0 0 0 0 0 0 0 2 0 3 0 0 0 0 0三、编程1、 分别用矩形公式、梯形公式、辛普森公式、Gauss-Lobatto公式及随机模拟方法计算数值积分,并

13、与符号运算计算的结果进展比较。format long*=0:0.01:pi/2;y=e*p(3*).*sin(2*);s1=sum(y)*0.01;s2=trapz(*,y);s3=quad(e*p(3*).*sin(2*),0,pi/2);s4=quadl(e*p(3*).*sin(2*),0,pi/2);n=10000;*=unifrnd(0,pi/2,1,n);y=unifrnd(0,e*p(5.5),1,n);k=0;for i=1:nif y(i)=e*p(3*(i).*sin(2*(i) k=k+1;endends5=k/n*pi/2*e*p(5.5);syms *s=int(e*p(3*).*sin(2*),0,pi/2);s6=double(s);s1,s2,s3,s4,s5,s6输出结果:ans = Columns 1 through 3 17.7868 17.6092 17.7587 Columns 4 through 6 17.7087 17.4841 17.86502、

展开阅读全文
相关资源
正为您匹配相似的精品文档
相关搜索

最新文档


当前位置:首页 > 办公文档 > 工作计划

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号