第2次上机作业我自己留的.docxVIP

第二次上机答案 1.1做三维图形 ①解 n=100; x1=linspace(0,1,n); y1=linspace(1,2,n); [x,y]=meshgrid(x1,y1); z=cos(x)+sin(y); mesh(x1,y1,z); ②解 n=100; x1=linspace(0,1,n); y1=linspace(1,2,n); [x,y]=meshgrid(x1,y1); mesh(cos(x)+sin(y)); 得图如下 使用surf作图 n=100; x1=linspace(0,1,n); y1=linspace(1,2,n); [x,y]=meshgrid(x1,y1); z=cos(x)+sin(y); surf(z); 得图如下 使用waterfall作图 n=100; x1=linspace(0,1,n); y1=linspace(1,2,n); [x,y]=meshgrid(x1,y1); z=cos(x)+sin(y); waterfall(z); 得图如下 使用plot3作图 n=100; x1=linspace(0,1,n); y1=linspace(1,2,n); [x,y]=meshgrid(x1,y1); z=cos(x)+sin(y); plot3(x1,y1,z); 得图如下 使用meshz作图 n=100; x1=linspace(0,1,n); y1=linspace(1,2,n); [x,y]=meshgrid(x1,y1); z=cos(x)+sin(y); meshz(z); 得图如下 1.2进行数值积分 ①解 dblquad(@(x,y)(cos(x)+sin(y)),0,1,1,2) ans = 1.7979 我用的是matlab7.0，所以这个可以实现 ②解 fun=@(x,y)(cos(x)+sin(y)); dblquad(fun,0,1,1,2) ans = 1.7979 ③解 fun=cos(x)+sin(y); dblquad(@(x,y)fun,0,1,1,2) ans = 1.8444 PS：不知道答案为什么会不一样 ④解 fun=inline(cos(x)+sin(y),x,y) fun = Inline function: fun(x,y) = cos(x)+sin(y) jf=dblquad(fun,0,1,1,2) 这个是抄老师的，但是我实现不了，不知道是不是抄错了 ⑤解 jf=dblquad(inline(cos(x)+sin(y),x,y),0,1,1,2) 这个也是，实现不了 ⑥解 用int ③ f=cos(x)+sin(y) f = cos(x) + sin(y) g=int(f,x,0,1) g = sin(1) + sin(y) h=int(g,y,1,2) h = cos(1) - cos(2) + sin(1) 这个是上机的时候实现的，但是不知道为什么重新做的时候就实现不了了。。。。。。 1.3 使用MC求积分 n=1000; x=unifrnd(0,1,n,1); y=unifrnd(1,2,n,1); z=unifrnd(0,2,n,n); m=0; for i=1:n for j=1:n z1=cos(x(i))+sin(y(j)); if z(i,j)z1 m=m+1; end end end p=m/n^2 p = 0.9015 f=2*p f = 1.8030 两种方法比较得出：使用MC方法要比常规方法精确度小 2．使用MC方法得出投掷三枚骰子点数之和大于10的概率 n=10000; m=0; for i=1:n x=randperm(6); y=randperm(6); z=randperm(6); f=x(1)+y(1)+z(1); if f10 m=m+1; end end p=m/n p = 0.4986 3.使用MC方法求x^2*sin(x)cos(x)在0x2π内的最大值和最小值 n=10000; x=unifrnd(0,2*pi,n,1); y=x.^2.*sin(x).*cos(x); max(y) ans = 7.9544 min(y) ans = -15.3595 给该题函数作图 plot(x,y) set(gca,ytick,[-15,0,8]) ------在-15,0,8处分别给y画线 grid on --------------------