资源描述:
《北京科技大学matlab材料 2015春季学期》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、一、【实验任务】4.编写函数文件,计算,并求出当n=20时表达式的值。1.用不同的数据格式显示自然底数e的值,并分析各个数据格式之间有什么相同与不同之处。2.计算矩阵A*B,A.*B,并比较两者的区别。3.做矩阵A和B的简单运算和逻辑运算。二、【实验程序】4.1.23.一、【实验结果】4.1.结果在【实验程序】截图中分析:第一个数据格式是保留小数点后面4位有效数字;第二个数据格式是5位有效数字的科学计数法表示;第三个数据格式是15位有效数字的科学计数法表示;第四个数据格式是近似有理数表示。2.区别:第
2、一个答案是矩阵的乘法,按矩阵乘法规则,第二个答案是数组运算中的乘法,是元素对元素的运算。3.一、【实验任务】2015.5.7练习交换积分次序,分别画出两部分积分区域D1和D2。P341、π/4=1-1/3+1/5-1/7+....公式求π的近似值,直到最后一项的绝对值小于10^-6为止。3、编写函数,计算1!+2!+…+50!P791、绘制的图像,要求用蓝色的星号画图;并且画出其包络线的图像,用红色的点划线画图。3、在同一图形窗口画三个子图,要求使用指令gtext、axis、legend、title、
3、xlabel、ylabel:(1)y=xcosx,(2)y=xtan,(3)y=一、【实验程序】2015.5.7练习:D1:x=-0.5:0.1:1.5;y1=-(2*x+1).^(1/2);y2=(2*x+1).^(1/2);plot(x,y1,'g',x,y2,'g')xlabel('x轴'),ylabel('y轴')gtext('y1=-(2*x+1).^(1/2)'),gtext('y2=(2*x+1).^(1/2);')holdonplot([-0.5,-0.5],[-2,2],'g-')g
4、text('x=-0.5')plot([1.5,1.5],[-2,2],'g-')gtext('x=1.5')D2:x=-0.5:0.1:1.5;y1=-(2*x+1).^(1/2);y2=(2*x+1).^(1/2);plot(x,y1,'g',x,y2,'g')xlabel('x轴'),ylabel('y轴')gtext('y1=-(2*x+1).^(1/2)'),gtext('y2=(2*x+1).^(1/2);')holdonx=1.5:0.1:2;y3=-(16-8*x).^(1/2);y4
5、=(16-8*x).^(1/2);plot(x,y3,'r',x,y4,'r')gtext('y3=-(16-8*x).^(1/2)'),gtext('y4=(16-8*x).^(1/2)')plot([-0.5,-0.5],[-2,2],'g-')gtext('x=-0.5')plot([1.5,1.5],[-2,2],'r-')gtext('x=1.5')plot([2,2],[-2,2],'r-')gtext('x=2')P391、sum=0n=1while1/(2*n-1)>1e-6sum=s
6、um+1/(2*n-1)*((-1)^(n-1));n=n+1;endnsum*43、functions=jiecheng(n)s=1;1=1;forii=1:ns=s*ii;endP791、x=0:pi/10:4*pi;y1=exp(x/3).*sin(3*x);y2=exp(x/3);y3=-exp(x/3);plot(x,y1,'b*')holdonplot(x,y2,'r-.',x,y3,'r-.')3、x=-pi:pi/10:pi;y1=x.*cos(x);plot(x,y1,'b*-'),
7、gridonholdonx=pi:pi/100:4*pi;y2=x.*tan(1./x).*sin(x.^3);plot(x,y2,'rs-')x=1:0.1:8;y3=exp(1./x).*sin(x);plot(x,y3,'yp--')title('曲线y1=xcosx,y2=xtan(1/x)*sinx^3与y3=exp(1/x)*sinx')xlabel('x轴'),ylabel('y轴')gtext('y1=x.*cos(x)'),gtext('y2=x.*tan(1./x).*sin(x.
8、^3)')gtext('y3=exp(1./x).*sin(x)')legend('y1=x.*cos(x)','y2=x.*tan(1./x).*sin(x.^3)','y3=exp(1./x).*sin(x)')axis([-pi4*pi-44])一、【实验结果】2015.5.7练习:D1:D2:P391、>>Untitled3n=500001ans=3.14163、>>s=jiecheng(50)s=3.0414e+64P791、3、一、【实验任务