《北京科技大学matlab材料 2015春季学期》由会员分享，可在线阅读，更多相关《北京科技大学matlab材料 2015春季学期（41页珍藏版）》请在金锄头文库上搜索。
1、一、 【实验任务】4.编写函数文件，计算 ，并求出当 n=20 时表达式的值。nk1!1.用不同的数据格式显示自然底数 e 的值，并分析各个数据格式之间有什么相同与不同之处。2.计算矩阵 A*B，A.*B，并比较两者的区别。3.做矩阵 A 和 B 的简单运算和逻辑运算。二、 【实验程序】4.1.23.三、 【实验结果】4.1.结果在【实验程序】截图中分析：第一个数据格式是保留小数点后面 4 位有效数字；第二个数据格式是 5 位有效数字的科学计数法表示；第三个数据格式是 15 位有效数字的科学计数法表示；第四个数据格式是近似有理数表示。2.区别：第一个答案是矩阵的乘法，按矩阵乘法规则，第二个答案
2、是数组运算中的乘法，是元素对元素的运算。3.四、 【实验任务】2015.5.7 练习交换 积分次序，分别画出两部分32x12168x2 31df(,y)df(,y)d 积分区域 D1 和 D2。P34 1、/4=1-1/3+1/5-1/7+.公式求 的近似值，直到最后一项的绝对值小于 10-6 为止。3、编写函数，计算 1！+2！+50 ！P791、绘制 的图像，要求用蓝色的星号画图；并)3sin(xey）（ 4,0且画出其包络线 的图像，用红色的点划线画图。3x3、在同一图形窗口画三个子图，要求使用指令gtext、 axis、legend、 title、xlabel、ylabel:（1） y
3、=xcosx, ),(x（2） y=xtan ,13sin4（3） y= 8,1si1xex五、 【实验程序】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);)hold onplot(-0.5,-0.5,-2,2,g-)gtext(x=-0.5)plot(1.5,1.5,-2,2,g-)gtext(x=1.5)D2：x=-0.5:0.1:
4、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);)hold onx=1.5:0.1:2;y3=-(16-8*x).(1/2);y4=(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
5、,-2,2,r-)gtext(x=1.5)plot(2,2,-2,2,r-)gtext(x=2)P391、sum=0n=1while 1/(2*n-1)1e-6sum=sum+1/(2*n-1)*(-1)(n-1);n=n+1;endnsum*43、function s=jiecheng(n)s=1;1=1;for ii=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*)hold onplot(x,y2,r-.,x,y3,r-.)3、x=-pi:pi/10:p
6、i;y1=x.*cos(x);plot(x,y1,b*-),grid onhold onx=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)*sinx3 与 y3=exp(1/x)*sinx)xlabel(x 轴),ylabel(y 轴)gtext(y1=x.*cos(x),gtext(y2=x.*tan(1./x).*sin(x.3)gtext(y3=exp(1./x).*sin(x
7、)legend(y1=x.*cos(x),y2=x.*tan(1./x).*sin(x.3),y3=exp(1./x).*sin(x)axis(-pi 4*pi -4 4)六、 【实验结果】2015.5.7 练习：D1：D2：P391、 Untitled3n =500001ans =3.14163、 s=jiecheng(50)s =3.0414e+64P791、3、七、 【实验任务】画图练习：画出积分区域面 1,2，交线和在 xoy 平面上的投影P79 第七题第八题第九题八、 【实验程序】积分区域面 1：t=-1:.1:1;x,y=meshgrid(t);z1=x.2+2*y.2;subpl
8、ot(1,2,1),mesh(x,y,z1),title(曲面 z=x2+2y2),grid on积分区域面 2：t=-1:.1:1;x,y=meshgrid(t);z1=x.2+2*y.2;subplot(1,2,1),mesh(x,y,z1),title(曲面 z=x2+2y2)hold onz2=2-x.2;subplot(1,2,1),mesh(x,y,z2),title(曲面 z=x2+2y2 和曲面 z=2-x2)交线和在 xoy 平面上的投影t=-1:.1:1;x,y=meshgrid(t);z1=x.2+2*y.2;subplot(1,2,1),mesh(x,y,z1)hold
9、 onz2=2-x.2;subplot(1,2,1),mesh(x,y,z2),title(曲面 z=x2+2y2 和曲面 z=2-x2)r0=abs(z1-z2) A=1 2 3 4 5 6;2 3 4 5 6 7;3 4 5 6 7 8;4 5 6 7 8 9;1 4 5 6 3 4;3 5 4 6 7 5; AA =1 2 3 4 5 62 3 4 5 6 73 4 5 6 7 84 5 6 7 8 91 4 5 6 3 43 5 4 6 7 5 Aans =1 2 3 4 1 32 3 4 5 4 53 4 5 6 5 44 5 6 7 6 65 6 7 8 3 76 7 8 9 4
10、5 det(A)ans =0 rank(A)ans =4 rref(A)ans =1.0000 0 0 0 2.0000 2.00000 1.0000 0 0 -2.3333 -3.66670 0 1.0000 0 -2.3333 -0.66670 0 0 1.0000 3.6667 3.33330 0 0 0 0 00 0 0 0 0 0第 14 题 A=2 1 1;1 2 1;1 1 2; AA =2 1 11 2 11 1 2 p=poly(A)p =1.0000 -6.0000 9.0000 -4.0000 poly2str(p,x)ans =x3 - 6 x2 + 9 x - 4 d
11、=eig(A)d =1.00001.00004.0000 V,D=eig(A)V =0.4082 0.7071 0.57740.4082 -0.7071 0.5774-0.8165 0 0.5774D =1.0000 0 00 1.0000 0 0 0 4.0000第 21 题（1） A=1 1 2 -1;-1 1 3 0;2 -3 4 -1; rref(A)ans =1.0000 0 0 -0.56000 1.0000 0 -0.20000 0 1.0000 -0.1200A 的秩为 3，未知量个数 4，所以有无穷多个解。原方程对应的同解方程组为：取 =1得方程组的基础解系为：4321.05
12、6.x，所以方程组的通解为：12.12.056.432kx（2） A=1 -1 -1 1 0;1 -1 1 -3 1;1 -1 -2 3 -0.5A =1.0000 -1.0000 -1.0000 1.0000 01.0000 -1.0000 1.0000 -3.0000 1.00001.0000 -1.0000 -2.0000 3.0000 -0.5000 rref(A)ans =1.0000 -1.0000 0 -1.0000 0.50000 0 1.0000 -2.0000 0.50000 0 0 0 0A 的秩为 2，未知量个数 4，所以有无穷多个解。原方程对应的同解方程组为：可找到其
13、中一个特解为：21431x再求解对应的齐次线性方程组 ，可得到一个基础解系：0214321x，012因此，此方程组的通解为： Rccx2121432 ,00第 17 题function y=jifen(x)y=x.*sin(x)./(1+cos(x).2); x=0:pi/20:pi;y=x.*sin(x)./(1+cos(x).2);t=length(x);s1=sum(y(1:(t-1)*0.01s2=sum(y(2:t)*0.01s3=trapz(x,y)s4=quad(jifen,0,1)第 18 题： h=pi/20;x=0:h:pi/4;y=1./(1-sin(x);format longt
